Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs

Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

A novel, simplified method to prepare and preserve freeze-dried mouse sperm in plastic microtubes

Although freeze-drying sperm can save space, reduce maintenance costs, and facilitate the transportation of genetic samples, the current method requires breakable, custom-made, and expensive glass ampoules. In the present study, we developed a simple and economical method for collecting freeze-dried (FD) sperm using commercially available plastic microtubes. Mouse epididymal sperm suspensions were placed in 1.5 ml polypropylene tubes, frozen in liquid nitrogen, and dried in an acrylic freeze-drying chamber, after which they were closed under a vacuum. The drying duration did not differ between the microtube and glass ampoule methods (control); however, the sperm recovery rate was higher using the microtube method, and the physical damage to the sperm after rehydration was also reduced. Intracytoplasmic sperm injection (ICSI) using FD sperm stored in microtubes at -30°C yielded healthy offspring without reducing the success rate, even after 9 months of storage. Air infiltration into all microtubes stored at room temperature (RT) within 2 weeks of storage caused a drastic decrease in the fertilization rate of FD sperm; underwater storage did not prevent air infiltration. RT storage of FD sperm in microtubes for 1 week resulted in healthy offspring after ICSI (5-18%), but the addition of silica gel or CaCl2 did not improve the success rate. Our novel microtube method is currently the simplest and most effective method for treating FD sperm, contributing to the development of alternative low-cost approaches for preserving and transporting genetic resources.

Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs

Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

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.

Steric Hindrance On-Off Mass-Tagged Probe Set Enables Detection of Protein Homodimer in Living Cells

The major challenge in the detection of protein homodimers is that the identical monomers in a homodimer are indistinguishable using most conventional methods and cannot be sequentially recognized. In this study, a steric hindrance on-off mass-tagged probe set strategy was developed for the quantification of HER2 homodimer in living cells. The probe set contained a hindrance probe and a detection probe. The hindrance probe had a DNA dendrimer as a hindrance group to achieve the steric hindrance on-off function and thus the assignment of monomer identity. The detection probe contained a mass tag released for mass spectrometric quantification. Using the steric hindrance on-off mass-tagged probe set, the level of HER2 homodimer in various breast cancer cell lines was quantified. This is the first report to determine the quantity of protein homodimers, and the steric hindrance on-off probe set developed herein can facilitate the illustration of protein function in cancer.

Effect of Molecular Crowding on the Stability of RNA G-Quadruplexes with Various Numbers of Quartets and Lengths of Loops

G-Quadruplexes are noncanonical structures formed by guanine-rich regions of not only DNA but also RNA. RNA G-quadruplexes are widely present in the transcriptome as mRNAs and noncoding RNAs and take part in various essential functions in cells. Furthermore, stable RNA G-quadruplexes control the extent of biological functions, such as mRNA translation and antigen presentation. To understand and regulate the functions controlled by RNA G-quadruplexes in cellular environments, which are molecularly crowded, we would be required to investigate the stability of G-quadruplexes in molecular crowding. Here, we systematically investigated the thermodynamic stability of RNA G-quadruplexes with different numbers of G-quartets and lengths of loops. The molecular crowding conditions of polyethylene glycol with an average molecular weight of 200 (PEG200) were found to stabilize RNA G-quadruplexes with three and four G-quartets, while G-quadruplexes with two G-quartets did not exhibit any stabilization upon addition of PEG200. On the other hand, no difference in stabilization by PEG200 was observed among the G-quadruplexes with different loop lengths. Thermodynamic analysis of the RNA G-quadruplexes revealed more appropriate motifs for identifying G-quadruplex-forming sequences. The informatics analysis with new motifs demonstrated that the distributions of G-quadruplexes in human noncoding RNAs differed depending on the number of G-quartets. Therefore, RNA G-quadruplexes with different numbers of G-quartets may play different roles in response to environmental changes in cells.

Origin of Slow Solvation Dynamics in DNA: DAPI in Minor Groove of Dickerson-Drew DNA

