Probing the dynamics of the P1 helix within the Tetrahymena group I intron

Thoughts on how to think (and talk) about RNA structure

Recent events have pushed RNA research into the spotlight. Continued discoveries of RNA with unexpected diverse functions in healthy and diseased cells, such as the role of RNA as both the source and countermeasure to a severe acute respiratory syndrome coronavirus 2 infection, are igniting a new passion for understanding this functionally and structurally versatile molecule. Although RNA structure is key to function, many foundational characteristics of RNA structure are misunderstood, and the default state of RNA is often thought of and depicted as a single floppy strand. The purpose of this perspective is to help adjust mental models, equipping the community to better use the fundamental aspects of RNA structural information in new mechanistic models, enhance experimental design to test these models, and refine data interpretation. We discuss six core observations focused on the inherent nature of RNA structure and how to incorporate these characteristics to better understand RNA structure. We also offer some ideas for future efforts to make validated RNA structural information available and readily used by all researchers.

Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution

Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution^1-3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure-function relationships⁴, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure-function relationships in ribozymes.

Fluorescence-based investigations of RNA-small molecule interactions

Interrogating non-coding RNA structures and functions with small molecules is an area of rapidly increasing interest among biochemists and chemical biologists. However, many biochemical approaches to monitoring RNA structures are time-consuming and low-throughput, and thereby are only of limited utility for RNA-small molecule studies. Fluorescence-based techniques are powerful tools for rapid investigation of RNA conformations, dynamics, and interactions with small molecules. Many fluorescence methods are amenable to high-throughput analysis, enabling library screening for small molecule binders. In this review, we summarize numerous fluorescence-based approaches for identifying and characterizing RNA-small molecule interactions. We describe in detail a high-information content dual-reporter FRET assay we developed to characterize small molecule-induced conformational and stability changes. Our assay is uniquely suited as a platform for both small molecule discovery and thorough characterization of RNA-small molecule binding mechanisms. Given the growing recognition of non-coding RNAs as attractive targets for therapeutic intervention, we anticipate our FRET assay and other fluorescence-based techniques will be indispensable for the development of potent and specific small molecule inhibitors targeting RNA.

Design and Experimental Evolution of trans-Splicing Group I Intron Ribozymes

Group I intron ribozymes occur naturally as cis-splicing ribozymes, in the form of introns that do not require the spliceosome for their removal. Instead, they catalyze two consecutive trans-phosphorylation reactions to remove themselves from a primary transcript, and join the two flanking exons. Designed, trans-splicing variants of these ribozymes replace the 3′-portion of a substrate with the ribozyme’s 3′-exon, replace the 5′-portion with the ribozyme’s 5′-exon, or insert/remove an internal sequence of the substrate. Two of these designs have been evolved experimentally in cells, leading to variants of group I intron ribozymes that splice more efficiently, recruit a cellular protein to modify the substrate’s gene expression, or elucidate evolutionary pathways of ribozymes in cells. Some of the artificial, trans-splicing ribozymes are promising as tools in therapy, and as model systems for RNA evolution in cells. This review provides an overview of the different types of trans-splicing group I intron ribozymes that have been generated, and the experimental evolution systems that have been used to improve them.

Photophysical Characterization of Enhanced 6-Methylisoxanthopterin Fluorescence in Duplex DNA

The structure and dynamic motions of bases in DNA duplexes and other constructs are important for understanding mechanisms of selectivity and recognition of DNA-binding proteins. The fluorescent guanine analogue, 6-methylisoxanthopterin 6-MI, is well suited to this purpose as it exhibits an unexpected 3- to 4-fold increase in relative quantum yield upon duplex formation when incorporated into the following sequences: ATFAA, AAFTA, or ATFTA (where F represents 6-MI). To better understand some of the factors leading to the 6-MI fluorescence increase upon duplex formation, we characterized the effect of local sequence and structural perturbations on 6-MI photophysics through temperature melts, quantum yield measurements, fluorescence quenching assays, and fluorescence lifetime measurements. By examining 21 sequences we have determined that the duplex-enhanced fluorescence (DEF) depends on the composition of bases adjacent to 6-MI and the presence of adenines at locations n ± 2 from the probe. Investigation of duplex stability and local solvent accessibility measurements support a model in which the DEF arises from a constrained geometry of 6-MI in the duplex, which remains H-bonded to cytosine, stacked with adjacent bases and inaccessible to quenchers. Perturbation of DNA structure through the introduction of an unpaired base 3' to 6-MI or a mismatched basepair increases 6-MI dynamic motion leading to fluorescence quenching and a reduction in quantum yield. Molecular dynamics simulations suggest the enhanced fluorescence results from a greater degree of twist at the X-F step relative to the quenched duplexes examined. These results point to a model where adenine residues located at n ± 2 from 6-MI induce a structural geometry with greater twist in the duplex that hinders local motion reducing dynamic quenching and producing an increase in 6-MI fluorescence.

