The Jahn-Teller effect of the Cr(2+) ion in aqueous solution: Ab initio QM/MM molecular dynamics simulations

The hydration structure of Cr(2+) has been studied using molecular dynamics (MD) simulations including three-body corrections and combined ab initio quantum mechanical/molecular mechanical (QM/MM) MD simulations at the Hartree-Fock level. The structural properties are determined in terms of radial distribution functions, coordination numbers, and several angle distributions. The mean residence time was evaluated for describing ligand exchange processes in the second hydration shell. The Jahn-Teller distorted octahedral [Cr(H(2)O)(6)](2+) complex was pronounced in the QM/MM MD simulation. The first-shell distances of Cr(2+) are in the range of 1.9-2.8 A, which are slightly larger than those observed in the cases of Cu(2+) and Ti(3+). No first-shell water exchange occurred during the simulation time of 35 ps. Several water-exchange processes were observed in the second hydration shell with a mean residence time of 7.3 ps.     

NMR solvent shifts of adenine in aqueous solution from hybrid QM/MM molecular dynamics simulations

We present first principles calculations of the NMR solvent shift of adenine in aqueous solution. The calculations are based on snapshots sampled from a molecular dynamics simulation, which were obtained via a hybrid quantum-mechanical/mechanical modeling approach, using an all-atom force field (TIP3P). We find that the solvation via the strongly fluctuating hydrogen bond network of water leads to nontrivial changes in the NMR spectra of the solutes regarding the ordering of the resonance lines. Although there are still sizable deviations from experiment, the overall agreement is satisfactory for the 1H and 15N NMR shifts. Our work is another step toward a realistic first-principles prediction of NMR chemical shifts in complex chemical environments.     

Structure and dynamics of Au+ ion in aqueous solution: ab initio QM/MM MD simulations

Structure and dynamics of hydrated Au(+) have been investigated by means of molecular dynamics simulations based on ab initio quantum mechanical molecular mechanical forces at Hartree-Fock level for the treatment of the first hydration shell. The outer region of the system was described using a newly constructed classical three-body corrected potential. The structure was evaluated in terms of radial and angular distribution functions and coordination number distributions. Water exchange processes between coordination shells and bulk indicate a very labile structure of the first hydration shell whose average coordination number of 4.7 is a mixture of 3-, 4-, 5-, 6-, and 7-coordinated species. Fast water exchange reactions between first and second hydration shell occur, and the second hydration shell is exceptionally large. Therefore, the mean residence time of water molecules in the first hydration shell (5.6 ps/7.5 ps for t*= 0.5 ps/2.0 ps) is shorter than that in the second shell (9.4 ps/21.2 ps for t*= 0.5 ps/2.0 ps), leading to a quite specific picture of a "structure-breaking" effect.     

Structure and dynamics of high-spin Ru in aqueous solution: Ab initio QM/MM molecular dynamics simulation

The structural and dynamical properties of high-spin Ru²⁺ in aqueous solution have been theoretically studied using molecular dynamics (MD) simulations. The conventional MD simulation based on pair potentials gives the overestimated average first shell coordination number of 9, whereas the value of 5.9 was observed when the three-body corrected function was included. A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to take into account the many-body effects on the hydration shell structure of Ru²⁺. The most important region, the first hydration shell, was treated by ab initio quantum mechanics at UHF level using the SBKJC VDZ ECP basis set for Ru²⁺ and the 6-31G^∗ basis sets for water. An exact coordination number of 6 for the first hydration shell was obtained from the QM/MM simulation. The QM/MM simulation predicts the average Ru²⁺–O distance of 2.42 Å for the first hydration shell, whereas the values of 2.34 and 2.46 Å are resulted from the pair potentials without and with the three-body corrected simulations, respectively. Several other structural properties representing position and orientation of the solvate molecules were evaluated for describing the hydration shell structure of the Ru²⁺ ion in dilute aqueous solution. A mean residence time of 7.1 ps was obtained for water ligands residing in the second hydration shell.     

Structure and dynamics of the CrIII ion in aqueous solution: Ab initio QM/MM molecular dynamics simulation