The measurement and understanding of collective solvation dynamics in DNA have vital biological implications, as protein and ligand binding to DNA can be directly controlled by complex electrostatic interactions of anionic DNA and surrounding dipolar water, and ions. Time-resolved fluorescence Stokes shift (TRFSS) experiments revealed anomalously slow solvation dynamics in DNA much beyond 100 ps that follow either power-law or slow multiexponential decay over several nanoseconds. The origin of such dispersed dynamics remains difficult to understand. Here we compare results of TRFSS experiments to molecular dynamics (MD) simulations of well-known 4',6-diamidino-2-phenylindole (DAPI)/Dickerson-Drew DNA complex over five decades of time from 100 fs to 10 ns to understand the origin of such dispersed dynamics. We show that the solvation time-correlation function (TCF) calculated from 200 ns simulation trajectory (total 800 ns) captures most features of slow dynamics as measured in TRFSS experiments. Decomposition of TCF into individual components unravels that slow dynamics originating from dynamically coupled DNA-water motion, although contribution from coupled water-Na⁺ motion is non-negligible. The analysis of residence time of water molecules around the probe (DAPI) reveals broad distribution from ∼6 ps to ∼3.5 ns: Several (49 nos.) water molecules show residences time greater than 500 ps, of which at least 14 water molecules show residence times of more than 1 ns in the first solvation shell of DAPI. Most of these slow water molecules are found to occupy two hydration sites in the minor groove near DAPI binding site. The residence time of Na⁺, however, is found to vary within ∼17-120 ps. Remarkably, we find that freezing the DNA fluctuations in simulation eliminates slower dynamics beyond ∼100 ps, where water and Na⁺ dynamics become faster, although strong anticorrelation exists between them. These results indicate that primary origin of slow dynamics lies within the slow fluctuations of DNA parts that couple with nearby slow water and ions to control the dispersed collective solvation dynamics in DNA minor groove.

A New Generation of Minor-Groove-Binding-Heterocyclic Diamidines That Recognize G·C Base Pairs in an AT Sequence Context

We review the preparation of new compounds with good solution and cell uptake properties that can selectively recognize mixed A·T and G·C bp sequences of DNA. Our underlying aim is to show that these new compounds provide important new biotechnology reagents as well as a new class of therapeutic candidates with better properties and development potential than other currently available agents. In this review, entirely different ways to recognize mixed sequences of DNA by modifying AT selective heterocyclic cations are described. To selectively recognize a G·C base pair an H-bond acceptor must be incorporated with AT recognizing groups as with netropsin. We have used pyridine, azabenzimidazole and thiophene-N-methylbenzimidazole GC recognition units in modules crafted with both rational design and empirical optimization. These modules can selectively and strongly recognize a single G·C base pair in an AT sequence context. In some cases, a relatively simple change in substituents can convert a heterocyclic module from AT to GC recognition selectivity. Synthesis and DNA interaction results for initial example lead modules are described for single G·C base pair recognition compounds. The review concludes with a description of the initial efforts to prepare larger compounds to recognize sequences of DNA with more than one G·C base pairs. The challenges and initial successes are described along with future directions.

Microscopic picture of water-ethylene glycol interaction near a model DNA by computer simulation: Concentration dependence, structure, and localized thermodynamics

It is known that crowded molecular environment affects the structure, thermodynamics, and dynamics of macromolecules. Most of the previous works on molecular crowding have majorly focused on the behavior of the macromolecule with less emphasis on the behavior of the crowder and water molecules. In the current study, we have precisely focused on the behavior of the crowder, (ethylene glycol (EG)), salt ions, and water in the presence of a DNA with the increase of the EG concentration. We have probed the behavior of water and crowder using molecular dynamics (MD) simulation and by calculating localized thermodynamic properties. Our results show an interesting competition between EG and water molecules to make hydrogen bonds (H-bond) with DNA. Although the total number of H-bonds involving DNA with both EG and water remains essentially same irrespective of the increase in EG concentration, there is a proportional change in the H-bonding pattern between water-water, EG-EG, and EG-water near DNA and in bulk. At low concentrations of EG, the displacement of water molecules near DNA is relatively easy. However, the displacement of water becomes more difficult as the concentration of EG increases. The density of Na⁺ (Cl^-) near DNA increases (decreases) as the concentration of EG is increased. The density of Cl^- near Na⁺ increases with the increase in EG concentration. It was also found that the average free energy per water in the first solvation shell increases with the increase in EG concentration. Putting all these together, a microscopic picture of EG, water, salt interaction in the presence of DNA, as a function of EG concentration, has emerged.

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6-41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA. Communicated by Ramaswamy H. Sarma.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Drastic stability change of X-X mismatch in d(CXG) trinucleotide repeat disorders under molecular crowding condition

