Effects of non-nearest neighbors on the thermodynamic stability of RNA GNRA hairpin tetraloops

Nanopore molecular trajectories of a eukaryotic reverse transcriptase reveal a long-range RNA structure sensing mechanism

Eukaryotic reverse transcriptases (RTs) can have essential or deleterious roles in normal human physiology and disease. Compared to well-studied helicases, it remains unclear how RTs overcome the ubiquitous RNA structural barriers during reverse transcription. Herein, we describe the development of a Mycobacterium smegmatis porin A (MspA) nanopore technique to sequence RNA to quantify the single-molecule kinetics of an RT from Bombyx mori with single-nucleotide resolution. By establishing a quadromer map that correlates RNA sequence and MspA ion current, we were able to quantify the RT's dwell time at every single nucleotide step along its RNA template. By challenging the enzyme with various RNA structures, we found that during cDNA synthesis the RT can sense and actively destabilize RNA structures 11-12 nt downstream of its front boundary. The ability to sequence single molecules of RNA with nanopores paves the way to investigate the single-nucleotide activity of other processive RNA translocases.

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

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.

Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA

In contrast to GNRA tetraloop receptors that are common in RNA, receptors for the more thermostable UNCG loops have remained elusive for almost three decades. An analysis of all RNA structures with resolution ≤3.0 Å from the PDB allowed us to identify three previously unnoticed receptors for UNCG and GNRA tetraloops that adopt a common UNCG fold, named ‘Z-turn’ in agreement with our previously published nomenclature. These receptors recognize the solvent accessible second Z-turn nucleotide in different but specific ways. Two receptors participating in a complex network of tertiary interactions are associated with the rRNA UUCG and GAAA Z-turns capping helices H62 and H35a in rRNA large subunits. Structural comparison of fully assembled ribosomes and comparative sequence analysis of >6500 rRNA sequences helped us recognize that these motifs are almost universally conserved in rRNA, where they may contribute to organize the large subunit around the subdomain-IV four-way junction. The third UCCG receptor was identified in a rRNA/protein construct crystallized at acidic pH. These three non-redundant Z-turn receptors are relevant for our understanding of the assembly of rRNA and other long-non-coding RNAs, as well as for the design of novel folding motifs for synthetic biology.

Analysis of RNA nearest neighbor parameters reveals interdependencies and quantifies the uncertainty in RNA secondary structure prediction

RNA secondary structure prediction is often used to develop hypotheses about structure-function relationships for newly discovered RNA sequences, to identify unknown functional RNAs, and to design sequences. Secondary structure prediction methods typically use a thermodynamic model that estimates the free energy change of possible structures based on a set of nearest neighbor parameters. These parameters were derived from optical melting experiments of small model oligonucleotides. This work aims to better understand the precision of structure prediction. Here, the experimental errors in optical melting experiments were propagated to errors in the derived nearest neighbor parameter values and then to errors in RNA secondary structure prediction. To perform this analysis, the optical melting experimental values were systematically perturbed within the estimates of experimental error and alternative sets of nearest neighbor parameters were then derived from these error-bounded values. Secondary structure predictions using either the perturbed or reference parameter sets were then compared. This work demonstrated that the precision of RNA secondary structure prediction is more robust than suggested by previous work based on perturbation of the nearest neighbor parameters. This robustness is due to correlations between parameters. Additionally, this work identified weaknesses in the parameter derivation that makes accurate assessment of parameter uncertainty difficult. Considerations for experimental design are provided to mitigate these weaknesses are provided.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

DNA binding of Pd(TC3), a conformable cationic porphyrin with a long-lived triplet state

