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.     

QM/MM Studies on Cyclodextrin-Alcohol Interaction

The main objective of this study is to provide an insight into the interactions involved during adsorption of the alcohols on β-CD composite nanostructured membrane. Interactions between β-cyclodextrin (β-CD) and alcohols (methanol, ethanol and butanol) are studied using the QM/MM method. Magnitude of interaction energies show that the alcohols are adsorbed on the membrane. In addition, the thermochemical analysis suggests that the formation of these host-guest complexes is enthalpy driven.     

Interfacing Q-Chem and CHARMM to perform QM/MM reaction path calculations

A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with Hartree-Fock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O...H2O, H2O...Na+, and H2O...Cl- complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is -23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.     

Spectral Features of Canthaxanthin in HCP2. A QM/MM Approach

The increased interest in sequencing cyanobacterial genomes has allowed the identification of new homologs to both the N-terminal domain (NTD) and C-terminal domain (CTD) of the Orange Carotenoid Protein (OCP). The N-terminal domain homologs are known as Helical Carotenoid Proteins (HCPs). Although some of these paralogs have been reported to act as singlet oxygen quenchers, their distinct functional roles remain unclear. One of these paralogs (HCP2) exclusively binds canthaxanthin (CAN) and its crystal structure has been recently characterized. Its absorption spectrum is significantly red-shifted, in comparison to the protein in solution, due to a dimerization where the two carotenoids are closely placed, favoring an electronic coupling interaction. Both the crystal and solution spectra are red-shifted by more than 50 nm when compared to canthaxanthin in solution. Using molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) studies of HCP2, we aim to simulate these shifts as well as obtain insight into the environmental and coupling effects of carotenoid–protein interactions.     

Infrared spectroscopy of N-methylacetamide in water from high-level QM/MM calculations简

The infrared(IR) spectra of the N-methylacetamide molecule in water are calculated by using the MD simulation with high-level QM/MM corrections. The B3LYP and MP2 levels with 6-311++G** basis set are used for the QM region, respectively. Our results show all IR spectra at the B3LYP level are well consistent with the corresponding MP2 results. A dynamical charge fluctuation is observed for each atom along the simulation trajectories due to the electrostatic polarization(EP) effects from surrounding solvent environment. We find that the QM/MM corrected IR spectra satisfactorily reproduce the experimental vibrational features of amide I–III modes.     

Pseudobond ab initio QM/MM approach and its applications to enzyme reactions

This perspective article mainly focuses on the development and applications of a pseudobond ab initio QM/MM approach to study enzyme reactions. The following aspects of methodology development are discussed: the approaches for the QM/MM covalent boundary problem, an efficient iterative optimization procedure, the methods to determine enzyme reaction paths, and the approaches to calculate free energy change in enzyme reactions. Several applications are described to illustrate the capability of the methods. Finally, future directions are discussed.     

QM/MM study of the taxadiene synthase mechanism

Combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the reaction mechanism of taxadiene synthase (TXS). TXS catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to taxadiene (T) and four minor cyclic products. All these products originate from the deprotonation of carbocation intermediates. The reaction profiles for the conversion of GGPP to T as well as to minor products were calculated for different configurations of relevant TXS carbocation complexes. The QM region was treated at the M06-2X/TZVP level, while the CHARMM27 force field was used to describe the MM region. The QM/MM calculations suggest a reaction pathway for the conversion of GGPP to T, which slightly differs from previous proposals regarding the number of reaction steps and the conformation of the carbocations. The QM/MM results also indicate that the formation of minor products via water-assisted deprotonation of the carbocations is highly exothermic, by about −7 to −23 kcal/mol. Curiously, however, the computed barriers and reaction energies indicate that the formation of some of the minor products is more facile than the formation of T. Thus, the present QM/MM calculations provide detailed insights into possible reaction pathways and into the origin of the promiscuity of TXS, but they do not reproduce the product distribution observed experimentally. © 2019 Wiley Periodicals, Inc.     

Producing DFT/MM enzyme reaction trajectories from SCC‐DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction

The SCC‐DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC‐DFTB/MIO/CHARMM reaction surface to drive B3LYP/6‐31G/MM and B3LYP/6‐31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6‐31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6‐31G/MM‐computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi‐empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.     

QM/MM method for metal–organic interfaces

We present a QM/MM method for modeling metal/organic interfaces, which incorporates contributions from long‐range electron correlation, characteristic to metals and non‐bonded interactions in organic systems. This method can be used to study structurally irregular systems. We apply the method to model finite size domains of self‐assembled monolayers on the gold (111) surface and discuss the influence of boundary effects on the electrostatic and electronic properties of these systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

