Predictions and analyses of RNA nearest neighbor parameters for modified nucleotides

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.

Diaminopurine in Nonenzymatic RNA Template Copying

In the RNA World before the emergence of an RNA polymerase, nonenzymatic template copying would have been essential for the transmission of genetic information. However, the products of chemical copying with the canonical nucleotides (A, U, C, and G) are heavily biased toward the incorporation of G and C, which form a more stable base pair than A and U. We therefore asked whether replacing adenine (A) with diaminopurine (D) might lead to more efficient and less biased nonenzymatic template copying by making a stronger version of the A:U pair. As expected, primer extension substrates containing D bound to U in the template more tightly than substrates containing A. However, primer extension with D exhibited elevated reaction rates on a C template, leading to concerns about fidelity. Our crystallographic studies revealed the nature of the D:C mismatch by showing that D can form a wobble-type base pair with C. We then asked whether competition with G would decrease the mismatched primer extension. We performed nonenzymatic primer extension with all four activated nucleotides on randomized RNA templates containing all four letters and used deep sequencing to analyze the products. We found that the DUCG genetic system exhibited a more even product distribution and a lower mismatch frequency than the canonical AUCG system. Furthermore, primer extension is greatly reduced following all mismatches, including the D:C mismatch. Our study suggests that D deserves further attention for its possible role in the RNA World and as a potentially useful component of artificial nonenzymatic RNA replication systems.

Isotope-Edited Variable Temperature Infrared Spectroscopy for Measuring Transition Temperatures of Single A-T Watson-Crick Base Pairs in DNA Duplexes

Experimental methods to determine transition temperatures for individual base pair melting events in DNA duplexes are lacking despite intense interest in these thermodynamic parameters. Here, we determine the dimensions of the thymine (T) C2═O stretching vibration when it is within the DNA duplex via isotopic substitutions at other atomic positions in the structure. First, we determined that this stretching state was localized enough to specific atoms in the molecule to make submolecular scale measurements of local structure and stability in high molecular weight complexes. Next, we develop a new isotope-edited variable temperature infrared method to measure melting transitions at various locations in a DNA structure. As an initial test of this "sub-molecular scale thermometer", we applied our T13C2 difference infrared signal to measure location-dependent melting temperatures (TmL) in a DNA duplex via variable temperature attenuated total reflectance Fourier transform infrared (VT-ATR-FTIR) spectroscopy. We report that the TmL of a single Watson-Crick A-T base pair near the end of an A-T rich sequence (poly T) is ∼34.9 ± 0.7°C. This is slightly lower than the TmL of a single base pair near the middle position of the poly T sequence (TmL ∼35.6±0.2°C). In addition, we also report that the TmL of a single Watson-Crick A-T base pair near the end of a 50% G-C sequence (12-mer) is ∼52.5 ± 0.3°C, which is slightly lower than the global melting Tm of the 12-mer sequence (TmL ∼54.0±0.9°C). Our results provide direct physical evidence for end fraying in DNA sequences with our novel spectroscopic methods.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Transient states during the annealing of mismatched and bulged oligonucleotides

Oligonucleotide hybridization is crucial in various biological, prebiotic and nanotechnological processes, including gene regulation, non-enzymatic primer extension and DNA nanodevice assembly. Although extensive research has focused on the thermodynamics and kinetics of nucleic acid hybridization, the behavior of complex mixtures and the outcome of competition for target binding remain less well understood. In this study, we investigate the impact of mismatches and bulges in a 12 bp DNA or RNA duplex on its association (kon) and dissociation (koff) kinetics. We find that such defects have relatively small effects on the association kinetics, while the dissociation kinetics vary in a position-dependent manner by up to 6 orders of magnitude. Building upon this observation, we explored a competition scenario involving multiple oligonucleotides, and observed a transient low specificity of probe hybridization to fully versus partially complementary targets in solution. We characterize these long-lived metastable states and their evolution toward equilibrium, and show that sufficiently long-lived mis-paired duplexes can serve as substrates for prebiotically relevant chemical copying reactions. Our results suggest that transient low accuracy states may spontaneously emerge within all complex nucleic acid systems comprising a large enough number of competing strands, with potential repercussions for gene regulation in the realm of modern biology and the prebiotic preservation of genetic information.

Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.

Predicting nearest neighbor free energies of modified RNA with LIE: results for pseudouridine and N1-methylpseudouridine within RNA duplexes

