Changes In RNA Structure Research Articles

The role of structure in regulatory RNA elements.

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e. RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Recruitment of multi-segment genomic RNAs by Bluetongue virus requires a preformed RNA network.

How do segmented RNA viruses correctly recruit their genome has yet to be clarified. Bluetongue virus is a double-stranded RNA virus with 10 segments of different sizes, but it assembles its genome in single-stranded form through a series of specific RNA-RNA interactions prior to packaging. In this study, we determined the structure of each BTV transcript, individually and in different combinations, using 2'-hydroxyl acylation analysed by primer extension and mutational profiling (SHAPE-MaP). SHAPE-MaP identified RNA structural changes during complex formation and putative RNA-RNA interaction sites. Our data also revealed a core RNA-complex of smaller segments which serves as the foundation ('anchor') for the assembly of a complete network composed of ten ssRNA segments. The same order of core RNA complex formation was identified in cells transfected with viral RNAs. No viral protein was required for these assembly reactions. Further, substitution mutations in the interacting bases within the core assemblies, altered subsequent segment addition and affected virus replication. These data identify a wholly RNA driven reaction that may offer novel opportunities for designed attenuation or antiviral therapeutics.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression.

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Monitoring Real-Time Temperature Dynamics of a Short RNA Hairpin Using Förster Resonance Energy Transfer and Circular Dichroism.

RNA molecules play crucial roles in gene expression regulation and cellular signaling, and these functions are governed by the formation of RNA secondary and tertiary structures. These structures are highly dynamic and subject to rapid changes in response to environmental cues, temperature in particular. Thermosensitive RNA secondary structures have been harnessed by multiple organisms to survey their temperature environment and to adjust gene expression accordingly. It is thus highly desirable to observe RNA structural changes in real time over a range of temperatures. Multiple approaches have been developed to study structural dynamics, but many of these require extensive processing of the RNA, large amounts of RNA input, and/or cannot be applied under physiological conditions. Here, we describe the use of a dually fluorescently labeled RNA oligonucleotide (containing a predicted hairpin structure) to monitor subtle RNA structural dynamics in vitro by Förster resonance energy transfer (FRET) and circular dichroism (CD) spectroscopy. These approaches can be employed under physiologically relevant conditions over a range of temperatures and with RNA concentrations as low as 200nM; they enable us to observe RNA structural dynamics in real time and to correlate these dynamics with changes in biological processes such as translation.

Ribocentre-switch: a database of riboswitches.

Riboswitches are regulatory elements found in the untranslated regions (UTRs) of certain mRNA molecules. They typically comprise two distinct domains: an aptamer domain that can bind to specific small molecules, and an expression platform that controls gene expression. Riboswitches work by undergoing a conformational change upon binding to their specific ligand, thus activating or repressing the genes downstream. This mechanism allows gene expression regulation in response to metabolites or small molecules. To systematically summarise riboswitch structures and their related ligand binding functions, we present Ribocentre-switch, a comprehensive database of riboswitches, including the information as follows: sequences, structures, functions, ligand binding pocketsand biological applications. It encompasses 56 riboswitches and 26 orphan riboswitches from over 430 references, with a total of 89591 sequences. It serves as a good resource for comparing different riboswitches and facilitating the identification of potential riboswitch candidates. Therefore, it may facilitate the understanding of RNA structural conformational changes in response to ligand signaling. The database is publicly available at https://riboswitch.ribocentre.org.

Zα domain proteins mediate the immune response.

The Zα domain has a compact α/β architecture containing a three-helix bundle flanked on one side by a twisted antiparallel β sheet. This domain displays a specific affinity for double-stranded nucleic acids that adopt a left-handed helical conformation. Currently, only three Zα-domain proteins have been identified in eukaryotes, specifically ADAR1, ZBP1, and PKZ. ADAR1 is a double-stranded RNA (dsRNA) binding protein that catalyzes the conversion of adenosine residues to inosine, resulting in changes in RNA structure, function, and expression. In addition to its editing function, ADAR1 has been shown to play a role in antiviral defense, gene regulation, and cellular differentiation. Dysregulation of ADAR1 expression and activity has been associated with various disease states, including cancer, autoimmune disorders, and neurological disorders. As a sensing molecule, ZBP1 exhibits the ability to recognize nucleic acids with a left-handed conformation. ZBP1 harbors a RIP homotypic interaction motif (RHIM), composed of a highly charged surface region and a leucine-rich hydrophobic core, enabling the formation of homotypic interactions between proteins with similar structure. Upon activation, ZBP1 initiates a downstream signaling cascade leading to programmed cell death, a process mediated by RIPK3 via the RHIM motif. PKZ was identified in fish, and contains two Zα domains at the N-terminus. PKZ is essential for normal growth and development and may contribute to the regulation of immune system function in fish. Interestingly, some pathogenic microorganisms also encode Zα domain proteins, such as, Vaccinia virus and Cyprinid Herpesvirus. Zα domain proteins derived from pathogenic microorganisms have been demonstrated to be pivotal contributors in impeding the host immune response and promoting virus replication and spread. This review focuses on the mammalian Zα domain proteins: ADAR1 and ZBP1, and thoroughly elucidates their functions in the immune response.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis.

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

