Understanding the Relative Flexibility of RNA and DNA Duplexes: Stretching and Twist-Stretch Coupling

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Enhancing Neuronal Cell Uptake of Therapeutic Nucleic Acids with Tetrahedral DNA Nanostructures

The successful translation of therapeutic nucleic acids (NAs) for the treatment of neurological disorders depends on their safe and efficient delivery to neural cells, in particular neurons. DNA nanostructures can be a promising NAs delivery vehicle. Nonetheless, the potential of DNA nanostructures for neuronal cell delivery of therapeutic NAs is unexplored. Here, tetrahedral DNA nanostructures (TDN) as siRNA delivery scaffolds to neuronal cells, exploring the influence of functionalization with two different reported neuronal targeting ligands: C4-3 RNA aptamer and Tet1 peptide are investigated. Nanostructures are characterized in vitro, as well as in silico using molecular dynamic simulations to better understand the overall TDN structural stability. Enhancement of neuronal cell uptake of TDN functionalized with the C4-3 Aptamer (TDN-Apt), not only in neuronal cell lines but also in primary neuronal cell cultures is demonstrated. Additionally, TDN and TDN-Apt nanostructures carrying siRNA are shown to promote silencing in a process aided by chloroquine-induced endosomal disruption. This work presents a thorough workflow for the structural and functional characterization of the proposed TDN as a nano-scaffold for neuronal delivery of therapeutic NAs and for targeting ligands evaluation, contributing to the future development of new neuronal drug delivery systems based on DNA nanostructures.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study

We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.

5-Methyl-cytosine stabilizes DNA but hinders DNA hybridization revealed by magnetic tweezers and simulations

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.

Gold nanoparticle design for RNA compaction

RNA-based therapeutics hold a great promise in treating a variety of diseases. However, double-stranded RNAs (dsRNAs) are inherently unstable, highly charged, and stiff macromolecules that require a delivery vehicle. Cationic ligand functionalized gold nanoparticles (AuNPs) are able to compact nucleic acids and assist in RNA delivery. Here, we use large-scale all-atom molecular dynamics simulations to show that correlations between ligand length, metal core size, and ligand excess free volume control the ability of nanoparticles to bend dsRNA far below its persistence length. The analysis of ammonium binding sites showed that longer ligands that bind deep within the major groove did not cause bending. By limiting ligand length and, thus, excess free volume, we have designed nanoparticles with controlled internal binding to RNA's major groove. NPs that are able to induce RNA bending cause a periodic variation in RNA's major groove width. Density functional theory studies on smaller models support large-scale simulations. Our results are expected to have significant implications in packaging of nucleic acids for their applications in nanotechnology and gene delivery.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Accurate Sequence-Dependent Coarse-Grained Model for Conformational and Elastic Properties of Double-Stranded DNA

We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar, at the base moiety, and at the phosphate group, respectively. The sequence dependence is included by considering a step-dependent parametrization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from the literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist–stretch and twist–bend couplings, the persistence length, and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of single-molecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines. Particularly, an implementation of the model in LAMMPS is made available on an online repository to ease its usage within the DNA research community.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.

Systematic mapping of DNAzymes research from 1995 to 2019

DNAzymes (catalytic DNA) have gained significant diagnostic and therapeutic applications with increasing research output over the years. Functional oligonucleotides are used as molecular recognition elements within biosensors for detection of analytes and viral infections such as SARS-CoV-2. DNAzymes are also applied for silencing and regulating cancer specific genes. However, there has not been any report on systematic analysis to track research status, reveal hotspots, and map knowledge in this field. Therefore, in the present study, research articles on DNAzymes from 1995 to 2019 were extracted from Web of Science (SCI-Expanded) after which, 1037 articles were imported into Rstudio (version 3.6.2) and analysed accordingly. The highest number of articles was published in 2019 (n = 138), while the least was in 1995 (n = 1). The articles were published across 216 journals by 2344 authors with 2337 multi-author and 7 single authors. The most prolific authors were Li Y (n = 47), Liu J (n = 46), Wang L (n = 33), Willner I (n = 33) and Zhang L (n = 33). The top three most productive countries were China (n = 2018), USA (n = 447) and Canada (n = 251). The most productive institutions were Hunan University, China (n = 141), University of Illinois, USA (n = 139) and Fuzhou University, China (n = 101). Despite the increasing interest in this field, international collaborations between institutions were very low which requires immediate attention to mitigate challenges such as limited funding, access to facilities, and existing knowledge gap.

