The influence of quantum forces on molecular dynamics simulation results for hydrated aluminium(III)

Anomalous Water Penetration in Al3+ Dissolution

The physicochemical characterization of trivalent ions is limited due to a lack of accurate force fields. By leveraging the latest machine learning force field to model aqueous AlCl3, we discover that upon dissolution of Al3+, water molecules beyond the second hydration shell are involved in the hydration process. A combination of scissoring of coordinating water is followed by synchronized secondary motion of water in the second solvation shell due to hydrogen bonding. Consequently, the water beyond the second solvation penetrates through the second solvation shell and coordinates to the Al3+. Our study reveals a novel microscopic understanding of solvation dynamics for the trivalent ion.

The Jahn-Teller effect in mixed aqueous solution: the solvation of Cu2+ in 18.6% aqueous ammonia obtained from ab initio quantum mechanical charge field molecular dynamics

The solvation structure and dynamics of Cu2+ in 18.6 % aqueous ammonia have been investigated using an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach at the Hartree–Fock (HF) level of theory applying the LANL2DZ ECP and Dunning DZP basis sets for Cu2+, ammonia and water, respectively. During a simulation time of 20 ps, only NH3 molecules are observed within the first solvation shell of Cu2+, resulting in the formation of an octahedral [Cu(NH3)6]2+ complex. While no exchange of these ligands with the second solvation shell are observed along the simulation, the monitoring of the associated N-Ntrans distances highlight the dynamics of the associated Jahn-Teller distortions, showing on average 2 elongated axial (2.19 Å) and 4 equatorial Cu–N bonds (2.39 Å). The observed structural properties are found in excellent agreement with experimental studies. In addition, an NBO analysis was carried out, confirming the strong electrostatic character of the Cu2+–NH3 interaction.

Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li+ Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations

A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δ_s G_M⁺,wat^◦ of the Li⁺ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χ_wat^◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol^-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol^-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li⁺ are also compared to those for Na⁺ and K⁺ from a previous study in terms of relative hydration free energies ΔΔ_s G_M⁺,wat^◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔ_s G_M⁺,wat^• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li⁺ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.

Ultrasensitivity of Water Exchange Kinetics to the Size of Metal Ion

Metal ions play a vital role in many biological processes. An important factor in these processes is the dynamics of exchange between ion bound-water molecules and the bulk. Although structural and dynamical properties of labile waters bound to metal ions, such as Na⁺ and Ca²⁺, can be elucidated using molecular dynamics simulations, direct evaluation of rates of exchange of waters rigidly bound to high charge density Mg²⁺, has been elusive. Here, we report a universal relationship, allowing us to determine the water exchange time on metal ions as a function of valence and hydration radius. The proposed relationship, which covers times spanning 14 orders of magnitude, highlights the ultrasensitivity of water lifetime to the ion size, as exemplified by divalent ions, Ca²⁺ (∼100 ps) and Mg²⁺ (∼1.5 μs). We show that even when structures, characterized by radial distributions are similar, a small difference in hydration radius leads to a qualitatively different (associative or dissociative) mechanism of water exchange. Our work provides a theoretical basis for determination of hydration radius, which is critical for accurately modeling the water dynamics around multivalent ions, and hence in describing all electrostatically driven events such as ribozyme folding and catalysis.

Urea in aqueous solution studied by quantum mechanical charge field-molecular dynamics (QMCF-MD)

This work presents a quantum mechanical charge field-molecular dynamics (QMCF-MD) simulation of urea in dilute aqueous solution. Detailed data for structure and dynamics are provided and compared to previous works of other groups. Radial and angular distributions are employed, as well as higher degree spatial investigations, two-dimensional particle mapping, volume maps and the previously proposed SLICE formalism. Information on dynamical properties are presented in the form of hydrogen bond correlation functions and mean lifetime analysis based on weighted Voronoi decomposition. Dihedral and tilt/theta angle distributions substantiate the previous findings of other groups, that urea is far from being planar within aqueous solution. In addition to the analysis of the complete hydration shell, several specific regions of hydration have been identified, for which individual analysis has been performed in terms of hydrogen bond lifetime correlation functions and re-orientational times. A decomposition study based on Laguerre tessellation further investigates the structure and dynamics of the individual hydration layers. It is found that urea does not show properties found in the case of typical structure breaking agents, such as Rb(+) or Cs(+), which is in accordance with spectroscopic data of Rezus and Bakker.

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.

