Wrapping Up Hydrophobic Hydration: Locality Matters.

Water, being the universal solvent, acts as a competing agent in fundamental processes, such as folding, aggregation or biomolecular recognition. A molecular understanding of hydrophobic hydration is of central importance to understanding the subtle free energy differences, which dictate function. Ab initio and classical molecular dynamics simulations yield two distinct hydration water populations in the hydration shell of solvated tert-butanol noted as "HB-wrap" and "HB-hydration2bulk". The experimentally observed hydration water spectrum can be dissected into two modes, centered at 164 and 195 cm-1. By comparison to the simulations, these two bands are attributed to the "HB-wrap" and "HB-hydration2bulk" populations, respectively. We derive a quantitative correlation between the population in each of these two local water coordination motifs and the temperature dependence of the solvation entropy. The crossover from entropy to enthalpy dominated solvation at elevated temperatures, as predicted by theory and observed experimentally, can be rationalized in terms of the distinct temperature stability and thermodynamic signatures of "HB-wrap" and "HB-hydration2bulk".

Graph theory for automatic structural recognition in molecular dynamics simulations.

Graph theory algorithms have been proposed in order to identify, follow in time, and statistically analyze the changes in conformations that occur along molecular dynamics (MD) simulations. The atomistic granularity level of the MD simulations is maintained within the graph theoric algorithms proposed here, isomorphism is a key component together with keeping the chemical nature of the atoms. Isomorphism is used to recognize conformations and construct the graphs of transitions, and the reduction in complexity of the isomorphism has been achieved by the introduction of "orbits" and "reference snapshots." The proposed algorithms are applied to MD trajectories of gas phase molecules and clusters as well as condensed matter. The changes in conformations followed over time are hydrogen bond(s), proton transfer(s), coordination number(s), covalent bond(s), multiple fragmentation(s), and H-bonded membered rings. The algorithms provide an automatic analysis of multiple trajectories in parallel, and can be applied to ab initio and classical MD trajectories alike, and to more coarse grain representations.

Infrared predissociation vibrational spectroscopy of Li⁺(H2O)(3-4)Ar(0,1) reanalyzed using density functional theory molecular dynamics.

The experimental IR-PD (infrared predissociation) spectra of Li(+)(H2O)(3-4)Ar and Li(+)(H2O)(3-4) clusters, monitoring two different loss channels and thus different temperatures, have been reanalyzed using DFT-MD (density functional theory based molecular dynamics) simulations for finite temperature and anharmonic theoretical spectroscopy. The use of DFT-MD to calculate IR-PD spectra at low and elevated temperatures was found remarkably accurate and useful in precise structural characterization. The dynamical spectra have in particular provided the opportunity to estimate the clusters temperatures in the IR-PD experiments. The temperatures for Li(+)(H2O)(3-4)Ar are estimated at 50-60 K whereas Li(+)(H2O)3 and Li(+)(H2O)4 have been estimated at around 500-600 and 400 K, respectively.

Modeling proton-induced damage on 2-deoxy-D-ribose. Conformational analysis.

Modeling proton-induced damage in biological systems, in particular in DNA building blocks, is of major concern in studies on cancer proton therapy. This is indeed an extremely complex process and analysis of the mechanism at the molecular level is of crucial interest. Such collision reactions of protons on biological targets induce different reactions: excitation and ionization of the biomolecule, fragmentation of the ionized species, and charge transfer from the projectile ion toward the biomolecular target. In order to have an insight into such mechanisms, we have performed a theoretical approach of two of the most important steps, the fragmentation and the charge transfer processes. For that purpose, we have considered collision of protons with isolated 2-deoxy-D-ribose by means of ab-initio molecular dynamics and quantum chemistry molecular methods. The conformation of the sugar moiety has been analyzed and appears to induce important effects, in particular different fragmentation patterns have been pointed out with regard to the conformation, and significant variations of the charge transfer cross sections have been exhibited.

