Melting studies of RNA unfolding and RNA-ligand interactions

Tradeoffs in the design of RNA thermometers

The synthesis of RNA thermometers is aimed at achieving temperature responses with desired thresholds and sensitivities. Although previous works have generated thermometers with a variety of thresholds and sensitivities as well as guidelines for design, possible constraints in the achievable thresholds and sensitivities remain unclear. We addressed this issue using a two-state model and its variants, as well as melt profiles generated from thermodynamic computations. In the two-state model, we found that the threshold was inversely proportional to the sensitivity, in the case of a fixed energy difference between the two states. Notably, this constraint could persist in variations of the two-state model with sequentially unfolding states and branched parallel pathways. Furthermore, the melt profiles generated from a library of thermometers exhibited a similar constraint. These results should inform the design of RNA thermometers as well as other responses that are mediated in a similar fashion.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p

The s2m, a highly conserved 41-nt hairpin structure in the SARS-CoV-2 genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and the subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a >99% correlation of a single nucleotide mutation at the 15th position (G15U) in Delta SARS-CoV-2. Based on 1H NMR spectroscopy assignments comparing the imino proton resonance region of s2m and the s2m G15U at 19°C, we show that the U15-A29 base pair closes, resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we show that the s2m G15U mutation drastically impacts the binding of host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on immune responses.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p

The stem loop 2 motif (s2m), a highly conserved 41-nucleotide hairpin structure in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a greater than 99% correlation of a single nucleotide mutation at the 15 th position (G15U) in Delta SARS-CoV-2. Based on 1 H NMR assignments comparing the imino proton resonance region of s2m and the G15U at 19°C, we find that the U15-A29 base pair closes resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we find that the s2m G15U mutation drastically reduces the binding affinity of the host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on the immune response.

A low-cost 3D-printable differential scanning fluorometer for protein and RNA melting experiments

Differential scanning fluorimetry (DSF) is a widely used biophysical technique with applications to drug discovery and protein biochemistry. DSF experiments are commonly performed in commercial real-time polymerase chain reaction (qPCR) thermal cyclers or nanoDSF instruments. Here, we report the construction, validation, and example applications of an open-source DSF system for 176 €, which, in addition to protein-DSF experiments, also proved to be a versatile biophysical instrument for less conventional RNA-DSF experiments. Using 3D-printed parts made of polyoxymethylene, we were able to fabricate a thermostable machine chassis for protein-melting experiments. The combination of blue high-power LEDs as the light source and stage light foil as filter components was proven to be a reliable and affordable alternative to conventional optics equipment for the detection of SYPRO Orange or Sybr Gold fluorescence. The ESP32 microcontroller is the core piece of this openDSF instrument, while the in-built I²S interface was found to be a powerful analog-to-digital converter for fast acquisition of fluorescence and temperature data. Airflow heating and inline temperature control by thermistors enabled high-accuracy temperature management in PCR tubes (±0.1 °C) allowing us to perform high-resolution thermal shift assays (TSA) from exemplary biological applications.

Effect of Mutations on the Binding of Kanamycin-B to RNA Hairpins Derived from the Mycobacterium tuberculosis Ribosomal A-Site

Kanamycin is an aminoglycoside antibiotic used in the treatment of drug-resistant tuberculosis. Mutations at the rRNA A-site have been associated with kanamycin resistance in Mycobacterium tuberculosis clinical isolates. Understanding the effect of these mutations on the conformation of the M. tuberculosis A-site is critical for understanding the mechanisms of antibiotic resistance in M. tuberculosis. In this work, we have studied RNA hairpins derived from the M. tuberculosis A-site, the wild type and three mutants at the following positions (M. tuberculosis/Escherichia coli numbering): A1400/1408 → G, C1401/1409 → U, and the double mutant G1483/1491 C1401/1409 → UA. Specifically, we used circular dichroism, ultraviolet spectroscopy, and fluorescence spectroscopy to characterize the conformation, stability, and binding affinity of kanamycin-B and other aminoglycoside antibiotics for these RNA hairpins. Our results show that the mutations affect the conformation of the decoding site, with the mutations at position 1401/1409 resulting in significant destabilizations. Interestingly, the mutants bind paromomycin with weaker affinity than the wild type, but they bind kanamycin-B with similar affinity than the wild type. The results suggest that the presence of mutations does not prevent kanamycin-B from binding. Instead, kanamycin may promote different interactions with a third partner in the mutants compared to the wild type. Furthermore, our results with longer and shorter hairpins suggest that the region of the A-site that varies among organisms may have modulating effects on the binding and interactions of the A-site.

