Céline Houriez

Prediction of nitroxide hyperfine coupling constants in solution from combined nanosecond scale simulations and quantum computations

We present a combined theoretical approach based on analyzing molecular dynamics trajectories (at the nanosecond scale) generated by use of classical polarizable force fields and on quantum calculations to compute averaged hyperfine coupling constants. That method is used to estimate the constant of a prototypical nitroxide: the dimethylnitroxide. The molecule is embedded during the simulations in a cubic box containing about 500 water molecules and the molecular dynamics is generated using periodic conditions. Once the trajectories are achieved, the nitroxide and its first hydration shell molecules are extracted, and the coupling constants are computed by considering the latter aggregates by means of quantum computations. However, all the water molecules of the bulk are also accounted for during those computations by means of the electrostatic potential fitted method. Our results exhibit that in order to predict accurate and reliable coupling constants, one needs to describe carefully the out-of-plane motion of the nitroxide nitrogen and to sample trajectories with a time interval of 400 fs at least to generate an uncorrelated large set of nitroxide structures. Compared to Car-Parrinello molecular dynamics techniques, our approach can be used readily to compute hyperfine coupling constants of large systems, such as nitroxides of great size interacting with macromolecules such as proteins or polymers.

Further insights into the environmental effects on the computed hyperfine coupling constants of nitroxides in aqueous solution.

We investigated the main two factors influencing the mean hyperfine coupling constants of small nitroxide radicals in aqueous solution, i.e., the out-of-plane displacement of their nitrogen atom and the environmental effects (solvent effects), by means of the approach we previously developed and fine-tuned to study the solvation of the dimethyl nitroxide radical. Our methodology efficiently combines classical molecular dynamics based on a polarizable force field at the nanosecond scale and quantum mechanics/molecular mechanics (QM/MM) computations to account for the bulk instantaneous electrostatic environmental effect. Our method has been applied to five small nitroxides, namely methyl nitroxide, ethyl nitroxide, dimethyl nitroxide, di-tert-butyl nitroxide, and PROXYL. The theoretical nitrogen hyperfine coupling constant values for the five nitroxides in solution are in good agreement with experiment (difference of 0.3 G on average). Our approach showed that the solvent shift in nitroxide hyperfine coupling constants is almost constant whatever the nitroxide, and, particularly, whatever the nitroxide NO moietys accessibility to the solvent. This result contrasts with earlier results derived from 10 ps scale trajectories based on Car-Parrinello molecular dynamics approach. However, we show that if we consider on average these latter results, they are in agreement with our conclusion. We also present an attempt to identify the origin of this result by analyzing the solvent contributions in terms of effects of the nitroxide first hydration shell and of the bulk, and by investigating the relation between these two contributions and the solvent structure at the vicinity of the NO moiety.

Structure and spectromagnetic properties of the superoxide radical adduct of DMPO in water: elucidation by theoretical investigations.

In the field of spin trapping chemistry, the design of more efficient radical traps can be assisted by the development of theoretical methods able to give a quantitative evaluation of the electron paramagnetic resonance (EPR) spectrum features of the spin-adduct radical, even before initiating the experimental work. The superoxide radical adduct of the 5,5-dimethyl-1-pyrroline-N-oxide nitrone (DMPO-OOH) has been reported in a huge number of papers devoted to the study of the oxidative stress. Here, we present for the first time the theoretical study of DMPO-OOH in an explicit water solution, based on the combined QM/MM//MD protocol we recently proposed, featuring a full coupling between the solute and all the explicit water molecules. Our results show that the DMPO-OOH EPR spectrum, whose interpretation is still debated, can be explained in the light of two sites in chemical exchange, in agreement with the most recent experimental data. Moreover, we demonstrate that each site consists of an equilibrium between the two main 5-membered ring conformations of DMPO-OOH. We provide also an analysis of the solvent contribution to the hyperfine coupling constants (hccs) as well as an exhaustive study of the possible relationship between the hccs and the main structural characteristics of DMPO-OOH. Our QM/MM//MD protocol appears thus to be an accurate theoretical tool allowing the investigation of the magnetic properties of large nitroxide spin adducts in complex environments.

Simulated Solvation of Organic Ions II: Study of Linear Alkylated Carboxylate Ions in Water Nanodrops and in Liquid Water. Propensity for Air/Water Interface and Convergence to Bulk Solvation Properties.

