Hydration of the phosphate group in double-helical DNA

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.

A mammalian model reveals inorganic polyphosphate channeling into the nucleolus and induction of a hyper-condensate state

Inorganic polyphosphate (polyP) is a ubiquitous polymer that controls fundamental processes. To overcome the absence of a genetically tractable mammalian model, we developed an inducible mammalian cell line expressing Escherichia coli polyphosphate kinase 1 (EcPPK1). Inducing EcPPK1 expression prompted polyP synthesis, enabling validation of polyP analytical methods. Virtually all newly synthesized polyP accumulates within the nucleus, mainly in the nucleolus. The channeled polyP within the nucleolus results in the redistribution of its markers, leading to altered rRNA processing. Ultrastructural analysis reveals electron-dense polyP structures associated with a hyper-condensed nucleolus resulting from an exacerbation of the liquid-liquid phase separation (LLPS) phenomena controlling this membraneless organelle. The selective accumulation of polyP in the nucleoli could be interpreted as an amplification of polyP channeling to where its physiological function takes place. Indeed, quantitative analysis of several mammalian cell lines confirms that endogenous polyP accumulates within the nucleolus.

Structure and Formation of Z-DNA and Z-RNA

Despite structural differences between the right-handed conformations of A-RNA and B-DNA, both nucleic acids adopt very similar, left-handed Z-conformations. In contrast to their structural similarities and sequence preferences, RNA and DNA exhibit differences in their ability to adopt the Z-conformation regarding their hydration shells, the chemical modifications that promote the Z-conformation, and the structure of junctions connecting them to right-handed segments. In this review, we highlight the structural and chemical properties of both Z-DNA and Z-RNA and delve into the potential factors that contribute to both their similarities and differences. While Z-DNA has been extensively studied, there is a gap of knowledge when it comes to Z-RNA. Where such information is lacking, we try and extend the principles of Z-DNA stability and formation to Z-RNA, considering the inherent differences of the nucleic acids.

Detection of Trace Amounts of Water in Organic Solvents by DNA-Based Nanomechanical Sensors

The detection of trace amounts of water in organic solvents is of great importance in the field of chemistry and in the industry. Karl Fischer titration is known as a classic method and is widely used for detecting trace amounts of water; however, it has some limitations in terms of rapid and direct detection because of its time-consuming sample preparation and specific equipment requirements. Here, we found that a DNA-based nanomechanical sensor exhibits high sensitivity and selectivity to water vapor, leading to the detection and quantification of trace amounts of water in organic solvents as low as 12 ppm in THF, with a ppb level of LoD through their vapors. Since the present method is simple and rapid, it can be an alternative technique to the conventional Karl Fischer titration.

Zeptoliter DNA Origami Reactor to Reveal Cosolute Effects on Nanoconfined G-Quadruplexes