Identification of Spermidine Binding Site in T-box Riboswitch Antiterminator RNA

The T-box transcription antitermination riboswitch controls bacterial gene expression by structurally responding to uncharged, cognate tRNA. Previous studies indicated that cofactors, such as the polyamine spermidine, might serve a specific functional role in enhancing riboswitch efficacy. As riboswitch function depends on key RNA structural changes involving the antiterminator element, the interaction of spermidine with the T-box riboswitch antiterminator element was investigated. Spermidine binds antiterminator model RNA with high affinity (micromolar Kd ) based on isothermal titration calorimetry and fluorescence-monitored binding assays. NMR titration studies, molecular modeling, and inline and enzymatic probing studies indicate that spermidine binds at the 3' portion of the highly conserved seven-nucleotide bulge in the antiterminator. Together, these results support the conclusion that spermidine binds the T-box antiterminator RNA preferentially in a location important for antiterminator function. The implications of these findings are significant both for better understanding of the T-box riboswitch mechanism and for antiterminator-targeted drug discovery efforts.

Ligand-induced changes in T box antiterminator RNA stability

The T box antiterminator RNA element is an important component of the T box riboswitch that controls the transcription of vital genes in many Gram-positive bacteria. A series of 1,4-disubstituted 1,2,3-triazoles was screened in a fluorescence-monitored thermal denaturation assay to identify ligands that altered the stability of antiterminator model RNA. Several ligands were identified that significantly increased or decreased the melting temperature (T(m) ) of the RNA. The results indicate that this series of triazole ligands can alter the stability of antiterminator model RNA in a structure-dependent manner.

Characterization of a 1,4-disubstituted 1,2,3-triazole binding to T box antiterminator RNA

The T box riboswitch regulates the transcription of many bacterial genes by structurally responding to cognate non-aminoacylated (uncharged) tRNA. The riboswitch contains multiple conserved RNA elements including a key structural element, the antiterminator, which binds the tRNA acceptor end nucleotides. Previous studies identified a lead 1,4-disubstituted 1,2,3-triazole, GHB-7, that disrupted formation of a tRNA-antiterminator RNA model complex. The affinity and molecular interactions of GHB-7 binding to antiterminator model RNA were characterized as part of a comprehensive T box antiterminator RNA-targeted drug discovery project. In-line probing, UV-monitored thermal denaturation and docking studies all consistently indicated that GHB-7 likely binds to the bulge region of the antiterminator, reduces the flexibility of the bulge nucleotides and, overall, stabilizes the RNA secondary structure. These results begin to elucidate possible mechanisms for ligand-induced inhibition of tRNA binding to T box antiterminator RNA and contribute to the knowledge of how small molecules bind relatively simple RNA structural elements such as bulges.

Effect of salt concentration on the conformation of TAR RNA and its association with aminoglycoside antibiotics

