Ab initio QM/MM dynamics of anion–water hydrogen bonds in aqueous solution

Relative cooperativity in neutral and charged molecular clusters using QM/MM calculations

QM/MM methods have been used to study electronic structure properties and chemical reactivity in complex molecular systems where direct electronic structure calculations are not feasible. In our previous work, we showed that non-polarizable force fields, by design, describe intermolecular interactions through pairwise interactions, overlooking many-body interactions involving three or more particles. In contrast, polarizable force fields account partially for many-body effects through polarization, but still handle van der Waals and permanent electrostatic interactions pairwise. We showed that despite those limitations, polarizable and non-polarizable force fields can reproduce relative cooperativity achieved using density functional theory due to error compensation mechanisms. In this contribution, we assess the performance of QM/MM methods in reproducing these phenomena. Our study highlights the significance of the QM region size and force field choice in QM/MM calculations, emphasizing the importance of parameter validation to obtain accurate interaction energy predictions.

Cooperativity and Frustration Effects (or Lack Thereof) in Polarizable and Non-polarizable Force Fields

Understanding cooperativity and frustration is crucial for studying biological processes such as molecular recognition and protein aggregation. Force fields have been extensively utilized to explore cooperativity in the formation of protein secondary structures and self-assembled systems. Multiple studies have demonstrated that polarizable force fields provide more accurate descriptions of this phenomenon compared to fixed-charge pairwise nonpolarizable force fields, thanks to the incorporation of polarization effects. In this study, we assess the performance of the AMOEBA polarizable force field and the AMBER and OPLS nonpolarizable pairwise force fields in capturing positive and negative cooperativity recently explored in neutral and charged molecular clusters using density functional theory. Our findings show that polarizable and nonpolarizable force fields qualitatively reproduce the relative cooperativity observed in electron structure calculations. However, AMBER and OPLS fail to describe absolute cooperativity. In contrast, AMOEBA accounts for the absolute cooperativity by considering interactions beyond pairwise interactions. According to the energy decomposition analysis, it is observed that the electrostatic interactions calculated with the AMBER and OPLS force fields seem to play an important and counterintuitive role in reproducing the adiabatic interaction energies calculated with density functional theory. However, it is important to note that these force fields, due to their nature, do not explicitly incorporate many-body effects, which limits their ability to accurately describe cooperativity. On the other hand, frustration in polarizable and nonpolarizable force fields is caused by changes in bond stretching and angle bending terms of the building blocks when they are forming a complex.

Why Are Saccharides Dehydrated in the Presence of Electrolytes? Insights from Molecular Modeling and Thermodynamic Measurements

The mechanisms governing the interactions of neutral polar solutes with ions in aqueous solutions are still poorly understood, despite the importance of this phenomenon in many fields (chemistry, physicochemistry, biology, food industries). In order to go further through the understanding of the molecular mechanisms governing the ions' specific effects, this paper presents a generic method dealing with the characterization and understanding of interactions between saccharides and ions in aqueous systems. For that, an original approach combining a computational technique and experimental measurements (thermodynamic properties) is proposed to explain and rationalize the relationship between the solute hydration and the physical chemistry of the ions in solution (cation/anion, charge, size, and hydration). These relationships make it possible to evaluate the hydration state of a saccharide, a polar neutral molecule, according to the ionic composition, from the knowledge of the ions' hydration properties. This work proposes new insight into molecular mechanisms governing the polar neutral solute/ion interactions and a new understanding of the hydration phenomenon in electrolytic solutions.

Dynamics of actinyl ions in water: a molecular dynamics simulation study

The dynamics of actinyl ions (AnO₂ⁿ⁺) in aqueous solutions is important not only for the design of advanced separation processes but also for understanding the fate of actinides in the environment.

Structure of an accurate ab initio model of the aqueous Cl- ion at high temperatures

