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.     

Structural and dynamical properties of the Hg2+ aqua ion: a molecular dynamics study

Molecular dynamics simulations of the Hg2+ ion in aqueous solution have been carried out using an effective two-body potential derived from quantum mechanical calculations. A stable heptacoordinated structure of the Hg2+ first hydration shell has been observed and confirmed by extended X-ray absorption fine structure (EXAFS) experimental data. The structural properties of the Hg2+ hydration shells have been investigated using radial and angular distribution functions, while the dynamical behavior has been discussed in terms of reorientational correlation functions, mean residence times of water molecules in the first and second hydration shells, and self-diffusion coefficients. The effect of water-water interactions on the Hg2+ hydration properties has been evaluated using the SPC/E and TIP5P water models.     

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.     

Sr(II) in water: A labile hydrate with a highly mobile structure

Despite the large number of experimental as well as theoretical investigations available in the literature, some properties of the hydration structure of Sr(II), for example, the coordination number, are still ambiguous. The presented molecular dynamics study based on a most suitable ab initio QM/MM protocol allowed a detailed investigation of structural and dynamical properties of this hydrate, which shows a considerable degree of internal flexebility as well as ligand mobility within the first shell. Despite the high computational effort an exceptionally long QM/MM simulation had to be carried out to obtain sufficient information to investigate first shell ligand exchange reactions.     

New Insights into the Jahn–Teller Effect through ab initio Quantum-Mechanical/Molecular-Mechanical Molecular Dynamics Simulations of CuII in Water

The Cu^II hydration shell structure has been studied by means of classical molecular dynamics (MD) simulations including three-body corrections and hybrid quantum-mechanical/molecular-mechanical (QM/MM) molecular dynamics (MD) simulations at the Hartree–Fock level. The copper(II ) ion is found to be six-fold coordinated and [Cu(H₂O)₆]²⁺ exhibits a distorted octahedral structure. The QM/MM MD approach reproduces correctly the experimentally observed Jahn–Teller effect but exhibits faster inversions (<200 fs) and a more complex behaviour than expected from experimental data. The dynamic Jahn–Teller effect causes the high lability of [Cu(H₂O)₆]²⁺ with a ligand-exchange rate constant some orders of magnitude higher than its neighbouring ions Ni^II and Zn^II. Nevertheless, no first-shell water exchange occurred during a 30-ps simulation. The structure of the hydrated ion is discussed in terms of radial distribution functions, coordination numbers, and various angular distributions and the dynamical properties as librational and vibrational motions and reorientational times were evaluated, which lead to detailed information about the first hydration shell. Second-shell water-exchange processes could be observed within the simulation time scale and yielded a mean ligand residence time of ≈20 ps.     

An extended ab initio QM/MM MD approach to structure and dynamics of Zn(II) in aqueous solution

Structural and dynamical properties of Zn(II) in aqueous solution were investigated, based on an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation at double-zeta Hartree-Fock quantum mechanical level including the first and second hydration shells into the QM region. The inclusion of the second shell in the QM region resulted in significant changes in the properties of the hydrate. The first shell coordination number was found to be 6, the second shell consists of approximately 14 water molecules. The structural properties were determined in terms of RDF, ADF, tilt and theta angle distributions, while dynamics were characterized by mean ligand residence times, ion-ligand stretching frequencies and the vibrational and librational motions of water ligands.     

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.     

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.     

Quantum mechanical/molecular mechanical simulations of the Tl(III) ion in water

Classical molecular dynamics (MD) and combined quantum mechanical/molecular mechanical (QM/MM) MD simulations have been performed to investigate the structural and dynamical properties of the Tl(III) ion in water. A six-coordinate hydration structure with a maximum probability of the Tl-O distance at 2.21 A was observed, which is in good agreement with X-ray data. The librational and vibrational spectra of water molecules in the first hydration shell are blue-shifted compared with those of pure liquid water, and the Tl-O stretching force constant was evaluated as 148 Nm(-1). Both structural and dynamical properties show a distortion of the first solvation shell structure. The second shell ligands' mean residence time was determined as 12.8 ps. The Tl(III) ion can be classified as "structure forming" ion; the calculated hydration energy of -986 +/- 9 kcal mol agrees well with the experimental value of -986 kcal mol.     

Quantum mechanical charge field molecular dynamics simulation of the TiO2+ ion in aqueous solution

