Phosphate-Magnesium Ion Interactions in Water Probed by Ultrafast Two-Dimensional Infrared Spectroscopy

Structure and Dynamics of ATP and the ATP-Zn2+ Complex in Solution

Despite the crucial role of ATP in life and artificial life-like applications, fundamental aspects relevant to its function, such as its conformational properties and its interaction with water and ions, remain unclear. Here, by integrating linear and two-dimensional infrared spectroscopy with ab initio molecular dynamics, we provide a detailed characterization of the vibrational spectra of the phosphate groups in ATP and in its complex with Zn2+ in water. Our study highlights the role of conformational disorder and solvation dynamics, beyond the harmonic normal-mode analysis, and reveals a complex scenario in which electronic and environmental effects tune the coupling between phosphate vibrations. We identify βγ-bidentate and αβγ-tridentate modes as the preferential coordination modes of Zn2+, as was proposed in the literature for Mg2+, although this conclusion is reached by a different spectral interpretation.

Hydration structure and dynamics of phosphoric acid and its anions-Ultrafast 2D-IR spectroscopy and ab initio molecular dynamics simulations

The hydration shells of phosphate ions and phosphate groups of nucleotides and phospholipid membranes display markedly different structures and hydrogen-bond strengths. Understanding phosphate hydration requires insight into the spatial arrangements of water molecules around phosphates and in thermally activated structure fluctuations on ultrafast time scales. Femtosecond two-dimensional infrared spectroscopy of phosphate vibrations, particularly asymmetric stretching vibrations between 1000 and 1200 cm-1, and ab initio molecular dynamics (AIMD) simulations are combined to map and characterize dynamic local hydration structures and phosphate-water interactions. Phosphoric acid H3PO4 and its anions H2PO4-, HPO42-, and PO43- are studied in aqueous environments of different pH value. The hydration shells of phosphates providing OH donor groups in hydrogen bonds with the first water layer undergo ultrafast structural fluctuations, which induce a pronounced spectral diffusion of vibrational excitations on a sub-300 fs time scale. With a decreasing number of phosphate OH groups, the hydration shell becomes more ordered and rigid. The 2D-IR line shapes observed with hydrated PO43- ions display a pronounced inhomogeneous broadening, reflecting a distribution of hydration geometries without fast equilibration. The AIMD simulations allow for an in-depth characterization of the hydration geometries with different numbers of water molecules in the first hydration layer and different correlation functions of the fluctuating electric field that the water environment exerts on the vibrational phosphate oscillators.

Magnesium: properties and rich chemistry for new material synthesis and energy applications

Magnesium (Mg) has many unique properties suitable for applications in the fields of energy conversion and storage. These fields presently rely on noble metals for efficient performance. However, among other challenges, noble metals have low natural abundance, which undermines their sustainability. Mg has a high negative standard reduction potential and a unique crystal structure, and its low melting point at 650 °C makes it a good candidate to replace or supplement numerous other metals in various energy applications. These attractive features are particularly helpful for improving the properties and limits of materials in energy systems. However, knowledge of Mg and its practical uses is still limited, despite recent studies which have reported Mg's key roles in synthesizing new structures and modifying the chemical properties of materials. At present, information about Mg chemistry has been rather scattered without any organized report. The present review highlights the chemistry of Mg and its uses in energy applications such as electrocatalysis, photocatalysis, and secondary batteries, among others. Future perspectives on the development of Mg-based materials are further discussed to identify the challenges that need to be addressed.

Recycling of Cofactors in Crude Enzyme Hydrogels as Co-immobilized Heterogeneous Biocatalysts for Continuous-Flow Asymmetric Reduction of Ketones

Flow biocatalysis involving oxidoreductase is limited by the difficulty in recycling expensive cofactors. In this study, an enzyme-rich hydrogel monolithic microreactor was developed via in situ microfluidic assembly of inexpensive crude enzymes. This porous gel biocatalyst exhibited good tethering functions to nicotinamide cofactors; thus, they were retained by the hydrogel to controllably form a novel heterogeneous biocatalyst with self-sufficient cofactors. The flow asymmetric production of a chiral alcohol in this cofactor-entrapped gel microreactor achieved >99 % enantioselectivity and a high space-time yield of 46.3 g L^-1  h^-1 at 94.8 % conversion. Moreover, the turnover number of cofactors reached as high as 4800 after continuous operation of 160 reactor volumes, realizing significantly higher utilization of the cofactors compared with many reported strategies. Furthermore, this engineered heterogeneous biocatalyst exhibited improved performance in terms of product tolerance and storage stability, paving the way for a green, cost-effective, and sustainable continuous-flow production of enantiopure alcohols.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg²⁺ and SO₄ ^2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO₄ ^2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO₄ ^2- ions around 1100 cm^-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO₄ and for Na₂SO₄ in the presence of Mg²⁺ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na₂SO₄ solutions in the absence of Mg²⁺ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO₄ ^2--Mg²⁺ ion pairs.