The trinucleotide repeat d(CXG) (X = A, C, G or T) is the most common sequence causing repeat expansion disorders. The formation of non-canonical structures, such as hairpin structures with X-X mismatches, has been proposed to affect gene expression and regulation, which are important in pathological studies of these devastating neurological diseases. However, little information is available regarding the thermodynamics of the repeat sequence under crowded cellular conditions where many non-canonical structures such as G-quadruplexes are highly stabilized, while duplexes are destabilised. In this study, we investigated the different stabilities of X-X mismatches in the context of internal d(CXG) self-complementary sequences in an environment with a high concentration of cosolutes to mimic the crowding conditions in cells. The stabilities of full-matched duplexes and duplexes with A-A, G-G, and T-T mismatched base pairs under molecular crowding conditions were notably decreased compared to under dilute conditions. However, the stability of the DNA duplex with a C-C mismatch base pair was only slightly destabilised. Investigating different stabilities of X-X mismatches in d(CXG) sequences is important for improving our understanding of the formation and transition of multiple non-canonical structures in trinucleotide repeat diseases, and may provide insights for pathological studies and drug development.

Unzipping of A-Form DNA-RNA, A-Form DNA-PNA, and B-Form DNA-DNA in the α-Hemolysin Nanopore

Unzipping of double-stranded nucleic acids by an electric field applied across a wild-type α-hemolysin (αHL) nanopore provides structural information about different duplex forms. In this work, comparative studies on A-form DNA-RNA duplexes and B-form DNA-DNA duplexes with a single-stranded tail identified significant differences in the blockage current and the unzipping duration between the two helical forms. We observed that the B-form duplex blocks the channel 1.9 ± 0.2 pA more and unzips ∼15-fold more slowly than an A-form duplex at 120 mV. We developed a model to describe the dependence of duplex unzipping on structure. We demonstrate that the wider A-form duplex (d = 2.4 nm) is unable to enter the vestibule opening of αHL on the cis side, leading to unzipping outside of the nanopore with higher residual current and faster unzipping times. In contrast, the smaller B-form duplexes (d = 2.0 nm) enter the vestibule of αHL, resulting in decreased current blockages and slower unzipping. We investigated the effects of varying the length of the single-stranded overhang, and studied A-form DNA-PNA duplexes to provide additional support for the proposed model. This study identifies key differences between A- and B-form duplex unzipping that will be important in the design of future probe-based methods for detecting DNA or RNA.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue?

Realistic all-atom simulation of biological systems requires accurate modeling of both the biomolecules and their ionic environment. Recently, ion nucleation phenomena leading to the rapid growth of KCl or NaCl clusters in the vicinity of biomolecular systems have been reported. To better understand this phenomenon, molecular dynamics simulations of KCl aqueous solutions at three (1.0, 0.25, and 0.10 M) concentrations were performed. Two popular water models (TIP3P and SPC/E) and two Lennard-Jones parameter sets (AMBER and Dang) were combined to produce a total of 80 ns of molecular dynamics trajectories. Results suggest that the use of the Dang cation Lennard-Jones parameters instead of those adopted by the AMBER force-field produces a more accurate description of the ionic solution. In the later case, formation of salt aggregates is probably indicative of an artifact resulting from misbalanced force-field parameters. Because similar results were obtained with two different water parameter sets, the simulations exclude a water model dependency in the formation of anomalous ionic clusters. Overall, the results strongly suggest that for accurate modeling of ions in biomolecular systems, great care should be taken in choosing balanced ionic parameters even when using the most popular force-fields. These results invite a reexamination of older data obtained using available force-fields and a thorough check of the quality of current parameters sets by performing simulations at finite (>0.25 M) instead of minimal salt conditions.

Competitive interaction of monovalent cations with DNA from 3D-RISM

The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.

Mechanistic insights into temperature-dependent regulation of the simple cyanobacterial hsp17 RNA thermometer at base-pair resolution

