Universal cold RNA phase transitions.

RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.

Phase Separation Modulates the Thermodynamics and Kinetics of RNA Hybridization.

Biomolecular condensates can influence cellular function in a number of ways, including by changing the structural dynamics and conformational equilibria of the molecules partitioned within them. Here we use methyl transverse relaxation optimized spectroscopy (methyl-TROSY) NMR in conjunction with 2'-O-methyl labeling of RNA to characterize the thermodynamics and kinetics of RNA-RNA base pairing in condensates formed by the C-terminal intrinsically disordered region of CAPRIN1, an RNA-binding protein involved in RNA transport, translation, and stability. CAPRIN1 condensates destabilize RNA-RNA base pairing, resulting from a ∼270-fold decrease and a concomitant ∼15-fold increase in the on- and off-rates for duplex formation, respectively. The ∼30-fold slower diffusion of RNA single strands within the condensed phase partially accounts for the reduced on-rate, but the further ∼9-fold reduction likely reflects shedding of CAPRIN1 chains that are interacting with the RNA prior to hybridization. Our study emphasizes the important role of protein solvation in modulating nucleic acid recognition processes inside condensates.

Predicting the effect of binding molecules on the shape and mechanical properties of structured DNA assemblies.

Chemo-mechanical deformation of structured DNA assemblies driven by DNA-binding ligands has offered promising avenues for biological and therapeutic applications. However, it remains elusive how to effectively model and predict their effects on the deformation and mechanical properties of DNA structures. Here, we present a computational framework for simulating chemo-mechanical change of structured DNA assemblies. We particularly quantify the effects of ethidium bromide (EtBr) intercalation on the geometry and mechanical properties of DNA base-pairs through molecular dynamics simulations and integrated them into finite-element-based structural analysis to predict the shape and properties of DNA objects. The proposed model captures various structural changes induced by EtBr-binding such as shape variation, flexibility modulation, and supercoiling instability. It enables a rational design of structured DNA assemblies with tunable shapes and mechanical properties by binding molecules.

Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field.

Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.

Structural basis of water-mediated cis Watson-Crick/Hoogsteen base-pair formation in non-CpG methylation.

Non-CpG methylation is associated with several cellular processes, especially neuronal development and cancer, while its effect on DNA structure remains unclear. We have determined the crystal structures of DNA duplexes containing -CGCCG- regions as CCG repeat motifs that comprise a non-CpG site with or without cytosine methylation. Crystal structure analyses have revealed that the mC:G base-pair can simultaneously form two alternative conformations arising from non-CpG methylation, including a unique water-mediated cis Watson-Crick/Hoogsteen, (w)cWH, and Watson-Crick (WC) geometries, with partial occupancies of 0.1 and 0.9, respectively. NMR studies showed that an alternative conformation of methylated mC:G base-pair at non-CpG step exhibits characteristics of cWH with a syn-guanosine conformation in solution. DNA duplexes complexed with the DNA binding drug echinomycin result in increased occupancy of the (w)cWH geometry in the methylated base-pair (from 0.1 to 0.3). Our structural results demonstrated that cytosine methylation at a non-CpG step leads to an anti→syntransition of its complementary guanosine residue toward the (w)cWH geometry as a partial population of WC, in both drug-bound and naked mC:G base pairs. This particular geometry is specific to non-CpG methylated dinucleotide sites in B-form DNA. Overall, the current study provides new insights into DNA conformation during epigenetic regulation.

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.

7-Deazapurine and Pyrimidine Nucleoside and Oligonucleotide Cycloadducts Formed by Inverse Diels-Alder Reactions with 3,6-Di(pyrid-2-yl)-1,2,4,5-tetrazine: Ethynylated and Vinylated Nucleobases for Functionalization and Impact of Pyridazine Adducts on DNA Base Pair Stability and Mismatch Discrimination.

The manuscript reports on 7-deazapurine and pyrimidine nucleoside and oligonucleotide cycloadducts formed by the inverse electron demand Diels-Alder (iEDDA) reaction with 3,6-di(pyrid-2-yl)-1,2,4,5-tetrazine. Cycloadducts were constructed from ethynylated and vinylated nucleobases. Oligonucleotides were synthesized containing iEDDA modifications, and the impact on duplex stability was investigated. iEDDA reactions were performed on nucleoside triple bond side chains. Oxidation was not required in these cases as dihydropyridazine intermediates are not formed. In contrast, oxidation is necessary for reactions performed on alkenyl compounds. This was verified on 5-vinyl-2'-deoxyuridine. A diastereomeric mixture of 1,2-dihydropyridazine cycloadduct intermediates was isolated, characterized, and later oxidized. 12-mer oligonucleotides containing 1,2-pyridazine inverse Diels-Alder cycloadducts and their precursors were hybridized to short DNA duplexes. For that, a series of phosphoramidites was prepared. DNA duplexes with 7-functionalized 7-deazaadenines and 5-functionalized pyrimidines display high duplex stability when spacer units are present between nucleobases and pyridazine cycloadducts. A direct connectivity of the pyridazine moiety to nucleobases as reported for metabolic labeling of vinyl nucleosides reduced duplex stability strongly. Oligonucleotides bearing linkers with and without pyridazine cycloadducts attached to the 7-deazaadenine nucleobase significantly reduced mismatch formation with dC and dG.

