Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

One of the critical variables that determine the rate of any reaction is temperature. For biological systems, the effects of temperature are convoluted with myriad (and often opposing) contributions from enzyme catalysis, protein stability, and temperature-dependent regulation, for example. We have coined the phrase "macromolecular rate theory (MMRT)" to describe the temperature dependence of enzyme-catalyzed rates independent of stability or regulatory processes. Central to MMRT is the observation that enzyme-catalyzed reactions occur with significant values of ΔCp(‡) that are in general negative. That is, the heat capacity (Cp) for the enzyme-substrate complex is generally larger than the Cp for the enzyme-transition state complex. Consistent with a classical description of enzyme catalysis, a negative value for ΔCp(‡) is the result of the enzyme binding relatively weakly to the substrate and very tightly to the transition state. This observation of negative ΔCp(‡) has important implications for the temperature dependence of enzyme-catalyzed rates. Here, we lay out the fundamentals of MMRT. We present a number of hypotheses that arise directly from MMRT including a theoretical justification for the large size of enzymes and the basis for their optimum temperatures. We rationalize the behavior of psychrophilic enzymes and describe a "psychrophilic trap" which places limits on the evolution of enzymes in low temperature environments. One of the defining characteristics of biology is catalysis of chemical reactions by enzymes, and enzymes drive much of metabolism. Therefore, we also expect to see characteristics of MMRT at the level of cells, whole organisms, and even ecosystems.

The Diels-Alder reaction, a [4 + 2] cycloaddition of a conjugated diene to a dienophile, is one of the most powerful reactions in synthetic chemistry. Biocatalysts capable of unlocking new and efficient Diels-Alder reactions would have major impact. Here we present a molecular-level description of the reaction mechanism of the spirotetronate cyclase AbyU, an enzyme shown here to be a bona fide natural Diels-Alderase. Using enzyme assays, X-ray crystal structures, and simulations of the reaction in the enzyme, we reveal how linear substrate chains are contorted within the AbyU active site to facilitate a transannular pericyclic reaction. This study provides compelling evidence for the existence of a natural enzyme evolved to catalyze a Diels-Alder reaction and shows how catalysis is achieved.

Combined quantum mechanics/molecular mechanics (QM/MM) calculations with high levels of correlated ab initio theory can now provide benchmarks for enzyme-catalyzed reactions. Here, we use such methods to test various QM/MM methods and the sensitivity of the results to details of the models for an important enzyme reaction, proton abstraction from acetyl-coenzyme A in citrate synthase. We calculate multiple QM/MM potential energy surfaces up to the local coupled cluster theory (LCCSD(T0)) level, with structures optimized at hybrid density functional theory and Hartree-Fock levels. The influence of QM methods, basis sets, and QM region size is shown to be significant. Correlated ab initio QM/MM calculations give barriers in agreement with experiment for formation of the acetyl-CoA enolate intermediate. In contrast, B3LYP fails to identify the enolate as an intermediate, whereas BH&HLYP does. The results indicate that QM/MM methods and setup should be tested, ideally using high-level calculations, to draw reliable mechanistic conclusions.

