Role of conformational dynamics in the evolution of novel enzyme function

Unlocking the potential of enzyme engineering via rational computational design strategies

Enzymes play a pivotal role in various industries by enabling efficient, eco-friendly, and sustainable chemical processes. However, the low turnover rates and poor substrate selectivity of enzymes limit their large-scale applications. Rational computational enzyme design, facilitated by computational algorithms, offers a more targeted and less labor-intensive approach. There has been notable advancement in employing rational computational protein engineering strategies to overcome these issues, it has not been comprehensively reviewed so far. This article reviews recent developments in rational computational enzyme design, categorizing them into three types: structure-based, sequence-based, and data-driven machine learning computational design. Case studies are presented to demonstrate successful enhancements in catalytic activity, stability, and substrate selectivity. Lastly, the article provides a thorough analysis of these approaches, highlights existing challenges and potential solutions, and offers insights into future development directions.

Computation-guided engineering of distal mutations in an artificial enzyme

Artificial enzymes are valuable biocatalysts able to perform new-to-nature transformations with the precision and (enantio-)selectivity of natural enzymes. Although they are highly engineered biocatalysts, they often cannot reach catalytic rates akin those of their natural counterparts, slowing down their application in real-world industrial processes. Typically, their designs only optimise the chemistry inside the active site, while overlooking the role of protein dynamics on catalysis. In this work, we show how the catalytic performance of an already engineered artificial enzyme can be further improved by distal mutations that affect the conformational equilibrium of the protein. To this end, we subjected a specialised artificial enzyme based on the lactococcal multidrug resistance regulator (LmrR) to an innovative algorithm that quickly inspects the whole protein sequence space for hotpots which affect the protein dynamics. From an initial predicted selection of 73 variants, two variants with mutations distant by more than 11 Å from the catalytic pAF residue showed increased catalytic activity towards the new-to-nature hydrazone formation reaction. Their recombination displayed a 66% higher turnover number and 14 °C higher thermostability. Microsecond time scale molecular dynamics simulations evidenced a shift in the distribution of productive enzyme conformations, which are the result of a cascade of interactions initiated by the introduced mutations.

Kinetic View of Enzyme Catalysis from Enhanced Sampling QM/MM Simulations

The rate constants of enzyme-catalyzed reactions (kcat) are often approximated from the barrier height of the reactive step. We introduce an enhanced sampling QM/MM approach that directly calculates the kinetics of enzymatic reactions, without introducing the transition-state theory assumptions, and takes into account the dynamical equilibrium between the reactive and non-reactive conformations of the enzyme/substrate complex. Our computed kcat values are in order-of-magnitude agreement with the experimental data for two representative enzymatic reactions.

Catalytic Redundancies and Conformational Plasticity Drives Selectivity and Promiscuity in Quorum Quenching Lactonases

Several enzymes from the metallo-β-lactamase-like family of lactonases (MLLs) degrade N- acyl-L-homoserine lactones (AHLs). In doing so, they play a role in a microbial communication system, quorum sensing, which contributes to pathogenicity and biofilm formation. There is currently great interest in designing quorum quenching ( QQ ) enzymes that can interfere with this communication and be used in a range of industrial and biomedical applications. However, tailoring these enzymes for specific targets requires a thorough understanding of their mechanisms and the physicochemical properties that determine their substrate specificities. We present here a detailed biochemical, computational, and structural study of the MLL GcL, which is highly proficient, thermostable, and has broad substrate specificity. Strikingly, we show that GcL does not only accept a broad range of substrates but is also capable of utilizing different reaction mechanisms that are differentially used in function of the substrate structure or the remodeling of the active site via mutations. Comparison of GcL to other lactonases such as AiiA and AaL demonstrates similar mechanistic promiscuity, suggesting this is a shared feature across lactonases in this enzyme family. Mechanistic promiscuity has previously been observed in the lactonase/paraoxonase PON1, as well as with protein tyrosine phosphatases that operate via a dual general-acid mechanism. The apparent prevalence of this phenomenon is significant from both a biochemical and an engineering perspective: in addition to optimizing for specific substrates, it is possible to optimize for specific mechanisms, opening new doors not just for the design of novel quorum quenching enzymes, but also of other mechanistically promiscuous enzymes.

