Ultrafast phosphate hydration dynamics in bulk H2O

Impact of molecular symmetry on crystallization pathways in highly supersaturated KH2PO4 solutions

Solute structure and its evolution in supersaturated aqueous solutions are key clues to understand Ostwald's step rule. Here, we measure the structural evolution of solute molecules in highly supersaturated solutions of KH2PO4 (KDP) and NH4H2PO4 (ADP) using a combination of electrostatic levitation and synchrotron X-ray scattering. The measurement reveals the existence of a solution-solution transition in KDP solution, caused by changing molecular symmetries and structural evolution of the solution with supersaturation. Moreover, we find that the molecular symmetry of H2PO4- impacts on phase selection. These findings manifest that molecular symmetry and its structural evolution can govern the crystallization pathways in aqueous solutions, explaining the microscopic origin of Ostwald's step rule.

Dioxygen and glucose force motion of the electron-transfer switch in the iron(III) flavohemoglobin-type nitric oxide dioxygenase

Kinetic and structural investigations of the flavohemoglobin-type NO dioxygenase have suggested critical roles for transient Fe(III)O2 complex formation and O2-forced movements affecting hydride transfer to the FAD cofactor and electron-transfer to the Fe(III)O2 complex. Stark-effect theory together with structural models and dipole and internal electrostatic field determinations provided a semi-quantitative spectroscopic method for investigating the proposed Fe(III)O2 complex and O2-forced movements. Deoxygenation of the enzyme causes Stark effects on the ferric heme Soret and charge-transfer bands revealing the Fe(III)O2 complex. Deoxygenation also elicits Stark effects on the FAD that expose forces and motions that create a more restricted NADH access to FAD for hydride transfer and switch electron-transfer off. Glucose also forces the enzyme toward an off state. Amino acid substitutions at the B10, E7, E11, G8, D5, and F7 positions influence the Stark effects of O2 on resting heme spin states and FAD consistent with the proposed roles of the side chains in the enzyme mechanism. Deoxygenation of ferric myoglobin and hemoglobin A also induces Stark effects on the hemes suggesting a common 'oxy-met' state. The ferric myoglobin and hemoglobin heme spectra are also glucose-responsive. A conserved glucose or glucose-6-phosphate binding site is found bridging the BC-corner and G-helix in flavohemoglobin and myoglobin suggesting novel allosteric effector roles for glucose or glucose-6-phosphate in the NO dioxygenase and O2 storage functions. The results support the proposed roles of a ferric O2 intermediate and protein motions in regulating electron-transfer during NO dioxygenase turnover.

Physical insight into the entropy-driven ion association

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na⁺ , Cs⁺ , Ca²⁺ , F^- , I^- , CO₃ ^2- , and HPO₄ ^2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

Bursting the bubble: A molecular understanding of surfactant-water interfaces

Surfactant science has historically emphasized bulk, thermodynamic measurements to understand the microemulsion properties of greatest industrial significance, such as interfacial tensions, phase behavior, and thermal stability. Recently, interest in the molecular properties of surfactants has grown among the physical chemistry community. This has led to the application of cutting-edge spectroscopic methods and advanced simulations to understand the specific interactions that give rise to the previously studied bulk characteristics. In this Perspective, we catalog key findings that describe the surfactant-oil and surfactant-water interfaces in molecular detail. We emphasize the role of ultrafast spectroscopic methods, including two-dimensional infrared spectroscopy and sum-frequency-generation spectroscopy, in conjunction with molecular dynamics simulations, and the role these techniques have played in advancing our understanding of interfacial properties in surfactant microemulsions.

Aqueous Contact Ion Pairs of Phosphate Groups with Na+, Ca2+ and Mg2+ – Structural Discrimination by Femtosecond Infrared Spectroscopy and Molecular Dynamics Simulations

The extent of contact and solvent shared ion pairs of phosphate groups with Na+, Ca2+ and Mg2+ ions in aqueous environment and their relevance for the stability of polyanionic DNA and RNA structures is highly debated. Employing the asymmetric phosphate stretching vibration of dimethyl phosphate (DMP), a model system of the sugar-phosphate backbone of DNA and RNA, we present linear infrared, femtosecond infrared pump-probe and absorptive 2D-IR spectra that report on contact ion pair formation via the presence of blue shifted spectral signatures. Compared to the linear infrared spectra, the nonlinear spectra reveal contact ion pairs with increased sensitivity because the spectra accentuate differences in peak frequency, transition dipole moment strength, and excited state lifetime. The experimental results are corroborated by long time scale MD simulations, benchmarked by density functional simulations on phosphate-ion-water clusters. The microscopic interpretation reveals subtle structural differences of ion pairs formed by the phosphate group and the ions Na+, Ca2+ and Mg2+. Intricate properties of the solvation shell around the phosphate group and the ion are essential to explain the experimental observations. The present work addresses a challenging to probe topic with the help of a model system and establishes new experimental data of contact ion pair formation, thereby underlining the potential of nonlinear 2D-IR spectroscopy as an analytical probe of phosphate-ion interactions in complex biological systems.

Water Librations in the Hydration Shell of Phospholipids