The cyanobacterial hsp17 ribonucleicacid thermometer (RNAT) is one of the smallest naturally occurring RNAT. It forms a single hairpin with an internal 1×3-bulge separating the start codon in stem I from the ribosome binding site (RBS) in stem II. We investigated the temperature-dependent regulation of hsp17 by mapping individual base-pair stabilities from solvent exchange nuclear magnetic resonance (NMR) spectroscopy. The wild-type RNAT was found to be stabilized by two critical CG base pairs (C14-G27 and C13-G28). Replacing the internal 1×3 bulge by a stable CG base pair in hsp17^rep significantly increased the global stability and unfolding cooperativity as evidenced by circular dichroism spectroscopy. From the NMR analysis, remote stabilization and non-nearest neighbour effects exist at the base-pair level, in particular for nucleotide G28 (five nucleotides apart from the side of mutation). Individual base-pair stabilities are coupled to the stability of the entire thermometer within both the natural and the stabilized RNATs by enthalpy–entropy compensation presumably mediated by the hydration shell. At the melting point the Gibbs energies of the individual nucleobases are equalized suggesting a consecutive zipper-type unfolding mechanism of the RBS leading to a dimmer-like function of hsp17 and switch-like regulation behaviour of hsp17^rep. The data show how minor changes in the nucleotide sequence not only offset the melting temperature but also alter the mode of temperature sensing. The cyanobacterial thermosensor demonstrates the remarkable adjustment of natural RNATs to execute precise temperature control.

Analyzing ion distributions around DNA: sequence-dependence of potassium ion distributions from microsecond molecular dynamics

Microsecond molecular dynamics simulations of B-DNA oligomers carried out in an aqueous environment with a physiological salt concentration enable us to perform a detailed analysis of how potassium ions interact with the double helix. The oligomers studied contain all 136 distinct tetranucleotides and we are thus able to make a comprehensive analysis of base sequence effects. Using a recently developed curvilinear helicoidal coordinate method we are able to analyze the details of ion populations and densities within the major and minor grooves and in the space surrounding DNA. The results show higher ion populations than have typically been observed in earlier studies and sequence effects that go beyond the nature of individual base pairs or base pair steps. We also show that, in some special cases, ion distributions converge very slowly and, on a microsecond timescale, do not reflect the symmetry of the corresponding base sequence.

Dangling ends perturb the stability of RNA duplexes responsive to surrounding conditions

Unpaired terminal nucleotides (dangling ends) occur in various biologically important RNA structures. We studied the thermal stability of RNA duplexes with dangling ends under conditions that mimic those in cells. Dangling ends of one or two nucleotides stabilized a duplex up to approximately 2.7 kcal mol(-1) in the absence of cosolutes. RNA duplexes with dangling purine nucleotides were more stable than those with pyrimidine nucleotides. Interestingly, in the presence of various cosolutes, RNA duplexes with purine dangling ends were significantly destabilized, although those with pyrimidine dangling ends were destabilized slightly. For example, in 30 wt % poly(ethylene glycol), stabilization resulting from adenine dangling ends was reduced by 1.4 kcal mol(-1) . Our quantitative analyses also showed that the number of water molecules bound to the dangling ends in an aqueous solution was independent of the nucleotide type but dependent on the stability of the dangling-end region. It has been considered that dangling ends stabilize helices; however, our results suggest that the stabilization is responsive to the surrounding conditions.

DNA conformational flexibility study using phosphate backbone neutralization model

Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.

Ion counting from explicit-solvent simulations and 3D-RISM

The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.

MD and NMR analyses of choline and TMA binding to duplex DNA: on the origins of aberrant sequence-dependent stability by alkyl cations in aqueous and water-free solvents