Engineering of Halide Methyltransferase BxHMT through Dynamic Cross-Correlation Network Analysis

Halide methyltransferases (HMTs) provide an effective way to regenerate S-adenosyl methionine (SAM) from S-adenosyl homocysteine and reactive electrophiles, such as methyl iodide (MeI) and methyl toluene sulfonate (MeOTs). As compared with MeI, the cost-effective unnatural substrate MeOTs can be accessed directly from cheap and abundant alcohols, but shows only limited reactivity in SAM production. In this study, we developed a dynamic cross-correlation network analysis (DCCNA) strategy for quickly identifying hot spots influencing the catalytic efficiency of the enzyme, and applied it to the evolution of HMT from Paraburkholderia xenovorans. Finally, the optimal mutant, M4 (V55T/C125S/L127T/L129P), exhibited remarkable improvement, with a specific activity of 4.08 U/mg towards MeOTs, representing an 82-fold increase as compared to the wild-type (WT) enzyme. Notably, M4 also demonstrated a positive impact on the catalytic ability with other methyl donors. The structural mechanism behind the enhanced enzyme activity was uncovered by molecular dynamics simulations. Our work not only contributes a promising biocatalyst for the regeneration of SAM, but also offers a strategy for efficient enzyme engineering.

Design of Efficient Artificial Enzymes Using Crystallographically Enhanced Conformational Sampling

The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing the catalytic efficiency of de novo enzymes by several orders of magnitude without relying on directed evolution and high-throughput screening. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models, with catalytic contacts present as designed and transition-state root-mean-square deviations of ≤0.65 Å. Our work shows how ensemble-based design can generate efficient artificial enzymes by exploiting the true conformational ensemble to design improved active sites.

Epistasis arises from shifting the rate-limiting step during enzyme evolution of a β-lactamase

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a β-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

The shortest path method (SPM) webserver for computational enzyme design

SPMweb is the online webserver of the Shortest Path Map (SPM) tool for identifying the key conformationally-relevant positions of a given enzyme structure and dynamics. The server is built on top of the DynaComm.py code and enables the calculation and visualization of the SPM pathways. SPMweb is easy-to-use as it only requires three input files: the three-dimensional structure of the protein of interest, and the two matrices (distance and correlation) previously computed from a Molecular Dynamics simulation. We provide in this publication information on how to generate the files for SPM construction even for non-expert users and discuss the most relevant parameters that can be modified. The tool is extremely fast (it takes less than one minute per job), thus allowing the rapid identification of distal positions connected to the active site pocket of the enzyme. SPM applications expand from computational enzyme design, especially if combined with other tools to identify the preferred substitution at the identified position, but also to rationalizing allosteric regulation, and even cryptic pocket identification for drug discovery. The simple user interface and setup make the SPM tool accessible to the whole scientific community. SPMweb is freely available for academia at http://spmosuna.com/.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Binding of Cholesterol to the N-Terminal Domain of the NPC1L1 Transporter: Analysis of the Epimerization-Related Binding Selectivity and Loop Mutations