Multivalent Cations Reverse the Twist-Stretch Coupling of RNA

When stretched, both DNA and RNA duplexes change their twist angles through twist-stretch coupling. The coupling is negative for DNA but positive for RNA, which is not yet completely understood. Here, our magnetic tweezers experiments show that the coupling of RNA reverses from positive to negative by multivalent cations. Combining with the previously reported tension-induced negative-to-positive coupling reversal of DNA, we propose a unified mechanism of the couplings of both RNA and DNA based on molecular dynamics simulations. Two deformation pathways are competing when stretched: shrinking the radius causes positive couplings but widening the major groove causes negative couplings. For RNA whose major groove is clamped by multivalent cations and canonical DNA, their radii shrink when stretched, thus exhibiting positive couplings. For elongated DNA whose radius already shrinks to the minimum and canonical RNA, their major grooves are widened when stretched, thus exhibiting negative couplings.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

Nucleic Acid Identity, Structure, and Flexibility Affect the Electrochemical Signal of Tethered Redox Molecules upon Biopolymer Collapse

We demonstrate that cation condensation can induce the collapse of surface-bound nucleic acids and that the electrochemical signal from a tethered redox molecule (methylene blue) upon collapse reports on nucleic acid identity, structure, and flexibility. Furthermore, the correlation of the electrochemical signal and structure is consistent with theoretical considerations of nucleic acid collapse. Changes in solution dielectric permittivity or the concentration of trivalent cations cause the structure of nucleic acids to become more compact due to an increase in attractive electrostatic interactions between the charged biopolymer backbone and multivalent ions in the solution. Consequently, the compaction of nucleic acids results in a change in the dynamics and location of the terminally appended redox marker, which is reflected in the faradaic current measured using cyclic voltammetry. In comparison to ssDNA, nucleic acid duplexes (dsDNA, DNA/peptide nucleic acid, and dsRNA) require nucleic-acid-composition-specific solution conditions for the collapse to occur. Moreover, the magnitude of current increase observed after the collapse is different for each nucleic structure, and we find here that these changes are dictated by physical parameters of the nucleic acids including the axial charge spacing and the periodicity of the helix. The work here aims to provide quantitative and predicative measures of the effects of the nucleic acid structure on the electrochemical signal produced from distal-end appended redox markers. This architecture is commonly employed in functional nucleic acid sensors and a better understanding of structure-to-signal correlations will enable the rational design of sensitive sensing architectures.

A molecular view of DNA flexibility

DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity

RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Multiscale modelling reveals higher charge transport efficiencies of DNA relative to RNA independent of mechanism

In this study, we compare the charge transport properties of multiple double-stranded (ds)RNA sequences with corresponding dsDNA sequences. Recent studies have presented a contradictory picture of relative charge transport efficiencies in A-form DNA : RNA hybrids and dsDNA. Using a multiscale modelling framework, we compute conductance of dsDNA and dsRNA using Landauer formalism in the coherent limit and Marcus-Hush theory in the incoherent limit. We find that dsDNA conducts better than dsRNA in both the charge transport regimes. Our analysis shows that the structural differences in the twist angle and slide of dsDNA and dsRNA are the main reasons behind the higher conductance of dsDNA in the incoherent hopping regime. In the coherent limit however, for the same base pair length, the conductance of dsRNA is higher than that of dsDNA for the morphologies where dsRNA has a smaller end-to-end length relative to that of dsDNA.