Exploring the mechanisms behind the E. coli pH-responsive riboswitch.

Riboswitches are 5’-untranslated regions of mRNA that change their conformation in response to ligand binding, dictating whether the downstream gene is transcribed and/or translated into protein. A “pH-responsive element” (PRE) riboswitch discovered in E. coli, however, does not require a ligand binding event in a traditional sense to rearrange its structure and instead senses the environmental H+ concentration. Specifically, at neutral pH, the PRE folds into a translationally inactive structure with an occluded ribosome binding sequence (RBS), whereas at alkaline pH, the PRE adopts a translationally active structure with an exposed RBS. The precise mechanism that controls pH-responsive riboswitching has remained unclear, with two proposed models that are yet to be distinguished; pH either regulates the transcription rate of RNA polymerase (RNAP) or acts on the RNA itself. Previous work suggested that alkaline pH enhances RNAP pause longevity during PRE synthesis to drive folding of structural motifs in the active conformer. I investigated this proposal by performing the first rigorous kinetic study of E. coli RNAP pausing at key sites for PRE folding under different pH conditions. I find that pH influences RNAP pausing, but not in the proposed manner; instead, my experiments uncovered that alkaline pH either decreases or has no effect on RNAP pause longevity. These results indicate that modulation of RNAP pausing is not the sole mechanism by which pH affects PRE folding and invite the possibility that the riboswitch RNA is actively involved in sensing the cytosolic pH (Stephen and Mishanina, J. Biol. Chem. 2022, 298:102302). I am currently conducting chemical probing and cryo-EM studies of PRE structure to explore the pH-induced RNA structural changes in rigorous detail, which will ultimately illuminate the mechanism for how pH-dependent riboswitching occurs.

Snapshots of the first-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM.

Tetrahymena ribozyme is a group I intron, whose self-splicing is the result of two sequential ester-transfer reactions. To understand how it facilitates catalysis in the first self-splicing reaction, we used cryogenic electron microscopy (cryo-EM) to resolve the structures of L-16 Tetrahymena ribozyme complexed with a 11-nucleotide 5'-splice site analog substrate. Four conformations were achieved to 4.14, 3.18, 3.09 and 2.98 Å resolutions, respectively, corresponding to different splicing intermediates during the first enzymatic reaction. Comparison of these structures reveals structural alterations, including large conformational changes in IGS/IGSext (P1-P1ext duplex) and J5/4, as well as subtle local rearrangements in the G-binding site. These structural changes are required for the enzymatic activity of the Tetrahymena ribozyme. Our study demonstrates the ability of cryo-EM to capture dynamic RNA structural changes, ushering in a new era in the analysis of RNA structure-function by cryo-EM.

The HIV-1 Integrase C-Terminal Domain Induces TAR RNA Structural Changes Promoting Tat Binding.

Recent evidence indicates that the HIV-1 Integrase (IN) binds the viral genomic RNA (gRNA), playing a critical role in the morphogenesis of the viral particle and in the stability of the gRNA once in the host cell. By combining biophysical, molecular biology, and biochemical approaches, we found that the 18-residues flexible C-terminal tail of IN acts as a sensor of the peculiar apical structure of the trans-activation response element RNA (TAR), interacting with its hexaloop. We show that the binding of the whole IN C-terminal domain modifies TAR structure, exposing critical nucleotides. These modifications favour the subsequent binding of the HIV transcriptional trans-activator Tat to TAR, finally displacing IN from TAR. Based on these results, we propose that IN assists the binding of Tat to TAR RNA. This working model provides a mechanistic sketch accounting for the emerging role of IN in the early stages of proviral transcription and could help in the design of anti-HIV-1 therapeutics against this new target of the viral infectious cycle.

Global 5'-UTR RNA structure regulates translation of a SERPINA1 mRNA.

