Search

Many well-established classical biomolecular force fields, fitted on the solvation properties of single ions, do not necessarily describe all the details of ion pairing accurately, especially for complex polyatomic ions. Depending on the target application, it might not be sufficient to reproduce the thermodynamics of ion pairing, but it may also be necessary to correctly capture structural details, such as the coordination mode. In this work, we analyzed how classical force fields can be optimized to yield a realistic description of these different aspects of ion pairing. Given the prominent role of the interactions of negatively charged amino-acid side chains and divalent cations in many biomolecular systems, we chose calcium acetate as a benchmark system to devise a general optimization strategy that we applied to two popular force fields, namely, GROMOS and OPLS-AA. Using experimental association constants and first-principles molecular dynamics simulations as a reference, we found that small modifications of the van der Waals ion-ion interaction parameters allow a systematic improvement of the essential thermodynamic and structural properties of ion pairing.

Although they have been widely used as models for oxide surfaces, the deprotonation behaviors of the Keggin-ions (MeAl127+) and typical oxide surfaces are very different. On Keggin-ions, the deprotonation occurs over a very narrow pH range at odds with the broad charging curve of larger oxide surfaces. Depending on the Me concentration, the deprotonation curve levels off sooner (high Me concentration) or later (for low Me concentration). The leveling off shows the onset of aggregation before which the Keggin-ions are present as individual units. We show that the atypical titration data previously observed for some GaAl12 solutions in comparison to the originally reported data can be explained by the presence of Ga2Al11 ions. The pKa value of aquo-groups bound to octahedral Ga was determined from ab initio molecular dynamics simulations relative to the pure GaAl12 ions. Using these results within a surface complexation model, the onset of deprotonation of the crude solution is surprisingly well predicted and the ratio between the different species is estimated to be in the proportion 20 (Ga2Al11) : 20 (Al13) : 60 (GaAl12). [ABSTRACT FROM AUTHOR]

Room temperature ionic liquids (RTILs) have been shown to exhibit spatial heterogeneity or structural heterogeneity in the sense that they form hydrophobic and ionic domains. Yet studies of the relationship between this structural heterogeneity and the ∼picosecond motion of the molecular constituents remain limited. In order to obtain insight into the time scales relevant to this structural heterogeneity, we perform molecular dynamics simulations of a series of RTILs. To investigate the relationship between the structures, i.e., the presence of hydrophobic and ionic domains, and the dynamics, we gradually increase the size of the hydrophobic part of the cation from ethylammonium nitrate (EAN), via propylammonium nitrate (PAN), to butylammonium nitrate (BAN). The two ends of the organic cation, namely, the charged Nhead–H group and the hydrophobic Ctail–H group, exhibit rotational dynamics on different time scales, evidencing dynamical heterogeneity. The dynamics of the Nhead–H group is slower because of the strong coulombic interaction with the nitrate counter-ionic anions, while the dynamics of the Ctail–H group is faster because of the weaker van der Waals interaction with the surrounding atoms. In particular, the rotation of the Nhead–H group slows down with increasing cationic chain length, while the rotation of the Ctail–H group shows little dependence on the cationic chain length, manifesting that the dynamical heterogeneity is enhanced with a longer cationic chain. The slowdown of the Nhead–H group with increasing cationic chain length is associated with a lower number of nitrate anions near the Nhead–H group, which presumably results in the increase of the energy barrier for the rotation. The sensitivity of the Nhead–H rotation to the number of surrounding nitrate anions, in conjunction with the varying number of nitrate anions, gives rise to a broad distribution of Nhead–H reorientation times. Our results suggest that the asymmetry of the cations and the larger excluded volume for longer cationic chain are important for both the structural heterogeneity and the dynamical heterogeneities. The observed dynamical heterogeneities may affect the rates of chemical reactions depending on where the reactants are solvated in ionic liquids and provide an additional guideline for the design of RTILs as solvents. [ABSTRACT FROM AUTHOR]

A combination of density functional theory (DFT) and classical molecular dynamics is used to quantitatively estimate the key redox properties of two variants of rubredoxin, a small and comparatively simple iron-sulfur protein. The results have confirmed that in the case of rubredoxin the redox potential modulation is obtained by the short-range contribution, while the reorganization free energy is dominated by long-range effects.