It has been known for decades that alkylammonium ions, such as tetramethyl ammonium (TMA), alter the usual correlation between DNA GC-content and duplex stability. In some cases it is even possible for an AT-rich duplex to be more stable than a GC-rich duplex of the same length. There has been much speculation regarding the origin of this aberration in sequence-dependent DNA duplex stability, but no clear resolution. Using a combination of molecular dynamics simulations and NMR spectroscopy we demonstrate that choline (2-hydroxy-N,N,N-trimethylethanaminium) and TMA are preferentially localized in the minor groove of DNA duplexes at A·T base pairs and these same ions show less pronounced localization in the major groove compared to what has been demonstrated for alkali and alkali earth metal ions. Furthermore, free energy calculations show that single-stranded GC-rich sequences exhibit more favorable solvation by choline than single-stranded AT-rich sequences. The sequence-specific nature of choline and TMA binding provides a rationale for the enhanced stability of AT-rich sequences when alkyl-ammonium ions are used as the counterions of DNA. Our combined theoretical and experimental study provides one of the most detailed pictures to date of cations localized along DNA in the solution state, and provides insights that go beyond understanding alkyl-ammonium ion binding to DNA. In particular, because choline and TMA bind to DNA in a manner that is found to be distinct from that previously reported for Na(+), K(+), Mg(2+), and Ca(2+), our results reveal the important but underappreciated role that most other cations play in sequence-specific duplex stability.

Fast-SAXS-pro: a unified approach to computing SAXS profiles of DNA, RNA, protein, and their complexes

A generalized method, termed Fast-SAXS-pro, for computing small angle x-ray scattering (SAXS) profiles of proteins, nucleic acids, and their complexes is presented. First, effective coarse-grained structure factors of DNA nucleotides are derived using a simplified two-particle-per-nucleotide representation. Second, SAXS data of a 18-bp double-stranded DNA are measured and used for the calibration of the scattering contribution from excess electron density in the DNA solvation layer. Additional test on a 25-bp DNA duplex validates this SAXS computational method and suggests that DNA has a different contribution from its hydration surface to the total scattering compared to RNA and protein. To account for such a difference, a sigmoidal function is implemented for the treatment of non-uniform electron density across the surface of a protein/nucleic-acid complex. This treatment allows differential scattering from the solvation layer surrounding protein/nucleic-acid complexes. Finally, the applications of this Fast-SAXS-pro method are demonstrated for protein/DNA and protein/RNA complexes.

Quantitative analyses of nucleic Acid stability under the molecular crowding condition induced by cosolutes

A variety of biomolecules, including nucleic acids, proteins, polysaccharides, and other soluble and insoluble low-molecular weight components, are present in living cells. These molecules occupy a significant fraction of the cellular volume (up to 40%), resulting in a highly crowded intracellular environment. This situation is referred to as molecular crowding. Although the thermodynamic stabilities of DNA structures are known to be altered in a crowded environment, less is known about the behavior of nucleic acids and their interactions with cations and water molecules under such conditions. This unit describes methods that can be used to quantitatively analyze the molecular crowding effects caused by cosolutes on the thermodynamic stability, hydration, and cation binding of nucleic acid structures.

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Hydration changes upon DNA folding studied by osmotic stress experiments

The thermal stability of nucleic acid structures is perturbed under the conditions that mimic the intracellular environment, typically rich in inert components and under osmotic stress. We now describe the thermodynamic stability of DNA oligonucleotide structures in the presence of high background concentrations of neutral cosolutes. Small cosolutes destabilize the basepair structures, and the DNA structures consisting of the same nearest-neighbor composition show similar thermodynamic parameters in the presence of various types of cosolutes. The osmotic stress experiments reveal that water binding to flexible loops, unstable mismatches, and an abasic site upon DNA folding are almost negligible, whereas the binding to stable mismatch pairs is significant. The studies using the basepair-mimic nucleosides and the peptide nucleic acid suggest that the sugar-phosphate backbone and the integrity of the basepair conformation make important contributions to the binding of water molecules to the DNA bases and helical grooves. The study of the DNA hydration provides the basis for understanding and predicting nucleic acid structures in nonaqueous solvent systems.

Use of nucleic Acid analogs for the study of nucleic Acid interactions

Unnatural nucleosides have been explored to expand the properties and the applications of oligonucleotides. This paper briefly summarizes nucleic acid analogs in which the base is modified or replaced by an unnatural stacking group for the study of nucleic acid interactions. We also describe the nucleoside analogs of a base pair-mimic structure that we have examined. Although the base pair-mimic nucleosides possess a simplified stacking moiety of a phenyl or naphthyl group, they can be used as a structural analog of Watson-Crick base pairs. Remarkably, they can adopt two different conformations responding to their interaction energies, and one of them is the stacking conformation of the nonpolar aromatic group causing the site-selective flipping of the opposite base in a DNA double helix. The base pair-mimic nucleosides can be used to study the mechanism responsible for the base stacking and the flipping of bases out of a nucleic acid duplex.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.

