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

Angiogenin-catalyzed cleavage within tRNA anticodon-loops identified by cP-RNA-seq

Angiogenin (Ang), an endoribonuclease belonging to the RNase A superfamily, cleaves the anticodon-loops of tRNAs to produce tRNA-half molecules. Although previous studies have demonstrated the involvement of Ang in the pathobiology of neurodegenerative disorders, the characterization of Ang-generated tRNA halves in neuronal cells remains limited. This is partly due to the technical limitations of standard RNA-seq methods, which cannot capture Ang-generated RNAs containing a 2',3'-cyclic phosphate (cP). In this report, we established an Ang-treatment model using SH-SY5Y, a human neuroblastoma cell line, and demonstrated Ang-dependent accumulation of tRNA halves. By performing cP-RNA-seq, which selectively captures cP-containing RNAs, we identified Ang-generated tRNA halves and the specific cleavage positions within tRNA anticodon-loops responsible for their generation. Our results provide insights into the anticodon-loop cleavage and the selective production of a specific subset of tRNA halves by Ang.

RNA Order Regulates Its Interactions with Zwitterionic Lipid Bilayers

RNA-lipid interactions directly influence RNA activity, which plays a crucial role in the development of new applications in medicine and biotechnology. However, while specific preferential behaviors between RNA and lipid bilayers have been identified experimentally, their molecular origin remains unexplored. Here we use molecular dynamics simulations to investigate the interaction between RNA and membranes composed of zwitterionic lipids at the atomistic level. Our data reproduce and rationalize previous experimental observations, including that short-chain RNAs rich in guanine have a higher affinity for gel-phase membranes compared to RNA sequences rich in other nucleotides and that RNA prefers gel-phase membranes to fluid bilayers. Our simulations reveal that RNA order is a key molecular determinant of RNA-zwitterionic phospholipid interactions. Our data provide a wealth of information at the atomic level that will help accelerate research on RNA-lipid assemblies for task-specific applications such as designing lipid-based nanocarriers for RNA delivery.

Modeling Infrared Spectroscopy of Nucleic Acids: Integrating Vibrational Non-Condon Effects with Machine Learning Schemes

Vibrational non-Condon effects, which describe how molecular vibrational transitions are influenced by a system's rotational and translational degrees of freedom, are often overlooked in spectroscopy studies of biological macromolecules. In this work, we explore these effects in the modeling of infrared (IR) spectra for nucleic acids in the 1600-1800 cm-1 region. Through electronic structure calculations, we reveal that the transition dipole moments of the C═O and C═C stretching modes in nucleobases are highly sensitive to solvation, hydrogen bonding, and base stacking conditions. To incorporate vibrational non-Condon effects into spectroscopy modeling, we use local electric fields on chromophore atoms as collective coordinates and leverage experimental IR spectra of oligonucleotides to develop deep neural network-based transition dipole strength (TDS) maps for the C═O and C═C chromophores. By integrating molecular dynamics simulations with a mixed quantum/classical treatment of the line shape theory, we apply the TDS maps to calculate the IR spectra of nucleoside 5'-monophosphates, DNA double helices and yeast phenylalanine tRNA. The resulting theoretical spectra show quantitative agreement with experimental measurements. While the predictions for nucleoside 5'-monophosphates are comparable to baseline performance, the TDS maps yield significantly improved IR peak intensities across all oligonucleotides. This theoretical framework effectively bridges atomistic simulations and IR spectroscopy experiments, offering molecular insights into how vibrational non-Condon effects impact the observed spectral features.

Structure of Essential RNA Regulatory Elements in the West Nile Virus 3'-Terminal Stem Loop

Flaviviruses, such as West Nile and Dengue Virus, pose a significant and growing threat to global health. Central to the flavivirus life cycle are highly structured 5'- and 3'-untranslated regions (UTRs), which harbor conserved cis-acting RNA elements critical for viral replication and host adaptation. Despite their essential roles, detailed molecular insights into these RNA elements have been limited. By employing nuclear magnetic resonance (NMR) spectroscopy in conjunction with SAXS experiments, we determined the three-dimensional structure of the West Nile Virus (WNV) 3'-terminal stem-loop core, a highly conserved element critical for viral genome cyclization and replication. Single nucleotide mutations at several sites within this RNA abolish the ability of the virus to replicate. These critical sites are located within a short 18-nucleotide hairpin stem, a substructure notable for its conformational flexibility, while the adjoining main stem-loop adopts a well-defined extended helix interrupted by three non-Watson-Crick pairs. This study enhances our understanding of several metastable RNA structures that play key roles in regulating the flavivirus lifecycle, and thereby also opens up potential new avenues for the development of antivirals targeting these conserved RNA structures. In particular, the structure we observe suggests that the plastic junction between the small hairpin and the tail of the longer stem-loop could provide a binding pocket for small molecules, for example potentially stabilizing the RNA in a conformation which hinders the conformational rearrangements critical for viral replication.

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.

Developing MiR-133a Zipper Nanoparticles for Targeted Enhancement of Thermogenic Adipocyte Generation