Structural and dynamical properties of the TiO(2+) ion in aqueous solution have been investigated by using the new ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism, which does not require any other potential functions except those for solvent-solvent interactions. Both first and second hydration shell have been treated at Hartree-Fock (HF) quantum mechanical level. A Ti-O bond distance of 1.5 A was observed for the [Ti=O](2+) ion. The first hydration shell of the ion shows a varying coordination number ranging from 5 to 7, five being the dominant one and representing one axial and four equatorial water molecules directly coordinated to Ti, which are located at 2.3 A and 2.1 A, respectively. The flexibility in the coordination number reflects the fast exchange processes, which occur only at the oxo atom, where water ligands are weakly bound through hydrogen bonds. Considering the first shell hydration, the composition of the TiO(2+) hydrate can be characterized as [(H(2)O)(0.7)(H(2)O)(4) (eq)(H(2)O)(ax)](2+). The second shell consists in average of 12 water molecules located at a mean distance of 4.4 A. Several other structural parameters such as radial and angular distribution functions and coordination number distributions were analyzed to fully characterize the hydration structure of the TiO(2+) ion in aqueous solution. For the dynamics of the TiO(2+) ion, different sets of dynamical parameters such as Ti=O, Ti-O(eq), and Ti-O(ax) stretching frequencies and ligands' mean residence times were evaluated. During the simulation time of 15 ps, 3 water exchange processes in the first shell were observed at the oxo atom, corresponding to a mean residence time of 3.6 ps. The ligands' mean residence time for the second shell was determined as 3.5 ps.     

Characterization of structure and dynamics of an aqueous scandium(iii) ion by an extended ab initio QM/MM molecular dynamics simulation

Hydration structure and dynamics of an aqueous Sc(iii) solution were characterized by means of an extended ab initio quantum mechanical/molecular dynamical (QM/MM) molecular dynamics simulation at Hartree-Fock level. A monocapped trigonal prismatic structure composed of seven water molecules surrounding scandium(iii) ion was proposed by the QM/MM simulation including the quantum mechanical effects for the first and second hydration shells. The mean Sc(iii)-O bond length of 2.14 ? was identified for six prism water molecules with one capping water located at around 2.26 ?, reproducing well the X-ray diffraction data. The Sc(iii)-O stretching frequency of 432 cm(-1) corresponding to a force constant of 130 N m(-1), evaluated from the enlarged QM/MM simulation, is in good agreement with the experimentally determined value of 430 cm(-1) (128 N m(-1)). Various water exchange processes in the second hydration shell of the hydrated Sc(iii) ion predict a mean ligand residence time of 7.3 ps.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

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.     

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.     

Ab initio QM/MM dynamics of H3O+ in water

A molecular dynamics (MD) simulation based on a combined ab initio quantum mechanics/molecular mechanics (QM/MM) method has been performed to investigate the solvation structure and dynamics of H3O+ in water. The QM region is a sphere around the central H3O+ ion, and contains about 6-8 water molecules. It is treated at the Hartree-Fock (HF) level, while the rest of the system is described by means of classical pair potentials. The Eigen complex (H9O4+) is found to be the most prevalent species in the aqueous solution, partly due to the selection scheme of the center of the QM region. The QM/MM results show that the Eigen complex frequently converts back and forth into the Zundel (H5O2+) structure. Besides the three nearest-neighbor water molecules directly hydrogen-bonded to H3O+, other neighbor waters, such as a fourth water molecule which interacts preferentially with the oxygen atom of the hydronium ion, are found occasionally near the ion. Analyses of the water exchange processes and the mean residence times of water molecules in the ion's hydration shell indicate that such next-nearest neighbor water molecules participate in the rearrangement of the hydrogen bond network during fluctuative formation of the Zundel ion and, thus, contribute to the Grotthuss transport of the proton.     

Be(II) in aqueous solution--an extended ab initio QM/MM MD study

An ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation at double-zeta restricted Hartree-Fock (RHF) level was performed at 293.15 K, including first and second hydration shell in the QM region to study the structural and dynamical properties of the Be(II)-hydrate in aqueous solution. The first tetrahedrally arranged hydration shell, with the four water molecules located at a mean Be-O distance of 1.61 A, is highly inert with respect to ligand exchange processes. The second shell, however, consisting in average of approximately 9.2 water ligands at a mean Be-O distance of 3.7 A and the third shell at a mean Be-O distance of 5.4 A with approximately 19 ligands rapidly exchange water molecules between them and with the bulk, respectively. Other structural parameters such as radial and angular distribution functions (RDF and ADF) and tilt- and theta-angle distributions were also evaluated. The dynamics of the hydrate were studied in terms of ligand mean residence times (MRTs) and librational and vibrational frequencies. The mean residence times for second shell and third shell ligands were determined as 4.8 and 3.2 ps, respectively. The Be-O stretching frequency of 658 cm(-1), associated with a force constant of 147 N m(-1) could be overestimated but it certainly reflects the exceptional stability of the ion-ligand bond in the first hydration shell.