应用QM与ABEEM/MM结合的方法研究NAMI-A的结构性质

应用量子力学(QM)与ABEEM浮动电荷力场(ABEEM/MM)相结合的方法研究了抗癌药物NAMI-A在水溶液中的结构性质. 所有的结构优化都是在DFT的B3LYP方法下采用6-31G(d,p)和LanL2DZ基组完成的, 没有加入任何限制性条件. 结果表明, 优化得到的NAMI-A构型受不同环境及方法的影响均有变化. 与气相中得到的构型相比, QM/MM迭代优化得到构型要比PCM的构型变化更明显. QM/MM (ABEEM/MM)迭代优化得到的NAMI-A构型比QM/MM (OPLS-AA)的变化要小. 总之, 溶剂通过极化效应对NAMI-A结构、电荷分布及径向分布函数等性质均有影响, 客观地处理极化效应才能正确地反映QM区的性质.     

Molecular insight into the interaction mechanisms of amino‐2H‐imidazole derivatives with BACE1 protease: A QM/MM and QTAIM study

In this study, we described quantitatively the interactions between two new amino‐2H‐imidazole inhibitors ((R)‐1t and (S)‐1m) and BACE1 using a hybrid quantum mechanics‐molecular mechanical (QM/MM) method together with a quantum theory of atoms In Molecules (QTAIM) analysis. Our computational calculations revealed that the binding affinity of these compounds is mostly related to the amino‐2H‐imidazole core, which interact tightly with the aspartate dyad of the active site. The interactions were stronger when the inhibitors presented a bulky substituent with a hydrogen bond acceptor motif pointing toward Trp76, such as the 3,5‐dimethyl‐4‐methoxyphenyl group of compound (S)‐1m. Furthermore, the QTAIM analysis revealed that many hydrophobic interactions complement cooperatively the hydrogen bond which is not present when compound (R)‐1t is bound to the enzyme. The combined QM/MM‐QTAIM analysis allows identifying the interactions that account for the activity difference between compounds, even at a nanomolar range.     

Flexible effective fragment QM/MM method: validation through the challenging tests

A new version of the QM/MM method, which is based on the effective fragment potential (EFP) methodology [Gordon, M. et al., J Phys Chem A 2001, 105, 293] but allows flexible fragments, is verified through calculations of model molecular systems suggested by different authors as challenging tests for QM/MM approaches. For each example, the results of QM/MM calculations for a partitioned system are compared to the results of an all-electron ab initio quantum chemical study of the entire system. In each case we were able to achieve approximately similar or better accuracy of the QM/MM results compared to those described in original publications. In all calculations we kept the same set of parameters of our QM/MM scheme. A new test example is considered when calculating the potential of internal rotation in the histidine dipeptide around the C(alpha)bond;C(beta) side chain bond.     

MM and QM/MM Modeling of Threonyl-tRNA Synthetase: Model Testing and Simulations

Aminoacyl-tRNA synthetases are centrally important enzymes in protein synthesis. We have investigated threonyl-tRNA synthetase from E. coli, complexed with reactants, using molecular mechanics and combined quantum mechanical/molecular mechanical (QM/MM) techniques. These modeling methods have the potential to provide molecular level understanding of enzyme catalytic processes. Modeling of this enzyme presents a number of challenges. The procedure of system preparation and testing is described in detail. For example, the number of metal ions at the active site, and their positions, were investigated. Molecular dynamics simulations suggest that the system is most stable when it contains only one magnesium ion, and the zinc ion is removed. Two different QM/MM methods were tested in models based on the findings of MM molecular dynamics simulations. AM1/CHARMM calculations resulted in unrealistic structures for the phosphates in this system. This is apparently due to an error of AM1. PM3/CHARMM calculations proved to be more suitable for this enzyme system. These results will provide a useful basis for future modeling investigations of the enzyme mechanism and dynamics.     

QM/MM study of aqueous solvation of the uranyl fluoride [UO2F4(2-)] complex

The aqueous solvation of the uranylfluoride complex [UO(2)F(4) (2-)] was studied using full quantum mechanical (QM) and hybrid QM/molecular mechanics (MM) methods. Inclusion of a complete first solvation shell was found necessary to reproduce the experimentally observed heptacoordination of uranium. An efficient and accurate computational model is proposed that consists of structure optimization of the coordinated uranium complex as QM region, followed by single-point full QM calculations to compute relative energies. This method is proven feasible for studies of large solvated actinide complexes.