Cholesterol is a fat-like substance with a pivotal physiological relevance in humans, and its homeostasis is tightly regulated by various cellular processes, including the import in the small intestine and the reabsorption in the biliary ducts by the Niemann-Pick C1 Like 1 (NPC1L1) importer. NPC1L1 can mediate the absorption of a variety of sterols but strikingly exhibits a large sensitivity to cholesterol epimerization. This study examines the molecular basis of the epimerization-related selective binding of cholesterol by combining extended unbiased molecular dynamics simulations of the apo and holo species of the N-terminal domain of wild-type NPC1L1, in conjunction with relative binding free energy, umbrella sampling, and well-tempered metadynamics calculations. The analysis of the results discloses the existence of two distinct binding modes for cholesterol and epi-cholesterol. The former binds deeper in the cavity, forming key hydrogen-bond interactions with Q95, S56, and a water molecule. In contrast, epi-cholesterol is shifted ca. 3 Å to the mouth of the cavity and the transition to the Q95 site is prevented by an energetic barrier of 4.1 kcal·mol-1. Thus, the configuration of the hydroxyl group of cholesterol, together with the presence of a structural water molecule, is a key feature for effective absorption. Finally, whereas these findings may seemingly be challenged by single-point mutations that impair cholesterol transport but have a mild impact on the binding of cholesterol to the Q95 binding site, our results reveal that they have a drastic influence on the conformational landscape of the α8/β7 loop in the apo species compared to the wild-type protein. Overall, the results give support to the functional role played by the α8/β7 loop in regulating the access of ligands to NPC1L1, and hence to interpreting the impact of these mutations on diseases related to disruption of sterol absorption, paving the way to understanding certain physiological dysfunctions.

Pervasive epistasis exposes intramolecular networks in adaptive enzyme evolution

Enzyme evolution is characterized by constant alterations of the intramolecular residue networks supporting their functions. The rewiring of these network interactions can give rise to epistasis. As mutations accumulate, the epistasis observed across diverse genotypes may appear idiosyncratic, that is, exhibit unique effects in different genetic backgrounds. Here, we unveil a quantitative picture of the prevalence and patterns of epistasis in enzyme evolution by analyzing 41 fitness landscapes generated from seven enzymes. We show that >94% of all mutational and epistatic effects appear highly idiosyncratic, which greatly distorted the functional prediction of the evolved enzymes. By examining seemingly idiosyncratic changes in epistasis along adaptive trajectories, we expose several instances of higher-order, intramolecular rewiring. Using complementary structural data, we outline putative molecular mechanisms explaining higher-order epistasis along two enzyme trajectories. Our work emphasizes the prevalence of epistasis and provides an approach to exploring this phenomenon through a molecular lens.

Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability

Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.

Design of efficient artificial enzymes using crystallographically-enhanced conformational sampling

The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing catalytic efficiency of de novo enzymes to a level comparable to their natural counterparts without relying on directed evolution. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models. Our work shows how computational design can generate efficient artificial enzymes by exploiting the true conformational ensemble to more effectively stabilize the transition state.

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Glycosylation increases active site rigidity leading to improved enzyme stability and turnover

Glycosylation is the most prevalent protein post-translational modification, with a quarter of glycosylated proteins having enzymatic properties. Yet, the full impact of glycosylation on the protein structure-function relationship, especially in enzymes, is still limited. Here, we show that glycosylation rigidifies the important commercial enzyme horseradish peroxidase (HRP), which in turn increases its turnover and stability. Circular dichroism spectroscopy revealed that glycosylation increased holo-HRP's thermal stability and promoted significant helical structure in the absence of haem (apo-HRP). Glycosylation also resulted in a 10-fold increase in enzymatic turnover towards o-phenylenediamine dihydrochloride when compared to its nonglycosylated form. Utilising a naturally occurring site-specific probe of active site flexibility (Trp117) in combination with red-edge excitation shift fluorescence spectroscopy, we found that glycosylation significantly rigidified the enzyme. In silico simulations confirmed that glycosylation largely decreased protein backbone flexibility, especially in regions close to the active site and the substrate access channel. Thus, our data show that glycosylation does not just have a passive effect on HRP stability but can exert long-range effects that mediate the 'native' enzyme's activity and stability through changes in inherent dynamics.

AlphaFold2 and Deep Learning for Elucidating Enzyme Conformational Flexibility and Its Application for Design

