Javier Iglesias-Fernández

Hydroxynitrile lyase from rubber tree ( Hb HNL) shares 45% identical amino acid residues with the homologous esterase from tobacco, SABP2, but the two enzymes catalyze different reactions. The x-ray structures reveal a serine-histidine-aspartate catalytic triad in both enzymes along with several differing amino acid residues within the active site. Previous exchange of three amino acid residues in the active site of Hb HNL with the corresponding amino acid residue in SABP2 (T11G-E79H-K236M) created variant HNL3, which showed low esterase activity toward p-nitrophenyl acetate. Further structure comparison reveals additional differences surrounding the active site. Hb HNL contains an improperly positioned oxyanion hole residue and differing solvation of the catalytic aspartate. We hypothesized that correcting these structural differences would impart good esterase activity on the corresponding HbHNL variant. To predict the amino acid substitutions needed to correct the structure, we calculated shortest path maps for both Hb HNL and SABP2, which reveal correlated movements of amino acids in the two enzymes. Replacing four amino acid residues (C81L-N104T-V106F-G176S) whose movements are connected to the movements of the catalytic residues yielded variant HNL7TV (stabilizing substitution H103V was also added), which showed an esterase catalytic efficiency comparable to that of SABP2. The x-ray structure of an intermediate variant, HNL6V, showed an altered solvation of the catalytic aspartate and a partially corrected oxyanion hole. This dramatic increase in catalytic efficiency demonstrates the ability of shortest path maps to predict which residues outside the active site contribute to catalytic activity.

Halohydrin dehalogenases (HHDH) are industrially relevant biocatalysts exhibiting a promiscuous epoxide-ring opening reactivity in the presence of small nucleophiles, thus giving access to novel carbon–carbon, carbon–oxygen, carbon–nitrogen, and carbon–sulfur bonds. Recently, the repertoire of HHDH has been expanded, providing access to some novel HHDH subclasses exhibiting a broader epoxide substrate scope. In this work, we develop a computational approach based on the application of linear and non-linear dimensionality reduction techniques to long time-scale Molecular Dynamics (MD) simulations to study the HHDH conformational landscapes. We couple the analysis of the conformational landscapes to CAVER calculations to assess their impact on the active site tunnels and potential ability towards bulky epoxide ring opening reaction. Our study indicates that the analyzed HHDHs subclasses share a common breathing motion of the halide binding pocket, but present large deviations in the loops adjacent to the active site pocket and N-terminal regions. Such conformational differences affect the available tunnels for epoxide binding to the active site. The superior activity of the HHDH G subclass towards bulkier substrates is explained by the additional structural elements delimiting the active site region, its rich conformational heterogeneity, and the substantially wider and frequently observed active site tunnels. This study therefore provides key information for HHDH promiscuity and engineering.

Multimeric enzyme complexes are ubiquitous in nature and catalyze a broad range of useful biological transformations. They are often characterized by a tight allosteric coupling between subunits, making them highly inefficient when isolated. A good example is Tryptophan synthase (TrpS), an allosteric heterodimeric enzyme in the form of an αββα complex that catalyzes the biosynthesis of L-tryptophan. In this study, we decipher the allosteric regulation existing in TrpS from Pyrococcus furiosus (PfTrpS), and how the allosteric conformational ensemble is recovered in laboratory-evolved stand-alone β-subunit variants. We find that recovering the conformational ensemble of a subdomain of TrpS affecting the relative stabilities of open, partially closed, and closed conformations is a prerequisite for enhancing the catalytic efficiency of the β-subunit in the absence of its binding partner. The distal mutations resuscitate the allosterically driven conformational regulation and alter the populations and rates of exchange between these multiple conformational states, which are essential for the multistep reaction pathway of the enzyme. Interestingly, these distal mutations can be a priori predicted by careful analysis of the conformational ensemble of the TrpS enzyme through computational methods. Our study provides the enzyme design field with a rational approach for evolving allosteric enzymes toward improved stand-alone function for biosynthetic applications.