Monomorphic RNA G-quadruplex and polymorphic DNA G-quadruplex structures responding to cellular environmental factors

We systematically and quantitatively investigated the structure and thermodynamics of G-quadruplexes of RNAs and corresponding DNAs of the same sequences under molecular crowding conditions that mimic the high osmotic stress induced by the numerous molecules inside of living cells. Structural analyses demonstrated that various telomere RNA sequences folded into parallel-stranded G-quadruplexes in a manner independent of the surrounding conditions with different cations under both dilute and molecular crowding conditions. In contrast, DNA G-quadruplexes showed structural polymorphism. Moreover, we demonstrated that the G-quadruplexes of the RNA sequences were more stable than those of the same DNA sequences. These results show that a single and robust RNA G-quadruplex structure can exist in a manner independent of the sequence and surrounding conditions. To confirm this, we studied a guanine-rich sequence located in the 5'-untranslated region of human bcl-2 mRNA that is thought to play a role in translation. The results revealed a stable parallel G-quadruplex that formed under all conditions tested. For example, a bcl-RNA G-quadruplex in the presence of 5 mM KCl [free energy change at 25 degrees C (DeltaG degrees (25)) of -5.42 kcal/mol] was more stable than its corresponding DNA G-quadruplex (DeltaG degrees (25) = -2.31 kcal/mol). Our results further indicated that water molecules binding to the 2'-OH group of RNA G-quadruplexes play a critical role in their formation and stability.

Effect of high hydrostatic pressure on hydration and activity of ribozymes

Formation and stabilization of RNA structure in the cell depends on its interaction with solvent and metal ions. High hydrostatic pressure (HHP) is a convenient tool in an analysis of the role of small molecules in the structure stabilization of biological macromolecules. Analysis of HHP effect and various concentrations of ions showed that water induce formation of the active ribozyme structure. So, it is clear that water is the driving force of conformational changes of nucleic acid.

Assessing the performance of implicit solvation models at a nucleic acid surface

Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.

Comparison of DNA hydration patterns obtained using two distinct computational methods, molecular dynamics simulation and three-dimensional reference interaction site model theory

Because proteins and DNA interact with each other and with various small molecules in the presence of water molecules, we cannot ignore their hydration when discussing their structural and energetic properties. Although high-resolution crystal structure analyses have given us a view of tightly bound water molecules on their surface, the structural data are still insufficient to capture the detailed configurations of water molecules around the surface of these biomolecules. Thanks to the invention of various computational algorithms, computer simulations can now provide an atomic view of hydration. Here, we describe the apparent patterns of DNA hydration calculated by using two different computational methods: Molecular dynamics (MD) simulation and three-dimensional reference interaction site model (3D-RISM) theory. Both methods are promising for obtaining hydration properties, but until now there have been no thorough comparisons of the calculated three-dimensional distributions of hydrating water. This rigorous comparison showed that MD and 3D-RISM provide essentially similar hydration patterns when there is sufficient sampling time for MD and a sufficient number of conformations to describe molecular flexibility for 3D-RISM. This suggests that these two computational methods can be used to complement one another when evaluating the reliability of the calculated hydration patterns.

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.

Observation of a dynamic crossover in RNA hydration water which triggers a dynamic transition in the biopolymer

High-resolution quasielastic neutron scattering spectroscopy was used to measure H2O and D2O hydrated RNA samples. The contribution of scattering from RNA was subtracted out by taking the difference of the signals between the two samples. The measurements were made at a series of temperatures from 270 K down to 180 K. The relaxing-cage model was used to analyze the difference quasielastic spectra. We observed clear evidence of a fragile-to-strong dynamic crossover (FSC) at TL=220 K in RNA hydration water. We further show that the mean-square displacements of the hydrogen atoms in both RNA and its hydration water exhibit a sharp change in slope at approximately the same temperature 220 K. This latter fact suggests that the dynamic transition in RNA is triggered by the abrupt change of mobility of the hydration water at its FSC temperature.