Existing delivery methods for RNAi therapeutics encounter challenges, including stability, specificity, and off-target effects, which restrict their clinical effectiveness. In this study, a novel miR-133a zipper nanoparticle (NP) system that integrates miRNA zipper technology with rolling circle transcription (RCT) to achieve targeted delivery and specific regulation of miR-133a in adipocytes, is presented. This innovative approach can greatly enhance the delivery and release of miR-133a zippers, increasing the expression of thermogenic genes and mitochondrial biogenesis. he miR-133a zipper NP is utilized for the delivery of miRNA zipper-blocking miR-133a, an endogenous inhibitor of Prdm16 expression, to enhance the thermogenic activity of adipocytes by modulating their transcriptional program. Inhibition of miR-133a through the miR-133a zipper NP leads to more significant upregulation of thermogenic gene expression (Prdm16 and Ucp1) than with the free miR-133a zipper strand. Furthermore, miR-133a zipper NPs increase the number of mitochondria and induce heat production, reducing the size of 3D adipose spheroids. In short, this study emphasizes the role of RNA NPs in improving RNAi stability and specificity and paves the way for broader applications in gene therapy. Moreover, this research represents a significant advancement in RNAi-based treatments, pointing toward a promising direction for future therapeutic strategies.

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.

Unraveling the Structure-Spectrum Relationship of Yeast Phenylalanine Transfer RNA: Insights from Theoretical Modeling of Infrared Spectroscopy

Yeast phenylalanine tRNA (tRNAphe) is a paradigmatic model in structural biology. In this work, we combine molecular dynamics simulations and spectroscopy modeling to establish a direct link between its structure, conformational dynamics, and infrared (IR) spectra. Employing recently developed vibrational frequency maps and coupling models, we apply a mixed quantum/classical treatment of the line shape theory to simulate the IR spectra of tRNAphe in the 1600-1800 cm-1 region across its folded and unfolded conformations and under varying concentrations of Mg2+ ions. The predicted IR spectra of folded and unfolded tRNAphe are in good agreement with experimental measurements, validating our theoretical framework. We then elucidate how the characteristic L-shaped tertiary structure of the tRNA and its modulation in response to diverse chemical environments give rise to distinct IR absorption peaks and line shapes. These calculations effectively bridge IR spectroscopy experiments and atomistic molecular simulations, unraveling the molecular origins of the observed IR spectra of tRNAphe. This work presents a robust theoretical protocol for modeling the IR spectroscopy of nucleic acids, which will facilitate its application as a sensitive probe for detecting the fluctuating secondary and tertiary structures of these essential biological macromolecules.

Insights on the comparative affinity of ribonucleic acids with plant-based beta carboline alkaloid, harmine: Spectroscopic, calorimetric and computational evaluation

Small molecules as ligands target multifunctional ribonucleic acids (RNA) for therapeutic engagement. This study explores how the anticancer DNA intercalator harmine interacts various motifs of RNAs, including the single-stranded A-form poly (rA), the clover leaf tRNAphe, and the double-stranded A-form poly (rC)-poly (rG). Harmine showed the affinity to the polynucleotides in the order, poly (rA) > tRNAphe > poly (rC)·poly (rG). While no induced circular dichroism change was detected with poly (rC)poly (rG), significant structural alterations of poly (rA) followed by tRNAphe and occurrence of concurrent initiation of optical activity in the attached achiral molecule of alkaloid was reported. At 25 °C, the affinity further showed exothermic and entropy-driven binding. The interaction also highlighted heat capacity (ΔC o p ) and Gibbs energy contribution from the hydrophobic transfer (ΔG hyd) of binding with harmine. Molecular docking calculations indicated that harmine exhibits higher affinity for poly (rA) compared to tRNAphe and poly (rC)·poly (rG). Subsequent molecular dynamics simulations were conducted to investigate the binding mode and stability of harmine with poly(A), tRNAphe, and poly (rC)·poly (rG). The results revealed that harmine adopts a partial intercalative binding with poly (rA) and tRNAphe, characterized by pronounced stacking forces and stronger binding free energy observed with poly (rA), while a comparatively weaker binding free energy was observed with tRNAphe. In contrast, the stacking forces with poly (rC)·poly (rG) were comparatively less pronounced and adopts a groove binding mode. It was also supported by ferrocyanide quenching analysis. All these findings univocally provide detailed insight into the binding specificity of harmine, to single stranded poly (rA) over other RNA motifs, probably suggesting a self-structure formation in poly (rA) with harmine and its potential as a lead compound for RNA based drug targeting.

Elucidating the influence of RNA modifications and Magnesium ions on tRNA Phe conformational dynamics in S. cerevisiae : Insights from Replica Exchange Molecular Dynamics simulations