Increased efficiency of evolved group I intron spliceozymes by decreased side product formation

The group I intron ribozyme from Tetrahymena was recently reengineered into a trans-splicing variant that is able to remove 100-nt introns from pre-mRNA, analogous to the spliceosome. These spliceozymes were improved in this study by 10 rounds of evolution in Escherichia coli cells. One clone with increased activity in E. coli cells was analyzed in detail. Three of its 10 necessary mutations extended the substrate binding duplexes, which led to increased product formation and reduced cleavage at the 5'-splice site. One mutation in the conserved core of the spliceozyme led to a further reduction of cleavage at the 5'-splice site but an increase in cleavage side products at the 3'-splice site. The latter was partially reduced by six additional mutations. Together, the mutations increased product formation while reducing activity at the 5'-splice site and increasing activity at the 3'-splice site. These results show the adaptation of a ribozyme that evolved in nature for cis-splicing to trans-splicing, and they highlight the interdependent function of nucleotides within group I intron ribozymes. Implications for the possible use of spliceozymes as tools in research and therapy, and as a model for the evolution of the spliceosome, are discussed.

Unveiling Inherent Degeneracies in Determining Population-Weighted Ensembles of Interdomain Orientational Distributions Using NMR Residual Dipolar Couplings: Application to RNA Helix Junction Helix Motifs

A growing number of studies employ time-averaged experimental data to determine dynamic ensembles of biomolecules. While it is well-known that different ensembles can satisfy experimental data to within error, the extent and nature of these degeneracies, and their impact on the accuracy of the ensemble determination remains poorly understood. Here, we use simulations and a recently introduced metric for assessing ensemble similarity to explore degeneracies in determining ensembles using NMR residual dipolar couplings (RDCs) with specific application to A-form helices in RNA. Various target ensembles were constructed representing different domain-domain orientational distributions that are confined to a topologically restricted (<10%) conformational space. Five independent sets of ensemble averaged RDCs were then computed for each target ensemble and a "sample and select" scheme used to identify degenerate ensembles that satisfy RDCs to within experimental uncertainty. We find that ensembles with different ensemble sizes and that can differ significantly from the target ensemble (by as much as ∑Ω ∼ 0.4 where ∑Ω varies between 0 and 1 for maximum and minimum ensemble similarity, respectively) can satisfy the ensemble averaged RDCs. These deviations increase with the number of unique conformers and breadth of the target distribution, and result in significant uncertainty in determining conformational entropy (as large as 5 kcal/mol at T = 298 K). Nevertheless, the RDC-degenerate ensembles are biased toward populated regions of the target ensemble, and capture other essential features of the distribution, including the shape. Our results identify ensemble size as a major source of uncertainty in determining ensembles and suggest that NMR interactions such as RDCs and spin relaxation, on their own, do not carry the necessary information needed to determine conformational entropy at a useful level of precision. The framework introduced here provides a general approach for exploring degeneracies in ensemble determination for different types of experimental data.

Roles of long-range tertiary interactions in limiting dynamics of the Tetrahymena group I ribozyme