The recent success of AlphaFold2 (AF2) and other deep learning (DL) tools in accurately predicting the folded three-dimensional (3D) structure of proteins and enzymes has revolutionized the structural biology and protein design fields. The 3D structure indeed reveals key information on the arrangement of the catalytic machinery of enzymes and which structural elements gate the active site pocket. However, comprehending enzymatic activity requires a detailed knowledge of the chemical steps involved along the catalytic cycle and the exploration of the multiple thermally accessible conformations that enzymes adopt when in solution. In this Perspective, some of the recent studies showing the potential of AF2 in elucidating the conformational landscape of enzymes are provided. Selected examples of the key developments of AF2-based and DL methods for protein design are discussed, as well as a few enzyme design cases. These studies show the potential of AF2 and DL for allowing the routine computational design of efficient enzymes.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

Tuning Selectivity in CalA Lipase: Beyond Tunnel Engineering

Engineering studies of Candida (Pseudozyma) antarctica lipase A (CalA) have demonstrated the potential of this enzyme in the selective hydrolysis of fatty acid esters of different chain lengths. CalA has been shown to bind substrates preferentially through an acyl-chain binding tunnel accessed via the hydrolytic active site; it has also been shown that selectivity for substrates of longer or shorter chain length can be tuned, for instance by modulating steric hindrance within the tunnel. Here we demonstrate that, whereas the tunnel region is certainly of paramount importance for substrate recognition, residues in distal regions of the enzyme can also modulate substrate selectivity. To this end, we investigate variants that carry one or more substitutions within the substrate tunnel as well as in distal regions. Combining experimental determination of the substrate selectivity using natural and synthetic substrates with computational characterization of protein dynamics and of tunnels, we deconvolute the effect of key substitutions and demonstrate that epistatic interactions contribute to procuring selectivity toward either long-chain or short/medium-chain fatty acid esters. We demonstrate that various mechanisms contribute to the diverse selectivity profiles, ranging from reshaping tunnel morphology and tunnel stabilization to obstructing the main substrate-binding tunnel, highlighting the dynamic nature of the substrate-binding region. This work provides important insights into the versatility of this robust lipase toward diverse applications.

Mechanism-Guided Computational Design of ω-Transaminase by Reprograming of High-Energy-Barrier Steps

ω-Transaminases (ω-TAs) show considerable potential for the synthesis of chiral amines. However, their low catalytic efficiency towards bulky substrates limits their application, and complicated catalytic mechanisms prevent precise enzyme design. Herein, we address this challenge using a mechanism-guided computational enzyme design strategy by reprograming the transition and ground states in key reaction steps. The common features among the three high-energy-barrier steps responsible for the low catalytic efficiency were revealed using quantum mechanics (QM). Five key residues were simultaneously tailored to stabilize the rate-limiting transition state with the aid of the Rosetta design. The 14 top-ranked variants showed 16.9-143-fold improved catalytic activity. The catalytic efficiency of the best variant, M9 (Q25F/M60W/W64F/I266A), was significantly increased, with a 1660-fold increase in k_cat /K_m and a 1.5-26.8-fold increase in turnover number (TON) towards various indanone derivatives.

Stochastic Lag Time Parameterization for Markov State Models of Protein Dynamics

Markov state models (MSMs) play a key role in studying protein conformational dynamics. A sliding count window with a fixed lag time is widely used to sample sub-trajectories for transition counting and MSM construction. However, sub-trajectories sampled with a fixed lag time may not perform well under different selections of lag time, which requires strong prior practice and leads to less robust estimation. To alleviate it, we propose a novel stochastic method from a Poisson process to generate perturbative lag time for sub-trajectory sampling and utilize it to construct a Markov chain. Comprehensive evaluations on the double-well system, WW domain, BPTI, and RBD-ACE2 complex of SARS-CoV-2 reveal that our algorithm significantly increases the robustness and power of a constructed MSM without disturbing the Markovian properties. Furthermore, the superiority of our algorithm is amplified for slow dynamic modes in complex biological processes.