RNA is an important biological target because it plays essential roles in many pathogenic and normal cellular processes. The design of inhibitors that target RNA involves optimization of noncovalent interactions, including van der Waals, hydrogen bond, and electrostatic interactions. Although sometimes regarded as nonspecific, electrostatic interactions are important in this optimization because the specific position of the phosphates may allow for specific charge-charge interactions with bound ligands. In this work, we have investigated the contribution of electrostatic interactions to the binding affinity of aminoglycoside antibiotics for TAR RNA. Because the charges in aminoglycoside antibiotics are provided by protonated amino groups, it is difficult to separate the contribution of hydrogen bonds and electrostatics to their binding specificity. Hence, we have investigated the dependence of the binding affinity on salt concentration, which should affect only the electrostatic contributions. Our results show that four aminoglycoside antibiotics (paromomycin, kanamycin-B, gentamycin, and tobramycin) bind TAR RNA with different affinities. Furthermore, the dependence of the binding affinity on salt concentration is different for kanamycin-B and paromomycin, with kanamycin-B showing a stronger dependence. Because all these antibiotics contain five positive charges, the results suggest that each antibiotic orients its charges in different ways when bound to TAR RNA. Our overall results support the idea that charge-charge interactions can contribute significantly to the specific binding of antibiotics to TAR RNA. Hence, the exact position of the charges should be considered in the design of any inhibitor of the interactions of TAR RNA.

Thermal methods for the analysis of RNA folding pathways

Once a model of the secondary structure of an RNA has been deduced, thermal melting analysis can be used to determine whether the model accounts for all intramolecular interactions of the RNA, or whether noncanonical and tertiary interactions make the structure more stable than predicted, or link parts of the structure in unexpected ways. It is also useful to determine the pH, salt, and temperature ranges under which the RNA adopts a stably folded structure, or to analyze unfolding pathways. This unit discusses sample preparation, instrumentation, and theoretical background. It also provide a sample analysis of tRNA unfolding.

Site-specific variations in RNA folding thermodynamics visualized by 2-aminopurine fluorescence

The fluorescent base analogue 2-aminopurine (2-AP) is commonly used to study specific conformational and protein binding events involving nucleic acids. Here, combinations of steady-state and time-resolved fluorescence spectroscopy of 2-AP were employed to monitor conformational transitions within a model hairpin RNA from diverse structural perspectives. RNA substrates adopting stable, unambiguous secondary structures were labeled with 2-AP at an unpaired base, within the loop, or inside the base-paired stem. Steady-state fluorescence was monitored as the RNA hairpins made the transitions between folded and unfolded conformations using thermal denaturation, urea titration, and cation-mediated folding. Unstructured control RNA substrates permitted the effects of higher-order RNA structures on 2-AP fluorescence to be distinguished from stimulus-dependent changes in intrinsic 2-AP photophysics and/or interactions with adjacent residues. Thermodynamic parameters describing local conformational changes were thus resolved from multiple perspectives within the model RNA hairpin. These data provided energetic bases for construction of folding mechanisms, which varied among different folding-unfolding stimuli. Time-resolved fluorescence studies further revealed that 2-AP exhibits characteristic signatures of component fluorescence lifetimes and respective fractional contributions in different RNA structural contexts. Together, these studies demonstrate localized conformational events contributing to RNA folding and unfolding that could not be observed by approaches monitoring only global structural transitions.

Interactions of the N-terminal domain of ribosomal protein L11 with thiostrepton and rRNA

Ribosomal protein L11 has two domains: the C-terminal domain (L11-C76) binds rRNA, whereas the N-terminal domain (L11-NTD) may variously interact with elongation factor G, the antibiotic thiostrepton, and rRNA. To begin to quantitate these interactions, L11 from Bacillus stearothermophilus has been overexpressed and its properties compared with those of L11-C76 alone in a fluorescence assay for protein-rRNA binding. The assay relies on 2'-amino-butyryl-pyrene-uridine incorporated in a 58-nucleotide rRNA fragment, which gives approximately 15-fold enhancement when L11 or L11-C76 is bound. Although the pyrene tag weakens protein binding, unbiased protein-RNA association constants were obtained in competition experiments with untagged RNA. It was found that (i) intact B. stearothermophilus L11 binds rRNA with K approximately 1.2 x 10(9) m(-1) in buffers with 0.2 m KCl, about 100-fold tighter than Escherichia coli L11; (ii) the N-terminal domain makes a small, salt-dependent contribution to the overall L11-RNA binding affinity (approximately 8-fold enhancement at 0.2 m KCl), (iii) L11 stimulates thiostrepton binding by 2.3 +/- 0.6 x 10(3)-fold, predicting an overall thiostrepton affinity for the ribosome of approximately 10(9) m(-1), and (iv) the yeast homolog of L11 shows no stimulation of thiostrepton binding. The latter observation resolves the question of why eukaryotes are insensitive to the antibiotic. These measurements also show that it is plausible for thiostrepton to compete directly with EF-G.GDP for binding to the L11-RNA complex, and provide a quantitative basis for further studies of L11 function and thiostrepton mechanism.