We determined the effects of mutating the long-range tertiary contacts of the Tetrahymena group I ribozyme on the dynamics of its substrate helix (referred to as P1) and on catalytic activity. Dynamics were assayed by fluorescence anisotropy of the fluorescent base analogue, 6-methyl isoxanthopterin, incorporated into the P1 helix, and fluorescence anisotropy and catalytic activity were measured for wild type and mutant ribozymes over a range of conditions. Remarkably, catalytic activity correlated with P1 anisotropy over 5 orders of magnitude of activity, with a correlation coefficient of 0.94. The functional and dynamic effects from simultaneous mutation of the two long-range contacts that weaken P1 docking are cumulative and, based on this RNA's topology, suggest distinct underlying origins for the mutant effects. Tests of mechanistic predictions via single molecule FRET measurements of rate constants for P1 docking and undocking suggest that ablation of the P14 tertiary interaction frees P2 and thereby enhances the conformational space explored by the undocked attached P1 helix. In contrast, mutation of the metal core tertiary interaction disrupts the conserved core into which the P1 helix docks. Thus, despite following a single correlation, the two long-range tertiary contacts facilitate P1 helix docking by distinct mechanisms. These results also demonstrate that a fluorescence anisotropy probe incorporated into a specific helix within a larger RNA can report on changes in local helical motions as well as differences in more global dynamics. This ability will help uncover the physical properties and behaviors that underlie the function of RNAs and RNA/protein complexes.

In vivo evolution of a catalytic RNA couples trans-splicing to translation

How does a non-coding RNA evolve in cells? To address this question experimentally we evolved a trans-splicing variant of the group I intron ribozyme from Tetrahymena over 21 cycles of evolution in E.coli cells. Sequence variation was introduced during the evolution by mutagenic and recombinative PCR, and increasingly active ribozymes were selected by their repair of an mRNA mediating antibiotic resistance. The most efficient ribozyme contained four clustered mutations that were necessary and sufficient for maximum activity in cells. Surprisingly, these mutations did not increase the trans-splicing activity of the ribozyme. Instead, they appear to have recruited a cellular protein, the transcription termination factor Rho, and facilitated more efficient translation of the ribozyme’s trans-splicing product. In addition, these mutations affected the expression of several other, unrelated genes. These results suggest that during RNA evolution in cells, four mutations can be sufficient to evolve new protein interactions, and four mutations in an RNA molecule can generate a large effect on gene regulation in the cell.

A single-stranded junction modulates nanosecond motional ordering of the substrate recognition duplex of a group I ribozyme

Rigid spinning: Site-directed spin-labeling studies using a rigid nitroxide spin label (Ç) reveal that both length and sequence of a single-stranded junction (J1/2) modulate nanosecond motional ordering of the substrate-recognition duplex (P1) of the 120 kD group I ribozyme. The studies demonstrate an approach for experimental measurements of nanosecond dynamics in high-molecular-weight RNA complexes.

Applying 6-methylisoxanthopterin-enhanced fluorescence to examine protein-DNA interactions in the picomolar range

Incorporation of fluorescent nucleoside analogues into duplex DNA usually leads to a reduction in quantum yield, which significantly limits their potential use and application. We have identified two pentamer DNA sequences containing 6-methylisoxanthopterin (6-MI) (ATFAA and AAFTA, where F is 6-MI) that exhibit significant enhancement of fluorescence upon formation of duplex DNA with quantum yields close to that of monomeric 6-MI. The enhanced fluorescence dramatically increases the utility and sensitivity of the probe and is used to study protein-DNA interactions of nanomolar specificity in this work. The increased sensitivity of 6-MI allows anisotropy binding measurements to be performed at DNA concentrations of 1 nM and fluorescence intensity measurements at 50 pM DNA. The ATFAA sequence was incorporated into DNA constructs to measure the binding affinity of four different protein-DNA interactions that exhibit sequence-specific and non-sequence-specific recognition. In all cases, the K(d) values obtained were consistent with previously reported values measured by other methods. Time-resolved and steady-state fluorescence measurements demonstrate that 6-MI fluorescence is very sensitive to local distortion and reports on different degrees of protein-induced perturbations with single-base resolution, where the largest changes occur at the site of protein binding.

Computational prediction of efficient splice sites for trans-splicing ribozymes