High energy conformers of M(+)(APE)(H2O)(0-1)Ar(0-1) clusters revealed by combined IR-PD and DFT-MD anharmonic vibrational spectroscopy.

IR-PD vibrational spectroscopy and DFT-based molecular dynamics simulations are combined in order to unravel the structures of M(+)(APE)(H2O)0-1 ionic clusters (M = Na, K), where APE (2-amino-1-phenyl ethanol) is commonly used as an analogue for the noradrenaline neurotransmitter. The strength of the synergy between experiments and simulations presented here is that DFT-MD provides anharmonic vibrational spectra that unambiguously help assign the ionic clusters structures. Depending on the interacting cation, we have found that the lowest energy conformers of K(+)(APE)(H2O)0-1 clusters are formed, while the lowest energy conformers of Na(+)(APE)(H2O)0-1 clusters can only be observed through water loss channel (i.e. without argon tagged to the clusters). Trapping of higher energy conformers is observed when the argon loss channel is recorded in the experiment. This has been rationalized by transition state energies. The dynamical anharmonic vibrational spectra unambiguously provide the prominent OH stretch due to the OH···NH2 H-bond, within 10 cm(-1) of the experiment, hence reproducing the 240-300 cm(-1) red-shift (depending on the interacting cation) from bare neutral APE. When this H-bond is not present, the dynamical anharmonic spectra provide the water O-H stretches as well as the rotational motion of the water molecule at finite temperature, as observed in the experiment.

Formation of the OOH(·) radical at steps of the boehmite surface and its inhibition by gallic acid: A theoretical study including DFT-based dynamics.

The ability of boehmite surface to stabilize reactive oxygen species (ROS) was investigated with density functional theory (DFT) calculations. We provide evidence that the O2(·-)radical is stabilized at the AlOOH boehmite terrace defect-free surface, and that step favors the formation of the OOH(·) radical without activation energy required. These tendencies are confirmed when considering the explicit presence of water solvent. We propose that gallic acid (GA), a non-expensive, non-toxic, natural product and radical scavenger, may form a full layer on the (101) stepped boehmite surface. The quenching of O2- to OOH(·) reaction on the GA-functionalized boehmite surface is investigated. It is shown that gallic acid passivates the surface and that the formation of the OOH(·) radical is consequently inhibited. The interaction of the GA-functionalized boehmite surface with water mimicking the body fluid is also investigated by means of DFT-based molecular dynamics simulations at room temperature. Hydrogen bonds between the functionalized GA layer and the immediate interfacial water are characterized.

Ultrafast nonadiabatic fragmentation dynamics of doubly charged uracil in a gas phase.

A combination of time-dependent density functional theory and Born-Oppenheimer molecular dynamics methods is used to investigate fragmentation of doubly charged gas-phase uracil in collisions with 100 keV protons. The results are in good agreement with ion-ion coincidence measurements. Orbitals of similar energy and/or localized in similar bonds lead to very different fragmentation patterns, thus showing the importance of intramolecular chemical environment. In general, the observed fragments do not correspond to the energetically most favorable dissociation path, which is due to dynamical effects occurring in the first few femtoseconds after electron removal.

Theoretical investigation of the ultrafast dissociation of ionised biomolecules immersed in water: direct and indirect effects.

