Efficient Computation of Free Energy Surfaces of Diels?Alder Reactions in Explicit Solvent at Ab Initio QM/MM Level.
Ontology highlight
ABSTRACT:
SUBMITTER: Li P
PROVIDER: S-EPMC6222833 | biostudies-other | 2018 Sep
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We performed a computational investigation of the mechanism by which cyclodextrins (CDs) catalyze Diels-Alder reactions between 9-anthracenemethanol and N-cyclohexylmaleimide. Hydrogen bonds (Hbonds) between N-cyclohexylmaleimide and the hydroxyl groups of cyclodextrins were suggested to play an important role in this catalytic process. However, our free energy calculations and molecular dynamics simulations showed that these Hbonds are not stable, and quantum mechanical calculations suggested that the reaction is not promoted by these Hbonds. The binding of 9-anthracenemethanol and N-cyclohexylmaleimide to cyclodextrins was the key to the catalytic process. Cyclodextrins act as a container to hold the two reactants in the cavity, pre-organize them for the reactions, and thus reduce the entropy penalty to the activation free energy. Dimethyl-?-CD was a better catalyst for this specific reaction than ?-CD because of its stronger van der Waals interaction with the pre-organized reactants and its better performance in reducing the activation energy. This computational work sheds light on the mechanism of the catalytic reaction by cyclodextrins and introduces new perspectives of supramolecular catalysis.
| S-EPMC6011829 | biostudies-literature
Project description:Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).
| S-EPMC3797183 | biostudies-literature
Project description:Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the covalent inhibition mechanism between AChE and the nerve toxin soman and determined its free energy profile for the first time. Our results indicate that phosphonylation of the catalytic serine by soman employs an addition-elimination mechanism, which is highly associative and stepwise: in the initial addition step, which is also rate-limiting, His440 acts as a general base to facilitate the nucleophilic attack of Ser200 on the soman's phosphorus atom to form a trigonal bipyrimidal pentacovalent intermediate; in the subsequent elimination step, Try121 of the catalytic gorge stabilizes the leaving fluorine atom prior to its dissociation from the active site. Together with our previous characterization of the aging mechanism of soman inhibited AChE, our simulations have revealed detailed molecular mechanistic insights into the damaging function of the nerve agent soman.
| S-EPMC4183371 | biostudies-literature
Project description:In QM/MM calculations, it is essential to handle electrostatic interactions between the QM and MM subsystems accurately and efficiently. To achieve maximal efficiency, it is convenient to adopt a hybrid scheme, where the QM electron density is used explicitly in the evaluation of short-range QM/MM electrostatic interactions, while a multipolar representation for the QM electron density is employed to account for the long-range QM/MM electrostatic interactions. In order to avoid energy discontinuity at the cutoffs, which separate the short- and long-range QM/MM electrostatic interactions, a switching function should be utilized to ensure a smooth potential energy surface. In this study, we benchmarked the accuracy of such hybrid embedding schemes for QM/MM electrostatic interactions using different multipolar representations, switching functions and cutoff distances. For test systems (neutral and anionic oxyluciferin in MM (aqueous and enzyme) environments), the best accuracy was acquired with a combination of QM electrostatic potential (ESP) charges and dipoles and two switching functions (long-range electrostatic corrections (LREC) and Switch) in the treatment of long-range QM/MM electrostatics. It allowed us to apply a 10Å distance cutoff and still obtain QM/MM electrostatics/polarization energies within 0.1 kcal/mol and time-dependent density functional theory (TDDFT)/MM vertical excitation energies within 10-3 eV from theoretical reference values.
| S-EPMC6222767 | biostudies-literature
Project description:Although quantum mechanical/molecular mechanics (QM/MM) methods are now routinely applied to the studies of chemical reactions in condensed phases and enzymatic reactions, they may experience technical difficulties when the reactive region is varying over time. For instance, when the solvent molecules are directly participating in the reaction, the exchange of water molecules between the QM and MM regions may occur on a time scale comparable to the reaction time. To cope with this situation, several adaptive QM/MM schemes have been proposed. However, these methods either add significantly to the computational cost or introduce artificial restraints to the system. In this work, we developed a novel adaptive QM/MM scheme and applied it to the study of a nucleophilic addition reaction. In this scheme, the configuration sampling was performed with a small QM region (without solvent molecules), and the thermodynamic properties under another potential energy function with a larger QM region (with a certain number of solvent molecules and/or different levels of QM theory) are computed via extrapolation using the reference-potential method. Our simulation results show that this adaptive QM/MM scheme is numerically stable, at least for the case studied in this work. Furthermore, this method also offers an inexpensive way to examine the convergence of the QM/MM calculation with respect to the size of the QM region.
| S-EPMC8335528 | biostudies-literature
Project description:Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1-2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.
| S-EPMC2919153 | biostudies-literature
Project description:AlkB is the title enzyme of a family of DNA dealkylases that catalyze the direct oxidative dealkylation of nucleobases. The conventional mechanism for the dealkylation of N(1)-methyl adenine (1-meA) catalyzed by AlkB after the formation of Fe(IV)-oxo is comprised by a reorientation of the oxo moiety, hydrogen abstraction, OH rebound from the Fe atom to the methyl adduct, and the dissociation of the resulting methoxide to obtain the repaired adenine base and formaldehyde. An alternative pathway with hydroxide as a ligand bound to the iron atom is proposed and investigated by QM/MM simulations. The results show OH(-) has a small impact on the barriers for the hydrogen abstraction and OH rebound steps. The effects of the enzyme and the OH(-) ligand on the hydrogen abstraction by the Fe(IV)-oxo moiety are discussed in detail. The new OH rebound step is coupled with a proton transfer to the OH(-) ligand and results in a novel zwitterion intermediate. This zwitterion structure can also be characterized as Fe-O-C complex and facilitates the formation of formaldehyde. In contrast, for the pathway with H2O bound to iron, the hydroxyl product of the OH rebound step first needs to unbind from the metal center before transferring a proton to Glu136 or other residue/substrate. The consistency between our theoretical results and experimental findings is discussed. This study provides new insights into the oxidative repair mechanism of DNA repair by nonheme Fe(II) and ?-ketoglutarate (?-KG) dependent dioxygenases and a possible explanation for the substrate preference of AlkB.
| S-EPMC4230374 | biostudies-literature
Project description:Ten-eleven translocation 2 (TET2) is an Fe/?-ketoglutarate (?-KG) dependent enzyme that dealkylates 5-methylcytosine (5mC). The reaction mechanism involves a series of three sequential oxidations that convert 5mC to 5-hydroxy-methylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). Our previous biochemical and computational studies uncovered an active site scaffold that is required for wild-type (WT) stepwise oxidation (Nat. Chem. Bio., 13, 181). We showed that the mutation of a single residue, T1372 to some amino acids, such as Glu, can impact the iterative oxidation steps and stop the oxidation of 5hmC to 5fC/caC. However, the source of the stalling at the first oxidation step by some mutant TET proteins still remains unclear. Here, we studied the catalytic mechanism of oxidation of 5hmC to 5fC by WT and T1372E TET2 using an ab initio quantum mechanical/molecular mechanical (QM/MM) approach. Our results suggest that the rate limiting step for WT TET2 involves a hydrogen atom abstraction from the hydroxyl group of 5hmC by the ferryl moiety in the WT. By contrast, our calculations for the T1372E mutant indicate that the rate limiting step for this variant corresponds to a second proton abstraction and the calculated barrier is almost twice as large as for WT TET2. Our results suggest that the large barrier for the 5hmC to 5fC oxidation in this mutant is due (at least in part) to the unfavorable orientation of the substrate in the active site. Combined electron localization function (ELF) and non-covalent interaction (NCI) analyses provide a qualitative description of the evolution of the electronic structure of the active site along the reaction path. Energy decomposition analysis (EDA) has been performed on the WT to investigate the impact of each MM residue on catalytic activity.
| S-EPMC6244454 | biostudies-literature
Project description:Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.
| S-EPMC3124314 | biostudies-literature
Project description:We describe a new interface of the GPU parallelized TeraChem electronic structure package and the Amber molecular dynamics package for quantum mechanical (QM) and mixed QM and molecular mechanical (MM) molecular dynamics simulations. This QM/MM interface is used for computation of the absorption spectra of the photoactive yellow protein (PYP) chromophore in vacuum, aqueous solution, and protein environments. The computed excitation energies of PYP require a very large QM region (hundreds of atoms) covalently bonded to the chromophore in order to achieve agreement with calculations that treat the entire protein quantum mechanically. We also show that 40 or more surrounding water molecules must be included in the QM region in order to obtain converged excitation energies of the solvated PYP chromophore. These results indicate that large QM regions (with hundreds of atoms) are a necessity in QM/MM calculations.
| S-EPMC3590007 | biostudies-literature