Post-transcriptional modifications in RNA can significantly impact their structure and function. In particular, transfer RNAs (tRNAs) are heavily modified, with around 100 different naturally occurring nucleotide modifications contributing to codon bias and decoding efficiency. Here, we describe our efforts to investigate the impact of RNA modifications on the structure and stability of tRNA Phenylalanine (tRNA Phe ) from S. cerevisiae using molecular dynamics (MD) simulations. Through temperature replica exchange MD (T-REMD) studies, we explored the unfolding pathway to understand how RNA modifications influence the conformational dynamics of tRNA Phe , both in the presence and absence of magnesium ions (Mg 2+ ). We observe that modified nucleotides in key regions of the tRNA establish a complex network of hydrogen bonds and stacking interactions which is essential for tertiary structure stability of the tRNA. Furthermore, our simulations show that modifications facilitate the formation of ion binding sites on the tRNA. However, high concentrations of Mg 2+ ions can stabilize the tRNA tertiary structure in the absence of modifications. Our findings illuminate the intricate interactions between modifications, magnesium ions, and RNA structural stability.

Collision-Induced Unfolding Reveals Disease-Associated Stability Shifts in Mitochondrial Transfer Ribonucleic Acids

Ribonucleic acids (RNAs) remain challenging targets for structural biology, creating barriers to understanding their vast functions in cellular biology and fully realizing their applications in biotechnology. The inherent dynamism of RNAs creates numerous obstacles in capturing their biologically relevant higher-order structures (HOSs), and as a result, many RNA functions remain unknown. In this study, we describe the development of native ion mobility-mass spectrometry and collision-induced unfolding (CIU) for the structural characterization of a variety of RNAs. We evaluate the ability of these techniques to preserve native structural features in the gas phase across a wide range of functional RNAs. Finally, we apply these tools to study the elusive mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes-associated A3243G mutation. Our data demonstrate that our experimentally determined conditions preserve some solution-state memory of RNAs via the correlated complexity of CIU fingerprints and RNA HOS, the observation of predicted stability shifts in the control RNA samples, and the retention of predicted magnesium binding events in gas-phase RNA ions. Significant differences in collision cross section and stability are observed as a function of the A3243G mutation across a subset of the mitochondrial tRNA maturation pathway. We conclude by discussing the potential application of CIU for the development of RNA-based biotherapeutics and, more broadly, transcriptomic characterization.

Direct sequencing of total Saccharomyces cerevisiae tRNAs by LC-MS/MS

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

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.

Investigation of the ionic conditions in SiRNA-mediated delivery through its carriers in the cell membrane: a molecular dynamic simulation

SiRNA is a new generation of drug molecules and a new approach for treating a variety of diseases such as cancer and viral infections. SiRNA delivery to cells and translocation into cytoplasm are the main challenges in the clinical application of siRNA. Lipid carriers are one of the most successful carriers for siRNA delivery. In this study, we investigated the interaction of siRNA with a zwitterionic bilayer and how ion concentration and lipid conjugation can affect it. The divalent cation such as Mg²⁺ ions could promote the siRNA adsorption on the bilayer surface. The cation ions can bind to the head groups of lipids and the grooves of siRNA molecules and form bridges between the siRNA and bilayer surface. Our findings demonstrated the bridges formed by divalent ions could facilitate the attachment of siRNA to the membrane surface. We showed that the divalent cations can regulate the bridging-driven membrane attachment and it seems the result of this modulation can be used for designing biomimetic devices. In the following, we examined the effect of cations on the interaction between siRNA modified by cholesterol and the membrane surface. Our MD simulations showed that in the presence of Mg²⁺, the electrostatic and vdW energy between the membrane and siRNA were higher compared to those in the presence of NA⁺. We showed that the electrostatic interaction between membrane and siRNA cannot be facilitated only by cholesterol conjugated. Indeed, cations are essential to create coulomb repulsion and enable membrane attachment. This study provides important insight into liposome carriers for siRNA delivery and could help us in the development of siRNA-based therapeutics. Due to the coronavirus pandemic outbreak, these results may shed light on the new approach for treating these diseases and their molecular details.

NMR assignment of non-modified tRNAIle from Escherichia coli

tRNAs are L-shaped RNA molecules of ~ 80 nucleotides that are responsible for decoding the mRNA and for the incorporation of the correct amino acid into the growing peptidyl-chain at the ribosome. They occur in all kingdoms of life and both their functions, and their structure are highly conserved. The L-shaped tertiary structure is based on a cloverleaf-like secondary structure that consists of four base paired stems connected by three to four loops. The anticodon base triplet, which is complementary to the sequence of the mRNA, resides in the anticodon loop whereas the amino acid is attached to the sequence CCA at the 3'-terminus of the molecule. tRNAs exhibit very stable secondary and tertiary structures and contain up to 10% modified nucleotides. However, their structure and function can also be maintained in the absence of nucleotide modifications. Here, we present the assignments of nucleobase resonances of the non-modified 77 nt tRNA^Ile from the gram-negative bacterium Escherichia coli. We obtained assignments for all imino resonances visible in the spectra of the tRNA as well as for additional exchangeable and non-exchangeable protons and for heteronuclei of the nucleobases. Based on these assignments we could determine the chemical shift differences between modified and non-modified tRNA^Ile as a first step towards the analysis of the effect of nucleotide modifications on tRNA's structure and dynamics.

RNA modifications stabilize the tertiary structure of tRNAfMet by locally increasing conformational dynamics

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA's secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

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.

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.