The goal of this work has been to synthesize and investigate Pd(TC3), an intercalating porphyrin that has conformable substituents capable of groove binding to B-form DNA. (TC3 denotes the doubly deprotonated form of 5,10,15,20-tetra[3-(3'-methylimidazolium-1'-yl)prop-1-yl]porphyrin.) Palladium(ii) is an apt choice for the central metal ion because it remains strictly four-coordinate and provides for a luminescent triplet excited state with a long lifetime. The DNA hosts are hairpin-forming sequences programmed to differ in base composition. Luminescence, absorbance, and circular dichroism results are consistent with the idea that congruent structural reorganization takes place at the host and ligand during uptake. Photoexcitation of DNA-bound Pd(TC3) generates a comparatively modest steady state concentration of singlet oxygen, due to a relatively slow reaction with molecular oxygen in solution. The sheer size of the substituent groups disfavors quenching, but groove-binding interactions compound the problem by inhibiting mobility. The results show how ligand design affects adduct structure as well as function.

Blind tests of RNA nearest-neighbor energy prediction

The predictive modeling and design of biologically active RNA molecules requires understanding the energetic balance among their basic components. Rapid developments in computer simulation promise increasingly accurate recovery of RNA's nearest-neighbor (NN) free-energy parameters, but these methods have not been tested in predictive trials or on nonstandard nucleotides. Here, we present, to our knowledge, the first such tests through a RECCES-Rosetta (reweighting of energy-function collection with conformational ensemble sampling in Rosetta) framework that rigorously models conformational entropy, predicts previously unmeasured NN parameters, and estimates these values' systematic uncertainties. RECCES-Rosetta recovers the 10 NN parameters for Watson-Crick stacked base pairs and 32 single-nucleotide dangling-end parameters with unprecedented accuracies: rmsd of 0.28 kcal/mol and 0.41 kcal/mol, respectively. For set-aside test sets, RECCES-Rosetta gives rmsd values of 0.32 kcal/mol on eight stacked pairs involving G-U wobble pairs and 0.99 kcal/mol on seven stacked pairs involving nonstandard isocytidine-isoguanosine pairs. To more rigorously assess RECCES-Rosetta, we carried out four blind predictions for stacked pairs involving 2,6-diaminopurine-U pairs, which achieved 0.64 kcal/mol rmsd accuracy when tested by subsequent experiments. Overall, these results establish that computational methods can now blindly predict energetics of basic RNA motifs, including chemically modified variants, with consistently better than 1 kcal/mol accuracy. Systematic tests indicate that resolving the remaining discrepancies will require energy function improvements beyond simply reweighting component terms, and we propose further blind trials to test such efforts.

Systematic reconstruction of binding and stability landscapes of the fluorogenic aptamer spinach

Fluorogenic RNAs that are based on the complex formed by 3,5-difluoro-4-hydroxybenzylidene imidazolinone (DFHBI) derivatives and the RNA aptamer named Spinach were used to engineer a new generation of in vitro and in vivo sensors for bioanalytics. With the resolved crystal structure of the RNA/small molecule complex, the engineering map becomes available, but comprehensive information regarding the thermodynamic profile of the molecule is missing. Here, we reconstructed the full thermodynamic binding and stability landscapes between DFHBI and a truncated sequence of first-generation Spinach. For this purpose, we established a systematic screening procedure for single- and double-point mutations on a microfluidic large-scale integrated chip platform for 87-nt long RNAs. The thermodynamic profile with single base resolution was used to engineer an improved fluorogenic spinach generation via a directed rather than evolutional approach.

Improved Model for Predicting the Free Energy Contribution of Dinucleotide Bulges to RNA Duplex Stability