We investigated the solvation of carboxylate ions from formate to hexanoate, in droplets of 50 to 1000 water molecules and neat water, by computations using standard molecular dynamics and sophisticated polarizable models. The carboxylate ions from methanoate to hexanoate show strong propensity for the air/water interface in small droplets. Only the ions larger than propanoate retain propensity for the interface in larger droplets, where their enthalpic stabilization by ion/water dispersion is reduced there by 3 kcal mol(-1) per CH2 group. This is compensated by entropy effects over +3.3 cal mol(-1) K(-1) per CH2 group. On the surface, the anionic headgroups are strongly oriented toward the aqueous core, while the hydrophobic alkyl chains are repelled into air and lose their structure-making effects. These results reproduce the structure-making effects of alkyl groups in solution, and suggest that the hydrocarbon chains of ionic headgroups and alkyl substituents solvate independently. Extrapolation to bulk solution using standard extrapolation schemes yields absolute carboxylate solvation energies. The results for formate and acetate yield a proton solvation enthalpy of about 270 kcal mol(-1), close to the experiment-based value. The largest carboxylate ions yield a value smaller by about 10 kcal mol(-1), which requires studies in much larger droplets.

Simulated Solvation of Organic Ions: Protonated Methylamines in Water Nanodroplets. Convergence toward Bulk Properties and the Absolute Proton Solvation Enthalpy

We applied an alternative, purely theoretical route to estimate thermodynamical properties of organic ions in bulk solution. The method performs a large ensemble of simulations of ions solvated in water nanodroplets of different sizes, using a polarizable molecular dynamics approach. We consider protonated ammonia and methylamines, and K(+) for comparison, solvated in droplets of 50-1000 water molecules. The parameters of the model are assigned from high level quantum computations of small clusters. All the bulk phase results extrapolated from droplet simulations match, and confirm independently, the relative and absolute experiment-based ion solvation energies. Without using experiment-based parameters or assumptions, the results confirm independently the solvation enthalpy of the proton, as -270.3 ± 1.1 kcal mol(-1). The calculated relative solvation enthalpies of these ions are constant from small water clusters, where only the ionic headgroups are solvated, up to bulk solution. This agrees with experimental thermochemistry, that the relative solvation energies of alkylammonium ions by only four H2O molecules reproduce the relative bulk solvation energies, although the small clusters lack major bulk solvation factors. The droplet results also show a slow convergence of ion solvation properties toward their bulk limit, and predict that the stepwise solvation enthalpies of ion/water droplets are very close to those of pure neutral water droplets already after 50 water molecules. Both the ionic and neutral clusters approach the bulk condensation energy very gradually up to 10,000 water molecules, consistent with the macroscopic liquid drop model for pure water droplets. Compared to standard computational methods based on infinite periodic systems, our protocol represents a new purely theoretical approach to investigate the solvation properties of ions. It is applicable to the solvation of organic ions, which are pivotal in environmental, industrial, and biophysical chemistry but have been little investigated theoretically up to the present.

Structural and atoms-in-molecules analysis of hydrogen-bond network around nitroxides in liquid water.

In this study, we investigated the hydrogen-bond network patterns involving the NO moieties of five small nitroxides in liquid water by analyzing nanosecond scale molecular dynamics trajectories. To this end, we implemented two types of hydrogen-bond definitions, based on electronic structure, using Baders atoms-in-molecules analysis and based on geometric criteria. In each definition framework, the nitroxide/water hydrogen-bond networks appear very variable from a nitroxide to another. Moreover, each definition clearly leads to a different picture of nitroxide hydration. For instance, the electronic structure-based definition predicts a number of hydrogen bonds around the nitroxide NO moiety usually larger than geometric structure-based ones. One particularly interesting result is that the strength of a nitroxide/water hydrogen bond does not depend on its linearity, leading us to question the relevance of geometric definition based on angular cutoffs to study this type of hydrogen bond. Moreover, none of the hydrogen-bond definitions we consider in the present study is able to quantitatively correlate the strength of nitroxide/water hydrogen-bond networks with the aqueous nitroxide spin properties. This clearly exhibits that the hydrogen-bonding concept is not reliable enough to draw quantitative conclusions concerning such properties.

Extrapolating Single Organic Ion Solvation Thermochemistry from Simulated Water Nanodroplets