Analysis of thermal melting curves

T(m) is defined as Temperature of melting or, more accurately, as temperature of midtransition. This term is often used for nucleic acids (DNA and RNA, oligonucleotides and polynucleotides). A thermal denaturation experiment determines the stability of the secondary structure of a DNA or RNA and aids in the choice of the sequences for antisense oligomers or PCR primers. Beyond a simple numerical value (the T(m)), a thermal denaturation experiment, in which the folded fraction of a structure is plotted vs. temperature, yields important thermodynamic information. We present the classic problems encountered during these experiments and try to demonstrate that a number of useful pieces of information can be extracted from these experimental curves.

Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA that resides in the HDV genome and regulates its replication. The native fold of the ribozyme is complex, having two pseudoknots. Earlier work implicated four non-native pairings in slowing pseudoknot formation: Alt 1, Alt 2, Alt 3, and Alt P1. The goal of the present work was design of a kinetically simplified and maximally reactive construct for in vitro mechanistic and structural studies. The initial approach chosen was site-directed mutagenesis in which known alternative pairings were destabilized while leaving the catalytic core intact. Based on prior studies, the G11C/U27Delta double mutant was prepared. However, biphasic kinetics and antisense oligonucleotide response trends opposite those of the well-studied G11C mutant were observed suggesting that new alternative pairings with multiple registers, termed Alt X and Alt Y, had been created. Enzymatic structure mapping of oligonucleotide models supported this notion. This led to a model wherein Alt 2 and the phylogenetically conserved Alt 3 act as "folding guides", facilitating folding of the major population of the RNA molecules by hindering formation of the Alt X and Alt Y registers. Attempts to eliminate the strongest of the Alt X pairings by rational design of a quadruple mutant only resulted in more complex kinetic behavior. In an effort to simultaneously destabilize multiple alternative pairings, studies were carried out on G11C/U27Delta in the presence of urea or increased monovalent ion concentration. Inclusion of physiological ionic strength allowed the goal of monophasic, fast-folding (kobs approximately 60 min(-1)) kinetics to be realized. To account for this, a model is developed wherein Na+, which destabilizes secondary and tertiary structures in the presence of Mg2+, facilitates native folding by destabilizing the multiple alternative secondary structures with a higher-order dependence.

A conserved element in the yeast RNase MRP RNA subunit can participate in a long-range base-pairing interaction

RNase MRP is a ribonucleoprotein endoribonuclease involved in eukaryotic pre-rRNA processing. The enzyme possesses a putatively catalytic RNA subunit, structurally related to that of RNase P. A thorough structure analysis of Saccharomyces cerevisiae MRP RNA, entailing enzymatic and chemical probing, mutagenesis and thermal melting, identifies a previously unrecognised stem that occupies a position equivalent to the P7 stem of RNase P. Inclusion of this P7-like stem confers on yeast MRP RNA a greater degree of similarity to the core RNase P RNA structure than that described previously and better delimits domain 2, the proposed specificity domain. The additional stem is created by participation of a conserved sequence element (ymCR-II) in a long-range base-pairing interaction. There is potential for this base-pairing throughout the known yeast MRP RNA sequences. Formation of a P7-like stem is not required, however, for the pre-rRNA processing or essential function of RNase MRP. Mutants that can base-pair are nonetheless detrimental to RNase MRP function, indicating that the stem will form in vivo but that only the wild-type pairing is accommodated. Although the alternative MRP RNA structure described is clearly not part of the active RNase MRP enzyme, it would be the more stable structure in the absence of protein subunits and the probability that it represents a valid intermediate species in the process of yeast RNase MRP assembly is discussed.

A thermodynamic framework for the magnesium-dependent folding of RNA