Unravelling Enzymatic Features in a Supramolecular Iridium Catalyst by Computational Calculations

Non-biological catalysts following the governing principles of enzymes are attractive systems to disclose unprecedented reactivities. Most of those existing catalysts feature an adaptable molecular recognition site for substrate binding that are prone to undergo conformational selection pathways. Herein, we present a non-biological catalyst that is able to bind substrates via the induced fit model according to in-depth computational calculations. The system, which is constituted by an inflexible substrate-recognition site derived from a zinc-porphyrin in the second coordination sphere, features destabilization of ground states as well as stabilization of transition states for the relevant iridium-catalyzed C-H bond borylation of pyridine. In addition, this catalyst appears to be most suited to tightly bind the transition state rather than the substrate. Besides these features, which are reminiscent of the action modes of enzymes, new elementary catalytic steps (i. e. C-B bond formation and catalyst regeneration) have been disclosed owing to the unique distortions encountered in the different intermediates and transition states.

Estimating conformational heterogeneity of tryptophan synthase with a template‐based Alphafold2 approach

The three‐dimensional structure of the enzymes provides very relevant information on the arrangement of the catalytic machinery and structural elements gating the active site pocket. The recent success of the neural network Alphafold2 in predicting the folded structure of proteins from the primary sequence with high levels of accuracy has revolutionized the protein design field. However, the application of Alphafold2 for understanding and engineering function directly from the obtained single static picture is not straightforward. Indeed, understanding enzymatic function requires the exploration of the ensemble of thermally accessible conformations that enzymes adopt in solution. In the present study, we evaluate the potential of Alphafold2 in assessing the effect of the mutations on the conformational landscape of the beta subunit of tryptophan synthase (TrpB). Specifically, we develop a template‐based Alphafold2 approach for estimating the conformational heterogeneity of several TrpB enzymes, which is needed for enhanced stand‐alone activity. Our results show the potential of Alphafold2, especially if combined with molecular dynamics simulations, for elucidating the changes induced by mutation in the conformational landscapes at a rather reduced computational cost, thus revealing its plausible application in computational enzyme design.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation (TDHDX-MS) has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend TDHDX from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine or DAA) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different HDX-MS properties for the two ligands, as expected from their 10⁶-fold differences in strength of binding. By contrast, analogous patterns for TDHDX-MS emerge in the presence of both DAA and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended TDHDX to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility vs. site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

Distal Mutations Shape Substrate-Binding Sites during Evolution of a Metallo-Oxidase into a Laccase

Laccases are in increasing demand as innovative solutions in the biorefinery fields. Here, we combine mutagenesis with structural, kinetic, and in silico analyses to characterize the molecular features that cause the evolution of a hyperthermostable metallo-oxidase from the multicopper oxidase family into a laccase (k _cat 273 s^-1 for a bulky aromatic substrate). We show that six mutations scattered across the enzyme collectively modulate dynamics to improve the binding and catalysis of a bulky aromatic substrate. The replacement of residues during the early stages of evolution is a stepping stone for altering the shape and size of substrate-binding sites. Binding sites are then fine-tuned through high-order epistasis interactions by inserting distal mutations during later stages of evolution. Allosterically coupled, long-range dynamic networks favor catalytically competent conformational states that are more suitable for recognizing and stabilizing the aromatic substrate. This work provides mechanistic insight into enzymatic and evolutionary molecular mechanisms and spots the importance of iterative experimental and computational analyses to understand local-to-global changes.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Computational Enzyme Design at Zymvol

Directed evolution is the most recognized methodology for enzyme engineering. The main drawback resides in its random nature and in the limited sequence exploration; both require screening of thousands (if not millions) of variants to achieve a target function. Computer-driven approaches can limit laboratorial screening to a few hundred candidates, enabling and accelerating the development of industrial enzymes. In this book chapter, the technology adopted at Zymvol is described. An overview of the current development and future directions in the company is also provided.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