Pseudouridine (Ψ) and N1-methylpseudouridine (m1Ψ) are among the key modifications in the field of mRNA therapeutics and vaccine research. The accuracy of the design and development of therapeutic RNAs containing such modifications depends on the accuracy of the secondary structure prediction, which in turn depends on the nearest neighbor (NN) thermodynamic parameters for the standard and modified residues. Here, we propose a simple approach based on molecular dynamics simulations and linear interaction energy (LIE) approximation that is able to predict the NN free energy parameters for U-A, Ψ-A and m1Ψ-A pairs in reasonable agreement with the recent experimental reports. We report the NN thermodynamic parameters for different U, Ψ and m1Ψ base pairs, which might be helpful for a deeper understanding of the effect of these modifications in RNA. The predicted NN free energy parameters in this study are able to closely reproduce the folding free energies of duplexes containing internal Ψ for which the thermodynamic data were available. Additionally, we report the predicted folding free energies for the duplexes containing internal m1Ψ.

Modified Nucleotides and RNA Structure Prediction

Nucleotide modifications are occurrent in all types of RNA and play an important role in RNA structure formation and stability. Modified bases not only possess the ability to shift the RNA structure ensemble towards desired functional confirmations. By changes in the base pairing partner preference, they may even enlarge or reduce the conformational space, i.e., the number and types of structures the RNA molecule can adopt. However, most methods to predict RNA secondary structure do not provide the means to include the effect of modifications on the result. With the help of a heavily modified transfer RNA (tRNA) molecule, this chapter demonstrates how to include the effect of different base modifications into secondary structure prediction using the ViennaRNA Package. The constructive approach demonstrated here allows for the calculation of minimum free energy structure and suboptimal structures at different levels of modified base support. In particular we, show how to incorporate the isomerization of uridine to pseudouridine ( Ψ ) and the reduction of uridine to dihydrouridine (D).

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.

The diverse role of RNA methylation in esophageal cancer

Esophageal cancer is one of the major life-threatening diseases in the world. RNA methylation is the most common post-transcriptional modification and a wide-ranging regulatory system controlling gene expression. Numerous studies have revealed that dysregulation of RNA methylation is critical for cancer development and progression. However, the diverse role of RNA methylation and its regulators in esophageal cancer remains to be elucidated and summarized. In this review, we focus on the regulation of major RNA methylation, including m ⁶A, m ⁵C, and m ⁷G, as well as the expression patterns and clinical implications of its regulators in esophageal cancer. We systematically summarize how these RNA modifications affect the "life cycle" of target RNAs, including mRNA, microRNA, long non-coding RNA, and tRNA. The downstream signaling pathways associated with RNA methylation during the development and treatment of esophageal cancer are also discussed in detail. Further studies on how these modifications function together in the microenvironment of esophageal cancer will draw a clearer picture of the clinical application of novel and specific therapeutic strategies.

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.

Molecular Simulations Matching Denaturation Experiments for N6-Methyladenosine

Post-transcriptional modifications are crucial for RNA function and can affect its structure and dynamics. Force-field-based classical molecular dynamics simulations are a fundamental tool to characterize biomolecular dynamics, and their application to RNA is flourishing. Here, we show that the set of force-field parameters for N⁶-methyladenosine (m⁶A) developed for the commonly used AMBER force field does not reproduce duplex denaturation experiments and, specifically, cannot be used to describe both paired and unpaired states. Then, we use reweighting techniques to derive new parameters matching available experimental data. The resulting force field can be used to properly describe paired and unpaired m⁶A in both syn and anti conformation, which thus opens the way to the use of molecular simulations to investigate the effects of N6 methylations on RNA structural dynamics.

Role of main RNA modifications in cancer: N6-methyladenosine, 5-methylcytosine, and pseudouridine

Cancer is one of the major diseases threatening human life and health worldwide. Epigenetic modification refers to heritable changes in the genetic material without any changes in the nucleic acid sequence and results in heritable phenotypic changes. Epigenetic modifications regulate many biological processes, such as growth, aging, and various diseases, including cancer. With the advancement of next-generation sequencing technology, the role of RNA modifications in cancer progression has become increasingly prominent and is a hot spot in scientific research. This review studied several common RNA modifications, such as N⁶-methyladenosine, 5-methylcytosine, and pseudouridine. The deposition and roles of these modifications in coding and noncoding RNAs are summarized in detail. Based on the RNA modification background, this review summarized the expression, function, and underlying molecular mechanism of these modifications and their regulators in cancer and further discussed the role of some existing small-molecule inhibitors. More in-depth studies on RNA modification and cancer are needed to broaden the understanding of epigenetics and cancer diagnosis, treatment, and prognosis.