Unexpected large charge transfer rate mediated by adenine in twisted DNA structure.

DNA exhibits remarkable charge transfer ability, which is crucial for its biological functions and potential electronic applications. The charge transfer process in DNA is widely recognized as primarily mediated by guanine, while the contribution of other nucleobases is negligible. Using the tight-binding models in conjunction with first-principles calculations, we investigated the charge transfer behavior of homogeneous GC and AT pairs. We found that the charge transfer rate of adenine significantly changes. With overstretching, the charge transfer rate of adenine can even surpass that of guanine, by as much as five orders of magnitude at a twist angle of around 26°. Further analysis reveals that it is attributed to the turnover of the relative coupling strength between homogeneous GC and AT base pairs, which is caused by the symmetry exchange between the two highest occupied molecular orbitals of base pairs occurring at different twist angles. Given the high degree of flexibility of DNA in vivo and in vitro conditions, these findings prompt us to reconsider the mechanism of biological functions concerning the charge transfer in DNA molecules and further open the potential of DNA as a biomaterial for electronic applications.

GeF-seq: A Simple Procedure for Base-Pair Resolution ChIP-seq.

Nucleotide sequences recognized and bound by DNA-binding proteins (DBPs) are critical to controlling and maintaining gene expression, replication, chromosome segregation, cell division, and nucleoid structure in bacterial cells. Therefore, determination of the binding sequences of DBPs is important not only to study DBP recognition mechanisms but also to understand the fundamentals of cell homeostasis. While ChIP-seq analysis appears to be an effective way to determine DBP binding sites on the genome, the resolution is sometimes not sufficient to identify the sites precisely. Here we introduce a simple and effective method named Genome Footprinting with high-throughput sequencing (GeF-seq) to determine binding sites of DBPs with single base-pair resolution. GeF-seq detects binding sites of DBPs as sharp peaks and thus makes it possible to identify the recognition sequence in each "binding peak" more easily and accurately compared to the common ChIP-seq.

Knotify_V2.0: Deciphering RNA Secondary Structures with H-Type Pseudoknots and Hairpin Loops.

Accurately predicting the pairing order of bases in RNA molecules is essential for anticipating RNA secondary structures. Consequently, this task holds significant importance in unveiling previously unknown biological processes. The urgent need to comprehend RNA structures has been accentuated by the unprecedented impact of the widespread COVID-19 pandemic. This paper presents a framework, Knotify_V2.0, which makes use of syntactic pattern recognition techniques in order to predict RNA structures, with a specific emphasis on tackling the demanding task of predicting H-type pseudoknots that encompass bulges and hairpins. By leveraging the expressive capabilities of a Context-Free Grammar (CFG), the suggested framework integrates the inherent benefits of CFG and makes use of minimum free energy and maximum base pairing criteria. This integration enables the effective management of this inherently ambiguous task. The main contribution of Knotify_V2.0 compared to earlier versions lies in its capacity to identify additional motifs like bulges and hairpins within the internal loops of the pseudoknot. Notably, the proposed methodology, Knotify_V2.0, demonstrates superior accuracy in predicting core stems compared to state-of-the-art frameworks. Knotify_V2.0 exhibited exceptional performance by accurately identifying both core base pairing that form the ground truth pseudoknot in 70% of the examined sequences. Furthermore, Knotify_V2.0 narrowed the performance gap with Knotty, which had demonstrated better performance than Knotify and even surpassed it in Recall and F1-score metrics. Knotify_V2.0 achieved a higher count of true positives (tp) and a significantly lower count of false negatives (fn) compared to Knotify, highlighting improvements in Prediction and Recall metrics, respectively. Consequently, Knotify_V2.0 achieved a higher F1-score than any other platform. The source code and comprehensive implementation details of Knotify_V2.0 are publicly available on GitHub.

Structure of an RNA G-quadruplex from the West Nile virus genome.