Analytical and functional aspects of protein-ligand interactions: Beyond induced fit and conformational selection

Ligand-dependent changes in protein conformation are foundational to biology. Historical mechanistic models for substrate-specific proteins are induced fit (IF) and conformational selection (CS), which invoke a change in protein conformation after ligand binds or before ligand binds, respectively. These mechanisms have important, but rarely discussed, functional relevance because IF vs. CS can differentially affect a protein's substrate specificity or promiscuity, and its regulatory properties. The modern view of proteins as conformational ensembles in both ligand free and bound states, together with the realization that most proteins exhibit some substrate promiscuity, demands a deeper interpretation of the historical models and provides an opportunity to improve mechanistic analyses. Here we describe alternative analytical strategies for distinguishing the historical models, including the more complex expanded versions of IF and CS. Functional implications of the different models are described. We provide an alternative perspective based on protein ensembles interacting with ligand ensembles that clarifies how a single protein can 'apparently' exploit different mechanisms for different ligands. Mechanistic information about protein ensembles can be optimized when they are probed with multiple ligands.

In Silico Identification and Experimental Validation of Distal Activity-Enhancing Mutations in Tryptophan Synthase

Allostery is a central mechanism for the regulation of multi-enzyme complexes. The mechanistic basis that drives allosteric regulation is poorly understood but harbors key information for enzyme engineering. In the present study, we focus on the tryptophan synthase complex that is composed of TrpA and TrpB subunits, which allosterically activate each other. Specifically, we develop a rational approach for identifying key amino acid residues of TrpB distal from the active site. Those residues are predicted to be crucial for shifting the inefficient conformational ensemble of the isolated TrpB to a productive ensemble through intra-subunit allosteric effects. The experimental validation of the conformationally driven TrpB design demonstrates its superior stand-alone activity in the absence of TrpA, comparable to those enhancements obtained after multiple rounds of experimental laboratory evolution. Our work evidences that the current challenge of distal active site prediction for enhanced function in computational enzyme design has become within reach.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

A distal regulatory strategy of enzymes: from local to global conformational dynamics

Modulating the distribution of various states in protein ensembles through distal sites may be promising in the evolution of enzymes in desired directions. However, the prediction of distal mutation hotspots that stabilize the favoured states from a computational perspective remains challenging. Here, we presented a strategy based on molecular dynamics (MD) and Markov state models (MSM) to predict distal mutation sites. Extensive MD combined with MSM was applied to determine the principally distributed metastable states interconverting at a slow timescale. Then, molecular docking was used to classify these states into active states and inactive ones. Distal mutation hotspots were targeted based on comparing the conformational features between active and inactive states, where mutations destabilize the inactive states and show little influence on the active state. The proposed strategy was used to explore the highly dynamic MHETase, which shows a potential application in the biodegradation of poly(ethylene terephthalate) (PET). Seven principally populated interrelated metastable states were identified, and the atomistic picture of their conformational changes was unveiled. Several residues at distal positions were found to adopt more H-bond occupancies in inactive states than active states, making them potential mutation hotspots for stabilizing the favoured conformations. In addition, the detailed mechanism revealed the significance of calcium ions at a distance from the catalytic centre in reshaping the free energy landscape. This study deepens the understanding of the conformational dynamics of α/β hydrolases containing a lid domain and advances the study of enzymatic plastic degradation.

Molecular Mechanism of Methanol Inhibition in CALB-Catalyzed Alcoholysis: Analyzing Molecular Dynamics Simulations by a Markov State Model