The structure of an accurate ab initio model of aqueous chloride ion was calculated at two high-temperature state points (573 K, 0.725 g/cm(3) and 723 K, 0.0098 g/cm(3)) by a two-step procedure. First, the structure of an approximate model was calculated from a molecular dynamics simulation of the model. Then the difference between the structure of the ab initio model and the approximate model was calculated by non-Boltzmann weighting of a sample of configurations taken from the approximate model simulation. Radial distribution functions, average coordination numbers, the distribution of coordination numbers, an analysis of orientations of water in the first coordination shell, and the free energy of hydration of the chloride ion are reported for both state points. The most common water structure has one hydrogen close to the chloride ion and one pointing away (46% at 573 K and 57% at 723 K). Waters in the first coordination shell that are not strongly bound to the chloride ions are common. Several variations of the method were tested. Models in which the water-water interaction is calculated with ab initio methods predict only a slightly different structure than models in which water-water interactions are determined from the approximate models. Similarly, using the approximate model for solute-water interactions when the water is far from the chloride ion did not affect the results. Uncertainties due to the limited sample of configurations are estimated and found to be small. The results are in qualitative agreement with X-ray and neutron diffraction experiments and with simulations of approximate models.

Specific Anion Effects on the Optical Rotation of α-Amino Acids

Changes in optical rotation of some alpha-amino acids are induced by electrolytes. Such effects on l- and d-enantiomers of a range of amino acids are explored for sodium salts with varying anion. The amino acids studied were alanine, aspartic acid, glutamic acid, glutamine, proline, threonine, and tryptophan. The anion's polarizability in solution accounts for the change in [alpha] only for the halides. Self-association of amino acids in solution and pH changes due to the presence of the electrolytes do not account for the observed variations in optical activity. Specific interactions of anions with the chiral amino acids (Hofmeister effects) and salt-induced perturbations of the amino acid hydration shell appear to be responsible for the effects, and conformational changes in the chiral solutes due to the presence of ionic species are discussed.

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.

Structure and Dynamics of Sulfate Ion in Aqueous SolutionAn ab initio QMCF MD Simulation and Large Angle X-ray Scattering Study

The hydrated sulfate ion has been characterized in aqueous solution in structural and dynamic aspects using ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) simulation and large angle X-ray scattering (LAXS) methods. The LAXS data show an average coordination number of the sulfate ion of up to 12 water molecules bound through hydrogen bonding, while the QMCF MD simulation displays a wide range of coordination numbers between 8 and 14 with an average value of approximately 11. The Os...Ow distance cannot be distinguished from the Ow...Ow distance in the LAXS experiment; the weighted mean O...O distance is 2.880(10) A. In the simulation, the Os...Ow and Ow...Ow distances are found to be very similar, namely, 2.86 and 2.84 A, respectively. The S-Os bond and S...Ow distance have been determined by the LAXS experiment as 1.495(6) and 3.61(2) A, respectively, indicating an average nearly tetrahedral S-Os...Ow angle. The approximately 5% deviations of simulation distances (1.47 and 3.82 A) from the experimental ones can probably be ascribed to the neglect of correlation energy in the quantum mechanical method. The mean residence time of water ligands at O atoms, 2.57 ps, is longer than that in pure water, 1.7 ps, characterizing the sulfate ion as a weak structure maker.

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.

Simulations of Liquid Ammonia Based on the Combined Quantum Mechanical/Molecular Mechanical (QM/MM) Approach

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely, HF/MM and B3LYP/MM, have been performed to investigate the local structure and dynamics of liquid ammonia. The most interesting region, a sphere containing a central reference molecule and all its nearest surrounding molecules (first coordination shell), was treated by the Hartree-Fock (HF) and hybrid density functional B3LYP methods, whereas the rest of the system was described by the classical pair potentials. On the basis of both HF and B3LYP methods, it is observed that the hydrogen bonding in this peculiar liquid is weak. The structure and dynamics of this liquid are suggested to be determined by the steric packing effects, rather than by the directional hydrogen bonding interactions. Compared to previous empirical as well as Car-Parrinello (CP) molecular dynamics studies, our QM/MM simulations provide detailed information that is in better agreement with experimental data.