We propose a new force field for trimethylamine N-oxide (TMAO), which is designed to reproduce the long-lived and highly directional hydrogen bond between the TMAO oxygen (OTMAO) atom and surrounding water molecules. Based on the data obtained by ab initio molecular dynamics simulations, we introduce three dummy sites around OTMAO to mimic the OTMAO lone pairs and we migrate the negative charge on the OTMAO to the dummy sites. The force field model developed here improves both structural and dynamical properties of aqueous TMAO solutions. Moreover, it reproduces the experimentally observed dependence of viscosity upon increasing TMAO concentration quantitatively. The simple procedure of the force field construction makes it easy to implement in molecular dynamics simulation packages and makes it compatible with the existing biomolecular force fields. This paves the path for further investigation of protein-TMAO interaction in aqueous solutions. [ABSTRACT FROM AUTHOR]

We investigate the dynamics of water in contact with solid calcium fluoride, where at low pH, localized charges can develop upon fluorite dissolution. We use 2D surface‐specific vibrational spectroscopy to quantify the heterogeneity of the interfacial water (D2O) molecules and provide information about the sub‐picosecond vibrational‐energy‐relaxation dynamics at the buried solid/liquid interface. We find that strongly H‐bonded OD groups, with a vibrational frequency below 2500 cm−1, display very rapid spectral diffusion and vibrational relaxation; for weakly H‐bonded OD groups, above 2500 cm−1, the dynamics slows down substantially. Atomistic simulations based on electronic‐structure theory reveal the molecular origin of energy transport through the local H‐bond network. We conclude that strongly oriented H‐bonded water molecules in the adsorbed layer, whose orientation is pinned by the localized charge defects, can exchange vibrational energy very rapidly due to the strong collective dipole, compensating for a partially missing solvation shell. [ABSTRACT FROM AUTHOR]

Copyright of Angewandte Chemie is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The thermochemical constants for the oxidation of tyrosine and tryptophan through proton coupled electron transfer in aqueous solution have been computed applying a recently developed density functional theory (DFT) based molecular dynamics method for reversible elimination of protons and electrons. This method enables us to estimate the solvation free energy of a proton (H+) in a periodic model system from the free energy for the deprotonation of an aqueous hydronium ion (H3O+). Using the computed solvation free energy of H+ as reference, the deprotonation and oxidation free energies of an aqueous species can be converted to pKa and normal hydrogen electrode (NHE) potentials. This conversion requires certain thermochemical corrections which were first presented in a similar study of the oxidation of hydrobenzoquinone [J. Cheng, M. Sulpizi, and M. Sprik, J. Chem. Phys. 131, 154504 (2009)]. Taking a different view of the thermodynamic status of the hydronium ion, these thermochemical corrections are revised in the present work. The key difference with the previous scheme is that the hydronium is now treated as an intermediate in the transfer of the proton from solution to the gas-phase. The accuracy of the method is assessed by a detailed comparison of the computed pKa, NHE potentials and dehydrogenation free energies to experiment. As a further application of the technique, we have analyzed the role of the solvent in the oxidation of tyrosine by the tryptophan radical. The free energy change computed for this hydrogen atom transfer reaction is very similar to the gas-phase value, in agreement with experiment. The molecular dynamics results however, show that the minimal solvent effect on the reaction free energy is accompanied by a significant reorganization of the solvent. [ABSTRACT FROM AUTHOR]

The density functional theory based molecular dynamics (DFTMD) method for the computation of redox free energies presented in previous publications and the more recent modification for computation of acidity constants are reviewed. The method uses a half reaction scheme based on reversible insertion/removal of electrons and protons. The proton insertion is assisted by restraining potentials acting as chaperones. The procedure for relating the calculated deprotonation free energies to Bro\nsted acidities (pKa) and the oxidation free energies to electrode potentials with respect to the normal hydrogen electrode is discussed in some detail. The method is validated in an application to the reduction of aqueous 1,4-benzoquinone. The conversion of hydroquinone to quinone can take place via a number of alternative pathways consisting of combinations of acid dissociations, oxidations, or dehydrogenations. The free energy changes of all elementary steps (ten in total) are computed. The accuracy of the calculations is assessed by comparing the energies of different pathways for the same reaction (Hess’s law) and by comparison to experiment. This two-sided test enables us to separate the errors related with the restrictions on length and time scales accessible to DFTMD from the errors introduced by the DFT approximation. It is found that the DFT approximation is the main source of error for oxidation free energies. [ABSTRACT FROM AUTHOR]

Photoemission and electrochemical data are used for showing the difference between the vertical and adiabatic ionization energy of aqueous hydroxide anion to be 2.9 eV. The solvent response to ionization is nonlinear and the studies have predicted a 4.1 eV difference between the energy for vertical attachment of an electron to the aqueous hydroxyl radical and the corresponding adiabatic electron affinity.