Theoretical simulations are particularly well suited to investigate, at a molecular level, direct and indirect effects of ionising radiations in DNA, as in the particular case of irradiation by swift heavy ions such as those used in hadron therapy. In the past recent years, we have developed the modeling at the microscopic level of the early stages of the Coulomb explosion of DNA molecules immersed in liquid water that follows the irradiation by swift heavy ions. To that end, Time-Dependent Density Functional Theory molecular dynamics simulations (TD-DFT MD) have been developed where localised Wannier orbitals are propagated. This latter enables to separate molecular orbitals of each water molecule from the molecular orbitals of the biomolecule. Our main objective is to demonstrate that the double ionisation of one molecule of the liquid sample, either one water molecule from the solvent or the biomolecule, may be in some cases responsible for the formation of an atomic oxygen as a direct consequence of the molecule Coulomb explosion. Our hypothesis is that the molecular double ionisation arising from irradiation by swift heavy ions (about 10% of ionisation events by ions whose velocity is about the third of speed of light), as a primary event, though maybe less probable than other events resulting from the electronic cascading (for instance, electronic excitations, electron attachments), may be systematically more damageable (and more lethal), as supported by experiments that have been carried out in our group in the 1990s (in studies of damages created by K holes in DNA). The chemical reactivity of the produced atomic oxygen with other radicals present in the medium will ultimately lead to chemical products that are harmful to DNA. In the present paper, we review our theoretical methodology in an attempt that the community be familiar with our new approach. Results on the production of atomic oxygen as a result of the double ionisation of water or as a result of the double ionisation of the Uracil RNA base will be presented.

Unravelling the conformational dynamics of the aqueous alanine dipeptide with first-principle molecular dynamics.

First principle DFT-based molecular dynamics simulations are performed in order to bring new insights into the conformational dynamics of the alanine dipeptide analogue Ac-Ala-NHMe (with methyl group caps at the extremities) immersed in liquid water at ambient temperature. Two simulations have been run for a total of 100 ps, which allows for a relevant statistical sampling of the phase space, at the ab initio level. PII-beta equilibrium and (PII-beta)-alphaR conformational interconversions are obtained, without a free-energy barrier for the Phi angle and with a rather low barrier of 2-3 kcal/mol for the Psi angle, easily overcome from solute-solvent energy transfers. We furthermore give first insights into the rather strong zwitterionic character of the peptide bonds of the dipeptide when immersed in the liquid. The structural and zwitterionic properties extracted from first-principle dynamics in the liquid phase will be useful as benchmarks for force field developments.

Vibrational Spectra of Small Protonated Peptides from Finite Temperature MD Simulations and IRMPD Spectroscopy.

Finite temperature Born-Oppenheimer DFT-based molecular dynamics simulations are presented for the vibrational spectroscopy of the prototype gas-phase Ala2H(+) and Ala3H(+) protonated peptides. The dynamics and the vibrational signatures are used to interpret IR-MPD spectra recorded in the NH/OH stretch region. Molecular dynamics simulations are one way to go beyond the harmonic approximations commonly applied for the calculations of infrared spectra, naturally including all anharmonicities, i.e. mode couplings, vibrational and dipole anharmonicities. The dynamics of the peptides allows understanding of the evolution of the shape and width of the N-H bands when increasing the size of the peptide, as demonstrated here with the two small prototypes Ala2H(+) and Ala3H(+). Hence, the conformational dynamics of Ala2(+) at room temperature participates to the broadening of the IR active bands. The complex N-H broadband of Ala3H(+) is shown to result from the dynamics of the N-H groups in the different peptide families, with a special role from breaking/reforming of hydrogen bonds involving the N-H groups. Taking this dynamics into account is thus mandatory for the understanding of this band in the 300-400 K experimental spectrum.

Molecular Dynamics and Room Temperature Vibrational Properties of Deprotonated Phosphorylated Serine.

The local structure of phosphorylated residues in peptides and proteins may have a decisive role on their functional properties. Recent IRMPD experiments have started to provide spectroscopic signatures of such structural details; however, a proper modeling of these signatures beyond the harmonic approximation, taking into account temperature and entropic effects, is still lacking. In order to bridge this gap, DFT-based Car-Parrinello molecular dynamics simulations have been carried out for the first time on a phosphorylated amino acid, gaseous deprotonated phosphoserine. It is found that all vibrational signatures are successfully reproduced, and new deconvolution techniques enable the assignment of the vibrational spectrum directly from the dynamics results and the comparison of vibrational modes at several temperatures. The lowest energy structure is found to involve a strong hydrogen bond between the deprotonated phosphate and the acid with relatively small free energy barriers to proton transfer; however, we find that proton shuttling between the two sites does not occur frequently. Anharmonicities turn out to be important to reproduce the frequencies and shapes of several experimental bands. Comparison of room temperature and 13 K, effectively harmonic dynamics, allows insight to be obtained into vibrational anharmonicities. In particular, a significant blue-shift and broadening of the C═O stretching frequency from 13 to 300 K can be ascribed to intrinsic anharmonicity rather than to anharmonic coupling to other modes. On the other hand, significant couplings are found for the stretching motions of the hydrogen bonded P-O bond and of the free P-OH bond, mainly with modes within the phosphate group.