Structure and dynamics of the U4+ ion in aqueous solution: an ab initio quantum mechanical charge field molecular dynamics study

The structure and dynamics of the stable four-times positively charged uranium(IV) cation in aqueous solution have been investigated by ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) simulation at the Hartree-Fock double-zeta quantum mechanical level. The QMCF-MD approach enables investigations with the accuracy of a quantum mechanics/molecular mechanics approach without the need for the construction of solute-solvent potentials. Angular distribution functions; radial distribution functions; coordination numbers of the first, second, and third shell (9, 19, and 44, respectively); coordination number distribution functions; tilt- and Theta-angle distribution functions; as well as local density corrected triangle distribution functions have been employed for the evaluation of the hydrated ion's structure. Special attention was paid to the determination of the geometry of the first hydration layer, and the results were compared to experimental large-angle X-ray scattering and extended X-ray absorption fine structure data. The solvent dynamics around the ion were also investigated using mean ligand residence times and related data and, resulting from the unavailability of any experimental data, were compared to ions with similar properties.

Al(III) hydration revisited. An ab initio quantum mechanical charge field molecular dynamics study

To assess the novel quantum mechanical charge field (QMCF) molecular dynamics (MD) approach, two simulations of hydrated Al(III) have been carried out, as this system proved to be a well-suited test case for hybrid ab initio/molecular mechanics simulations. Two different population analysis schemes according to Mulliken and Lowdin have been applied to evaluate the atomic charges in the QM region. It is shown that the QMCF MD approach yields a substantially improved description of the system and that, due to the fact that solute-solvent potentials can be renounced, the QMCF MD framework is a more convenient approach to investigate solvated systems compared to conventional ab initio QM/MM MD approaches.

A theoretical study of the principles regulating the specificity for Al(III) against Mg(II) in protein cavities

Several toxic effects arise from Al's presence in living systems, one of them being the alteration of the natural role of enzymes and non-enzyme proteins. Al(III) is capable of entering protein active sites that in normal conditions should be occupied by other metals. Even if Mg(II) is an ubiquitous metal in biological systems, the interference of aluminium in magnesium metabolism is not clear yet. In this work, a systematic study of the affinity of Al(III) for different protein binding sites is presented, with special attention on structural parameters, the role of the charge and the presence of different ligands in the protein cavity. The specificity of the binding site for Al(III) against Mg(II) has been studied, and also the thermodynamical propensity of a Mg(II)/Al(III) exchange. Quantum mechanical methods that proved to be reliable in previous works have been used, namely, the density functional theory (DFT) and polarizable continuum model (PCM).

Adaptive Partitioning in Combined Quantum Mechanical and Molecular Mechanical Calculations of Potential Energy Functions for Multiscale Simulations

In many applications of multilevel/multiscale methods, an active zone must be modeled by a high-level electronic structure method, while a larger environmental zone can be safely modeled by a lower-level electronic structure method, molecular mechanics, or an analytic potential energy function. In some cases though, the active zone must be redefined as a function of simulation time. Examples include a reactive moiety diffusing through a liquid or solid, a dislocation propagating through a material, or solvent molecules in a second coordination sphere (which is environmental) exchanging with solvent molecules in an active first coordination shell. In this article, we present a procedure for combining the levels smoothly and efficiently in such systems in which atoms or groups of atoms move between high-level and low-level zones. The method dynamically partitions the system into the high-level and low-level zones and, unlike previous algorithms, removes all discontinuities in the potential energy and force whenever atoms or groups of atoms cross boundaries and change zones. The new adaptive partitioning (AP) method is compared to Rode's "hot spot" method and Morokuma's "ONIOM-XS" method that were designed for multilevel molecular dynamics (MD) simulations. MD simulations in the microcanonical ensemble show that the AP method conserves both total energy and momentum, while the ONIOM-XS method fails to conserve total energy and the hot spot method fails to conserve both total energy and momentum. Two versions of the AP method are presented, one scaling as O(2N) and one with linear scaling in N, where N is the number of groups in a buffer zone separating the active high-level zone from the environmental low-level zone. The AP method is also extended to systems with multiple high-level zones to allow, for example, the study of ions and counterions in solution using the multilevel approach.

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.