Structure and dynamics of La(III) in aqueous solution – An ab initio QM/MM MD approach

Ab initio QM/MM MD simulations have allowed to clarify some of the ambiguities arising from various studies on the hydrated La(III) ion. Both nine- and ten-coordinated hydrates co-exist and interchange in a dissociative process on the nano- or even subnanosecond scale, and thus much faster than any other trivalent main group or transition metal ions. The weak ion–ligand bond (53 N/m) supplies a reasonable explanation for it. The simulation results for La(III) are also compared to those for the isoelectronic ions Cs(I) and Ba(II) obtained by the same ab initio MD procedure, leading to conclusions on the influence of central ion charge on structural and dynamic properties of hydrate complexes.     

Revisiting the structure and dynamics of hydrated Cd2+ in aqueous solutions: Insights from the RI-SCS-MP2/MM molecular dynamics simulation

The spin component scale MP2/molecular mechanics molecular dynamics simulation investigated the hydration shell formation and hydrated Cd²⁺ dynamics in the water environment. At the first hydration shell, six water molecules with 2.27 Å for the average distance between water and Cd²⁺. Dynamical properties were analyzed by computing the water molecule's mean residence time (MRT) in its first and second hydration shells. The MRT of each shell was determined to be 31.8 and 1.92 ps, suggesting the strong influence of Cd²⁺ in the first hydration shell. The second shell was labile, with an average number of water molecules being 18. Despite the strong interaction between Cd²⁺ and water molecules in the first shell, the influence of ions in the second hydration shell remained weak.     

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.     

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.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

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.     

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.     

Reaching multi‐nanosecond timescales in combined QM/MM molecular dynamics simulations through parallel horsetail sampling

We report an enhanced sampling technique that allows to reach the multi‐nanosecond timescale in quantum mechanics/molecular mechanics molecular dynamics simulations. The proposed technique, called horsetail sampling, is a specific type of multiple molecular dynamics approach exhibiting high parallel efficiency. It couples a main simulation with a large number of shorter trajectories launched on independent processors at periodic time intervals. The technique is applied to study hydrogen peroxide at the water liquid–vapor interface, a system of considerable atmospheric relevance. A total simulation time of a little more than 6 ns has been attained for a total CPU time of 5.1 years representing only about 20 days of wall‐clock time. The discussion of the results highlights the strong influence of the solvation effects at the interface on the structure and the electronic properties of the solute. © 2017 Wiley Periodicals, Inc.     

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     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

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.     

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.     

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.     

Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach

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.     

应用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区的性质.     

Structure-breaking effects of solvated Rb(I) in dilute aqueous solution--an ab initio QM/MM MD approach

Structural properties of the hydrated Rb(I) ion have been investigated by ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations at the double-zeta HF quantum mechanical level. The first shell coordination number was found to be 7.1, and several other structural parameters such as angular distribution functions, radial distribution functions and tilt- and theta-angle distributions allowed the full characterization of the hydration structure of the Rb(I) ion in dilute aqueous solution. Velocity autocorrelation functions were used to calculate librational and vibrational motions, ion-ligand motions, as well as reorientation times. Different dynamical parameters such as water reorientation, mean ligand residence time, the number of ligand exchange processes, and rate constants were also analyzed. The mean ligand residence time for the first shell was determined as tau = 2.0 ps.