Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

In-Cell Stability Prediction of RNA/DNA Hybrid Duplexes for Designing Oligonucleotides Aimed at Therapeutics

In cells, the formation of RNA/DNA hybrid duplexes regulates gene expression and modification. The environment inside cellular organelles is heterogeneously crowded with high concentrations of biomolecules that affect the structure and stability of RNA/DNA hybrid duplexes. However, the detailed environmental effects remain unclear. Therefore, the mechanistic details of the effect of such molecular crowding were investigated at the molecular level by using thermodynamic and nuclear magnetic resonance analyses, revealing structure-dependent destabilization of the duplexes under crowded conditions. The transition from B- to A-like hybrid duplexes due to a change in conformation of the DNA strand guided by purine-pyrimidine asymmetry significantly increased the hydration number, which resulted in greater destabilization by the addition of cosolutes. By quantifying the individual contributions of environmental factors and the bulk structure of the duplex, we developed a set of parameters that predict the stability of hybrid duplexes with conformational dissimilarities under diverse crowding conditions. A comparison of the effects of environmental conditions in living cells and in vitro crowded solutions on hybrid duplex formation using the Förster resonance energy transfer technique established the applicability of our parameters to living cells. Moreover, our derived parameters can be used to estimate the efficiency of transcriptional inhibition, genome editing, and silencing techniques in cells. This supports the usefulness of our parameters for the visualization of cellular mechanisms of gene expression and the development of nucleic acid-based therapeutics targeting different cells.

Improving Precise Genome Editing Using Donor DNA/gRNA Hybrid Duplex Generated by Complementary Bases

In precise genome editing, site-specific DNA double-strand breaks (DSBs) induced by the CRISPR/Cas9 system are repaired via homology-directed repair (HDR) using exogenous donor DNA templates. However, the low efficiency of HDR-mediated genome editing is a barrier to widespread use. In this study, we created a donor DNA/guide RNA (gRNA) hybrid duplex (DGybrid) that was composed of sequence-extended gRNA and single-stranded oligodeoxynucleotide (ssODN) combined with complementary bases without chemical modifications to increase the concentration of donor DNA at the cleavage site. The efficiency of genome editing using DGybrid was evaluated in Saccharomyces cerevisiae. The results show a 1.8-fold (from 35% to 62%) improvement in HDR-mediated editing efficiency compared to genome editing in which gRNA and donor DNA were introduced separately. In addition, analysis of the nucleic acid introduction efficiency using flow cytometry indicated that both RNA and ssODNs are efficiently incorporated into cells together by using the DNA/RNA hybrid. Our technique would be preferred as a universal and concise tool for improving the efficiency of HDR-mediated genome editing.

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

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

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.

Why are Hoogsteen base pairs energetically disfavored in A-RNA compared to B-DNA?

A(syn)-U/T and G(syn)-C⁺ Hoogsteen (HG) base pairs (bps) are energetically more disfavored relative to Watson–Crick (WC) bps in A-RNA as compared to B-DNA by >1 kcal/mol for reasons that are not fully understood. Here, we used NMR spectroscopy, optical melting experiments, molecular dynamics simulations and modified nucleotides to identify factors that contribute to this destabilization of HG bps in A-RNA. Removing the 2′-hydroxyl at single purine nucleotides in A-RNA duplexes did not stabilize HG bps relative to WC. In contrast, loosening the A-form geometry using a bulge in A-RNA reduced the energy cost of forming HG bps at the flanking sites to B-DNA levels. A structural and thermodynamic analysis of purine-purine HG mismatches reveals that compared to B-DNA, the A-form geometry disfavors syn purines by 1.5–4 kcal/mol due to sugar-backbone rearrangements needed to sterically accommodate the syn base. Based on MD simulations, an additional penalty of 3–4 kcal/mol applies for purine-pyrimidine HG bps due to the higher energetic cost associated with moving the bases to form hydrogen bonds in A-RNA versus B-DNA. These results provide insights into a fundamental difference between A-RNA and B-DNA duplexes with important implications for how they respond to damage and post-transcriptional modifications.

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.

Thermostability Trends of TNA:DNA Duplexes Reveal Strong Purine Dependence

The development of high fidelity polymerases and streamlined synthesis of threose nucleic acid (TNA) triphosphates and phosphoramidites has made TNA accessible as a motif for generating nuclease-resistant high-affinity aptamers, antisense oligos, and synthetic genetic biopolymers. Little is known, however, about the thermostability trends of TNA:DNA duplexes. Here we investigate the thermostability of 14 TNA:DNA duplexes with the goal of elucidating the fundamental factors governing TNA:DNA duplex stability. We find that purine content in TNA significantly influences the stability and conformation of TNA:DNA duplexes. Low TNA purine content destabilizes duplexes, with T_m values often 5 °C lower than analogous DNA:DNA and RNA:DNA duplexes. By contrast, TNA:DNA duplexes having high TNA purine content display greater stability than DNA:DNA or RNA:DNA duplexes having the same sequences. High TNA purine content leads TNA:DNA duplexes to adopt conformations similar to RNA:RNA (A-form) configuration, whereas duplexes with low TNA purine content have conformations more similar to DNA:DNA (B-form) configuration. These insights provide a basis for understanding and predicting TNA:DNA duplex stability, which is anticipated to guide the practical use of TNA in biotechnology applications.

Comparison of the Binding Geometry of Free-Base and Hexacoordinated Cationic Porphyrins to A- and B-Form DNA

Although the transition from B-DNA to the A-form is essential for many biological concerns, the properties of this transition have not been resolved. The B to A equilibrium can be analyzed conveniently because of the significant changes in circular dichroism (CD) and absorption spectrum. CD and linear dichroism (LD) methods were used to examine the binding of water-soluble meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (TMPyP) and its derivatives, Co-TMPyP, with B- and A-calf thymus DNA. B- to A-transitions occurred when the physiological buffer was replaced with a water-ethanol mixture (∼80 v/v %), and the fluorescence emission spectra of TMPyP bound to DNA showed a different pattern under ethanol-water conditions and water alone. The featureless broad emission bands of TMPyP were split into two peaks near at 658 and 715 nm in the presence of DNA under an aqueous solution. In the case of an ethanol-water system, however, the emission bands are split in two peaks near at 648 and 708 nm and 656 and 715 nm with and without DNA, respectively. This may be due to a change in the solution polarity. On the basis of the CD and LD data, TMPyP interacts with B-DNA via intercalation at a low ratio under a low ionic strength, 1 mM sodium phosphate. On the other hand, the interaction with A-DNA (80 v/v % ethanol-water system) occurs in a nonintercalating manner. This difference might be because the structural conformations, such as the groove of A-DNA, are not as deep as in B-DNA and the bases are much more tilted. In the case of Co-TMPyP, porphyrin binds preferably via an outside self-stacking mode with B- and A-DNA.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.