Ab Initio Molecular Dynamics Approach to Quantitative Description of Ion Pairing in Water

Ion pairing of lithium and fluoride in water is described quantitatively using ab initio molecular dynamics simulations. We design a reliable computational protocol for evaluating the ion-ion potential of mean force using density functional based simulation methods. By comparison to classical molecular dynamics with empirical force fields we establish the statistical error of the procedure. We also check the accuracy of the electronic structure description by comparison to experimental structural data and to higher level calculations for model systems. The present approach not only points to deficiencies in force field calculations of potentials of mean force for difficult cases of high charge density ions like the aqueous lithium fluoride pair, but also allows extracting electronic information, such as the amount of charge transfer to solvent and its dependence on the ion-ion distance

Peptide salt bridge stability: From gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipep- tide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desol- vated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hy- drate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within similar to 30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O -> OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM/MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with similar to 40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O((aq))(+) species. (C) 2011 American Institute of Physics. doi:10.1063/1.3664746

DNA Lesion Can Facilitate Base Ionization Vertical Ionization Energies of Aqueous 8 Oxoguanine and its Nucleoside and Nucleotide

8 Oxoguanine is one of the key products of indirect radiation damage to DNA by reactive oxygen species. Here, we describe ionization of this damaged nucleobase and the corresponding nucleoside and nucleotide in aqueous phase, modeled by the nonequilibrium polarizable continuum model, establishing their lowest vertical ionization energies of 6.8 7.0 eV. We thus confirm that 8 oxoguanine has even lower ionization energy than the parental guanine, which is the canonical nucleobase with the lowest ionization energy. Therefore, it can act as a trap for the cationic hole formed by ionizing radiation and thus protect DNA from further radiation damage. We also model using time dependent density functional theory and measure by liquid jet photoelectron spectroscopy the valence photoelectron spectrum of 8 oxoguanine in water. We show that the calculated higher lying ionization states match well the experiment which, however, is not sensitive enough to capture the electron signal corresponding to the lowest ionization process due to the low solubility of 8 oxoguanine in wate

Aqueous Cation-Amide Binding: Free Energies and IR Spectral Signatures by Ab Initio Molecular Dynamics

Understanding specific ion effects on proteins remains a considerable challenge. N-methylacetamide serves as a useful proxy for the protein backbone that can be well characterized both experimentally and theoretically. The spectroscopic signatures in the amide I band reflecting the strength of the interaction of alkali cations and alkali earth dications with the carbonyl group remain difficult to assign and controversial to interpret. Herein, we directly compute the IR shifts corresponding to the binding of either sodium or calcium to aqueous N-methylacetamide using ab initio molecular dynamics simulations. We show that the two cations interact with aqueous N-methylacetamide with different affinities and in different geometries. Since sodium exhibits a weak interaction with the carbonyl group, the resulting amide I band is similar to an unperturbed carbonyl group undergoing aqueous solvation. In contrast, the stronger calcium binding results in a clear IR shift with respect to N-methylacetamide in pure water

Transforming Anion Instability into Stability Contrasting Photoionization of Three Protonation Forms of the Phosphate Ion upon Moving into Water

We use photoelectron emission spectroscopy with vacuum microjet technique and quantum chemistry calculations to investigate electronic structure and stability of aqueous phosphate anions. On the basis of the measured photoelectron spectra of sodium phosphates at different pH, we report the lowest vertical ionization energies of monobasic 9.5 eV , dibasic 8.9 eV , and tribasic 8.4 eV anions. Electron binding energies were in tandem modeled with ab initio methods, using a mixed dielectric solvation model together with up to 64 explicitly solvating water molecules. We demonstrate that two solvation layers of explicit water molecules are needed to obtain converged values of vertical ionization energies VIEs within this mixed solvation model, leading to very good agreement with experiment. We also show that the highly charged PO43 anion, which is electronically unstable in the gas phase, gains the electronic stability with about 16 water molecules, while only 2 3 water molecules are sufficient to stabilize the doubly charged phosphate anion. We also investigate the effect of ion pairing on the vertical ionization energy. In contrast to protonation leading to a formation of covalent O H bond , sodiation leading to an anion Na ion pair has only a weak effect on the electron binding energ