Structure-mechanics statistical learning unravels the linkage between local rigidity and global flexibility in nucleic acids

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression, but tracing their molecular origin has been difficult due to the structural and chemical complexity. We posit that concepts from machine learning can help to tackle this long-standing challenge. Here, we demonstrate the feasibility and advantage of this strategy through developing a structure-mechanics statistical learning scheme to elucidate how local rigidity in double-stranded (ds)DNA and dsRNA may lead to their global flexibility in bend, stretch, and twist. Specifically, the mechanical parameters in a heavy-atom elastic network model are computed from the trajectory data of all-atom molecular dynamics simulation. The results show that the inter-atomic springs for backbone and ribose puckering in dsRNA are stronger than those in dsDNA, but are similar in strengths for base-stacking and base-pairing. Our analysis shows that the experimental observation of dsDNA being easier to bend but harder to stretch than dsRNA comes mostly from the respective B- and A-form topologies. The computationally resolved composition of local rigidity indicates that the flexibility of both nucleic acids is mostly due to base-stacking. But for properties like twist-stretch coupling, backbone springs are shown to play a major role instead. The quantitative connection between local rigidity and global flexibility sets foundation for understanding how local binding and chemical modification of genetic materials effectuate longer-ranged regulatory signals.

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression. We devise structural mechanics statistical learning method to reveal their molecular origin in terms of chemical interactions.

Crystal structures of an unmodified bacterial tRNA reveal intrinsic structural flexibility and plasticity as general properties of unbound tRNAs

Ubiquitous across all domains of life, tRNAs constitute an essential component of cellular physiology, carry out an indispensable role in protein synthesis, and have been historically the subject of a wide range of biochemical and biophysical studies as prototypical folded RNA molecules. Although conformational flexibility is a well-established characteristic of tRNA structure, it is typically regarded as an adaptive property exhibited in response to an inducing event, such as the binding of a tRNA synthetase or the accommodation of an aminoacyl-tRNA into the ribosome. In this study, we present crystallographic data of a tRNA molecule to expand on this paradigm by showing that structural flexibility and plasticity are intrinsic properties of tRNAs, apparent even in the absence of other factors. Based on two closely related conformations observed within the same crystal, we posit that unbound tRNAs by themselves are flexible and dynamic molecules. Furthermore, we demonstrate that the formation of the T-loop conformation by the tRNA TΨC stem-loop, a well-characterized and classic RNA structural motif, is possible even in the absence of important interactions observed in fully folded tRNAs.

Opposite Effects of High-Valent Cations on the Elasticities of DNA and RNA Duplexes Revealed by Magnetic Tweezers

We report that trivalent cobalt hexammine cations decrease the persistence length, stretching modulus, helical density, and size of plectonemes formed under torque of DNA but increase those of RNA. Divalent magnesium cations, however, decrease the persistence lengths, contour lengths, and sizes of plectonemes while increasing the helical densities of both DNA and RNA. The experimental results are explained by different binding modes of the cations on DNA and RNA in our all-atom molecular dynamics simulations. The significant variations of the helical densities and structures of DNA and RNA duplexes induced by high-valent cations may affect interactions of the duplexes with proteins.

Molecular chains under tension: Thermal and mechanical activation of statistically interacting extension and contraction particles

This work introduces a methodology for the statistical mechanical analysis of polymeric chains under tension controlled by optical or magnetic tweezers at thermal equilibrium with an embedding fluid medium. The response of single bonds between monomers or of entire groups of monomers to tension is governed by the activation of statistically interacting particles representing quanta of extension or contraction. This method of analysis is capable of describing thermal unbending of the freely jointed or wormlike chain kind, linear or nonlinear contour elasticity, and structural transformations including effects of cooperativity. The versatility of this approach is demonstrated in an application to double-stranded DNA undergoing torsionally unconstrained stretching across three regimes of mechanical response including an overstretching transition. The three-regime force-extension characteristic, derived from a single free-energy expression, accurately matches empirical evidence.

Sequence-dependent mechanical properties of double-stranded RNA

The mechanical properties of double-stranded RNA (dsRNA) are involved in many of its biological functions and are relevant for future nanotechnology applications. DsRNA must tightly bend to fit inside viral capsids or deform upon the interaction with proteins that regulate gene silencing or the immune response against viral attacks. However, the question of how the nucleotide sequence affects the global mechanical properties of dsRNA has so far remained largely unexplored. Here, we have employed state-of-the-art atomistic molecular dynamics simulations to unveil the mechanical response of different RNA duplexes to an external force. Our results reveal that, similarly to dsDNA, the mechanical properties of dsRNA are highly sequence-dependent. However, we find that the nucleotide sequence affects in a strikingly different manner the stretching and twisting response of RNA and DNA duplexes under force. We find that the elastic response of dsRNA is dominated by the local high flexibility of pyrimidine-purine steps. Moreover, the flexibility of pyrimidine-purine steps is independent of the sequence context, and the global flexibility of the duplex reasonably scales with the number of this kind of base-pair dinucleotides. We conclude that disparities of the mechanical response of dinucleotides are responsible for the differences observed in the mechanical properties of RNA and DNA duplexes.