Lipases are widely used enzymes that catalyze hydrolysis and alcoholysis of fatty acid esters. At high concentrations of small alcohols such as methanol or ethanol, many lipases are inhibited by the substrate. The molecular basis of the inhibition of Candida antarctica lipase B (CALB) by methanol was investigated by unbiased molecular dynamics (MD) simulations, and the substrate binding kinetics was analyzed by Markov state models (MSMs). The modeled fluxes of productive methanol binding at concentrations between 50 mM and 5.5 M were in good agreement with the experimental activity profile of CALB, with a peak at 300 mM. The kinetic and structural analysis uncovered the molecular basis of CALB inhibition. Beyond 300 mM, the kinetic bottleneck results from crowding of methanol in the substrate access channel, which is caused by the gradual formation of methanol patches close to Leu140 (helix α5), Leu278, and Ile285 (helix α10) at a distance of 4-5 Å from the active site. Our findings demonstrate the usefulness of unbiased MD simulations to study enzyme-substrate interactions at realistic substrate concentrations and the feasibility of scale-bridging by an MSM analysis to derive kinetic information.

Evolution of dynamical networks enhances catalysis in a designer enzyme

Activation heat capacity is emerging as a crucial factor in enzyme thermoadaptation, as shown by the non-Arrhenius behaviour of many natural enzymes. However, its physical origin and relationship to the evolution of catalytic activity remain uncertain. Here we show that directed evolution of a computationally designed Kemp eliminase reshapes protein dynamics, which gives rise to an activation heat capacity absent in the original design. These changes buttress transition-state stabilization. Extensive molecular dynamics simulations show that evolution results in the closure of solvent-exposed loops and a better packing of the active site. Remarkably, this gives rise to a correlated dynamical network that involves the transition state and large parts of the protein. This network tightens the transition-state ensemble, which induces a negative activation heat capacity and non-linearity in the activity-temperature dependence. Our results have implications for understanding enzyme evolution and suggest that selectively targeting the conformational dynamics of the transition-state ensemble by design and evolution will expedite the creation of novel enzymes.

Rate-Perturbing Single Amino Acid Mutation for Hydrolases: A Statistical Profiling

Hydrolases are a critical component for modern chemical, pharmaceutical, and environmental sciences. Identifying mutations that enhance catalytic efficiency presents a roadblock to design and to discover new hydrolases for broad academic and industrial uses. Here, we report the statistical profiling for rate-perturbing mutant hydrolases with a single amino acid substitution. We constructed an integrated structure-kinetics database for hydrolases, IntEnzyDB, which contains 3907 k_cats, 4175 K_Ms, and 2715 Protein Data Bank IDs. IntEnzyDB adopts a relational architecture with a flattened data structure, enabling facile and efficient access to clean and tabulated data for machine learning uses. We conducted statistical analyses on how single amino acids mutations influence the turnover number (i.e., k_cat) and efficiency (i.e., k_cat/K_M), with a particular emphasis on profiling the features for rate-enhancing mutations. The results show that mutation to bulky nonpolar residues with a hydrocarbon chain involves a higher likelihood for rate acceleration than to other types of residues. Linear regression models reveal geometric descriptors of substrate and mutation residues that mediate rate-perturbing outcomes for hydrolases with bulky nonpolar mutations. On the basis of the analyses of the structure-kinetics relationship, we observe that the propensity for rate enhancement is independent of protein sizes. In addition, we observe that distal mutations (i.e., >10 Å from the active site) in hydrolases are significantly more prone to induce efficiency neutrality and avoid efficiency deletion but involve similar propensity for rate enhancement. The studies reveal the statistical features for identifying rate-enhancing mutations in hydrolases, which will potentially guide hydrolase discovery in biocatalysis.

Rigidifying a De Novo Enzyme Increases Activity and Induces a Negative Activation Heat Capacity