We compute the ion/water interaction energies of methylated ammonium cations and alkylated carboxylate anions solvated in large nanodroplets of 10 000 water molecules using 10 ns molecular dynamics simulations and an all-atom polarizable force-field approach. Together with our earlier results concerning the solvation of these organic ions in nanodroplets whose molecular sizes range from 50 to 1000, these new data allow us to discuss the reliability of extrapolating absolute single-ion bulk solvation energies from small ion/water droplets using common power-law functions of cluster size. We show that reliable estimates of these energies can be extrapolated from a small data set comprising the results of three droplets whose sizes are between 100 and 1000 using a basic power-law function of droplet size. This agrees with an earlier conclusion drawn from a model built within the mean spherical framework and paves the road toward a theoretical protocol to systematically compute the solvation energies of complex organic ions.

Physisorbed H2@Cu(100) surface: Potential and spectroscopy

Using an embedding approach, a 2-D potential energy function has been calculated to describe the physisorption interaction of H2 with a Cu(100) surface. For this purpose, a cluster model of the system calculated with highly correlated wavefunctions is combined with a periodic Density-Functional-Theory method using van der Waals-DF2 functional. Rotational and vibrational energy levels of physisorbed H2, as well as D2 and HD, are calculated using the 2D embedding corrected potential energy function. The calculated transitions are in a very good agreement with Electron-Energy-Loss-Spectroscopy observations.

Organic ion association in aqueous phase and ab initio-based force fields: The case of carboxylate/ammonium salts

We performed molecular dynamics simulations of carboxylate/methylated ammonium ion pairs solvated in bulk water and of carboxylate/methylated ammonium salt solutions at ambient conditions using an ab initio-based polarizable force field whose parameters are assigned to reproduce only high end quantum computations, at the Møller-Plesset second-order perturbation theory/complete basis set limit level, regarding single ions and ion pairs as isolated and micro-hydrated in gas phase. Our results agree with the available experimental results regarding carboxylate/ammonium salt solutions. For instance, our force field approach predicts the percentage of acetate associated with ammonium ions in CH3COO-/CH3NH3+ solutions at the 0.2-0.8M concentration scale to range from 14% to 35%, in line with the estimates computed from the experimental ion association constant in liquid water. Moreover our simulations predict the number of water molecules released from the ion first hydration shell to the bulk upon ion association to be about 2.0 ± 0.6 molecules for acetate/protonated amine ion pairs, 3.1 ± 1.5 molecules for the HCOO-/NH4+ pair and 3.3 ± 1.2 molecules for the CH3COO-/(CH3)4N+ pair. For protonated amine-based ion pairs, these values are in line with experiment for alkali/halide pairs solvated in bulk water. All these results demonstrate the promising feature of ab initio-based force fields, i.e., their capacity in accurately modeling chemical systems that cannot be readily investigated using available experimental techniques.

Solvation of the Guanidinium Ion in Pure Aqueous Environments: A Theoretical Study from an “Ab Initio”-Based Polarizable Force Field

We report simulation results regarding the hydration process of the guanidinium cation in water droplets and in bulk liquid water, at a low concentration of 0.03 M, performed using a polarizable approach to model both water/water and ion/water interactions. In line with earlier theoretical studies, our simulations show a preferential orientation of guanidinium at water-vacuum interfaces, i.e., a parallel orientation of the guanidinium plane to the aqueous surface. In an apparent contradiction with earlier simulation studies, we show also that guanidinium has a stronger propensity for the cores of aqueous systems than the ammonium cation. However, our bulk simulation conditions correspond to weaker cation concentrations than in earlier studies, by 2 orders of magnitude, and that the same simulations performed using a standard nonpolarizable force field leads to the same conclusion. From droplet data, we extrapolate the guanidinium single hydration enthalpy value to be -82.9 ± 2.2 kcal mol-1. That is about half as large as the sole experimental estimate reported to date, about -144 kcal mol-1. Our result yields a guanidinium absolute bulk hydration free energy at ambiant conditions to be -78.4 ± 2.6 kcal mol-1, a value smaller by 3 kcal mol-1 compared to ammonium. The relatively large magnitude of our guanidinium hydration free energy estimate suggests the Gdm+ protein denaturing properties to result from a competition between the cation hydration effects and the cation/protein interactions, a competition that can be modulated by weak differences in the protein or in the cation chemical environment.