Potential G-quadruplex sites have been identified in the genomes of DNA and RNA viruses and proposed as regulatory elements. The genus Orthoflavivirus contains arthropod-transmitted, positive-sense, single-stranded RNA viruses that cause significant human disease globally. Computational studies have identified multiple potential G-quadruplex sites that are conserved across members of this genus. Subsequent biophysical studies established that some G-quadruplexes predicted in Zika and tickborne encephalitis virus genomes can form and known quadruplex binders reduced viral yields from cells infected with these viruses. The susceptibility of RNA to degradation and the variability of loop regions have made structure determination challenging. Despite these difficulties, we report a high-resolution structure of the NS5-B quadruplex from the West Nile virus genome. Analysis reveals two stacked tetrads that are further stabilized by a stacked triad and transient noncanonical base pairing. This structure expands the landscape of solved RNA quadruplex structures and demonstrates the diversity and complexity of biological quadruplexes. We anticipate that the availability of this structure will assist in solving further viral RNA quadruplexes and provides a model for a conserved antiviral target in Orthoflavivirus genomes.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials.

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

DNA-based ForceChrono probes for deciphering single-molecule force dynamics in living cells.

Understanding cellular force transmission dynamics is crucial in mechanobiology. We developed the DNA-based ForceChrono probe to measure force magnitude, duration, and loading rates at the single-molecule level within living cells. The ForceChrono probe circumvents the limitations of in vitro single-molecule force spectroscopy by enabling direct measurements within the dynamic cellular environment. Our findings reveal integrin force loading rates of 0.5-2 pN/s and durations ranging from tens of seconds in nascent adhesions to approximately 100 s in mature focal adhesions. The probe's robust and reversible design allows for continuous monitoring of these dynamic changes as cells undergo morphological transformations. Additionally, by analyzing how mutations, deletions, or pharmacological interventions affect these parameters, we can deduce the functional roles of specific proteins or domains in cellular mechanotransduction. The ForceChrono probe provides detailed insights into the dynamics of mechanical forces, advancing our understanding of cellular mechanics and the molecular mechanisms of mechanotransduction.

The binding of diruthenium (II,III) and dirhodium (II,II) paddlewheel complexes at DNA/RNA nucleobases: Computational evidences of an appreciable selectivity toward the AU base pairs.

Multiple medicinal strategies involve modifications of the structure of DNA or RNA, which disrupt their correct functioning. Metal complexes with medicinal effects, also known as metallodrugs, are among the agents intended specifically for the attack onto nucleosides. The diruthenium (II,III) and dirhodium (II,II) paddlewheel complexes constitute promising dual acting drugs due to their ability to release the therapeutically active bridging ligands upon their substitution by endogenous ligands. In this paper, we study the structure and the stability of the complexes formed by the diruthenium (II,III) and dirhodium (II,II) paddlewheel complexes coordinated in axial positions with the DNA/RNA nucleobases or base pairs, assuming the attainable metalation at all the accessible pyridyl nitrogens. Dirhodium complexes coordinate at the pyridyl nitrogens more strongly than the diruthenium complexes. On the other hand, we found that the diruthenium scaffold binds more selectively to nucleobase targets. Furthermore, we reveal a tighter coordination of diruthenium complex at the adenine-uracil base pair, compared to adenine-thymine, hence constituting a scarce instance of RNA-selectivity. We envision that the here reported computational outcomes may pace future experiments addressing the binding of diruthenium and dirhodium paddlewheel complexes at either single nucleobases or DNA/RNA fragments.

RNA Complexes with Nicks and Gaps: Thermodynamic and Kinetic Effects of Coaxial Stacking and Dangling Ends.

Multiple RNA strands can interact in solution and assume a large variety of configurations dictated by their potential for base pairing. Although duplex formation from two complementary oligonucleotides has been studied in detail, we still lack a systematic characterization of the behavior of higher order complexes. Here, we focus on the thermodynamic and kinetic effects of an upstream oligonucleotide on the binding of a downstream oligonucleotide to a common template, as we vary the sequence and structure of the contact interface. We show that coaxial stacking in RNA is well correlated with but much more stabilizing than helix propagation over an analogous intact double helix step (median ΔΔG°37 °C ≈ 1.7 kcal/mol). Consequently, approximating coaxial stacking in RNA with the helix propagation term leads to large discrepancies between predictions and our experimentally determined melting temperatures, with an offset of ≈10 °C. Our kinetic study reveals that the hybridization of the downstream probe oligonucleotide is impaired (lower kon) by the presence of the upstream oligonucleotide, with the thermodynamic stabilization coming entirely from an extended lifetime (lower koff) of the bound downstream oligonucleotide, which can increase from seconds to months. Surprisingly, we show that the effect of nicks is dependent on the length of the stacking oligonucleotides, and we discuss the binding of ultrashort (1-4 nt) oligonucleotides that are relevant in the context of the origin of life. The thermodynamic and kinetic data obtained in this work allow for the prediction of the formation and stability of higher-order multistranded complexes.