Modulation of protein–protein interactions (PPIs) is essential for understanding and tuning biologically relevant processes. Although inhibitors for PPIs are widely used, the field still lacks the targeted design of stabilizers. Here, we report unnatural stabilizers based on the combination of multivalency effects and the artificial building block guanidiniocarbonylpyrrol (GCP), an arginine mimetic. Unlike other GCP‐based ligands that modulate PPIs in different protein targets, only a tetrameric design shows potent activity as stabilizer of the 14‐3‐3ζ/C‐Raf and 14‐3‐3ζ/Tau complexes in the low‐micromolar range. This evidences the role of multivalency for achieving higher specificity in the modulation of PPIs.

Enzymes exist as an ensemble of conformational states, whose populations can be shifted by substrate binding, allosteric interactions, but also by introducing mutations to their sequence. Tuning the populations of the enzyme conformational states through mutation enables evolution towards novel activity. Herein, Markov state models are used to unveil hidden conformational states of monoamine oxidase from Aspergillus niger (MAO‐N). These hidden conformations, not previously observed by any other technique, play a crucial role in substrate binding and enzyme activity. This reveals how distal mutations regulate MAO‐N activity by stabilizing these hidden, catalytically important conformational states, but also by modulating the communication pathway between both MAO‐N subunits.

Biosynthesis of glycogen, the essential glucose (and hence energy) storage molecule in humans, animals and fungi¹, is initiated by the glycosyltransferase enzyme, glycogenin (GYG). Deficiencies in glycogen formation cause neurodegenerative and metabolic disease^2,3,4, and mouse knockout⁵ and inherited human mutations⁶ of GYG impair glycogen synthesis. GYG acts as a ‘seed core’ for the formation of the glycogen particle by catalysing its own stepwise autoglucosylation to form a covalently bound gluco-oligosaccharide chain at initiation site Tyr 195. Precise mechanistic studies have so far been prevented by an inability to access homogeneous glycoforms of this protein, which unusually acts as both catalyst and substrate. Here we show that unprecedented direct access to different, homogeneously glucosylated states of GYG can be accomplished through a palladium-mediated enzyme activation ‘shunt’ process using on-protein C–C bond formation. Careful mimicry of GYG intermediates recapitulates catalytic activity at distinct stages, which in turn allows discovery of triphasic kinetics and substrate plasticity in GYG’s use of sugar substrates. This reveals a tolerant but ‘proof-read’ mechanism that underlies the precision of this metabolic process. The present demonstration of direct, chemically controlled access to intermediate states of active enzymes suggests that such ligation-dependent activation could be a powerful tool in the study of mechanism.

Enzymes exist as a dynamic ensemble of conformations, each potentially playing a key role in substrate binding, the chemical transformation, or product release. We discuss recent advances in the evaluation of the enzyme conformational dynamics and its evolution towards new functions or substrate preferences.

Directed evolution involves iterative rounds of genetic diversification followed by selection of variants with improved enzyme properties. A common trend in these studies is the introduction of multiple scattered mutations remote from the active site. Rationalizing how distal mutations influence the conformational states or ensemble of conformations formed by enzymes, and its relation to enzymatic catalysis is highly appealing for the improvement of current rational enzyme design strategies. Directed evolution was applied to convert d‐sialic acid aldolase into an efficient l‐3‐deoxy‐manno‐2‐octulosonic acid (l‐KDO) aldolase. In this study, we computationally evaluate the conformational dynamics of both enzymes to shed light on the specificity of the evolved variant. We further demonstrate the role of distal mutations on the modulation of enzyme conformational dynamics and its relation to substrate accessibility and selectivity. Mutations markedly altered the active site shape and substrate access tunnels in the evolved l‐KDO aldolase, thus affecting the enzyme specificity.