Alanine polypeptide structural fingerprints at room temperature: what can be gained from non-harmonic Car-Parrinello molecular dynamics simulations.

Structural infrared fingerprints of neutral gas phase alanine peptides of increasing size and complexity (dipeptide, octapeptide, and beta-strand peptide) are characterized through DFT-based Car-Parrinello molecular dynamics simulations. Harmonic and nonharmonic vibrational signatures are calculated from the time correlation of the dipole moment of the gas phase peptide in a direct way (without any approximation) respectively from low temperature (20 K) and room temperature (300 K) molecular dynamics. Our main purpose is to answer the two following questions: (i) Is the direct inclusion of temperature for the calculation of infrared spectra mandatory for the comprehension of the vibrational signatures experimentally recorded at room temperature? (ii) To what extent is the amide I, II, and III domain sensitive enough to the local structure of the peptides, to provide vibrational signatures that can be definitely used to assess the peptide conformation at 300 K?

Resonant infrared multiphoton dissociation spectroscopy of gas-phase protonated peptides. Experiments and Car-Parrinello dynamics at 300 K.

The gas-phase structures of protonated peptides are studied by means of resonant infrared multiphoton dissociation spectroscopy (R-IRMPD) performed with a free electron laser. The peptide structures and protonation sites are obtained through comparison between experimental IR spectra and their prediction from quantum chemistry calculations. Two different analyses are conducted. It is first supposed that only well-defined conformations, sufficiently populated according to a Boltzmann distribution, contribute to the observed spectra. On the contrary, DFT-based Car-Parrinello molecular dynamics simulations show that at 300 K protonated peptides no longer possess well-defined structures, but rather dynamically explore the set of conformations considered in the first conventional approach.

Extracting effective normal modes from equilibrium dynamics at finite temperature.

A general method for obtaining effective normal modes of a molecular system from molecular dynamics simulations is presented. The method is based on a localization criterion for the Fourier transformed velocity time-correlation functions of the effective modes. For a given choice of the localization function used, the method becomes equivalent to the principal mode analysis (PMA) based on covariance matrix diagonalization. On the other hand, a proper choice of the localization function leads to a novel method with a strong analogy with the usual normal mode analysis of equilibrium structures, where the Hessian system at the minimum energy structure is replaced by the thermal averaged Hessian, although the Hessian itself is never actually calculated. This method does not introduce any extra numerical cost during the simulation and bears the same simplicity as PMA itself. It can thus be readily applied to ab initio molecular dynamics simulations. Three such examples are provided here. First we recover effective normal modes of an isolated formaldehyde molecule computed at 20 K in very good agreement with the results of a normal mode analysis performed at its equilibrium structure. We then illustrate the applicability of the method for liquid phase studies. The effective normal modes of a water molecule in liquid water and of a uracil molecule in aqueous solution can be extracted from ab initio molecular dynamics simulations of these two systems at 300 K.

Proton transfers induced by lead(II) in a uracil nucleobase: a study based on quantum chemistry calculations.