Cellular environments such as nanoconfinement and molecular crowding can change biomolecular properties. However, in nanoconfinement, it is extremely challenging to investigate effects of crowding cosolutes on macromolecules. By using optical tweezers, here, we elucidated the effects of hexaethylene glycol (HEG) on the mechanical stability of a telomeric G-quadruplex (GQ) in a zeptoliter DNA origami reactor (zepto-reactor). When HEG molecules were introduced in the GQ-containing zepto-reactor at different positions, we found that the GQ species split into two equilibrated populations, reflecting diverse effects of the oligoethylene glycol on the GQ via either a long-range dehydration effect or direct interactions. When the number of HEG molecules was increased, the stability of the GQ unexpectedly decreased, suggesting that the direct destabilizing interaction between the GQ and HEG is dominating over the long-range stabilizing dehydration effects of the HEG in hydrophilic nanocavities. These findings indicate that a nanoconfined environment can alter regular effects of cosolutes on biomacromolecules.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the `spine of hydration' in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

New damage model for simulating radiation-induced direct damage to biomolecular systems and experimental validation using pBR322 plasmid

In this work, we proposed a new damage model for estimating radiation-induced direct damage to biomolecular systems and validated its the effectiveness for pBR322 plasmids. The proposed model estimates radiation-induced damage to biomolecular systems by: (1) simulation geometry modeling using the coarse-grained (CG) technique to replace the minimum repeating units of a molecule with a single bead, (2) approximation of the threshold energy for radiation damage through CG potential calculation, (3) calculation of cumulative absorption energy for each radiation event in microscopic regions of CG models using the Monte Carlo track structure (MCTS) code, and (4) estimation of direct radiation damage to biomolecular systems by comparing CG potentials and absorption energy. The proposed model replicated measured data with an average error of approximately 14.2% in the estimation of radiation damage to pBR322 plasmids using the common MCTS code Geant4-DNA. This is similar to the results of previous simulation studies. However, in existing damage models, parameters are adjusted based on experimental data to increase the reliability of simulation results, whereas in the proposed model, they can be determined without using empirical data. Because the proposed model proposed is applicable to DNA and various biomolecular systems with minimal experimental data, it provides a new method that is convenient and effective for predicting damage in living organisms caused by radiation exposure.

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.

Hydration Shells of DNA from the Point of View of Terahertz Time-Domain Spectroscopy

Hydration plays a fundamental role in DNA structure and functioning. However, the hydration shell has been studied only up to the scale of 10-20 water molecules per nucleotide. In the current work, hydration shells of DNA were studied in a solution by terahertz time-domain spectroscopy. The THz spectra of three DNA solutions (in water, 40 mm MgCl₂ and 150 mM KCl) were transformed using an effective medium model to obtain dielectric permittivities of the water phase of solutions. Then, the parameters of two relaxation bands related to bound and free water molecules, as well as to intermolecular oscillations, were calculated. The hydration shells of DNA differ from undisturbed water by the presence of strongly bound water molecules, a higher number of free molecules and an increased number of hydrogen bonds. The presence of 40 mM MgCl₂ in the solution almost does not alter the hydration shell parameters. At the same time, 150 mM KCl significantly attenuates all the found effects of hydration. Different effects of salts on hydration cannot be explained by the difference in ionic strength of solutions, they should be attributed to the specific action of Mg²⁺ and K⁺ ions. The obtained results significantly expand the existing knowledge about DNA hydration and demonstrate a high potential for using the THz time-domain spectroscopy method.

Abnormal Properties of Low-Dimensional Confined Water

Water molecules confined to low-dimensional spaces exhibit unusual properties compared to bulk water. For example, the alternating hydrophilic and hydrophobic nanodomains on flat silicon wafer can induce the abnormal spreading of water (contact angles near 0°) which is caused by the 2D capillary effect. Hence, exploring the physicochemical properties of confined water from the nanoscale is of great value for understanding the challenges in material science and promoting the applications of nanomaterials in the fields of mass transport, nanofluidic designing, and fuel cell. The knowledge framework of confined water can also help to better understand the complex functions of the hydration layer of biomolecules, and even trace the origin of life. In this review, the physical properties, abnormal behaviors, and functions of the confined water are mainly summarized through several common low-dimensional water formats in the fields of solid/air-water interface, nanochannel confinement, and biological hydration layer. These researches indicate that the unusual behaviors of the confined water depend strongly on the confinement size and the interaction between the molecules and confining surface. These diverse properties of confined water open a new door to materials science and may play an important role in the future development of biology.

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.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Effects of Phosphoryl Guanidine Modification of Phosphate Residues on the Structure and Hybridization of Oligodeoxyribonucleotides

Phosphoryl guanidine oligonucleotides (PGOs) are promising tools for biological research and development of biosensors and therapeutics. We performed structural and hybridization analyses of octa-, deca-, and dodecamers with all phosphate residues modified by 1,3-dimethylimidazolidine-2-imine moieties. Similarity of the B-form double helix between native and modified duplexes was noted. In PGO duplexes, we detected a decrease in the proportion of C2'-endo and an increased proportion of C1'-exo sugar conformations of the modified chain. Applicability of the two-state model to denaturation transition of all studied duplexes was proved for the first time. Sequence-dependent effects of this modification on hybridization properties were observed. The thermal stability of PGO complexes is almost native at 100 mM NaCl and slightly increases with decreasing ionic strength. An increase in water activity and dramatic changes in interaction with cations and in solvation of PGOs and their duplexes were noted, resulting in slight elevation of the melting temperature after an ionic-strength decrease from 1 M NaCl down to deionized water. Decreased binding of sodium ions and decreased water solvation were documented for PGOs and their duplexes. In contrast to DNA, the PGO duplex formation leads to a release of several cations. The water shell is significantly more disordered near PGOs and their complexes. Nevertheless, changes in solvation during the formation of native and PGO complexes are similar and indicate that it is possible to develop models for predictive calculations of the thermodynamic properties of phosphoryl guanidine oligomers. Our results may help devise an approach for the rational design of PGOs as novel improved molecular probes and tools for many modern methods involving oligonucleotides.

Detection of nucleotides in hydrated ssDNA via 2D h-BN nanopore with ionic-liquid/salt-water interface

Accomplishing slow translocation speed with high sensitivity has been the most critical mission for solid-state nanopore (SSN) device to electrically detect nucleobases in ssDNA. In this study, a method to detect nucleobases of ssDNA using a 2D SSN is introduced by considerably reducing the translocation speed and effectively increasing its sensitivity. The ultra-thin titanium dioxide coated hexagonal boron nitride nanopore was fabricated, along with an ionic-liquid 1-butyl-3-methylimidazolium hexafluorophosphate/2.0 M KCl aqueous (cis/trans) interface, for increasing both the spatial and the temporal resolutions. As the ssDNA molecules entered the nanopore, a brief surge of electrical conductivity occurred, which was followed by multiple resistive pulses from nucleobases during the translocation of ssDNA and another brief current surge flagging the exit of the molecule. The continuous detection of nucleobases using a 2D SSN device is a novel achievement: the water molecules bound to ssDNA increased the molecular conductivity and amplified electrical signals during the translocation. Along with the experiment, computational simulations using COMSOL Multiphysics are presented to explain the pivotal role of water molecules bound to ssDNA to detect nucleobases using a 2D SSN.

Ionic magnetic core-shell nanoparticles for DNA extraction

Magnetic nanoparticles with specific surface features are interesting materials for biomedical applications. The combination of molecular interactions on small particles with macroscopic cohesion forces offers unique opportunities. This work reports the synthesis of magnetic core-shell nanoparticles with alkylimidazolium coated surface for effective DNA extraction. A magnetic Fe2O3 core was coated with a silica shell and functionalized with an organic halide. This enabled a surface coating with organic cations to mediate effective molecular interactions with polyanionic DNA. The large surface area of the ∼20 nm small particles with a magnetization of 25 emu g-1 enabled high DNA particle loading of 1/30 m% with easy isolation based on an external magnetic field. Moreover, the coating of the particles stabilized DNA against ultrasound initiated fragmentation.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

Structural alphabets for conformational analysis of nucleic acids available at dnatco.datmos.org

A detailed description of the dnatco.datmos.org web server implementing the universal structural alphabet of nucleic acids is presented. It is capable of processing any mmCIF- or PDB-formatted files containing DNA or RNA molecules; these can either be uploaded by the user or supplied as the wwPDB or PDB-REDO structural database access code. The web server performs an assignment of the nucleic acid conformations and presents the results for the intuitive annotation, validation, modeling and refinement of nucleic acids.

Self-Assembling Glycopeptide Conjugate as a Versatile Platform for Mimicking Complex Polysaccharides

Polysaccharides are a class of carbohydrates that play pivotal roles in living systems such as being chemical messengers in many vital biological pathways. However, the complexity and heterogeneity of these natural structures have posed daunting challenges on their production, characterization, evaluation, and applications. While there have been various types of synthetic skeletons that could mimic some biological aspects of polysaccharides, a safer and more easily accessed system is still desired to avoid the unnatural components and difficulties in modifying the structures. In this work, conveniently accessible self-assembling glycopeptide conjugates are developed, where the natural O-glycosidic linkages and phosphoryl modifications assist the self-assembly and concurrently reduce the risk of toxicity. The generated nanoparticles in aqueous solution offer a multivalent display of structurally controllable carbohydrates as mimics of polysaccharides, among which a mannosylated version exhibits immunostimulatory effects in both cellular assays and vaccination of mice. The obtained results demonstrate the potential of this glycopeptide conjugate-derived platform in exploiting the intriguing properties of carbohydrates in a more structurally maneuverable fashion.

Hydrogen bonding of ionic liquids in the groove region of DNA controls the extent of its stabilization: synthesis, spectroscopic and simulation studies

Ionic liquids (ILs) have been extensively used for stabilization and long-term DNA storage. However, molecular level understanding of the role of the hydrogen bond of DNA with ILs in its stabilization is still inadequate. Two ILs, namely, 1,1,3,3-tetramethylguanidinium acetate (TMG) and 2,2-diethyl-1,1,3,3-tetramethylguanidinium acetate (DETMG), have been synthesized, of which TMG has a hydrogen bonding N-H group whereas DETMG does not contain any hydrogen bonding site. It has been found that both TMG and DETMG cations interact in the groove region of DNA; however, their mode of interaction is distinctly different, which causes the stabilization of DNA in the presence of TMG, whereas the effect is opposite in the case of DETMG. It is apparent from the data that only the accommodation of ILs in the groove region is not enough for the stabilization of DNA. MD simulation and spectroscopic studies combinedly indicate that the hydrogen bonding capability of the TMG cation enhances the hydrogen bonding between the Watson-Crick base pairs of DNA, resulting in its stabilization. In contrast, the bigger size as well as the absence of the hydrogen bonding site of the DETMG cation perturbs the minor groove width and base pair step parameters of DNA during its intrusion into the minor groove, which decreases the hydrogen bond between the Watson-Crick base pairs of DNA, leading to destabilization.

Effects of counterions and solvents on the geometrical and vibrational features of dinucleoside-monophosphate (dNMP): case of 3',5'-dideoxycytidine-monophosphate (dDCMP)

The effects of the interaction of the monovalent (Li⁺, Na⁺, K⁺) and divalent (Mg²⁺) counterions hexahydrated (6H₂O), with the PO₂^- group, on the geometrical and vibrational characteristics of 3', 5'-dDCMP, were studied using the DFT/B3LYP/6-31++G(d) method. These calculations were performed using the explicit (6H₂O) and hybrid (6H₂O/Continuum) solvation models. The optimizations reveal that in the conformation g^-g^- and in the explicit model of solvation, the small ions (Li⁺, Na⁺) deviate from the bisector plane of the angle O1-P-O2 and the large ions (K⁺ and Mg²⁺) remain in this plane, whereas in the hybrid model of solvation, the counterions deviate from this plane. However, when the conformer is g⁺g⁺, the monovalent counterions deviate and divide the remainder of the plane regardless of the type of solvation model. In addition, the g^-g^- conformer is the most stable in the presence of the explicit solvent, while the g⁺g⁺ conformer is the most stable in the presence of the hybrid solvent. Finally, the normal modes of the conformers g^-g^- and g⁺g⁺ in the presence of the counterions in the hybrid model show a better agreement with the available experimental data of the DNA forms A, B (g^-g^-), and Z (g⁺g⁺) relatively to the explicit model. This very good agreement is illustrated by the very small deviations ≤ 0.08% (g^-g^-) and ≤ 0.41% (g⁺g⁺) observed between the calculated and experimental data for the PO₂^- (asymmetric) stretching mode in the presence of the counterion K⁺ in the hybrid model. Graphical abstract.

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations

The water hydration shell has a decisive impact on the structural and functional properties of RNA. Changes of the RNA structure upon melting and in biochemical processes are accompanied by a change of hydration patterns, a process that is barely characterized. To discern hydration geometries around the backbone phosphate groups of an RNA double helix at the molecular level, we combine two-dimensional infrared spectroscopy of phosphate vibrations with theoretical simulations. There are three distinct coexisting hydration motifs of the RNA A-helix: an ordered chain-like arrangement of water molecules with links between neighboring phosphate groups, separate local hydration shells of up to six water molecules, and hydrated phosphate/counterion contact pairs. RNA disordering upon melting is connected with a transition from predominantly ordered water structures to local hydration shells around phosphate units. Structural fluctuations are dominated by librational water motions occurring on a 300 fs time scale, without exchange between hydration motifs.

Chromatin Compaction Multiscale Modeling: A Complex Synergy Between Theory, Simulation, and Experiment

Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.

Experimental Evidence of Solvent-Separated Ion Pairs as Metastable States in Electrostatic Interactions of Biological Macromolecules

Electrostatic interactions via ion pairs are vital for biological macromolecules. Regarding the free energy of each ion pair as a function of the interionic distance, continuum electrostatic models predict a single energy minimum corresponding to the contact ion-pair (CIP) state, whereas atomically detailed theoretical hydration studies predict multiple energy minima corresponding to the CIP and solvent-separated ion-pair (SIP) states. Through a statistical analysis of high-resolution crystal structures, we present experimental evidence of the SIP as a metastable state. The histogram of interionic distances between protein side-chain NH₃⁺ and DNA phosphate groups clearly shows two major peaks corresponding to the CIP and SIP states. The statistical data are consistent with the probability distribution of the CIP-SIP equilibria previously obtained with molecular dynamics simulations. Spatial distributions of NH₃⁺ ions and water molecules around phosphates reveal preferential sites for CIP and SIP formations and show how the ions compete with water molecules.

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.

Dynamical aspects of supercooled TIP3P-water in the grooves of DNA

We investigate by molecular dynamics simulations the mobility of the water located at the DNA minor and major grooves. We employ the TIP3P water model, and our system is analyzed for a range of temperatures 190-300 K. For high temperatures, the water at the grooves shows an Arrhenius behavior similar to that observed in the bulk water. At lower temperatures, a departure from the bulk behavior is observed. This slowing down in the dynamics is compared with the dynamics of the hydrogen of the DNA at the grooves and with the autocorrelation functions of the water hydrogen bonds. Our results indicate that the hydrogen bonds of the water at the minor grooves are highly correlated, which suggests that this is the mechanism for the slow dynamics at this high confinement.

Protein-Polymer Dynamics as Affected by Polymer Coating and Interactions

We investigate the relaxation dynamics of protein-polymer conjugates by neutron scattering spectroscopy to understand to which extent the coating of a protein by a polymer can replace water in promoting thermal structural fluctuations. For this purpose, we compare the dynamics of protein-polymer mixtures to that of conjugates with a variable number of polymers covalently attached to the protein. Results show that the flexibility of the protein is larger in protein-polymer mixtures than in native protein or in conjugates, even in the dry state. Upon hydration, both the native protein and the conjugate show equivalent dynamics, suggesting that the polymer grafted on the protein surface adsorbs all water molecules.

Hydration Structure of a Single DNA Molecule Revealed by Frequency-Modulation Atomic Force Microscopy

Hydration interaction shapes biomolecules and is a dominant intermolecular force. Mapping the hydration patterns of biomolecules is therefore essential for understanding molecular processes in biology. Numerous studies have been devoted to this challenge, but current methods cannot map the hydration of single biomolecules, let alone do so under physiological conditions. Here, we show that frequency-modulation atomic force microscopy (FM-AFM) can fill this gap and generate 3D hydration maps of single DNA molecules under near-physiological conditions. Additionally, we present real-space images of DNA in which the double helix is resolved with unprecedented resolution, clearly revealing individual phosphate groups along the DNA backbone. FM-AFM therefore emerges as a powerful enabling tool in the study of individual biomolecules and their hydration under physiological conditions.

DFT study of geometrical and vibrational features of a 3',5'-deoxydisugar-monophosphate (dDSMP) DNA model in the presence of counterions and solvent

The B3LYP/6-31++G* theoretical level was used to study the influence of various hexahydrated monovalent (Li⁺, Na⁺, K⁺) and divalent (Mg²⁺) metal counterions in interaction with the charged PO₂^- group, on the geometrical and vibrational characteristics of the DNA fragments of 3',5'-dDSMP, represented by four conformers (g⁺g⁺, g⁺t, g^-g^- and g^-t). All complexes were optimized through two solvation models [the explicit model (6H₂O) and the hybrid model (6H₂O/Continuum)]. The results obtained established that, in the hybrid model, counterions (Li⁺, Na⁺, K⁺, Mg²⁺) always remain in the bisector plane of the O1-P-O2 angle. When these counterions are explicitly hydrated, the smallest counterions (Li⁺, Na⁺) deviate from the bisector plane, while the largest counterions (K⁺ and Mg²⁺) always remain in the same plane. On the other hand, the present calculations reveal that the g⁺g⁺ conformer is the most stable in the presence of monovalent counterions, while conformers g⁺t and g^-t are the most stable in the presence of the divalent counterion Mg²⁺. Finally, the hybrid solvation model seems to be in better agreement with the available crystallographic and spectroscopic (Raman) experiments than the explicit model. Indeed, the six conformational torsions of the C4'-C3'-O3'-PO^-₂-O5'-C5'-C4' segment of all complexes of the g^-g^- conformer in 6H₂O/Continuum remain similar to the available experimental data of A- and B-DNA forms. The calculated wavenumbers of the g⁺g⁺ conformer in the presence of the monovalent counterion and of g^-t conformer in presence of the divalent counterion in the hybrid model are in good agreement with the Raman experimental data of A- and B-DNA forms. In addition, the maximum deviation between the calculated wavenumbers in the 6H₂O/Continuum for the g⁺g⁺ conformer and experimental value measured in an aqueous solution of the DMP-Na⁺ complex, is <1.07% for the PO₂^- (asymmetric and symmetric) stretching modes and <2.03% for the O5'-C5' and O3'-C3' stretching modes. Graphical abstract dDSMP-(OO)^- Mg²⁺/6W/Continuum.

Kinetics of the B-A transition of DNA: analysis of potential contributions to a reaction barrier

Because of open problems in the relation between results obtained by relaxation experiments and molecular dynamics simulations on the B-A transition of DNA, relaxation measurements of the B-A dynamics have been extended to a wider range of conditions. Field-induced reaction effects are measured selectively by the magic angle technique using a novel cell construction preventing perturbations from cell window anisotropy. The kinetics was recorded for the case of poly[d(AT)] up to the salt concentration limit of 4.4 mM, where aggregation does not yet interfere. Now experimental data on the B-A dynamics are available for poly[d(AT)] at salt concentrations of 0.18, 0.73, 2.44 and 4.4 mM. In all cases, a spectrum of time constants is found, ranging from ~ 10 μs up to components approaching ~ 1 ms. The relatively small dependence of these data on the salt concentration indicates that electrostatic effects on the kinetics are not as strong as may be expected. The ethanol content at the transition center is a linear function of the logarithm of the salt concentration, and the slope is close to that expected from polyelectrolyte theory. The B-A transition dynamics was also measured in D₂O at a salt concentration of 2.4 mM: the center of the transition is found at 20.0 mol/l H₂O and at 20.1 mol/l D₂O with an estimated accuracy of ± 0.1 mol/l; the spectrum of time constants at the respective transition centers is very similar. The experimental results are discussed regarding the data obtained by molecular dynamics simulations.

Effects of Near Infrared Radiation on DNA. DLS and ATR-FTIR Study

We presume that the primary effect of Near Infrared (NIR) radiation on aqueous solutions of biological molecules concerns modification of hydrogen bonded structures mainly the global and the hydration shell water molecules. Since water has a significant influence on the DNA structure, we expect that the thermal stability of DNA could be modified by NIR radiation. The herring sperm DNA was exposed to NIR radiation (700-1100nm) for 5, 10, and 20min periods. The temperature dependent infrared measurements were done for the thin films formed on the diamond ATR crystal from evaporated DNA solutions exposed and unexposed to NIR radiation. For the NIR-treated samples (at room temperature) the B form was better conserved than in the control sample independently of the irradiation period. Above 50°C a considerable increase in the A form was only observed for 10min NIR exposed samples. The hydrodynamic radius, (Rh), studied by the dynamic light scattering, showed drastic decrease with the increasing irradiation time. Principal components analysis (PCA) allowed to detect the spectral features correlated with the NIR effect and thermal stability of the DNA films. Obtained results strongly support the idea that the photoionization of water by NIR radiation in presence of DNA molecules is the main factor influencing on its physicochemical properties.

Reactivity of prehydrated electrons toward nucleobases and nucleotides in aqueous solution

DNA damage induced via dissociative attachment by low-energy electrons (0 to 20 eV) is well studied in both gas and condensed phases. However, the reactivity of ultrashort-lived prehydrated electrons ([Formula: see text]) with DNA components in a biologically relevant environment has not been fully explored to date. The electron transfer processes of [Formula: see text] to the DNA nucleobases G, A, C, and T and to nucleosides/nucleotides were investigated by using 7-ps electron pulse radiolysis coupled with pump-probe transient absorption spectroscopy in aqueous solutions. In contrast to previous results, obtained by using femtosecond laser pump-probe spectroscopy, we show that G and A cannot scavenge [Formula: see text] at concentrations of ≤50 mM. Observation of a substantial decrease of the initial yield of hydrated electrons ([Formula: see text]) and formation of nucleobase/nucleotide anion radicals at increasing nucleobase/nucleotide concentrations present direct evidence for the earliest step in reductive DNA damage by ionizing radiation. Our results show that [Formula: see text] is more reactive with pyrimidine than purine nucleobases/nucleotides with a reactivity order of T > C > A > G. In addition, analyses of transient signals show that the signal due to formation of the resulting anion radical directly correlates with the loss of the initial [Formula: see text] signal. Therefore, our results do not agree with the previously proposed dissociation of transient negative ions in nucleobase/nucleotide solutions within the timescale of these experiments. Moreover, in a molecularly crowded medium (for example, in the presence of 6 M phosphate), the scavenging efficiency of [Formula: see text] by G is significantly enhanced. This finding implies that reductive DNA damage by ionizing radiation depends on the microenvironment around [Formula: see text].

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.

Perspective: Structure and ultrafast dynamics of biomolecular hydration shells

The structure and function of biomolecules can be strongly influenced by their hydration shells. A key challenge is thus to determine the extent to which these shells differ from bulk water, since the structural fluctuations and molecular excitations of hydrating water molecules within these shells can cover a broad range in both space and time. Recent progress in theory, molecular dynamics simulations, and ultrafast vibrational spectroscopy has led to new and detailed insight into the fluctuations of water structure, elementary water motions, and electric fields at hydrated biointerfaces. Here, we discuss some central aspects of these advances, focusing on elementary molecular mechanisms and processes of hydration on a femto- to picosecond time scale, with some special attention given to several issues subject to debate.

An optically-gated AuNP-DNA protonic transistor

Bio-interface transistors, which manipulate the transportation of ions (i.e. protons), play an important role in bridging physical devices with biological functionalities, because electrical signals are carried by ions/protons in biological systems. All available ionic transistors use electrostatic gates to tune the ionic carrier density, which requires complicated interconnect wires. In contrast, an optical gate, which offers the advantages of remote control as well as multiple light wavelength selections, has never been explored for ionic devices. Here, we demonstrate a light-gated protonic transistor fabricated from an Au nanoparticle and DNA (AuNP-DNA) hybrid membrane. The device can be turned on and off completely by using light, with a high on/off current ratio of up to 2 orders of magnitude. Moreover, the device only responds to specific light wavelengths due to the plasmonic effect from the AuNPs, which enables the capability of wavelength selectivity. Our results open up new avenues for exploring remotely controlled ionic circuits, in vivo protonic switches, and other biomedical applications.

Prebiotic RNA polymerisation: energetics of nucleotide adsorption and polymerisation on clay mineral surfaces

A set of experimentally measured and internally consistent Gibbs free energies of binding between different model nucleotides and mineral surfaces is reported.

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.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.