共查询到20条相似文献,搜索用时 15 毫秒
1.
Letif Mones Andrew Jones Andreas W. Götz Teodoro Laino Ross C. Walker Ben Leimkuhler Gábor Csányi Noam Bernstein 《Journal of computational chemistry》2015,36(9):633-648
The implementation and validation of the adaptive buffered force (AdBF) quantum‐mechanics/molecular‐mechanics (QM/MM) method in two popular packages, CP2K and AMBER are presented. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM‐MM interface errors by discarding forces near the boundary according to the buffered force‐mixing approach. New adaptive thermostats, needed by force‐mixing methods, are also implemented. Different variants of the method are benchmarked by simulating the structure of bulk water, water autoprotolysis in the presence of zinc and dimethyl‐phosphate hydrolysis using various semiempirical Hamiltonians and density functional theory as the QM model. It is shown that with suitable parameters, based on force convergence tests, the AdBF QM/MM scheme can provide an accurate approximation of the structure in the dynamical QM region matching the corresponding fully QM simulations, as well as reproducing the correct energetics in all cases. Adaptive unbuffered force‐mixing and adaptive conventional QM/MM methods also provide reasonable results for some systems, but are more likely to suffer from instabilities and inaccuracies. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc. 相似文献
2.
Dr. Wesley F. Zandberg Ismael Compañón Yun Kong Dr. Francisco Corzana Prof. B. Mario Pinto Prof. Henrik Clausen Prof. Jesús M. Peregrina Prof. David J. Vocadlo Prof. Carme Rovira Dr. Ramon Hurtado‐Guerrero 《Angewandte Chemie (International ed. in English)》2014,53(31):8206-8210
The retaining glycosyltransferase GalNAc‐T2 is a member of a large family of human polypeptide GalNAc‐transferases that is responsible for the post‐translational modification of many cell‐surface proteins. By the use of combined structural and computational approaches, we provide the first set of structural snapshots of the enzyme during the catalytic cycle and combine these with quantum‐mechanics/molecular‐mechanics (QM/MM) metadynamics to unravel the catalytic mechanism of this retaining enzyme at the atomic‐electronic level of detail. Our study provides a detailed structural rationale for an ordered bi–bi kinetic mechanism and reveals critical aspects of substrate recognition, which dictate the specificity for acceptor Thr versus Ser residues and enforce a front‐face SNi‐type reaction in which the substrate N‐acetyl sugar substituent coordinates efficient glycosyl transfer. 相似文献
3.
Dr. João T. S. Coimbra Dr. Rui P. P. Neves Prof. Dr. Ana V. Cunha Prof. Dr. Maria J. Ramos Prof. Dr. Pedro A. Fernandes 《Chemistry (Weinheim an der Bergstrasse, Germany)》2022,28(42):e202201066
The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol−1. Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway. 相似文献
4.
《Journal of computational chemistry》2017,38(6):319-335
The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc. 相似文献
5.
6.
Martí S Andrés J Moliner V Silla E Tuñón I Bertrán J 《Chemistry (Weinheim an der Bergstrasse, Germany)》2008,14(2):596-602
The Diels-Alder reaction is one of the most important and versatile transformations available to organic chemists for the construction of complex natural products, therapeutics agents, and synthetic materials. Given the lack of efficient enzymes capable of catalyzing this kind of reaction, it is of interest to ask whether a biological catalyst could be designed from an antibody-combining site. In the present work, a theoretical study of the different behavior of a germline catalytic antibody (CA) and its matured form, 39 A-11, that catalyze a Diels-Alder reaction has been carried out. A free-energy perturbation technique based on a hybrid quantum-mechanics/molecular-mechanics scheme, together with internal energy minimizations, has allowed free-energy profiles to be obtained for both CAs. The profiles show a smaller barrier for the matured form, which is in agreement with the experimental observation. Free-energy profiles were obtained with this methodology, thereby avoiding the much more demanding two-dimensional calculations of the energy surfaces that are normally required to study this kind of reaction. Structural analysis and energy evaluations of substrate-protein interactions have been performed from averaged structures, which allows understanding of how the single mutations carried out during the maturation process can be responsible for the observed fourfold enhancement of the catalytic rate constant. The conclusion is that the mutation effect in this studied germline CA produces a complex indirect effect through coupled movements of the backbone of the protein and the substrate. 相似文献
7.
Alves CN Martí S Castillo R Andrés J Moliner V Tuñón I Silla E 《Chemistry (Weinheim an der Bergstrasse, Germany)》2007,13(27):7715-7724
Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-containing active sites. Furthermore, an energy decomposition analysis was performed to determine the contributions of individual residues to the enzyme-inhibitor interactions on averaged structures obtained from rather extensive conformational sampling. Analysis of the results reveals first that there is a correlation between protein-ligand interaction energy and experimental strand transfer into human chromosomes and secondly that the Asn-155, Lys-156 and Lys-159 residues and the Mg(2+) ion are crucial to anti-HIV IN activity. These results may explain the available experimental data. 相似文献
8.
Martin T. Stiebritz 《Journal of computational chemistry》2015,36(8):553-563
Metalloenzymes represent a particular challenge for any rational (re)design approach because the modeling of reaction events at their metallic cofactors requires time‐consuming quantum mechanical calculations, which cannot easily be reconciled with the fast, knowledge‐based approaches commonly applied in protein design studies. Here, an approach for the exploration of sequence‐reactivity relationships in metalloenzymes is presented (MetREx) that consists of force field‐based screening of mutants that lie energetically between a wild‐type sequence and the global minimum energy conformation and which should, therefore, be compatible with a given protein fold. Mutant candidates are subsequently evaluated with a fast and approximate quantum mechanical/molecular mechanical‐like procedure that models the influence of the protein environment on the active site by taking partial charges and van der Waals repulsions into account. The feasibility of the procedure is demonstrated for the active site of [FeFe] hydrogenase from Desulfovibrio desulfuricans. The method described allows for the identification of mutants with altered properties, such as inhibitor‐coordination energies, and the understanding of the robustness of enzymatic reaction steps with respect to variations in sequence space. © 2015 Wiley Periodicals, Inc. 相似文献
9.
György G. Ferenczy 《Journal of computational chemistry》2013,34(10):862-869
The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009 , 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self‐consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree–Fock–Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self‐consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave‐function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem. © 2013 Wiley Periodicals, Inc. 相似文献
10.
Alexandra T. P. Carvalho Ana F. S. Teixeira Maria J. Ramos 《Journal of computational chemistry》2013,34(18):1540-1548
Iron‐sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron‐sulfur proteins and applied them in the molecular dynamics simulations of three iron‐sulfur proteins. We used density‐functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H‐bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins. © 2013 Wiley Periodicals, Inc. 相似文献
11.
The approximate density‐functional tight‐binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle–mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free‐energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron‐induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization. © 2015 Wiley Periodicals, Inc. 相似文献
12.
Tatiana Vasilevskaya Maria G. Khrenova Alexander V. Nemukhin Walter Thiel 《Journal of computational chemistry》2015,36(21):1621-1630
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. 相似文献
13.
A modified CHARMM force‐field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). The ZHB potential is advantageous for quantum mechanics/molecular mechanics simulations as it minimizes the electron spill‐out problems. In the same spirit, here we propose a modified ZHB potential (MZHB) by reformulating its bonding potential, while retaining the nonbonding potential as in the ZHB force‐field. We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. Further, transferability of MZHB is also tested for reproducing the crystallographic structures of several polymorphs of silica. © 2015 Wiley Periodicals, Inc. 相似文献
14.
Quantum Mechanics/Molecular Mechanics Studies on the Mechanism of Action of Cofactor Pyridoxal 5′‐Phosphate in Ornithine 4,5‐Aminomutase
下载免费PDF全文