Structural and dynamical properties of the Cr(III) ion in aqueous solution have been investigated using a combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation. The hydration structure of Cr(III) was determined in terms of radial distribution functions, coordination numbers, and angular distributions. The QM/MM simulation gives coordination numbers of 6 and 15.4 for the first and second hydration shell, respectively. The first hydration shell is kinetically very inert but by no means rigid and variations of the first hydration shell geometry lead to distinct splitting in the vibrational spectra of Cr(H(2)O)(6) (3+). A mean residence time of 22 ps was obtained for water ligands residing in the second hydration shell, which is remarkably shorter than the experimentally estimated value. The hydration energy of -1108 +/- 7 kcal/mol, obtained from the QM/MM simulation, corresponds well to the experimental hydration enthalpy value.     

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.     

A combined QM/MM molecular dynamics simulations study of nitrate anion (NO3-) in aqueous solution

The structural and dynamical properties of NO3- in dilute aqueous solution have been investigated by means of two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely HF/MM and B3LYP/MM, in which the ion and its surrounding water molecules were treated at HF and B3LYP levels of accuracy, respectively, using the DZV+ basis set. On the basis of both HF and B3LYP methods, a well-defined first hydration shell of NO3- is obtainable, but the shell is quite flexible and the hydrogen-bond interactions between NO3- and water are rather weak. With respect to the detailed analysis of the geometrical arrangement and vibrations of NO3-, the experimentally observed solvent-induced symmetry breaking of the ion is well reflected. In addition, the dynamical information, i.e., the bond distortions and shifts in the corresponding bending and stretching frequencies as well as the mean residence time of water molecules surrounding the NO3- ion, clearly indicates the "structure-breaking" ability of this ion in aqueous solution. From a methodical point of view it seems that both the HF and B3LYP methods are not too different in describing this hydrated ion by means of a QM/MM simulation. However, the detailed analysis of the dynamics properties indicates a better suitability of the HF method compared to the B3LYP-DFT approach.     

Theoretical modeling of open-shell molecules in solution: a QM/MM molecular dynamics approach

In this work, the GLOB model, an effective and reliable computational approach well suited for ab initio and QM/MM molecular dynamics simulations of complex molecular systems in solution, has been applied to study two representative open-shell systems, the cobalt(II) ion and the glycine radical in aqueous solution, with special reference to their structural and magnetic properties. The main structural features of the solvent cage around the cobalt ion and the hydrogen bonding patterns around the neutral and zwitterionic forms of the glycine radical have been investigated in some detail. The general good agreement with experiments supports the use of the present model to investigate more challenging and biological/technological relevant open-shell systems.     

Insights into intrastrand cross-link lesions of DNA from QM/MM molecular dynamics simulations

DNA damages induced by oxidative intrastrand cross-links have been the subject of intense research during the past decade. Yet, the currently available experimental protocols used to isolate such lesions only allow to get structural information about linked dinucleotides. The detailed structure of the damaged DNA macromolecule has remained elusive. In this study we generated in silico the most frequent oxidative intrastrand cross-link adduct, G[8,5-Me]T, embedded in a solvated DNA dodecamer by means of quantum mechanics/molecular mechanics (QM/MM) Car-Parrinello simulations. The free energy of activation required to bring the reactant close together and to form the C-C covalent-bond is estimated to be ~10 kcal/mol. We observe that the G[8,5-Me]T tandem lesion is accommodated with almost no perturbation of the Watson-Crick hydrogen-bond network and induces bend and unwinding angles of ~20° and 8°, respectively. This rather small structural distortion of the DNA macromolecule compared to other well characterized intrastrand cross-links, such as cyclobutane pyrimidines dimers or cisplatin-DNA complex adduct, is a probable rationale for the known lack of efficient repair of oxidative damages.     

The Jahn-Teller effect of the TiIII ion in aqueous solution: extended ab initio QM/MM molecular dynamics simulations.

Combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, including only the first and the first and second hydration shells in the QM region, were performed for TiIII in aqueous solution. The hydration structure of TiIII is discussed in terms of radial distribution functions, coordination-number distributions and several angle distributions. Dynamical properties, such as librational and vibrational motions and TiIII-O vibrations, were evaluated. A fast dynamical Jahn-Teller effect of TiIII(aq) was observed in the QM/MM simulations, in particular when the second hydration shell was included into the QM region. The results justify the computational effort required for the inclusion of the second hydration shell into the QM region and show the importance of this effort for obtaining accurate hydration-shell geometries, dynamical properties, and details of the Jahn-Teller effect.     

Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations

Protein splicing is a post-translational process in which a biologically inactive protein is activated by the release of a segment denoted as an intein. The process involves four steps. In the third, the scission of the intein takes place after the cyclization of the last amino acid of the segment, an asparagine. Little is known about the chemical reaction necessary for this cyclization. Experiments demonstrate that two histidines (the penultimate amino acid of the intein, and a histidine located 10 amino acids upstream) are relevant in the cyclization of the asparagine. We have investigated the mechanism and determinants of reaction in the GyrA intein focusing on the requirements for asparagine activation for its cyclization. First, the influence that the protonation states of these two histidines have on the orientation of the asparagine side chain is investigated by means of molecular dynamics simulation. Molecular dynamics simulations using the CHARMM27 force field were carried out on the three possible protonation states for each of these two histidines. The results indicate that the only protonation state in which the conformation of the system is suitable for cyclization is when the penultimate histidine is fully protonated (positively charged), and the upstream histidine is in the His(ε) neutral tautomeric form. The free energy profile for the reaction in which the asparagine is activated by a proton transfer to the upstream histidine is presented, computed by hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics at the SCCDFTB/CHARMM27 level of theory. The calculated free energy barrier for the reaction is 19.0 kcal mol(-1). B3LYP/6-31+G(d) QM/MM single-point calculations give a qualitatively a similar energy profile, although with somewhat higher energy barriers, in good agreement with the value derived from experiment of 25 kcal mol(-1) at 60 °C. QM/MM molecular dynamics simulations of the reactant, activated reactant and intermediate states highlight the importance of the Arg181-Val182-Asp183 segment in catalysing the reaction. Overall, the results indicate that nucleophilic activation of the asparagine for its cyclization by the upstream histidine acting as the base is a plausible mechanism for the C-terminal cleavage in protein splicing.     

pKa calculations in solution and proteins with QM/MM free energy perturbation simulations: a quantitative test of QM/MM protocols

The accuracy of biological simulations depends, in large part, on the treatment of electrostatics. Due to the availability of accurate experimental values, calculation of pKa provides stringent evaluation of computational methods. The generalized solvent boundary potential (GSBP) and Ewald summation electrostatic treatments were recently implemented for combined quantum mechanical and molecular mechanics (QM/MM) simulations by our group. These approaches were tested by calculating pKa shifts due to differences in electronic structure and electrostatic environment; the shifts were determined for a series of small molecules in solution, using various electrostatic treatments, and two residues (His 31, Lys 102) in the M102K T4-lysozyme mutant with large pKa shifts, using the GSBP approach. The calculations utilized a free energy perturbation scheme with the QM/MM potential function involving the self-consistent charge density functional tight binding (SCC-DFTB) and CHARMM as the QM and MM methods, respectively. The study of small molecules demonstrated that inconsistent electrostatic models produced results that were difficult to correct in a robust manner; by contrast, extended electrostatics, GSBP, and Ewald simulations produced consistent results once a bulk solvation contribution was carefully chosen. In addition to the electrostatic treatment, the pKa shifts were also sensitive to the level of the QM method and the scheme of treating QM/MM Coulombic interactions; however, simple perturbative corrections based on SCC-DFTB/CHARMM trajectories and higher level single point energy calculations were found to give satisfactory results. Combining all factors gave a root-mean-square difference of 0.7 pKa units for the relative pKa values of the small molecules compared to experiment. For the residues in the lysozyme, an accurate pKa shift was obtained for His 31 with multiple nanosecond simulations. For Lys 102, however, the pKa shift was estimated to be too large, even after more than 10 nanosecond simulations for each lambda window; the difficulty was due to the significant, but slow, reorganization of the protein and water structure when Lys 102 was protonated. The simulations support that Lys 102 is deprotonated in the X-ray structure and the protein is highly destabilized when this residue is protonated.     

Structure and dynamics of hydrated NH(4) (+): an ab initio QM/MM molecular dynamics simulation