SERPINA1 mRNAs encode the protease inhibitor α-1-antitrypsin and are regulated through post-transcriptional mechanisms. α-1-antitrypsin deficiency leads to chronic obstructive pulmonary disease (COPD) and liver cirrhosis, and specific variants in the 5′-untranslated region (5′-UTR) are associated with COPD. The NM_000295.4 transcript is well expressed and translated in lung and blood and features an extended 5′-UTR that does not contain a competing upstream open reading frame (uORF). We show that the 5′-UTR of NM_000295.4 folds into a well-defined multi-helix structural domain. We systematically destabilized mRNA structure across the NM_000295.4 5′-UTR, and measured changes in (SHAPE quantified) RNA structure and cap-dependent translation relative to a native-sequence reporter. Surprisingly, despite destabilizing local RNA structure, most mutations either had no effect on or decreased translation. Most structure-destabilizing mutations retained native, global 5′-UTR structure. However, those mutations that disrupted the helix that anchors the 5′-UTR domain yielded three groups of non-native structures. Two of these non-native structure groups refolded to create a stable helix near the translation initiation site that decreases translation. Thus, in contrast to the conventional model that RNA structure in 5′-UTRs primarily inhibits translation, complex folding of the NM_000295.4 5′-UTR creates a translation-optimized message by promoting accessibility at the translation initiation site.

Using DMS-MaPseq to uncover the roles of DEAD-box proteins in ribosome assembly.

DMS (dimethylsulfate) is a time-tested chemical probe for nucleic acid secondary structure that has recently re-emerged as a powerful tool to study RNA structure and structural changes, by coupling it to high throughput sequencing techniques. This variant, termed DMS-MaPseq, allows for mapping of all RNAs in a cell at the same time. However, if an RNA adopts different structures, for example during the assembly of an RNA-protein complex, or as part of its functional cycle, then DMS-MaPseq cannot differentiate between these structures, and an ensemble average will be produced. This is especially challenging for long-lived RNAs, such as ribosomes, whose steady-state abundance far exceeds that of any assembly intermediates, rendering those inaccessible to DMS-MaPseq on total RNAs. These challenges can be overcome by purification of assembly intermediates stalled at specific assembly steps (or steps in the functional cycle), via a combination of affinity tags and mutants stalled at defined steps, and subsequent DMS probing of these intermediates. Interpretation of the differences in DMS accessibility is facilitated by additional structural information, e.g. from cryo-EM experiments, available for many functional RNAs. While this approach is generally useful for studying RNA folding or conformational changes within RNA-protein complexes, it can be particularly valuable for studying the role(s) of DEAD-box proteins, as these tend to lead to larger conformational rearrangements, often resulting from the release of an RNA-binding protein from a bound RNA. Here we provide an adaptation of the DMS-MaPseq protocol to study RNA conformational transitions during ribosome assembly, which addresses the challenges arising from the presence of many assembly intermediates, all at concentrations far below that of mature ribosomes. While this protocol was developed for the yeast S. cerevisiae, we anticipate that it should be readily transferable to other model organisms for which affinity purification has been established.

Inferring Potential Cancer Driving Synonymous Variants.

Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.

Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

RNA Secondary Structurome Revealed Distinct Thermoregulation in Plasmodium falciparum.

The development of Plasmodium parasites, a causative agent of malaria, requests two hosts and the completion of 11 different parasite stages during development. Therefore, an efficient and fast response of parasites to various complex environmental changes, such as ambient temperature, pH, ions, and nutrients, is essential for parasite development and survival. Among many of these environmental changes, temperature is a decisive factor for parasite development and pathogenesis, including the thermoregulation of rRNA expression, gametogenesis, and parasite sequestration in cerebral malaria. However, the exact mechanism of how Plasmodium parasites rapidly respond and adapt to temperature change remains elusive. As a fundamental and pervasive regulator of gene expression, RNA structure can be a specific mechanism for fine tuning various biological processes. For example, dynamic and temperature-dependent changes in RNA secondary structures can control the expression of different gene programs, as shown by RNA thermometers. In this study, we applied the in vitro and in vivo transcriptomic-wide secondary structurome approach icSHAPE to measure parasite RNA structure changes with temperature alteration at single-nucleotide resolution for ring and trophozoite stage parasites. Among 3,000 probed structures at different temperatures, our data showed structural changes in the global transcriptome, such as S-type rRNA, HRPII gene, and the erythrocyte membrane protein family. When the temperature drops from 37°C to 26°C, most of the genes in the trophozoite stage cause significantly more changes to the RNA structure than the genes in the ring stage. A multi-omics analysis of transcriptome data from RNA-seq and RNA structure data from icSHAPE reveals that the specific RNA secondary structure plays a significant role in the regulation of transcript expression for parasites in response to temperature changes. In addition, we identified several RNA thermometers (RNATs) that responded quickly to temperature changes. The possible thermo-responsive RNAs in Plasmodium falciparum were further mapped. To this end, we identified dynamic and temperature-dependent RNA structural changes in the P. falciparum transcriptome and performed a comprehensive characterization of RNA secondary structures over the course of temperature stress in blood stage development. These findings not only contribute to a better understanding of the function of the RNA secondary structure but may also provide novel targets for efficient vaccines or drugs.