Dr. Jiayun Pang Prof. Nigel S. Scrutton Prof. Michael J. Sutcliffe 《Chemistry (Weinheim an der Bergstrasse, Germany)》2014,20(36):11390-11401
A computational study was performed on the experimentally elusive cyclisation step in the cofactor pyridoxal 5′‐phosphate (PLP)‐dependent D ‐ornithine 4,5‐aminomutase (OAM)‐catalysed reaction. Calculations using both model systems and a combined quantum mechanics/molecular mechanics approach suggest that regulation of the cyclic radical intermediate is achieved through the synergy of the intrinsic catalytic power of cofactor PLP and the active site of the enzyme. The captodative effect of PLP is balanced by an enzyme active site that controls the deprotonation of both the pyridine nitrogen atom (N1) and the Schiff‐base nitrogen atom (N2). Furthermore, electrostatic interactions between the terminal carboxylate and amino groups of the substrate and Arg297 and Glu81 impose substantial “strain” energy on the orientation of the cyclic intermediate to control its trajectory. In addition the “strain” energy, which appears to be sensitive to both the number of carbon atoms in the substrate/analogue and the position of the radical intermediates, may play a key role in controlling the transition of the enzyme from the closed to the open state. Our results provide new insights into several aspects of the radical mechanism in aminomutase catalysis and broaden our understanding of cofactor PLP‐dependent reactions. 相似文献
15.
16.
The Quantum‐to‐molecular mechanics method (Q2MM) for converting quantum mechanical transition states (TSs) to molecular mechanical minima has been modified to allow a fit to the “natural” reaction mode eigenvalue. The resulting force field gives an improved representation of the energy curvature at the TS, but can potentially give false responses to steric interactions. Ways to address this problem while staying close to the “natural” TS force field are discussed. © 2014 Wiley Periodicals, Inc. 相似文献
17.
Hans Primas 《Theoretical chemistry accounts》1975,39(2):127-148
It is shown that the classical concept of an open system does not encompass quantal systems but has to be replaced by the non-Boolean notion of an entangled system. Molecular, chemical, or biological phenomena can be considered to be reduced to a fundamental theory like quantum mechanics only if the fundamental and the phenomenological theories are formally and interpretatively connected, and if the classifications used in the empirical sciences are shown to follow from a single set of fundamental dynamical laws. These conditions enforce a non-statistical and ontic interpretation of quantum mechanics, hence a non-Boolean calculus of propositions. In this interpretation the notion of a world state is well-defined, its Schmidt-decomposition defines a background-dependent model state for molecular systems and creates the phenomena we can observe. To any molecular system there is associated in an objective way a nonnegative number which we call the integrity. The integrity measures the inherent fuzziness of the system concept in a holistic theory, and is used to define recognizable molecular patterns. 相似文献
18.
The hydrogen‐capping method is one of the most popular and widely used coupling‐schemes for quantum mechanics/molecular mechanics (QM/MM)‐molecular dynamics simulations of macromolecular systems. This is mostly due to the fact that it is fairly convenient to implement and parametrize, thus providing an excellent compromise between accuracy and computational effort. In this work, a viable and straight‐forward approach to optimize the placing of the link atom on a suitable distance ratio between the frontier atoms is discussed. To further increase the accuracy, instead of global parameters for all amino acids, different parameter sets for each type of amino acid are derived. The dependency of the link bond parameters on the chemical environment and the used QM‐method is probed to assess the range of applicability of the parametrization. Suitable sets of parameters for RI‐MP2, B3LYP, (RI)‐B3LYP‐D3, and RI‐BLYP‐D3 at triple‐zeta level for all relevant proteinogenic amino acids are presented. Furthermore, the scope and range of the perturbation, stemming from the introduction of link bonds is evaluated through application of the presented QM/MM scheme in calculations of the active site of 15S‐lipoxygenase. © 2015 Wiley Periodicals, Inc. 相似文献
19.
We have run several molecular dynamics (MD) simulations on zinc-containing phosphotriesterase (PTE) with two bound substrates, sarin and paraoxon, and with the substrate analog diethyl 4-methylbenzylphosphonate. A standard nonbonded model was employed to treat the zinc ions with the commonly used charge of +2. In all the trajectories, we observed a tightly bound water (TBW) molecule in the active site that was coordinated to the less buried zinc ion. The phosphoryl oxygen of the substrate/inhibitor was found to be coordinated to the same zinc ion so that, considering all ligands, the less buried zinc was hexa-coordinated. The hexa-coordination of this zinc ion was not seen in the deposited X-ray pdb files for PTE. Several additional MD simulations were then performed using different charges (+1, +1.5) on the zinc ions, along with ab initio and density functional theory (DFT) calculations, to evaluate the following possibilities: the crystal diffraction data were not correctly interpreted; the hexa-coordinated zinc ion in PTE is only present in solution and not in the crystal; and the hexa-coordinated zinc ion in PTE is an artifact of the force field used. A charge of +1.5 leads to a coordination number (CN) of 5 on both zinc ions, which is consistent with the results from ab initio and DFT calculations and with the latest high resolution X-ray crystal structure. The commonly used charge of +2 produces a CN of 6 on the less buried zinc. The CN on the more buried zinc ion is 5 when the substrate/inhibitor is present in the simulation, and increases to 6 when the substrate/inhibitor is removed prior to the simulation. The results of both of the MD and quantum mechanical calculations lead to the conclusion that the zinc ions in the PTE active site are both penta-coordinated, and that the MD simulations performed with the charge of +2 overestimate the CN of the zinc ions in the PTE active site. The overall protein structures in the simulations remain unaffected by the change in zinc charge from +2 to +1.5. The results also suggest that the charge +1.5 is the most appropriate for the molecular dynamics simulations on zinc-containing PTE when a nonbonded model is used and no global thermodynamic conclusion is sought. We also show that the standard nonbonded model is not able to properly treat the CN and energy at the same time. A preliminary, promising charge-transfer model is discussed with the use of the zinc charge of +1.5. 相似文献
20.