Salt Dependence of A-Form RNA Duplexes: Structures and Implications

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

3D structure stability of the HIV-1 TAR RNA in ion solutions: A coarse-grained model study

As an extremely common structural motif, RNA hairpins with bulge loops [e.g., the human immunodeficiency virus type 1 (HIV-1) transactivation response (TAR) RNA] can play essential roles in normal cellular processes by binding to proteins and small ligands, which could be very dependent on their three-dimensional (3D) structures and stability. Although the structures and conformational dynamics of the HIV-1 TAR RNA have been extensively studied, there are few investigations on the thermodynamic stability of the TAR RNA, especially in ion solutions, and the existing studies also have some divergence on the unfolding process of the RNA. Here, we employed our previously developed coarse-grained model with implicit salt to predict the 3D structure, stability, and unfolding pathway for the HIV-1 TAR RNA over a wide range of ion concentrations. As compared with the extensive experimental/theoretical results, the present model can give reliable predictions on the 3D structure stability of the TAR RNA from the sequence. Based on the predictions, our further comprehensive analyses on the stability of the TAR RNA as well as its variants revealed that the unfolding pathway of an RNA hairpin with a bulge loop is mainly determined by the relative stability between different states (folded state, intermediate state, and unfolded state) and the strength of the coaxial stacking between two stems in folded structures, both of which can be apparently modulated by the ion concentrations as well as the sequences.

Atomic structures of RNA nanotubes and their comparison with DNA nanotubes

We present a computational framework to model RNA based nanostructures and study their microscopic structures. We model hexagonal nanotubes made of 6 dsRNA (RNTs) connected by double crossover (DX) at different positions. Using several hundred nano-second (ns) long all-atom molecular dynamics simulations, we study the atomic structure, conformational change and elastic properties of RNTs in the presence of explicit water and ions. Based on several structural quantities such as root mean square deviation (RMSD) and root mean square fluctuation (RMSF), we find that the RNTs are almost as stable as DNA nanotubes (DNTs). Although the central portion of the RNTs maintain its cylindrical shape, both the terminal regions open up to give rise to a gating like behavior which can play a crucial role in drug delivery. From the bending angle distribution, we observe that the RNTs are more flexible than DNTs. The calculated persistence length of the RNTs is in the micron range which is an order of magnitude higher than that of a single dsRNA. The stretch modulus of the RNTs from the contour length distribution is in the range of 4-7 nN depending on the sequence. The calculated persistence length and stretch modulus are in the same range of values as in the case of DNTs. To understand the structural properties of RNTs at the individual base-pair level we have also calculated all the helicoidal parameters and analyzed the relative flexibility and rigidity of RNTs having a different sequence. These findings emphasized the fascinating properties of RNTs which will expedite further theoretical and experimental studies in this field.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.

Identification of Myricetin as an Ebola Virus VP35-Double-Stranded RNA Interaction Inhibitor through a Novel Fluorescence-Based Assay

