Biocatalytic Routes to Lactone Monomers for Polymer Production

Monoterpenoids offer potential as biocatalytically derived monomer feedstocks for high-performance renewable polymers. We describe a biocatalytic route to lactone monomers menthide and dihydrocarvide employing Baeyer-Villiger monooxygenases (BVMOs) from Pseudomonas sp. HI-70 (CPDMO) and Rhodococcus sp. Phi1 (CHMOPhi1) as an alternative to organic synthesis. The regioselectivity of dihydrocarvide isomer formation was controlled by site-directed mutagenesis of three key active site residues in CHMOPhi1. A combination of crystal structure determination, molecular dynamics simulations, and mechanistic modeling using density functional theory on a range of models provides insight into the origins of the discrimination of the wild type and a variant CHMOPhi1 for producing different regioisomers of the lactone product. Ring-opening polymerizations of the resultant lactones using mild metal-organic catalysts demonstrate their utility in polymer production. This semisynthetic approach utilizing a biocatalytic step, non-petroleum feedstocks, and mild polymerization catalysts allows access to known and also to previously unreported and potentially novel lactone monomers and polymers.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

Dynamical origins of heat capacity changes in enzyme-catalysed reactions

Heat capacity changes are emerging as essential for explaining the temperature dependence of enzyme-catalysed reaction rates. This has important implications for enzyme kinetics, thermoadaptation and evolution, but the physical basis of these heat capacity changes is unknown. Here we show by a combination of experiment and simulation, for two quite distinct enzymes (dimeric ketosteroid isomerase and monomeric alpha-glucosidase), that the activation heat capacity change for the catalysed reaction can be predicted through atomistic molecular dynamics simulations. The simulations reveal subtle and surprising underlying dynamical changes: tightening of loops around the active site is observed, along with changes in energetic fluctuations across the whole enzyme including important contributions from oligomeric neighbours and domains distal to the active site. This has general implications for understanding enzyme catalysis and demonstrating a direct connection between functionally important microscopic dynamics and macroscopically measurable quantities.

Heat capacity changes affect the temperature dependence of enzyme catalysis, with implications for thermoadaptation, however their physical basis is unknown. Here the authors show that heat capacity changes are calculable by simulation, revealing distinct dynamical contributions from regions remote from the active site.

L718Q mutant EGFR escapes covalent inhibition by stabilizing a non-reactive conformation of the lung cancer drug osimertinib

Impact of L718Q mutation on the inhibitory activity of osimertinib on EGFR revealed by free-energy simulations.

Osimertinib is a third-generation inhibitor approved for the treatment of non-small cell lung cancer. It overcomes resistance to first-generation inhibitors by incorporating an acrylamide group which alkylates Cys797 of EGFR T790M. The mutation of a residue in the P-loop (L718Q) was shown to cause resistance to osimertinib, but the molecular mechanism of this process is unknown. Here, we investigated the inhibitory process for EGFR T790M (susceptible to osimertinib) and EGFR T790M/L718Q (resistant to osimertinib), by modelling the chemical step (i.e., alkylation of Cys797) using QM/MM simulations and the recognition step by MD simulations coupled with free-energy calculations. The calculations indicate that L718Q has a negligible impact on both the activation energy for Cys797 alkylation and the free-energy of binding for the formation of the non-covalent complex. The results show that Gln718 affects the conformational space of the EGFR–osimertinib complex, stabilizing a conformation of acrylamide which prevents reaction with Cys797.

Multiscale analysis of enantioselectivity in enzyme-catalysed 'lethal synthesis' using projector-based embedding

The action of fluoroacetate as a broad-spectrum mammalian pesticide depends on the 'lethal synthesis' of fluorocitrate by citrate synthase, through a subtle enantioselective enolization of fluoroacetyl-coenzyme A. In this work, we demonstrate how a projection-based embedding method can be applied to calculate coupled cluster (CCSD(T)) reaction profiles from quantum mechanics/molecular mechanics optimized pathways for this enzyme reaction. Comparison of pro-R and pro-S proton abstraction in citrate synthase at the CCSD(T)-in-DFT//MM level yields the correct enantioselectivity. We thus demonstrate the potential of projection-based embedding for determining stereoselectivity in enzymatic systems. We further show that the method is simple to apply, eliminates variability due to the choice of density functional theory functional and allows the efficient calculation of CCSD(T) quality enzyme reaction barriers.

Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms

The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.

Mechanistic Insights into the Reaction of Chlorination of Tryptophan Catalyzed by Tryptophan 7-Halogenase

Tryptophan 7-halogenase catalyzes chlorination of free tryptophan to 7-chlorotryptophan, which is the first step in the antibiotic pyrrolnitrin biosynthesis. Many biologically and pharmaceutically active natural products contain chlorine and thus, an understanding of the mechanism of its introduction into organic molecules is important. Whilst enzyme-catalyzed chlorination is accomplished with ease, it remains a difficult task for the chemists. Therefore, utilizing enzymes in the synthesis of chlorinated organic compounds is important, and providing atomistic mechanistic insights about the reaction mechanism of tryptophan 7-halogenase is vital and timely. In this work, we examined a mechanism for the reaction of tryptophan chlorination, performed by tryptophan 7-halogenase, by calculating potential energy and free energy surfaces using two different Combined Quantum Mechanical/Molecular Mechanical (QM/MM) methods both employing Density Functional Theory (DFT) for the QM region. Both computational strategies agree on the nature of the rate-limiting step and provided close results for the reaction barriers of the two reaction steps. The calculations for both the potential energy and the free energy profiles showed very similar geometric features and hydrogen bonding interactions for the characterized stationary points.

Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase

Creatininase catalyzes the conversion of creatinine (a biosensor for kidney function) to creatine via a two-step mechanism: water addition followed by ring opening. Water addition is common to other known cyclic amidohydrolases, but the precise mechanism for ring opening is still under debate. The proton donor in this step is either His178 or a water molecule bound to one of the metal ions, and the roles of His178 and Glu122 are unclear. Here, the two possible reaction pathways have been fully examined by means of combined quantum mechanics/molecular mechanics simulations at the SCC-DFTB/CHARMM22 level of theory. The results indicate that His178 is the main catalytic residue for the whole reaction and explain its role as proton shuttle during the ring-opening step. In the first step, His178 provides electrostatic stabilization to the gem-diolate tetrahedral intermediate. In the second step, His178 abstracts the hydroxyl proton of the intermediate and delivers it to the cyclic amide nitrogen, leading to ring opening. The latter is the rate-limiting step with a free energy barrier of 18.5 kcal/mol, in agreement with the experiment. We find that Glu122 must be protonated during the enzyme reaction, so that it can form a stable hydrogen bond with its neighboring water molecule. Simulations of the E122Q mutant showed that this replacement disrupts the H-bond network formed by three conserved residues (Glu34, Ser78, and Glu122) and water, increasing the energy barrier. Our computational studies provide a comprehensive explanation for previous structural and kinetic observations, including why the H178A mutation causes a complete loss of activity but the E122Q mutation does not.

Identification of the quinolinedione inhibitor binding site in Cdc25 phosphatase B through docking and molecular dynamics simulations

Cdc25 phosphatase B, a potential target for cancer therapy, is inhibited by a series of quinones. The binding site and mode of quinone inhibitors to Cdc25B remains unclear, whereas this information is important for structure-based drug design. We investigated the potential binding site of NSC663284 [DA3003-1 or 6-chloro-7-(2-morpholin-4-yl-ethylamino)-quinoline-5, 8-dione] through docking and molecular dynamics simulations. Of the two main binding sites suggested by docking, the molecular dynamics simulations only support one site for stable binding of the inhibitor. Binding sites in and near the Cdc25B catalytic site that have been suggested previously do not lead to stable binding in 50 ns molecular dynamics (MD) simulations. In contrast, a shallow pocket between the C-terminal helix and the catalytic site provides a favourable binding site that shows high stability. Two similar binding modes featuring protein-inhibitor interactions involving Tyr428, Arg482, Thr547 and Ser549 are identified by clustering analysis of all stable MD trajectories. The relatively flexible C-terminal region of Cdc25B contributes to inhibitor binding. The binding mode of NSC663284, identified through MD simulation, likely prevents the binding of protein substrates to Cdc25B. The present results provide useful information for the design of quinone inhibitors and their mechanism of inhibition.

Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128β. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

Structural Basis of Catalysis in the Bacterial Monoterpene Synthases Linalool Synthase and 1,8-Cineole Synthase

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.

Rapid Estimation of Catalytic Efficiency by Cumulative Atomic Multipole Moments: Application to Ketosteroid Isomerase Mutants

We propose a simple atomic multipole electrostatic model to rapidly evaluate the effects of mutation on enzyme activity and test its performance on wild-type and mutant ketosteroid isomerase. The predictions of our atomic multipole model are similar to those obtained with symmetry-adapted perturbation theory at a fraction of the computational cost. We further show that this approach is relatively insensitive to the precise amino acid side chain conformation in mutants and may thus be useful in computational enzyme (re)design.

Insights into the Mechanistic Basis of Plasmid-Mediated Colistin Resistance from Crystal Structures of the Catalytic Domain of MCR-1

The polymixin colistin is a "last line" antibiotic against extensively-resistant Gram-negative bacteria. Recently, the mcr-1 gene was identified as a plasmid-mediated resistance mechanism in human and animal Enterobacteriaceae, with a wide geographical distribution and many producer strains resistant to multiple other antibiotics. mcr-1 encodes a membrane-bound enzyme catalysing phosphoethanolamine transfer onto bacterial lipid A. Here we present crystal structures revealing the MCR-1 periplasmic, catalytic domain to be a zinc metalloprotein with an alkaline phosphatase/sulphatase fold containing three disulphide bonds. One structure captures a phosphorylated form representing the first intermediate in the transfer reaction. Mutation of residues implicated in zinc or phosphoethanolamine binding, or catalytic activity, restores colistin susceptibility of recombinant E. coli. Zinc deprivation reduces colistin MICs in MCR-1-producing laboratory, environmental, animal and human E. coli. Conversely, over-expression of the disulphide isomerase DsbA increases the colistin MIC of laboratory E. coli. Preliminary density functional theory calculations on cluster models suggest a single zinc ion may be sufficient to support phosphoethanolamine transfer. These data demonstrate the importance of zinc and disulphide bonds to MCR-1 activity, suggest that assays under zinc-limiting conditions represent a route to phenotypic identification of MCR-1 producing E. coli, and identify key features of the likely catalytic mechanism.