HIV-1 Packaging Visualised by In-Gel SHAPE.

HIV-1 packages two copies of its gRNA into virions via an interaction with the viral structural protein Gag. Both copies and their native RNA structure are essential for virion infectivity. The precise stepwise nature of the packaging process has not been resolved. This is largely due to a prior lack of structural techniques that follow RNA structural changes within an RNA–protein complex. Here, we apply the in-gel SHAPE (selective 2’OH acylation analysed by primer extension) technique to study the initiation of HIV-1 packaging, examining the interaction between the packaging signal RNA and the Gag polyprotein, and compare it with that of the NC domain of Gag alone. Our results imply interactions between Gag and monomeric packaging signal RNA in switching the RNA conformation into a dimerisation-competent structure, and show that the Gag–dimer complex then continues to stabilise. These data provide a novel insight into how HIV-1 regulates the translation and packaging of its genome.

Disease-associated human genetic variation through the lens of precursor and mature RNA structure.

Disease-associated variants (DAVs) are commonly considered either through a genomic lens that describes variant function at the DNA level, or at the protein function level if the variant is translated. Although the genomic and proteomic effects of variation are well-characterized, genetic variants disrupting post-transcriptional regulation is another mechanism of disease that remains understudied. Specific RNA sequence motifs mediate post-transcriptional regulation both in the nucleus and cytoplasm of eukaryotic cells, often by binding to RNA-binding proteins or other RNAs. However, many DAVs map far from these motifs, which suggests deeper layers of post-transcriptional mechanistic control. Here, we consider a transcriptomic framework to outline the importance of post-transcriptional regulation as a mechanism of disease-causing single-nucleotide variation in the human genome. We first describe the composition of the human transcriptome and the importance of abundant yet overlooked components such as introns and untranslated regions (UTRs) of messenger RNAs (mRNAs). We present an analysis of Human Gene Mutation Database variants mapping to mRNAs and examine the distribution of causative disease-associated variation across the transcriptome. Although our analysis confirms the importance of post-transcriptional regulatory motifs, a majority of DAVs do not directly map to known regulatory motifs. Therefore, we review evidence that regions outside these well-characterized motifs can regulate function by RNA structure-mediated mechanisms in all four elements of an mRNA: exons, introns, 5' and 3' UTRs. To this end, we review published examples of riboSNitches, which are single-nucleotide variants that result in a change in RNA structure that is causative of the disease phenotype. In this review, we present the current state of knowledge of how DAVs act at the transcriptome level, both through altering post-transcriptional regulatory motifs and by the effects of RNA structure.

Genome-wide RNA structure changes during human neurogenesis modulate gene regulatory networks

The distribution, dynamics, and function of RNA structures in human development are under-explored. Here, we systematically assayed RNA structural dynamics and their relationship with gene expression, translation, and decay during human neurogenesis. We observed that the human ESC transcriptome is globally more structurally accessible than differentiated cells and undergoes extensive RNA structure changes, particularly in the 3' UTR. Additionally, RNA structure changes during differentiation are associated with translation and decay. We observed that RBP and miRNA binding is associated with RNA structural changes during early neuronal differentiation, and splicing is associated during later neuronal differentiation. Furthermore, our analysis suggests that RBPs are major factors in structure remodeling and co-regulate additional RBPs and miRNAs through structure. We demonstrated an example of this by showing that PUM2-induced structure changes on LIN28A enable miR-30 binding. This study deepens our understanding of the widespread and complex role of RNA-based gene regulation during human development.

Vfold2D-MC: A Physics-Based Hybrid Model for Predicting RNA Secondary Structure Folding.

Accurate prediction of RNA structure and folding stability has a far-reaching impact on our understanding of RNA functions. Here we develop Vfold2D-MC, a new physics-based model, to predict RNA structure and folding thermodynamics from the sequence. The model employs virtual bond-based coarse-graining of RNA backbone conformation and generates RNA conformations through Monte Carlo sampling of the bond angles and torsional angles of the virtual bonds. Using a coarse-grained statistical potential derived from the known structures, we assign each conformation with a statistical weight. The weighted average over the conformational ensemble gives the entropy and free energy parameters for the hairpin, bulge, and internal loops, and multiway junctions. From the thermodynamic parameters, we predict RNA structures, melting curves, and structural changes from the sequence. Theory-experiment comparisons indicate that Vfold2D-MC not only gives improved structure predictions but also enables the interpretation of thermodynamic results for different RNA structures, including multibranched junctions. This new model sets a promising framework to treat more complicated RNA structures, such as pseudoknotted and intramolecular kissing loops, for which experimental thermodynamic parameters are often unavailable.