Ebola virus (EBOV) is a filovirus that causes a severe and rapidly progressing hemorrhagic syndrome; a recent epidemic illustrated the urgent need for novel therapeutic agents because no drugs have been approved for treatment of Ebola virus. A key contribution to the high lethality observed during EBOV outbreaks comes from viral evasion of the host antiviral innate immune response in which viral protein VP35 plays a crucial role, blocking interferon type I production, first by masking the viral double-stranded RNA (dsRNA) and preventing its detection by the pattern recognition receptor RIG-I. Aiming to identify inhibitors of the interaction of VP35 with the viral dsRNA, counteracting the VP35 viral innate immune evasion, we established a new methodology for high-yield recombinant VP35 (rVP35) expression and purification and a novel and robust fluorescence-based rVP35-RNA interaction assay ( Z' factor of 0.69). Taking advantage of such newly established methods, we screened a small library of Sardinian natural extracts, identifying Limonium morisianum as the most potent inhibitor extract. A bioguided fractionation led to the identification of myricetin as the component that can inhibit rVP35-dsRNA interaction with an IC₅₀ value of 2.7 μM. Molecular docking studies showed that myricetin interacts with the highly conserved region of the VP35 RNA binding domain, laying the basis for further structural optimization of potent inhibitors of VP35-dsRNA interaction.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Radial distribution function of semiflexible oligomers with stretching flexibility

The radial distribution of the end-to-end distance R_ee is crucial for quantifying the global size and flexibility of a linear polymer. For semiflexible polymers, several analytical formulas have been derived for the radial distribution of R_ee ignoring the stretching flexibility. However, for semiflexible oligomers, such as DNA or RNA, the stretching flexibility can be rather pronounced and can significantly affect the radial distribution of R_ee. In this study, we obtained an extended formula that includes the stretch modulus to describe the distribution of R_ee for semiflexible oligomers on the basis of previous formulas for semiflexible polymers without stretching flexibility. The extended formula was validated by extensive Monte Carlo simulations over wide ranges of the stretch modulus and persistence length, as well as all-atom molecular dynamics simulations of short DNAs and RNAs. Additionally, our analyses showed that the effect of stretching flexibility on the distribution of R_ee becomes negligible for DNAs longer than ∼130 base pairs and RNAs longer than ∼240 base pairs.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

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.

Precise Antibody-Independent m6A Identification via 4SedTTP-Involved and FTO-Assisted Strategy at Single-Nucleotide Resolution

Innovative detection techniques to achieve precise m6A distribution within mammalian transcriptome can advance our understanding of its biological functions. We specifically introduced the atom-specific replacement of oxygen with progressively larger atoms (sulfur and selenium) at 4-position of deoxythymidine triphosphate to weaken its ability to base pair with m6A, while maintaining A-T* base pair virtually the same as the natural one. 4SedTTP turned out to be an outstanding candidate that endowed m6A with a specific signature of RT truncation, thereby making this "RT-silent" modification detectable with the assistance of m6A demethylase FTO through next-generation sequencing. This antibody-independent, 4SedTTP-involved and FTO-assisted strategy is applicable in m6A identification, even for two closely gathered m6A sites, within an unknown region at single-nucleotide resolution.

Divalent Ion-Mediated DNA-DNA Interactions: A Comparative Study of Triplex and Duplex

Ion-mediated interaction between DNAs is essential for DNA condensation, and it is generally believed that monovalent and nonspecifically binding divalent cations cannot induce the aggregation of double-stranded (ds) DNAs. Interestingly, recent experiments found that alkaline earth metal ions such as Mg²⁺ can induce the aggregation of triple-stranded (ts) DNAs, although there is still a lack of deep understanding of the surprising findings at the microscopic level. In this work, we employed all-atom dynamic simulations to directly calculate the potentials of mean force (PMFs) between tsDNAs, between dsDNAs, and between tsDNA and dsDNA in Mg²⁺ solutions. Our calculations show that the PMF between tsDNAs is apparently attractive and becomes more strongly attractive at higher [Mg²⁺], although the PMF between dsDNAs cannot become apparently attractive even at high [Mg²⁺]. Our analyses show that Mg²⁺ internally binds into grooves and externally binds to phosphate groups for both tsDNA and dsDNA, whereas the external binding of Mg²⁺ is much stronger for tsDNA. Such stronger external binding of Mg²⁺ for tsDNA favors more apparent ion-bridging between helices than for dsDNA. Furthermore, our analyses illustrate that bridging ions, as a special part of external binding ions, are tightly and positively coupled to ion-mediated attraction between DNAs.

Resolving the Conformational Dynamics of DNA with Ångstrom Resolution by Pulsed Electron-Electron Double Resonance and Molecular Dynamics

Pulsed electron-electron double resonance (PELDOR/DEER) experiments of nucleic acids with rigid spin labels provide highly accurate distance and orientation information. Here we combine PELDOR experiments with molecular dynamics (MD) simulations to arrive at an atomistic view of the conformational dynamics of DNA. The MD simulations closely reproduce the PELDOR time traces, and demonstrate that bending, in addition to twist-stretch motions, underpin the sub-μs dynamics of DNA. PELDOR experiments correctly rank DNA force fields and resolve subtle differences in the conformational ensembles of nucleic acids, on the order of 1-2 Å. Long-range distance and angle measurements with rigid spin labels provide critical input for the refinement of computer models and the elucidation of the structure and dynamics of complex biomolecules.