A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase

Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.

Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions

Large-scale conformational change is a common feature in the catalytic cycles of enzymes. Many enzymes function as homodimers with active sites that contain elements from both chains. Symmetric and anti-symmetric cooperative motions in homodimers can potentially lead to correlated active site opening and/or closure, likely to be important for ligand binding and release. Here, we examine such motions in two different domain-swapped homodimeric enzymes: the DcpS scavenger decapping enzyme and citrate synthase. We use and compare two types of all-atom simulations: conventional molecular dynamics simulations to identify physically meaningful conformational ensembles, and rapid geometric simulations of flexible motion, biased along normal mode directions, to identify relevant motions encoded in the protein structure. The results indicate that the opening/closure motions are intrinsic features of both unliganded enzymes. In DcpS, conformational change is dominated by an anti-symmetric cooperative motion, causing one active site to close as the other opens; however a symmetric motion is also significant. In CS, we identify that both symmetric (suggested by crystallography) and asymmetric motions are features of the protein structure, and as a result the behaviour in solution is largely non-cooperative. The agreement between two modelling approaches using very different levels of theory indicates that the behaviours are indeed intrinsic to the protein structures. Geometric simulations correctly identify and explore large amplitudes of motion, while molecular dynamics simulations indicate the ranges of motion that are energetically feasible. Together, the simulation approaches are able to reveal unexpected functionally relevant motions, and highlight differences between enzymes.

High-level QM/MM calculations support the concerted mechanism for Michael addition and covalent complex formation in thymidylate synthase

Thymidylate synthase (TS) is a promising cancer target, due to its crucial function in thymine synthesis. It performs the reductive methylation of 2'-deoxyuridine-5'-phosphate (dUMP) to thymidine-5'-phosphate (dTMP), using N-5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) as a cofactor. After the formation of the dUMP/mTHF/TS noncovalent complex, and subsequent conformational activation, this complex has been proposed to react via nucleophilic attack (Michael addition) by Cys146, followed by methylene-bridge formation to generate the ternary covalent intermediate. Herein, QM/MM (B3LYP-D/6-31+G(d)-CHARMM27) methods are used to model the formation of the ternary covalent intermediate. A two-dimensional potential energy surface reveals that the methylene-bridged intermediate is formed via a concerted mechanism, as indicated by a single transition state on the minimum energy pathway and the absence of a stable enolate intermediate. A range of different QM methods (B3LYP, MP2 and SCS-MP2, and different basis sets) are tested for the calculation of the activation energy barrier for the formation of the methylene-bridged intermediate. We test convergence of the QM/MM results with respect to size of the QM region. Inclusion of Arg166, which interacts with the nucleophilic thiolate, in the QM region is important for reliable results; the MM model apparently does not reproduce energies for distortion of the guanidinium side chain correctly. The spin component scaled-Møller-Plessett perturbation theory (SCS-MP2) approach was shown to be in best agreement (within 1.1 kcal/mol) while the results obtained with MP2 and B3LYP also yielded acceptable values (deviating by less than 3 kcal/mol) compared with the barrier derived from experiment. Our results indicate that using a dispersion-corrected DFT method, or a QM method with an accurate treatment of electron correlation, increases the agreement between the calculated and experimental activation energy barriers, compared with the semiempirical AM1 method. These calculations provide important insight into the reaction mechanism of TS and may be useful in the design of new TS inhibitors.

Combined quantum mechanics/molecular mechanics (QM/MM) simulations for protein-ligand complexes: free energies of binding of water molecules in influenza neuraminidase

The applicability of combined quantum mechanics/molecular mechanics (QM/MM) methods for the calculation of absolute binding free energies of conserved water molecules in protein/ligand complexes is demonstrated. Here, we apply QM/MM Monte Carlo simulations to investigate binding of water molecules to influenza neuraminidase. We investigate five different complexes, including those with the drugs oseltamivir and peramivir. We investigate water molecules in two different environments, one more hydrophobic and one hydrophilic. We calculate the free-energy change for perturbation of a QM to MM representation of the bound water molecule. The calculations are performed at the BLYP/aVDZ (QM) and TIP4P (MM) levels of theory, which we have previously demonstrated to be consistent with one another for QM/MM modeling. The results show that the QM to MM perturbation is significant in both environments (greater than 1 kcal mol(-1)) and larger in the more hydrophilic site. Comparison with the same perturbation in bulk water shows that this makes a contribution to binding. The results quantify how electronic polarization differences in different environments affect binding affinity and also demonstrate that extensive, converged QM/MM free-energy simulations, with good levels of QM theory, are now practical for protein/ligand complexes.