Predicting the secondary structure of RNA is an intermediate in predicting RNA three-dimensional structure. Commonly, determining RNA secondary structure from sequence uses free energy minimization and nearest neighbor parameters. Current algorithms utilize a sequence-independent model to predict free energy contributions of dinucleotide bulges. To determine if a sequence-dependent model would be more accurate, short RNA duplexes containing dinucleotide bulges with different sequences and nearest neighbor combinations were optically melted to derive thermodynamic parameters. These data suggested energy contributions of dinucleotide bulges were sequence-dependent, and a sequence-dependent model was derived. This model assigns free energy penalties based on the identity of nucleotides in the bulge (3.06 kcal/mol for two purines, 2.93 kcal/mol for two pyrimidines, 2.71 kcal/mol for 5'-purine-pyrimidine-3', and 2.41 kcal/mol for 5'-pyrimidine-purine-3'). The predictive model also includes a 0.45 kcal/mol penalty for an A-U pair adjacent to the bulge and a -0.28 kcal/mol bonus for a G-U pair adjacent to the bulge. The new sequence-dependent model results in predicted values within, on average, 0.17 kcal/mol of experimental values, a significant improvement over the sequence-independent model. This model and new experimental values can be incorporated into algorithms that predict RNA stability and secondary structure from sequence.

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

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

Tuning riboswitch-mediated gene regulation by rational control of aptamer ligand binding properties

Riboswitch-mediated control of gene expression depends on ligand binding properties (kinetics and affinity) of its aptamer domain. A detailed analysis of interior regions of the aptamer, which affect the ligand binding properties, is important for both understanding natural riboswitch functions and for enabling rational design of tuneable artificial riboswitches. Kinetic analyses of binding reaction between flavin mononucleotide (FMN) and several natural and mutant aptamer domains of FMN-specific riboswitches were performed. The strong dependence of the dissociation rate (52.6-fold) and affinity (100-fold) on the identities of base pairs in the aptamer stem suggested that the stem region, which is conserved in length but variable in base-pair composition and context, is the tuning region of the FMN-specific aptamer. Synthetic riboswitches were constructed based on the same aptamer domain by rationally modifying the tuning regions. The observed 9.31-fold difference in the half-maximal effective concentration (EC50) corresponded to a 11.6-fold difference in the dissociation constant (K(D)) of the aptamer domains and suggested that the gene expression can be controlled by rationally adjusting the tuning regions.

Effect of sodium ions on RNA duplex stability

The standard sodium concentration for RNA optical melting experiments is 1.021 M. Algorithms that predict Tm, ΔG°37, and secondary structure from sequence generally rely on parameters derived from optical melting experiments performed in 1.021 M sodium. Physiological monovalent cation concentrations are much lower than 1.021 M. In fact, many molecular biology techniques require buffers containing monovalent cation concentrations other than 1.021 M. Predictions based on the 1.021 M Na(+) parameters may not be accurate when the monovalent cation concentration is not 1.021 M. Here, we report thermodynamic data from optical melting experiments for a set of 18 RNA duplexes, each melted over a wide range of sodium ion concentrations (71, 121, 221, and 621 mM). Using these data and previously published data for the same sequences melted in 1.021 M Na(+), we report Tm and ΔG°37 correction factors to scale the standard 1.021 M Na(+) RNA parameters to other sodium ion concentrations. The recommended Tm correction factor predicts the melting temperature within 0.7 °C, and the recommended ΔG°37 correction factor predicts the free energy within 0.14 kcal/mol. These correction factors can be incorporated into prediction algorithms that predict RNA secondary structure from sequence and provide Tm and ΔG°37 values for RNA duplexes.

Squaring theory with practice in RNA design

Ribonucleic acid (RNA) design offers unique opportunities for engineering genetic networks and nanostructures that self-assemble within living cells. Recent years have seen the creation of increasingly complex RNA devices, including proof-of-concept applications for in vivo three-dimensional scaffolding, imaging, computing, and control of biological behaviors. Expert intuition and simple design rules--the stability of double helices, the modularity of noncanonical RNA motifs, and geometric closure--have enabled these successful applications. Going beyond heuristics, emerging algorithms may enable automated design of RNAs with nucleotide-level accuracy but, as illustrated on a recent RNA square design, are not yet fully predictive. Looking ahead, technological advances in RNA synthesis and interrogation are poised to radically accelerate the discovery and stringent testing of design methods.