The goal of this review is to present a unified picture of the relationship between ion binding and RNA folding based on recent theoretical and computational advances. In particular, we present a model describing how the association of magnesium ions is coupled to the tertiary structure folding of several well-characterized RNA molecules. This model is developed in terms of the nonlinear Poisson-Boltzmann (NLPB) equation, which provides a rigorous electrostatic description of the interaction between Mg(2+) and specific RNA structures. In our description, most of the ions surrounding an RNA behave as a thermally fluctuating ensemble distributed according to a Boltzmann weighted average of the mean electrostatic potential around the RNA. In some cases, however, individual ions near the RNA may shed some of their surrounding waters to optimize their Coulombic interactions with the negatively charged ligands on the RNA. These chelated ions are energetically distinct from the surrounding ensemble and must be treated explicitly. This model is used to explore several different RNA systems that interact differently with Mg(2+). In each case, the NLPB equation accurately describes the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding without requiring any fitted parameters in the calculation. Based on this model, we present a physical description of how Mg(2+) binds and stabilizes specific RNA structures to promote the folding reaction.

Study of the influence of metal ions on tRNA(Phe) thermal unfolding equilibria by UV spectroscopy and multivariate curve resolution

The influence of metal ions (Na(+), Mg(2+) and Cd(2+)) on the thermal unfolding of phenylalanine transfer ribonucleic acid (tRNA(Phe)) was studied by UV spectroscopy-monitored melting experiments. Absorbance data were obtained during the unfolding process in the range 220-340 nm and later analyzed by a multivariate curve resolution approach (MCR-ALS) based on factor analysis. This procedure determines the number of spectroscopically distinct conformations present during the unfolding process and reveals their concentration profiles and pure spectra, without any initial assumption having to be made about the number of steps in the unfolding pathway. From the concentration profiles and pure spectra, information such as T(m) values can be recovered. The results were compared with those obtained previously in spectroscopic and calorimetric unfolding experiments, showing that the multivariate approach recovers information that complements that obtained in traditional spectroscopic melting experiments.

In vitro structure-function studies of the Bacillus subtilis tyrS mRNA antiterminator: evidence for factor-independent tRNA acceptor stem binding specificity

Expression of many aminoacyl-tRNA synthetase, amino acid biosynthesis and transport genes in Bacillus subtilis is controlled at the level of transcription termination using the T box system and requires the formation of specific secondary structures in the mRNA leader region. One structure functions as a transcriptional terminator, while an alternate form, the antiterminator, is necessary for transcription of the downstream coding regions. We have investigated the interaction of antiterminator model RNAs, based on the B.subtilis tyrS antiterminator with tRNA(Tyr) and tRNA acceptor stem models, using a gel shift assay. Binding of the antiterminator RNA to tRNA(Tyr) was dependent on complimentarity with the acceptor end of the tRNA or microhelix; affinity for the microhelix RNA was reduced relative to the tRNA. Alteration of a conserved position in the non-base pairing region of the bulge greatly reduced tRNA binding, consistent with in vivo studies. Therefore, it appears that some of the antiterminator-tRNA binding specificity is dependent on the structure of the antiterminator bulge alone and the complex it forms with tRNA in the absence of additional trans-acting factors. During the course of these studies we also discovered that the antiterminator can form a 'kissing' bulge complex, a unique RNA motif. The ease of formation of this RNA homodimer illustrates the propensity for the bulge of the antiterminator to bind RNA.

RNA folding pathways

A general overview of the questions and problems in RNA folding is presented. Topics include the differences in the folding problems/questions that apply to RNA versus proteins, methods for determination of final structures, folding versus unfolding, resolution of space and time, and conformational switching.

Stabilization of RNA tertiary structure by monovalent cations