The hydrophilic phosphate moiety in the headgroup of phospholipids forms strong hydrogen bonds with water molecules in the first hydration layer. Time-domain terahertz spectroscopy in a range from 100 to 1000 cm^-1 reveals the influence of such interactions on rotations of water molecules. We determine librational absorption spectra of water nanopools in phospholipid reverse micelles for a range from w₀ = 2 to 16 waters per phospholipid molecule. A pronounced absorption feature with maximum at 830 cm^-1 is superimposed on a broad absorption band between 300 and 1000 cm^-1. Molecular dynamics simulations of water in the reverse micelles suggest that the feature at 830 cm^-1 arises from water molecules forming one or two strong hydrogen bonds with phosphate groups, while the broad component comes from bulk-like environments. This behavior is markedly different from water interacting with less polar surfaces.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Molecular couplings and energy exchange between DNA and water mapped by femtosecond infrared spectroscopy of backbone vibrations

Molecular couplings between DNA and water together with the accompanying processes of energy exchange are mapped via the ultrafast response of DNA backbone vibrations after OH stretch excitation of the water shell. Native salmon testes DNA is studied in femtosecond pump-probe experiments under conditions of full hydration and at a reduced hydration level with two water layers around the double helix. Independent of their local hydration patterns, all backbone vibrations in the frequency range from 940 to 1120 cm^-1 display a quasi-instantaneous reshaping of the spectral envelopes of their fundamental absorption bands upon excitation of the water shell. The subsequent reshaping kinetics encompass a one-picosecond component, reflecting the formation of a hot ground state of the water shell, and a slower contribution on a time scale of tens of picoseconds. Such results are benchmarked by measurements with resonant excitation of the backbone modes, resulting in distinctly different absorption changes. We assign the fast changes of DNA absorption after OH stretch excitation to structural changes in the water shell which couple to DNA through the local electric fields. The second slower process is attributed to a flow of excess energy from the water shell into DNA, establishing a common heated ground state in the molecular ensemble. This interpretation is supported by theoretical calculations of the electric fields exerted by the water shell at different temperatures.

On the Structural and Dynamical Properties of DOPC Reverse Micelles

The structure and dynamics of phospholipid reverse micelles are studied by molecular dynamics. We report all-atom unconstrained simulations of 1,2-dioleoyl-sn-phosphatidylcholine (DOPC) reverse micelles in benzene of increasing sizes, with water-to-surfactant number ratios ranging from W₀ = 1 to 16. The aggregation number, i.e., the number of DOPC molecules per reverse micelle, is determined to fit experimental light-scattering measurements of the reverse micelle diameter. The simulated reverse micelles are found to be approximately spherical. Larger reverse micelles (W₀ > 4) exhibit a layered structure with a water core and the hydration structure of DOPC phosphate head groups is similar to that found in phospholipid membranes. In contrast, the structure of smaller reverse micelles (W₀ ≤ 4) cannot be described as a series of concentric layers successively containing water, surfactant head groups, and surfactant tails, and the head groups are only partly hydrated and frequently present in the core. The dynamics of water molecules within the phospholipid reverse micelles slow down as the reverse micelle size decreases, in agreement with prior studies on AOT and Igepal reverse micelles. However, the average water reorientation dynamics in DOPC reverse micelles is found to be much slower than in AOT and Igepal reverse micelles with the same W₀ ratio. This is explained by the smaller water pool and by the stronger interactions between water and the charged head groups, as confirmed by the red-shift of the computed infrared line shape with decreasing W₀.

Range, Magnitude, and Ultrafast Dynamics of Electric Fields at the Hydrated DNA Surface

Range and magnitude of electric fields at biomolecular interfaces and their fluctuations in a time window down to the subpicosecond regime have remained controversial, calling for electric-field mapping in space and time. Here, we trace fluctuating electric fields at the surface of native salmon DNA via their interactions with backbone vibrations in a wide range of hydration levels by building the water shell layer by layer. Femtosecond two-dimensional infrared spectroscopy and ab initio based theory establish water molecules in the first two layers as the predominant source of interfacial electric fields, which fluctuate on a 300 fs time scale with an amplitude of 25 MV/cm due to thermally excited water motions. The observed subnanometer range of these electric interactions is decisive for biochemical structure and function.

Ultrafast vibrational dynamics of the DNA backbone at different hydration levels mapped by two-dimensional infrared spectroscopy

DNA oligomers are studied at 0% and 92% relative humidity, corresponding to N < 2 and N > 20 water molecules per base pair. Two-dimensional (2D) infrared spectroscopy of DNA backbone modes between 920 and 1120 cm(-1) maps fluctuating interactions at the DNA surface. At both hydration levels, a frequency fluctuation correlation function with a 300 fs decay and a slow decay beyond 10 ps is derived from the 2D lineshapes. The fast component reflects motions of DNA helix, counterions, and water shell. Its higher amplitude at high hydration level reveals a significant contribution of water to the fluctuating forces. The slow component reflects disorder-induced inhomogeneous broadening.

Selectively Probing the Structures and Dynamics of β-Peptide Aggregates Using the Amide-A Vibrational Marker

The N-H stretching vibration in a β-peptide model compound, N-ethylpropionamide (NEPA), was characterized by one-dimensional infrared (1D IR) and two-dimensional (2D) IR experiments and ab initio anharmonic frequency computations. A narrowband pump-broadband probe 2D IR method was applied to selectively probe a subensemble of the N-H stretching vibrations from a mixture of different NEPA molecular aggregates that were formed via an intermolecular hydrogen bond. Vibrational lifetime and anharmonicity were found to be sensitive to the aggregation ensembles. In particular, diagonal anharmonicities were observed experimentally and confirmed computationally to be smaller for NEPA trimer than for dimer, which was explained by the presence of non-negligible off-diagonal anharmonicities in coupled N-H stretching modes.

Anharmonicities and coherent vibrational dynamics of phosphate ions in bulk H2O

2D IR spectroscopy reveals Fermi resonances and long lived quantum beats for phosphate ions in water.