Carbapenems are the most potent β-lactam antibiotics and key drugs for treating infections by Gram-negative bacteria. In such organisms, β-lactam resistance arises principally from β-lactamase production. Although carbapenems escape the activity of most β-lactamases, due in the class A enzymes to slow deacylation of the covalent acylenzyme intermediate, carbapenem-hydrolyzing class A β-lactamases are now disseminating in clinically relevant bacteria. The reasons why carbapenems are substrates for these enzymes, but inhibit other class A β-lactamases, remain to be fully established. Here, we present crystal structures of the class A carbapenemase SFC-1 from Serratia fonticola and of complexes of its Ser70 Ala (Michaelis) and Glu166 Ala (acylenzyme) mutants with the carbapenem meropenem. These are the first crystal structures of carbapenem complexes of a class A carbapenemase. Our data reveal that, in the SFC-1 acylenzyme complex, the meropenem 6α-1R-hydroxyethyl group interacts with Asn132, but not with the deacylating water molecule. Molecular dynamics simulations indicate that this mode of binding occurs in both the Michaelis and acylenzyme complexes of wild-type SFC-1. In carbapenem-inhibited class A β-lactamases, it is proposed that the deacylating water molecule is deactivated by interaction with the carbapenem 6α-1R-hydroxyethyl substituent. Structural comparisons with such enzymes suggest that in SFC-1 subtle repositioning of key residues (Ser70, Ser130, Asn132 and Asn170) enlarges the active site, permitting rotation of the carbapenem 6α-1R-hydroxyethyl group and abolishing this contact. Our data show that SFC-1, and by implication other such carbapenem-hydrolyzing enzymes, uses Asn132 to orient bound carbapenems for efficient deacylation and prevent their interaction with the deacylating water molecule.

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 “grips” the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple “open and shut” case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Terpenoids form the largest and stereochemically most diverse class of natural products, and there is considerable interest in producing these by biocatalysis with whole cells or purified enzymes, and by metabolic engineering. The monoterpenes are an important class of terpenes and are industrially important as flavors and fragrances. We report here structures for the recently discovered Streptomyces clavuligerus monoterpene synthases linalool synthase (bLinS) and 1,8-cineole synthase (bCinS), and we show that these are active biocatalysts for monoterpene production using biocatalysis and metabolic engineering platforms. In metabolically engineered monoterpene-producing E. coli strains, use of bLinS leads to 300-fold higher linalool production compared with the corresponding plant monoterpene synthase. With bCinS, 1,8-cineole is produced with 96% purity compared to 67% from plant species. Structures of bLinS and bCinS, and their complexes with fluorinated substrate analogues, show that these bacterial monoterpene synthases are similar to previously characterized sesquiterpene synthases. Molecular dynamics simulations suggest that these monoterpene synthases do not undergo large-scale conformational changes during the reaction cycle, making them attractive targets for structured-based protein engineering to expand the catalytic scope of these enzymes toward alternative monoterpene scaffolds. Comparison of the bLinS and bCinS structures indicates how their active sites steer reactive carbocation intermediates to the desired acyclic linalool (bLinS) or bicyclic 1,8-cineole (bCinS) products. The work reported here provides the analysis of structures for this important class of monoterpene synthase. This should now guide exploitation of the bacterial enzymes as gateway biocatalysts for the production of other monoterpenes and monoterpenoids.

development, combining different levels of theory to increase accuracy, aiming to connect chemical and molecular changes to macroscopic observables. In this review, we outline biomolecular simulation methods and highlight examples of its application to investigate questions in biology. enzyme, membrane, molecular dynamics, multiscale, protein, QM/MM 1 | INTRODUCTION Biomolecular simulations are now making significant contributions to a wide variety of problems in drug discovery, drug development, biocatalysis, biotechnology, nanotechnology, chemical biology, and medicine. Biomolecular simulation is a rapidly growing field in scale and impact, increasingly demonstrating its worth in understanding mechanisms and analyzing activities, and contributing to the design of drugs and biocatalysts. Physics-based simulations complement experiments in building a molecular-level understanding of biology: They can test hypotheses and interpret and analyze experimental data in terms of interactions at the atomic level. Different types of simulation techniques have been developed, which are applicable to a range of different problems in biomolecular science. Simulations have already shown their worth in helping to analyze how enzymes catalyze biochemical reactions, and how proteins adopt their functional structures, for example, within cell membranes. They contribute to the design of drugs and catalysts, and in understanding the molecular basis of disease. Simulations have played a key role in developing the conceptual framework now at the heart of biomolecular science: that the dynamics of biological molecules is central to their function. Developing methods from chemical physics and computational science will open exciting new opportunities in biomolecular science, including in drug design and development, biotechnology, and biocatalysis. With high-performance computing resources, large-scale atomistic simulations of biological machines such as the ribosome, proton pumps and motors, membrane receptor complexes, and even whole viruses have become possible. Useful simulations of smaller systems can be carried out with desktop resources, thanks to developments allowing, for example, graphics processing units (GPUs) to be used. A particular challenge across the field is the integration of simulations crossing the span of length-and timescales as different types of simulation method are required for different types of problems. 1 Biomolecular systems pose fundamental scientific challenges (e.g., protein folding, enzyme catalysis, gene regulation, disease mechanisms, and antimicrobial resistance) and are at the heart of many advanced technological developments (drug discovery, biotechnology, biocatalysis, biomaterials, and genetic engineering). Biomolecular systems are inherently complex and pose significant challenges in modeling. An essential underlying paradigm is the need to consider biomolecular ensembles and their dynamics, rather than simply static biomolecular structures to understand and predict their behavior and propertie...