Magnesium ions reversibly bind to DNA double stranded helix in thin films

Effects of magnesium (Mg²⁺) ions on the stability and structural properties of double-stranded DNA are vitally important for DNA folding and functional behavior. Complementing our previous study on highly hydrated thin films of DNA with sodium counterions, with no buffer (pH ≈ 6) and surrounded with Mg²⁺ cations, here we use Fourier transform infrared spectroscopy and band shape analysis to explore in detail the vibrational signatures of DNA-magnesium interaction in the case when DNA charges are neutralized solely by Mg²⁺ cations, hereafter called MgDNA. Ion atmosphere has been controlled by the magnesium to phosphate molar concentration ratio r which varied between 0.0067 and 10. For r = 0 we find that spectral features in the base region remain similar as in DNA, whereas changes in the backbone region indicate that the B conformation becomes fully stabilized. With increasing r a pronounced structural reshaping occurs in the phosphate backbone region indicating a blue shift of the asymmetric band, while the symmetric band does not show any displacement in frequency. The band shape analysis of overlapping peaks in the respective phosphate regions demonstrates that the number of constituent modes as well as their positions in frequency do not change, whereas their intensities and bandwidths display disparate changes. The results reflect a variety of local environments at the DNA backbone due to a heterogeneous ion atmosphere with randomly distributed magnesium ions and local patterns of hydrogen bonds which change with increasing r. Remarkably, after crowded r = 10 ion atmosphere is depleted, Mg induced spectral changes vanish and structural features of MgDNA (r ≈ 0) are fully restored. Overall results strongly suggest that in MgDNA on highly hydrated thin films the hydrogen-base pairing remains preserved and that Mg²⁺ ions, similar to sodium ions, retain their mobility and interact with double helix via water-mediated electrostatic forces.

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg2+) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg2+ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg2+ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg2+-excess concentrations and comparison with fluorescence titration measurements shows that on average 20-30% of the Mg2+ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.

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.

Multimolecular complexes of the phosphodiester anion with Zn(II) or Mg(II) and water molecules-Preliminary validations of a polarizable potential by ab initio quantum chemistry

Dimethyl phosphate (DMP^- ) is a model for the phosphodiester backbone of DNA, RNA, and phospholipids. It is central for the binding of divalent cations and water along the backbone of nucleic acids. Significant polarization and charge-transfer contributions and nonadditivity come into play in the multimolecular complexes organized around phosphate. Prior to large-scale molecular dynamics (MD) with advanced polarizable potentials, it is essential to evaluate how well the values and trends of intermolecular interaction energies (ΔE) from ab initio quantum chemistry (QC) and their individual contributions are reproduced in a diversity of such complexes. These differ by the starting binding modes of a divalent cation, Zn(II), namely direct, bi- or mono-dentate to anionic and/or ester oxygens, versus through-water binding. We present first the results from automated refinements of the individual contributions of the SIBFA potential with respect to their QC counterparts using a Zn(II) or a water probe. This is followed by validations on eight relaxed multimolecular complexes of DMP^- with Zn(II) or Mg(II) and seven waters, then on sixteen complexes of DMP^- with Zn(II) and eight waters in arrangements extracted from MD or energy-minimization on a droplet of sixty-four waters. This monitors the compared evolutions of SIBFA and QC ΔE and their individual contributions in the competing arrangements. Some waters, bridging Zn(II) and DMP^- , were found to have exceptionally large dipole moments, of up to 3.8 Debye. The perspectives of extension to a flexible phosphodiester backbone are discussed in the context of the SIBFA potential for DNA and RNA.

Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions

Electric interactions have a strong impact on the structure and dynamics of biomolecules in their native water environment. Given the variety of water arrangements in hydration shells and the femto- to subnanosecond time range of structural fluctuations, there is a strong quest for sensitive noninvasive probes of local electric fields. The stretching vibrations of phosphate groups, in particular the asymmetric (PO₂)⁻ stretching vibration ν_AS(PO₂)⁻, allow for a quantitative mapping of dynamic electric fields in aqueous environments via a field-induced redshift of their transition frequencies and concomitant changes of vibrational line shapes. We present a systematic study of ν_AS(PO₂)⁻ excitations in molecular systems of increasing complexity, including dimethyl phosphate (DMP), short DNA and RNA duplex structures, and transfer RNA (tRNA) in water. A combination of linear infrared absorption, two-dimensional infrared (2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives quantitative insight in electric-field tuning rates of vibrational frequencies, electric field and fluctuation amplitudes, and molecular interaction geometries. Beyond neat water environments, the formation of contact ion pairs of phosphate groups with Mg²⁺ ions is demonstrated via frequency upshifts of the ν_AS(PO₂)⁻ vibration, resulting in a distinct vibrational band. The frequency positions of contact geometries are determined by an interplay of attractive electric and repulsive exchange interactions.

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations

Ions interacting with hydrated RNA play a central role in defining its secondary and tertiary structure. While spatial arrangements of ions, water molecules, and phosphate groups have been inferred from X-ray studies, the role of electrostatic and other noncovalent interactions in stabilizing compact folded RNA structures is not fully understood at the molecular level. Here, we demonstrate that contact ion pairs of magnesium (Mg²⁺) and phosphate groups embedded in local water shells stabilize the tertiary equilibrium structure of transfer RNA (tRNA). Employing dialyzed tRNA^Phe from yeast and tRNA from Escherichia coli, we follow the population of Mg²⁺ sites close to phosphate groups of the ribose-phosphodiester backbone step by step, combining linear and nonlinear infrared spectroscopy of phosphate vibrations with molecular dynamics simulations and ab initio vibrational frequency calculations. The formation of up to six Mg²⁺/phosphate contact pairs per tRNA and local field-induced reorientations of water molecules balance the phosphate–phosphate repulsion in nonhelical parts of tRNA, thus stabilizing the folded structure electrostatically. Such geometries display limited sub-picosecond fluctuations in the arrangement of water molecules and ion residence times longer than 1 μs. At higher Mg²⁺ excess, the number of contact ion pairs per tRNA saturates around 6 and weakly interacting ions prevail. Our results suggest a predominance of contact ion pairs over long-range coupling of the ion atmosphere and the biomolecule in defining and stabilizing the tertiary structure of tRNA.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

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.

Contact Ion Pairs of Phosphate Groups in Water: Two-Dimensional Infrared Spectroscopy of Dimethyl Phosphate and ab Initio Simulations

The interaction of phosphate groups with ions in an aqueous environment has a strong impact on the structure and folding processes of DNA and RNA. The dynamic variety of ionic arrangements, including both contact pairs and water separated ions, and the molecular coupling mechanisms are far from being understood. In a combined experimental and theoretical approach, we address the properties of contact ion pairs of the prototypical system dimethyl phosphate with Na⁺, Ca²⁺, and Mg²⁺ ions in water. Linear and femtosecond two-dimensional infrared (2D-IR) spectroscopy of the asymmetric (PO₂)^- stretching vibration separates and characterizes the different species via their blue-shifted vibrational signatures and 2D-IR line shapes. Phosphate-magnesium contact pairs stand out as the most compact geometry while the contact pairs with Ca²⁺ and Na⁺ display a wider structural variation. Microscopic density functional theory simulations rationalize the observed frequency shifts and reveal distinct differences between the contact geometries.

Effect of Multivalent Cations on Intermolecular Association of Isotactic and Atactic Poly(Methacrylic Acid) Chains in Aqueous Solutions

The formation of nanoparticles of two poly(methacrylic acid) (PMA) isomers, atactic (aPMA) and isotactic (iPMA), was investigated in aqueous solutions in the presence of mono- (Na⁺) and multivalent cations (Mg²⁺ and La³⁺). Using dynamic (DLS) and static light scattering (SLS), we show that PMA nanoparticles have characteristics of microgel-like particles with a denser core and a swollen corona. iPMA aggregates are stable at a much higher degree of neutralization (α_N) than the aPMA ones, indicating a much stronger association between iPMA chains. This is explained by proposing segregation of ionized and unionized carboxyl groups within the iPMA aggregates and subsequent cooperative hydrogen-bonding between COOH groups. The calculated shape parameter (ρ) suggests different behavior of both isomers in the presence of Mg²⁺ ions on one hand and Na⁺ and La³⁺ on the other. The microgel-like particles formed in the presence of Mg²⁺ ions have a more even mass distribution (possibly a no core-shell structure) in comparison with those in the presence of Na⁺ and La³⁺ ions. Differences between the aggregate structures in the presence of different ions are reflected also in calorimetric experiments and supported by pH and fluorimetric measurements. Reasons for different behavior in the presence of Mg²⁺ ions lie in specific properties of this cation, in particular in its strong hydration and preference towards monodentate binding to carboxylate groups.