Group I introns have been engineered into trans-splicing ribozymes capable of replacing the 3'-terminal portion of an external mRNA with their own 3'-exon. Although this design makes trans-splicing ribozymes potentially useful for therapeutic application, their trans-splicing efficiency is usually too low for medical use. One factor that strongly influences trans-splicing efficiency is the position of the target splice site on the mRNA substrate. Viable splice sites are currently determined using a biochemical trans-tagging assay. Here, we propose a rapid and inexpensive alternative approach to identify efficient splice sites. This approach involves the computation of the binding free energies between ribozyme and mRNA substrate. We found that the computed binding free energies correlate well with the trans-splicing efficiency experimentally determined at 18 different splice sites on the mRNA of chloramphenicol acetyl transferase. In contrast, our results from the trans-tagging assay correlate less well with measured trans-splicing efficiency. The computed free energy components suggest that splice site efficiency depends on the following secondary structure rearrangements: hybridization of the ribozyme's internal guide sequence (IGS) with mRNA substrate (most important), unfolding of substrate proximal to the splice site, and release of the IGS from the 3'-exon (least important). The proposed computational approach can also be extended to fulfill additional design requirements of efficient trans-splicing ribozymes, such as the optimization of 3'-exon and extended guide sequences.

A role for a single-stranded junction in RNA binding and specificity by the Tetrahymena group I ribozyme

We have investigated the role of a single-stranded RNA junction, J1/2, that connects the substrate-containing P1 duplex to the remainder of the Tetrahymena group I ribozyme. Single-turnover kinetics, fluorescence anisotropy, and single-molecule fluorescence resonance energy transfer studies of a series of J1/2 mutants were used to probe the sequence dependence of the catalytic activity, the P1 dynamics, and the thermodynamics of docking of the P1 duplex into the ribozyme's catalytic core. We found that A29, the center A of three adenosine residues in J1/2, contributes 2 orders of magnitude to the overall ribozyme activity, and double-mutant cycles suggested that J1/2 stabilizes the docked state of P1 over the undocked state via a tertiary interaction involving A29 and the first base pair in helix P2 of the ribozyme, A31·U56. Comparative sequence analysis of this group I intron subclass suggests that the A29 interaction sets one end of a molecular ruler whose other end specifies the 5'-splice site and that this molecular ruler is conserved among a subclass of group I introns related to the Tetrahymena intron. Our results reveal substantial functional effects from a seemingly simple single-stranded RNA junction and suggest that junction sequences may evolve rapidly to provide important interactions in functional RNAs.

RNA dynamics: perspectives from spin labels

Dynamics are important and indispensible physical attributes that play essential roles in RNA function. RNA dynamics are complex, spanning vast timescales, and encompassing a large number of physical modes. The technique of site-directed spin labeling (SDSL), which derives information on local structural and dynamic features of a macromolecule by monitoring a chemically stable nitroxide radical using electron paramagnetic resonance spectroscopy, has been applied to monitor intrinsic dynamics at defined structural states as well as to probe conformational transition dynamics of RNAs. The current state of SDSL studies of RNA dynamics is summarized here. Further development and application of SDSL promise to open up many more opportunities for probing RNA dynamics and connecting dynamics to structure and function.

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.

Synthesis, photophysical characterization, and enzymatic incorporation of a microenvironment-sensitive fluorescent uridine analog

The synthesis of a microenvironment-sensitive base-modified fluorescent ribonucleoside analog based on a 5-(benzo[b]thiophen-2-yl)pyrimidine core, enzymatic incorporation of its corresponding triphosphate into RNA oligonucleotides, and photophysical characterization of fluorescently modified oligoribonucleotides are described.

Fluorescence polarization anisotropy to measure RNA dynamics

RNA requires helical motion to fold and carry out its function. As RNA helical motion occurs on the nanosecond timescale, the timescale probed by fluorescence dyes, fluorescence polarization anisotropy (FPA) is a simple, yet powerful, technique to study helical dynamics in RNA. With the recent development of several fluorescent base analogs that have a nanosecond timescale lifetime in a duplex, FPA has begun to be used for characterizing RNA dynamics. Using the probe 6-methylisoxanthopterin (6-MI) as an example, we describe the procedure for carrying out FPA experiments on model oligonucleotide systems and in a complex RNA, the Tetrahymena group I intron. For smaller RNA systems, isolating the motion of the target helix from the overall tumbling of the whole RNA system is necessary, and nucleic acids binding proteins can be incorporated into the RNA system to increase the overall size of the system, slow the overall tumbling, and thereby reduce the anisotropy contribution from the overall tumbling to negligible. The procedure for incorporating one such protein, the Lac Repressor, is given as an example.