Modelling aptamers with nucleic acid mimics (NAM): From sequence to three-dimensional docking

Aptamers are single-stranded oligonucleotides, formerly evolved by Systematic Evolution of Ligands by EXponential enrichment (SELEX), that fold into functional three-dimensional structures. Such conformation is crucial for aptamers' ability to bind to a target with high affinity and specificity. Unnatural nucleotides have been used to develop nucleic acid mimic (NAM) aptamers with increased performance, such as biological stability. Prior knowledge of aptamer-target interactions is critical for applying post-SELEX modifications with unnatural nucleotides since it can affect aptamers' structure and performance. Here, we describe an easy-to-apply in silico workflow using free available software / web servers to predict the tertiary conformation of NAM, DNA and RNA aptamers, as well as the docking with the target molecule. Representative 2'-O-methyl (2'OMe), locked nucleic acid (LNA), DNA and RNA aptamers, with experimental data deposited in Protein Data Bank, were selected to validate the workflow. All aptamers' tertiary structure and docking models were successfully predicted with good structural similarity to the experimental data. Thus, this workflow will boost the development of aptamers, particularly NAM aptamers, by assisting in the rational modification of specific nucleotides and avoiding trial-and-error approaches.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Secondary structure prediction for RNA sequences including N6-methyladenosine

There is increasing interest in the roles of covalently modified nucleotides in RNA. There has been, however, an inability to account for modifications in secondary structure prediction because of a lack of software and thermodynamic parameters. We report the solution for these issues for N⁶-methyladenosine (m⁶A), allowing secondary structure prediction for an alphabet of A, C, G, U, and m⁶A. The RNAstructure software now works with user-defined nucleotide alphabets of any size. We also report a set of nearest neighbor parameters for helices and loops containing m⁶A, using experiments. Interestingly, N⁶-methylation decreases folding stability for adenosines in the middle of a helix, has little effect on folding stability for adenosines at the ends of helices, and increases folding stability for unpaired adenosines stacked on a helix. We demonstrate predictions for an N⁶-methylation-activated protein recognition site from MALAT1 and human transcriptome-wide effects of N⁶-methylation on the probability of adenosine being buried in a helix.

RNA folding free energy nearest neighbor parameters were determined for sequences with the nucleotide m⁶A. The RNAstructure software package can accommodate modified nucleotides, enabling secondary structure prediction of sequences with m⁶A.

Recent advances in functional annotation and prediction of the epitranscriptome

RNA modifications, in particular N⁶-methyladenosine (m⁶A), participate in every stages of RNA metabolism and play diverse roles in essential biological processes and disease pathogenesis. Thanks to the advances in sequencing technology, tens of thousands of RNA modification sites can be identified in a typical high-throughput experiment; however, it remains a major challenge to decipher the functional relevance of these sites, such as, affecting alternative splicing, regulation circuit in essential biological processes or association to diseases. As the focus of RNA epigenetics gradually shifts from site discovery to functional studies, we review here recent progress in functional annotation and prediction of RNA modification sites from a bioinformatics perspective. The review covers naïve annotation with associated biological events, e.g., single nucleotide polymorphism (SNP), RNA binding protein (RBP) and alternative splicing, prediction of key sites and their regulatory functions, inference of disease association, and mining the diagnosis and prognosis value of RNA modification regulators. We further discussed the limitations of existing approaches and some future perspectives.

Deciphering nucleotide modification-induced structure and stability changes

Nucleotide modification in RNA controls a bevy of biological processes, including RNA degradation, gene expression, and gene editing. In turn, misregulation of modified nucleotides is associated with a host of chronic diseases and disorders. However, the molecular mechanisms driving these processes remain poorly understood. To partially address this knowledge gap, we used alchemical and temperature replica exchange molecular dynamics (TREMD) simulations on an RNA duplex and an analogous hairpin to probe the structural effects of modified and/or mutant nucleotides. The simulations successfully predict the modification/mutation-induced relative free energy change for complementary duplex formation, and structural analyses highlight mechanisms driving stability changes. Furthermore, TREMD simulations for a hairpin-forming RNA with and without modification provide reliable estimations of the energy landscape. Illuminating the impact of methylated and/or mutated nucleotides on the structure-function relationship and the folding energy landscape, the simulations provide insights into modification-induced alterations to the folding mechanics of the hairpin. The results here may be biologically significant as hairpins are widespread structure motifs that play critical roles in gene expression and regulation. Specifically, the tetraloop of the probed hairpin is phylogenetically abundant, and the stem mirrors a miRNA seed region whose modification has been implicated in epilepsy pathogenesis.