Mechanism of proteolysis in matrix metalloproteinase‐2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase‐2 (MMP‐2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all‐atom three‐dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace‐Gln‐Gly～Ile‐Ala‐Gly‐Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four‐step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and C? N bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate‐limiting step in MMP‐2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations. © 2015 Wiley Periodicals, Inc.     

Synthesis of Four Orthogonally Protected Rare l-Hexose Thioglycosides from d-Mannose by C-5 and C-4 Epimerization

l-Hexoses are important components of biologically relevant compounds and precursors of some therapeuticals. However, they typically cannot be obtained from natural sources and due to the complexity of their synthesis, their commercially available derivatives are also very expensive. Starting from one of the cheapest d-hexoses, d-mannose, using inexpensive and readily available chemicals, we developed a reaction pathway to obtain two orthogonally protected l-hexose thioglycoside derivatives, l-gulose and l-galactose, through the corresponding 5,6-unsaturated thioglycosides by C-5 epimerization. From these derivatives, the orthogonally protected thioglycosides of further two l-hexoses (l-allose and l-glucose) were synthesized by C-4 epimerization. The preparation of the key intermediates, the 5,6-unsaturated derivatives, was systematically studied using various protecting groups. By the method developed, we are able to produce highly functionalized l-gulose derivatives in 9 steps (total yields: 21–23%) and l-galactose derivatives in 12 steps (total yields: 6–8%) starting from d-mannose.     

A critical evaluation of different QM/MM frontier treatments with SCC-DFTB as the QM method

The performance of different link atom based frontier treatments in QM/MM simulations was evaluated critically with SCC-DFTB as the QM method. In addition to the analysis of gas-phase molecules as in previous studies, an important element of the present work is that chemical reactions in realistic enzyme systems were also examined. The schemes tested include all options available in the program CHARMM for SCC-DFTB/MM simulation, which treat electrostatic interactions due to the MM atoms close to the QM/MM boundary in different ways. In addition, a new approach, the divided frontier charge (DIV), has been implemented in which the partial charge associated with the frontier MM atom ("link host") is evenly distributed to the other MM atoms in the same group. The performance of these schemes was evaluated based on properties including proton affinities, deprotonation energies, dipole moments, and energetics of proton transfer reactions. Similar to previous work, it was found that calculated proton affinities and deprotonation energies of alcohols, carbonic acids, amino acids, and model DNA bases are very sensitive to the link atom scheme; the commonly used single link atom approach often gives error on the order of 15 to 20 kcal/mol. Other schemes give better and, on average, mutually comparable results. For proton transfer reactions, encouragingly, both activation barriers and reaction energies are fairly insensitive (within a typical range of 2-4 kcal/mol) to the link atom scheme due to error cancellation, and this was observed for both gas-phase and enzyme systems. Therefore, the effect of using different link atom schemes in QM/MM simulations is rather small for chemical reactions that conserve the total charge. Although the current study used an approximate DFT method as the QM level, the observed trends are expected to be applicable to QM/MM methods with use of other QM approaches. This observation does not mean to encourage QM/MM simulations without careful benchmark in the study of specific systems, rather it emphasizes that other technical details, such as the treatment of long-range electrostatics, tend to play a more important role and need to be handled carefully.     

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1–2 kcal ⋅ mol⁻¹. Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol⁻¹ and −6.3 to +3.8 kcal ⋅ mol⁻¹ depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.     

Mechanism of chemical reactions in the active site of aspartate N-acetyltransferase NAT8L revealed by molecular modeling

The results of a computational study of the synthesis of a key brain metabolite, N-acetyl-l-aspartate, catalyzed by aspartate N-acetyltransferase, encoded by the NAT8L gene, are reported. The reaction Gibbs energy profiles were computed using molecular dynamics simulations with interaction potentials estimated on-the-fly by the quantum mechanics/molecular mechanics QM(PBE0/6-31G**)/MM(CHARMM) approach. The revealed reaction mechanism includes four elementary steps with corresponding activation energies not exceeding 14 kcal mol^?1     

Electronic structure of compound I in human isoforms of cytochrome P450 from QM/MM modeling

Human cytochromes P450 play a vital role in drug metabolism. The key step in substrate oxidation involves hydrogen atom abstraction or C=C bond addition by the oxygen atom of the Compound I intermediate. The latter has three unpaired electrons, two on the Fe-O center and one shared between the porphyrin ring and the proximal cysteinyl sulfur atom. Changes in its electronic structure have been suggested to affect reactivity. The electronic and geometric structure of Compound I in three important human subfamilies of cytochrome P450 (P450, 2C, 2B, and 3A) that are major contributors to drug metabolism is characterized here using combined quantum mechanical/molecular mechanical (QM/MM) calculations at the B3LYP:CHARMM27 level. Compound I is remarkably similar in all isoforms, with the third unpaired electron located mainly on the porphyrin ring, and this prediction is not very sensitive to details of the QM/MM methodology, such as the DFT functional, the basis set, or the size of the QM region. The presence of substrate also has no effect. The main source of variability in spin density on the cysteinyl sulfur (from 26 to 50%) is the details of the system setup, such as the starting protein geometry used for QM/MM minimization. This conformational effect is larger than the differences between human isoforms, which are therefore not distinguishable on electronic grounds, so it is unlikely that observed large differences in substrate selectivity can be explained to a large extent in these terms.     

QM/MM simulations predict a covalent intermediate in the hen egg white lysozyme reaction with its natural substrate

Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations indicate that the reaction of native HEWL with its natural substrate involves a covalent intermediate, in contrast to the 'textbook' mechanism for this seminal enzyme.     

Methodological aspects of QM/MM calculations: A case study on matrix metalloproteinase‐2

We address methodological issues in quantum mechanics/molecular mechanics (QM/MM) calculations on a zinc‐dependent enzyme. We focus on the first stage of peptide bond cleavage by matrix metalloproteinase‐2 (MMP‐2), that is, the nucleophilic attack of the zinc‐coordinating water molecule on the carbonyl carbon atom of the scissile fragment of the substrate. This step is accompanied by significant charge redistribution around the zinc cation, bond cleavage, and bond formation. We vary the size and initial geometry of the model system as well as the computational protocol to demonstrate the influence of these choices on the results obtained. We present QM/MM potential energy profiles for a set of snapshots randomly selected from QM/MM‐based molecular dynamics simulations and analyze the differences in the computed profiles in structural terms. Since the substrate in MMP‐2 is located on the protein surface, we investigate the influence of the thickness of the water layer around the enzyme on the QM/MM energy profile. Thin water layers (0–2 Å) give unrealistic results because of structural reorganizations in the active‐site region at the protein surface. A 12 Å water layer appears to be sufficient to capture the effect of the solvent; the corresponding QM/MM energy profile is very close to that obtained from QM/MM/SMBP calculations using the solvent macromolecular boundary potential (SMBP). We apply the optimized computational protocol to explain the origin of the different catalytic activity of the Glu116Asp mutant: the energy barrier for the first step is higher, which is rationalized on structural grounds. © 2016 Wiley Periodicals, Inc.     

Total synthesis and biochemical characterization of mirror image barnase

In this study we synthesized and characterized mirror image barnase (B. amyloliquefaciens ribonuclease). d-Barnase was identical to l-barnase, when analyzed by liquid chromatography and mass-spectrometry. Proteolysis of the mirror image enzyme revealed that in contrast to its native counterpart, d-barnase was completely stable to digestive proteases. In enzymatic assays, d-barnase had the reciprocal chiral specificity and was fully active towards mirror image substrates. Interestingly, d-barnase also hydrolyzed the substrate of the native chirality, albeit 4000 times less efficiently. This effect was further confirmed by digesting a native 112-mer RNA with the enzyme. Additional studies revealed that barnase accommodates a range of substrates with various chiralities, but the prime requirement for guanosine remains. These studies point toward using mirror image enzymes as modern agents in biotechnology.     

An extensible interface for QM/MM molecular dynamics simulations with AMBER

We present an extensible interface between the AMBER molecular dynamics (MD) software package and electronic structure software packages for quantum mechanical (QM) and mixed QM and classical molecular mechanical (MM) MD simulations within both mechanical and electronic embedding schemes. With this interface, ab initio wave function theory and density functional theory methods, as available in the supported electronic structure software packages, become available for QM/MM MD simulations with AMBER. The interface has been written in a modular fashion that allows straight forward extensions to support additional QM software packages and can easily be ported to other MD software. Data exchange between the MD and QM software is implemented by means of files and system calls or the message passing interface standard. Based on extensive tests, default settings for the supported QM packages are provided such that energy is conserved for typical QM/MM MD simulations in the microcanonical ensemble. Results for the free energy of binding of calcium ions to aspartate in aqueous solution comparing semiempirical and density functional Hamiltonians are shown to demonstrate features of this interface. © 2013 Wiley Periodicals, Inc.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Molecular dynamics and free energy perturbation study of hydride-ion transfer step in dihydrofolate reductase using combined quantum and molecular mechanical model

We used molecular dynamics simulation and free energy perturbation (FEP) methods to investigate the hydride-ion transfer step in the mechanism for the nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reduction of a novel substrate by the enzyme dihydrofolate reductase (DHFR). The system is represented by a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1 semiempirical molecular orbital method for the reacting substrate and NADPH cofactor fragments, the AMBER force field for DHFR, and the TIP3P model for solvent water. The FEP calculations were performed for a number of choices for the QM system. The substrate, 8-methylpterin, was treated quantum mechanically in all the calculations, while the larger cofactor molecule was partitioned into various QM and MM regions with the addition of “link” atoms (F, CH₃, and H). Calculations were also carried out with the entire NADPH molecule treated by QM. The free energies of reaction and the net charges on the NADPH fragments were used to determine the most appropriate QM/MM model. The hydride-ion transfer was also carried out over several FEP pathways, and the QM and QM/MM component free energies thus calculated were found to be state functions (i.e., independent of pathway). A ca. 10 kcal/mol increase in free energy for the hydride-ion transfer with an activation barrier of ca. 30 kcal/mol was calculated. The increase in free energy on the hydride-ion transfer arose largely from the QM/MM component. Analysis of the QM/MM energy components suggests that, although a number of charged residues may contribute to the free energy change through long-range electrostatic interactions, the only interaction that can account for the 10 kcal/mol increase in free energy is the hydrogen bond between the carboxylate side chain of Glu30 (avian DHFR) and the activated (protonated) substrate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 977–988, 1998     

A crystal-structural study of Pauling–Corey rippled sheets

Following the seminal theoretical work on the pleated β-sheet published by Pauling and Corey in 1951, the rippled β-sheet was hypothesized by the same authors in 1953. In the pleated β-sheet the interacting β-strands have the same chirality, whereas in the rippled β-sheet the interacting β-strands are mirror-images. Unlike with the pleated β-sheet that is now common textbook knowledge, the rippled β-sheet has been much slower to evolve. Much of the experimental work on rippled sheets came from groups that study aggregating racemic peptide systems over the course of the past decade. This includes MAX1/DMAX hydrogels (Schneider), L/D-KFE8 aggregating systems (Nilsson), and racemic Amyloid β mixtures (Raskatov). Whether a racemic peptide mixture is “ripple-genic” (i.e., whether it forms a rippled sheet) or “pleat-genic” (i.e., whether it forms a pleated sheet) is likely governed by a complex interplay of thermodynamic and kinetic effects. Structural insights into rippled sheets remain limited to only a very few studies that combined sparse experimental structural constraints with molecular modeling. Crystal structures of rippled sheets are needed so we can rationally design rippled sheet architectures. Here we report a high-resolution crystal structure, in which (l,l,l)-triphenylalanine and (d,d,d)-triphenylalanine form dimeric antiparallel rippled sheets, which pack into herringbone layer structures. The arrangements of the tripeptides and their mirror-images in the individual dimers were in excellent agreement with the theoretical predictions by Pauling and Corey. A subsequent mining of the PDB identified three orphaned rippled sheets among racemic protein crystal structures.

Following the seminal theoretical work on the pleated β-sheet published by Pauling and Corey in 1951, the rippled β-sheet was hypothesized by the same authors in 1953.     

Hybrid quantum mechanics/molecular mechanics-based molecular dynamics simulation of acid-catalyzed dehydration of polyols in liquid water

We use the conversion of protonated glycerol to acrolein for a case study of the mechanism of acid-catalyzed dehydration of polyols in aqueous environments. We employ hybrid Quamtum Mechanics/Molecular Mechanics Molecular Dynamics (QM/MM MD) simulations with biased sampling and perform free energy calculations for the elementary steps of the reaction. We investigate the effects of solvent dynamics and in particular the role of quantum mechanical water in the dehydration mechanism. We present results supporting a mechanism that proceeds via water-mediated proton transfers and thus through an enol intermediate. We find that the first dehydration may take place by two, low-energy pathways requiring, respectively, 20.9 and 18.8 kcal/mol of activation free energy. The second dehydration requires 19.9 kcal/mol of activation free energy while for the overall reaction we compute a free energy change of -8 kcal/mol.     

Transition state structure invariance to model system size and calculation levels: a QM/MM study of the carboxylation step catalyzed by Rubisco

The present study elucidates structural features related to the molecular mechanism in the carboxylation step of the reaction catalyzed by Rubisco. Starting from the initial X-ray Protein Data Bank structure of a Rubisco monomer, the reactive subsystem in vacuo is subjected to quantum chemical semiempirical and ab initio studies, while the effects of the protein environments are included by means of a hybrid quantum mechanical/molecular mechanical (QM/MM) approach. The QM/MM is used to characterize the transition structure for carboxylation inside the protein. The calculations were made with the AM1/CHARMM/GRACE scheme. Comparisons between the in vacuo and in situ transition structures show remarkable invariance with respect to geometric parameters, index and transition vector amplitudes. The transition state couples the carbon dioxide attack to the C2 center of the substrate in its dienol form with a simultaneous intramolecular hydrogen transfer from the C2 atom to the hydroxyl group linked to the C3 center. This study suggests that carboxylation may be simultaneously coupled to the activation of the C3 center in the enzyme. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

A quantum mechanical/molecular mechanical study on the catalysis of the pyridoxal 5'-phosphate-dependent enzyme L-serine dehydratase

The catalytic mechanism of a pyridoxal 5'-phosphate-dependent enzyme, l-serine dehydratase, has been investigated using ab initio quantum mechanical/molecular mechanical (QM/MM) methods. New insights into the chemical steps have been obtained, including the chemical role of the substrate carboxyl group in the Schiff base formation step and a proton-relaying mechanism involving the phosphate of the cofactor in the beta-hydroxyl-leaving step. The latter step is of no barrier and follows sequentially after the elimination of the alpha-proton, leading to a single but sequential alpha, beta-elimination step. The rate-limiting transition state is specifically stabilized by the enzyme environment. At this transition state, charges are localized on the substrate carboxyl group, as well as on the amino group of Lys41. Specific interactions of the enzyme environment with these groups are able to lower the activation barrier significantly. One major difficulty associated with studies of complicated enzymatic reactions using ab initio QM/MM models is the appropriate choices of reaction coordinates. In this study, we have made use of efficient semiempirical models and pathway optimization techniques to overcome this difficulty.     

‘Reverse biomimetic’ synthesis of l-arogenate and its stabilized analogues from l-tyrosine

l-Arogenate (also known as l-pretyrosine) is a primary metabolite on a branch of the shikimate biosynthetic pathway to aromatic amino acids. It plays a key role in the synthesis of plant secondary metabolites including alkaloids and the phenylpropanoids that are the key to carbon fixation. Yet understanding the control of arogenate metabolism has been hampered by its extreme instability and the lack of a versatile synthetic route to arogenate and its analogues. We now report a practical synthesis of l-arogenate in seven steps from O-benzyl l-tyrosine methyl ester in an overall yield of 20%. The synthetic route also delivers the fungal metabolite spiroarogenate, as well as a range of stable saturated and substituted analogues of arogenate. The key step in the synthesis is a carboxylative dearomatization by intramolecular electrophilic capture of tyrosine''s phenolic ring using an N-chloroformylimidazolidinone moiety, generating a versatile, functionalizable spirodienone intermediate.

l-Tyrosine provides a precursor for a practical synthesis of the unstable primary metabolite l-arogenate and some stabilised arogenate analogues.