Versatile phosphoramidation reactions for nucleic acid conjugations with peptides, proteins, chromophores, and biotin derivatives

Chemical conjugations of nucleic acids with macromolecules or small molecules are common approaches to study nucleic acids in chemistry and biology and to exploit nucleic acids for medical applications. The conjugation of nucleic acids such as oligonucleotides with peptides is especially useful to circumvent cell delivery and specificity problems of oligonucleotides as therapeutic agents. However, current approaches are limited and inefficient in their ability to afford peptide-oligonucleotide conjugates (POCs). Here, we report an effective and reproducible approach to prepare POCs and other nucleic acid conjugates based on a newly developed nucleic acid phosphoramidation method. The development of a new nucleic acid phosphoramidation reaction was achieved by our successful synthesis of a novel amine-containing biotin derivative used to systematically optimize the reactions. The improved phosphoramidation reactions dramatically increased yields of nucleic acid-biotin conjugates up to 80% after 3 h reaction. Any nucleic acids with a terminal phosphate group are suitable reactants in phosphoramidation reactions to conjugate with amine-containing molecules such as biotin and fluorescein derivatives, proteins, and, most importantly, peptides to enable the synthesis of POCs for therapeutic applications. Polymerase chain reactions (PCRs) to study incorporation of biotin or fluorescein-tagged DNA primers into the reaction products demonstrated that appropriate controls of nucleic acid phosphoramidation reactions incur minimum adverse effects on inherited base-pairing characteristics of nucleotides in nucleic acids. The phosphoramidation approach preserves the integrity of hybridization specificity in nucleic acids when preparing POCs. By retaining integrity of the nucleic acids, their effectiveness as therapeutic reagents for gene silencing, gene therapy, and RNA interference is ensured. The potential for POC use was demonstrated by two-step phosphoramidation reactions to successfully synthesize nucleic acid-tetraglycine conjugates. In addition, phosphoramidation reactions provided a facile approach to prepare nucleic acid-BSA conjugates with good yields. In summary, the new approach to phosphoramidation reactions offers a universal method to prepare POCs and other nucleic acid conjugates with high yields in aqueous solutions. The methods can be easily adapted to typical chemistry or biology laboratory setups which will expedite the applications of POCs for basic research and medicine.

Photophysical and structural properties of the fluorescent nucleobase analogues of the tricyclic cytosine (tC) family

Fundamental insight into the unique fluorescence and nucleobase-mimicking properties of the fluorescent nucleobase analogues of the tC family is not only vital in explaining the behaviour of these probes in nucleic acid environments, but will also be profitable in the development of new and improved fluorescent base analogues. Here, temperature-dependent fluorescence quantum yield measurements are used to successfully separate and quantify the temperature-dependent and temperature-independent non-radiative excited-state decay processes of the three nucleobase analogues tC, tC(O) and tC(nitro); all of which are derivatives of a phenothiazine or phenoxazine tricyclic framework. These results strongly suggest that the non-radiative decay process dominating the fast deactivation of tC(nitro) is an internal conversion of a different origin than the decay pathways of tC and tC(O). tC(nitro) is reported to be fluorescent only in less dipolar solvents at room temperature, which is explained by an increase in excited-state dipole moment along the main non-radiative decay pathway, a suggestion that applies in the photophysical discussion of large polycyclic nitroaromatics in general. New insight into the ground and excited-state potential energy surfaces of the isolated tC bases is obtained by means of high level DFT and TDDFT calculations. The S(0) potential energy surfaces of tC and tC(nitro) possess two global minima corresponding to geometries folded along the middle sulfur-nitrogen axis separated by an energy barrier of 0.05 eV as calculated at the B3LYP/6-311+G(2d,p) level. The ground-state potential energy surface of tC(O) is also predicted to be shallow along the bending coordinate but with an equilibrium geometry corresponding to the planar conformation of the tricyclic framework, which may explain some of the dissimilar properties of tC and tC(O) in various confined (biological) environments. The S(1) equilibrium geometries of all three base analogues are predicted to be planar. These results are discussed in the context of the tC bases positioned in double-stranded DNA scenarios.