Simulating the Lyotropic Phase Behavior of a Partially Self-Complementary DNA Tetramer.

DNA oligomers in solution have been found to develop liquid crystal phases via a hierarchical process that involves Watson-Crick base pairing, supramolecular assembly into columns of duplexes, and long-range ordering. The multiscale nature of this phenomenon makes it difficult to quantitatively describe and assess the importance of the various contributions, particularly for very short strands. We performed molecular dynamics simulations based on the coarse-grained oxDNA model, aiming to depict all of the assembly processes involved and the phase behavior of solutions of the DNA GCCG tetramers. We find good quantitative matching to experimental data at both levels of molecular association (thermal melting) and collective ordering (phase diagram). We characterize the isotropic state and the low-density nematic and high-density columnar liquid crystal phases in terms of molecular order, size of aggregates, and structure, together with their effects on diffusivity processes. We observe a cooperative aggregation mechanism in which the formation of dimers is less thermodynamically favored than the formation of longer aggregates.

Bridge RNAs direct programmable recombination of target and donor DNA.

Genomic rearrangements, encompassing mutational changes in the genome such as insertions, deletions or inversions, are essential for genetic diversity. These rearrangements are typically orchestrated by enzymes that are involved in fundamental DNA repair processes, such as homologous recombination, or in the transposition of foreign genetic material by viruses and mobile genetic elements1,2. Here we report that IS110 insertion sequences, a family of minimal and autonomous mobile genetic elements, express a structured non-coding RNA that binds specifically to their encoded recombinase. This bridge RNA contains two internal loops encoding nucleotide stretches that base-pair with the target DNA and the donor DNA, which is the IS110 element itself. We demonstrate that the target-binding and donor-binding loops can be independently reprogrammed to direct sequence-specific recombination between two DNA molecules. This modularity enables the insertion of DNA into genomic target sites, as well as programmable DNA excision and inversion. The IS110 bridge recombination system expands the diversity of nucleic-acid-guided systems beyond CRISPR and RNA interference, offering a unified mechanism for the three fundamental DNA rearrangements-insertion, excision and inversion-that are required for genome design.

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH.

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

Efficient Silver(I)-Containing I-Motif DNA Hybrid Catalyst for Enantioselective Diels-Alder Reactions.

The inherent chiral structures of DNA serve as attractive scaffolds to construct DNA hybrid catalysts for valuable enantioselective transformations. Duplex and G-quadruplex DNA-based enantioselective catalysis has made great progress, yet novel design strategies of DNA hybrid catalysts are highly demanding and atomistic analysis of active centers is still challenging. DNA i-motif structures could be finely tuned by different cytosine-cytosine base pairs, providing a new platform to design DNA catalysts. Herein, we found that a human telomeric i-motif DNA containing cytosine-silver(I)-cytosine (C-Ag+-C) base pairs interacting with Cu(II) ions (i-motif DNA(Ag+)/Cu2+) could catalyze Diels-Alder reactions with full conversions and up to 95 % enantiomeric excess. As characterized by various physicochemical techniques, the presence of Ag+ is proved to replace the protons in hemiprotonated cytosine-cytosine (C : C+) base pairs and stabilize the DNA i-motif to allow the acceptance of Cu(II) ions. The i-motif DNA(Ag+)/Cu2+ catalyst shows about 8-fold rate acceleration compared with DNA and Cu2+. Based on DNA mutation experiments, thermodynamic studies and density function theory calculations, the catalytic center of Cu(II) ion is proposed to be located in a specific loop region via binding to one nitrogen-7 atom of an unpaired adenine and two phosphate-oxygen atoms from nearby deoxythymidine monophosphate and deoxyadenosine monophosphate, respectively.

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction.

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences. Evidence for specific base pairs turns out to be more informative than a position-wise profile for the conservation of the pairing status. A comparison with chemical probing data, furthermore, strongly suggests that phylogenetic base pairing data are more informative than position-specific data on (un)pairedness as obtained from chemical probing experiments. In this context we demonstrate, in addition, that the conversion of signal from probing data into pseudo-energies is possible using thermodynamic structure predictions as a reference instead of known RNA structures.