The effects of monovalent cations (Li(+), Na(+), K(+), Rb(+), Cs(+), and NH4(+)) on the thermal stability of RNA tertiary structure were investigated by UV melting. We show that with the RNA used here (nucleotides 1051-1108 of Escherichia coli 23 S rRNA with four base substitutions), monovalent cations and Mg(2+) compete in stabilizing the RNA tertiary structure, and that the competition takes place between two boundaries: one where Mg(2+) concentration is zero and the other where it is maximally stabilizing ("saturating"). The pattern of competition is the same for all monovalent cations and depends on the cation's ability to displace Mg(2+) from the RNA, its ability to stabilize tertiary structure in the absence of Mg(2+), and its ability to stabilize tertiary structure at saturating Mg(2+) concentrations. The stabilizing ability of a monovalent cation depends on its unhydrated ionic radius, and at a low monovalent cation concentration and saturating Mg(2+), there is a (calculated) net release of a single monovalent cation/RNA molecule when tertiary structure is denatured. The implications are that under these conditions there is at least one binding site for monovalent cations on the RNA, the site is specifically associated with formation of stable tertiary structure, K(+) is the most effective of the tested cations, and Mg(2+) appears ineffective at this site. At high ionic strength, and in the absence of Mg(2+), stabilization of tertiary structure is still monovalent-cation specific and ionic-radius dependent, but a larger number of cations ( approximately eight) are released upon RNA tertiary structure denaturation, and NH(4)(+) appears to be the most effective cation in stabilizing tertiary structure under these conditions. In the majority of the experiments, methanol was added as a cosolvent to the buffer. Its use allowed the examination of the behavior of monovalent ions under conditions where their effects would otherwise have been too weak to be observed. Methanol stabilizes tertiary but not secondary structure of the RNA. There was no evidence that it either causes qualitative changes in cation-binding properties of the RNA or a change in the pattern of monovalent cation/Mg(2+) competition.

Clinical severity and thermodynamic effects of iron-responsive element mutations in hereditary hyperferritinemia-cataract syndrome

Hereditary hyperferritinemia-cataract syndrome (HHCS) is a novel genetic disorder characterized by elevated serum ferritin and early onset cataract formation. The excessive ferritin production in HHCS patients arises from aberrant regulation of L-ferritin translation caused by mutations within the iron-responsive element (IRE) of the L-ferritin transcript. IREs serve as binding sites for iron regulatory proteins (IRPs), iron-sensing proteins that regulate ferritin translation. Previous observations suggested that each unique HHCS mutation conferred a characteristic degree of hyperferritinemia and cataract severity in affected individuals. Here we have measured the in vitro affinity of the IRPs for the mutant IREs and correlated decreases in binding affinity with clinical severity. Thermodynamic analysis of these IREs has also revealed that although some HHCS mutations lead to changes in the stability and secondary structure of the IRE, others appear to disrupt IRP-IRE recognition with minimal effect on IRE stability. HHCS is a noteworthy example of a human genetic disorder that arises from mutations within a protein-binding site of an mRNA cis-acting element. Analysis of the effects of these mutations on the energetics of the RNA-protein interaction explains the phenotypic variabilities of the disease state.

Thermodynamics of RNA unfolding: stabilization of a ribosomal RNA tertiary structure by thiostrepton and ammonium ion

RNAs with interesting secondary and tertiary structures tend to melt in several broad and overlapping transitions over a wide temperature range, and it has been consequently difficult to resolve the thermodynamics of individual unfolding steps. In the case that a ligand selectively binds a single folded state of the RNA, it is possible to obtain reliable thermodynamic parameters for both RNA unfolding and RNA-ligand binding simply from the hyperchromicity of RNA denaturation. The analysis procedure involves fitting a three-dimensional surface to absorbance data collected as a function of both temperature and ligand concentration. Analysis of the unfolding of a fragment of the large subunit ribosomal RNA (Escherichia coli sequence 1051 to 1109) is presented; both an antibiotic (thiostrepton) and ammonium ion specifically stabilize a tertiary structure within this RNA. A consistent set of thermodynamic parameters (delta H and tm) for the first two sequentially linked unfolding transitions is obtained from the experiments, and the binding constants obtained for the two ligands are consistent with other independent measurements. The approach is applicable to a variety of RNAs that specifically bind proteins, antibiotics, ions or other ligands.