A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to investigate solvation structure and dynamics of NH(4) (+) in water. The most interesting region, the sphere includes an ammonium ion and its first hydration shell, was treated at the Hartree-Fock level using DZV basis set, while the rest of the system was described by classical pair potentials. On the basis of detailed QM/MM simulation results, the solvation structure of NH(4) (+) is rather flexible, in which many water molecules are cooperatively involved in the solvation shell of the ion. Of particular interest, the QM/MM results show fast translation and rotation of NH(4) (+) in water. This phenomenon has resulted from multiple coordination, which drives the NH(4) (+) to translate and rotate quite freely within its surrounding water molecules. In addition, a "structure-breaking" behavior of the NH(4) (+) is well reflected by the detailed analysis on the water exchange process and the mean residence times of water molecules surrounding the ion.     

Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes

A new formalism for quantum mechanical / molecular mechanical (QM/MM) dynamics of chemical species in solution has been developed, which does not require the construction of any other potential functions except those for solvent–solvent interactions, maintains all the advantages of large simulation boxes and ensures the accuracy of ab initio quantum mechanics for all forces acting in the chemically most relevant region. Interactions between solute and more distant solvent molecules are incorporated by a dynamically adjusted force field corresponding to the actual molecular configuration of the simulated system and charges derived from the electron distribution in the solvate. The new formalism has been tested with some examples of hydrated ions, for which accurate conventional ab initio QM/MM simulations have been previously performed, and the comparison shows equivalence and in some aspects superiority of the new method. As this simulation procedure does not require any tedious construction of two-and three-body interaction potentials inherent to conventional QM/MM approaches, it opens the straightforward access to ab initio molecular dynamics simulations of any kind of solutes, such as metal complexes and other composite species in solution.     

Implementation of the IMOMM methodology for performing combined QM/MM molecular dynamics simulations and frequency calculations

A technique for implementing the integrated molecular orbital and molecular mechanics (IMOMM) methodology developed by Maseras and Morokuma that is used to perform combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, frequency calculations and simulations of macromolecules including explicit solvent is presented. Although the IMOMM methodology is generalized to any coordinate system, the implementation first described by Maseras and Morokuma requires that the QM and MM gradients be transformed into internal coordinates before they are added together. This coordinate transformation can be cumbersome for macromolecular systems and can become ill-defined during the course of a molecular dynamics simulation. We describe an implementation of the IMOMM method in which the QM and MM gradients are combined in the cartesian coordinate system, thereby avoiding potential problems associated with using the internal coordinate system. The implementation can be used to perform combined QM/MM molecular dynamics simulations and frequency calculations within the IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM frequency calculations. Received: 11 May 1998 / Accepted: 14 August 1998 / Published online: 16 November 1998     

Mechanism of triphosphate hydrolysis in aqueous solution: QM/MM simulations in water clusters

The mechanism of the hydrolysis reaction of the unprotonated methyl triphosphate (MTP) ester in water clusters has been modeled. The effective fragment potential based quantum mechanical-molecular mechanical (QM/MM) approach has been applied in the simulations. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step dissociative-type process similar to the case of the guanosine triphosphate (GTP) hydrolysis in the Ras-GAP protein complex (Grigorenko, B. L.; Nemukhin, A. V.; Topol, I. A.; Cachau, R. E.; Burt, S. K. Proteins: Struct., Funct., Bioinf. 2005, 60, 495). At the first stage, a unified action of environmental molecular groups and the catalytic water molecule leads to a substantial spatial separation of the gamma-phosphate group from the rest of the molecule. At the second stage, inorganic phosphate H2PO4- is formed from water and the metaphosphate anion PO3- through the chain of proton transfers along hydrogen bonds. The estimated activation barriers for MTP in aqueous solution at both stages (20 and 14 kcal/mol) are substantially higher than the corresponding barriers for the GTP hydrolysis in the protein.     

Implementation of the SCC-DFTB method for hybrid QM/MM simulations within the amber molecular dynamics package

Self-consistent charge density functional tight-binding (SCC-DFTB) is a semiempirical method based on density functional theory and has in many cases been shown to provide relative energies and geometries comparable in accuracy to full DFT or ab initio MP2 calculations using large basis sets. This article shows an implementation of the SCC-DFTB method as part of the new QM/MM support in the AMBER 9 molecular dynamics program suite. Details of the implementation and examples of applications are shown.