Within the context of metal biotoxicity, electrospray ionization mass spectrometry experiments (ESIMS) have recently been performed by us on the pyrimidine nucleobases (B) uracil and thymine complexed with lead(II) [Int. J. Mass. Spectrom. 2005, 243, 279]. Among the ions detected, [Pb(B)-H]+ complexes, where the base has been deprotonated, have been identified as producing intense signals. In the same study, quantum calculations based on density functional theory (DFT) have assessed the complexation sites and energies of [Pb(B)-H]+ ions. The present DFT investigations aim at giving an understanding on the energetics and mechanisms associated with uracil's loss of a proton. We specifically assess and quantify the role of lead binding in this process. For that purpose, intra- and intermolecular proton transfers have been considered. We have found that uracil (U) 1,3-tautomerization can be exergonic when uracil is complexed with Pb2+, in opposition to the situation without lead. The corresponding intramolecular processes were nonetheless found to occur at geological time scales. In contrast, the addition of a second body to [Pb(U)]2+ complexes, namely OH- or H2O (as found in the initial water droplet of ESIMS experiments), gives exergonic and fast uracil 1,3-proton transfers. Finally, we have shown that intermolecular proton transfers in uracil-H2O, uracil-OH-, or uracil-uracil complexes are able to explain the experimentally detected [Pb(U)-H]+ ions.

An interdisciplinary approach to investigate the impact of cobalt in a human keratinocyte cell line.

Since in nuclear power plants, risks of skin contact contamination by radiocobalt are significant, we focused on the impact of cobalt on a human cutaneous cell line, i.e. HaCaT keratinocytes. The present paper reports an interdisciplinary approach aimed at clarifying the biochemical mechanisms of metabolism and toxicity of cobalt in HaCaT cells. Firstly, a brief overview of the used instrumental techniques is reported. The following parts present description and discussion of results concerning: (i) toxicological studies concerning cobalt impact towards HaCaT cells (ii) structural and speciation fundamental studies of cobalt-bioligand systems, through X-ray absorption spectroscopy (XAS), ab initio and thermodynamic modelling (iii) preliminary results regarding intracellular cobalt speciation in HaCaT cells using size exclusion chromatography/inductively coupled plasma-atomic emission spectroscopy (SEC/ICP-AES) and direct in situ analysis by ion beam micropobe analytical techniques.

Ab initio molecular dynamics of protonated dialanine and comparison to infrared multiphoton dissociation experiments.

Finite temperature Car-Parrinello molecular dynamics simulations are performed for the protonated dialanine peptide in vacuo, in relation to infrared multiphoton dissociation experiments. The simulations emphasize the flexibility of the different torsional angles at room temperature and the dynamical exchange between different conformers which were previously identified as stable at 0 K. A proton transfer occurring spontaneously at the N-terminal side is also observed and characterized. The theoretical infrared absorption spectrum is computed from the dipole time correlation function, and, in contrast to traditional static electronic structure calculations, it accounts directly for anharmonic and finite temperature effects. The comparison to the experimental infrared multiphoton dissociation spectrum turns out very good in terms of both band positions and band shapes. It does help the identification of a predominant conformer and the attribution of the different bands. The synergy shown between the experimental and theoretical approaches opens the door to the study of the vibrational properties of complex and floppy biomolecules in the gas phase at finite temperature.

Infrared Spectroscopy of N-Methylacetamide Revisited by ab Initio Molecular Dynamics Simulations.

The density functional theory based molecular dynamics simulation method ("Car-Parrinello") was applied in a numerical study of the electronic properties, hydrogen bonding, and infrared spectroscopy of the trans and cis isomer of N-methylacetamide in aqueous solution. A detailed analysis of the electronic structure of the solvated molecules, in terms of localized Wannier functions and Born atomic charges, is presented. Two schemes for the computation of the solute infrared absorption spectrum are investigated: In the first method the spectrum is determined by Fourier transforming the time correlation function of the solute dipole as determined from the Wannier function analysis. The second method uses instead the molecular current-current correlation function computed from the Born charges and atomic velocities. The resulting spectral properties of trans- and cis-NMA are carefully compared to each other and to experimental results. We find that the two solvated isomers can be clearly distinguished by their infrared spectral profile in the 1000-2000 cm(-)(1) range.