Conformational sampling profoundly impacts the overall activity and temperature dependence of enzymes. Peroxidases have emerged as versatile platforms for high-value biocatalysis owing to their broad palette of potential biotransformations. Here, we explore the role of conformational sampling in mediating activity in the de novo peroxidase C45. We demonstrate that 2,2,2-triflouoroethanol (TFE) affects the equilibrium of enzyme conformational states, tending toward a more globally rigid structure. This is correlated with increases in both stability and activity. Notably, these effects are concomitant with the emergence of curvature in the temperature-activity profile, trading off activity gains at ambient temperature with losses at high temperatures. We apply macromolecular rate theory (MMRT) to understand enzyme temperature dependence data. These data point to an increase in protein rigidity associated with a difference in the distribution of protein dynamics between the ground and transition states. We compare the thermodynamics of the de novo enzyme activity to those of a natural peroxidase, horseradish peroxidase. We find that the native enzyme resembles the rigidified de novo enzyme in terms of the thermodynamics of enzyme catalysis and the putative distribution of protein dynamics between the ground and transition states. The addition of TFE apparently causes C45 to behave more like the natural enzyme. Our data suggest robust, generic strategies for improving biocatalytic activity by manipulating protein rigidity; for functional de novo protein catalysts in particular, this can provide more enzyme-like catalysts without further rational engineering, computational redesign, or directed evolution.

Exploiting enzyme evolution for computational protein design

Recent years have seen an explosion of interest in understanding the physicochemical parameters that shape enzyme evolution, as well as substantial advances in computational enzyme design. This review discusses three areas where evolutionary information can be used as part of the design process: (i) using ancestral sequence reconstruction (ASR) to generate new starting points for enzyme design efforts; (ii) learning from how nature uses conformational dynamics in enzyme evolution to mimic this process in silico; and (iii) modular design of enzymes from smaller fragments, again mimicking the process by which nature appears to create new protein folds. Using showcase examples, we highlight the importance of incorporating evolutionary information to continue to push forward the boundaries of enzyme design studies.

A simple approach for reconstruction of non-uniformly sampled pseudo-3D NMR data for accurate measurement of spin relaxation parameters

We explain how to conduct a pseudo-3D relaxation series NUS measurement so that it can be reconstructed by existing 3D NUS reconstruction methods to give accurate relaxation values. We demonstrate using reconstruction algorithms IST and SMILE that this 3D approach allows lower sampling densities than for independent 2D reconstructions. This is in keeping with the common finding that higher dimensionality increases signal sparsity, enabling lower sampling density. The approach treats the relaxation series as ordinary 3D time-domain data whose imaginary part in the pseudo-dimension is zero, and applies any suitably linear 3D NUS reconstruction method accordingly. Best results on measured and simulated data were achieved using acquisitions with 9 to 12 planes and exponential spacing in the pseudo-dimension out to ~ 2 times the inverse decay time. Given these criteria, in typical cases where 2D reconstructions require 50% sampling, the new 3D approach generates spectra reliably at sampling densities of 25%.

Reaction coupling in ADH1A alcohol dehydrogenase enzyme by exciplex formation with adenosine diphosphate moderated by low-energy electronic excited states

Two commonly accepted theories about enzymes were revisited. The first states that adenosine triphosphate (ATP)-stored energy is only released when the substrate is in place, because the substrate changes the enzyme structure when it is bound to the enzyme. In fact, as demonstrated and discussed presently, no structural changes are required, and ATP-stored energy is released when it can be used. The second states that ATP-released energy moves along the enzyme molecule in the form of molecular vibrations (Davydov's vibrational solitons). In fact, as reported presently, energy released upon ATP hydrolysis moves in the form of excited-state electrons (excitons), with no molecular vibrations involved. The relevant experimental evidence was obtained for the human ADH1A alcohol dehydrogenase enzyme. Spontaneous ATP hydrolysis in the absence of substrate was apparently prevented by electronically excited enzyme + adenosine diphosphate (ADP) + inorganic phosphate (P) complex (exciplex) formed upon ATP hydrolysis. This exciplex kept ADP + P bound and in place for the inverse reaction, until the excess energy was dissipated in the enzyme-catalyzed reaction or by energy transfer to a suitable acceptor. Additionally, and contrary to textbooks, ADH1A has required ATP, working orders of magnitude faster in its presence.