Structural specificity conferred by a group I RNA peripheral element

You get what you screen for: Standards for experimental design and data fitting in drug discovery

A common mantra in drug discovery is that "You get what you screen for." This is not a promise that you will always get an effective drug candidate, but rather a warning that inaccuracies in your protocol for screening will more likely produce a compound that fails to be an effective candidate because it matches the properties of your screen, not the desired features of an ideal lead compound. It is with this in mind that we examine the current protocols for evaluating drug candidates and highlight some deficiencies while pointing the way to better methods. Many of the errors in data fitting can be rectified by abandoning the traditional equation-based data fitting methods and adopting the more rigorous mechanism-based fitting afforded by computer simulation based on numerical integration of rate equations. Using these methods bypasses the errors in judgement in choosing the appropriate equation for data fitting and the approximations required to derive those equations. In this chapter we outline the limitations and systematic errors in conventional methods of data fitting and illustrate the advantages of computer simulation and introduce the methods of analysis.

History of advances in enzyme kinetic methods: From minutes to milliseconds

The last review on transient-state kinetic methods in The Enzymes was published three decades ago (Johnson, K.A., 1992. The Enzymes, XX, 1-61). In that review the foundations were laid out for the logic behind the design and interpretation of experiments. In the intervening years the instrumentation has improved mainly in providing better sample economy and shorter dead times. More significantly, in 1992 we were just introducing methods for fitting data based on numerical integration of rate equations, but the software was slow and difficult to use. Today, advances in numerical integration methods for data fitting have led to fast and dynamic software, making it easy to fit data without simplifying approximations. This approach overcomes multiple disadvantages of traditional data fitting based on equations derived by analytical integration of rate equations, requiring simplifying approximations. Mechanism-based fitting using computer simulation resolves mechanisms by accounting for the concentration dependence of the rates and amplitudes of the reaction to find a set of intrinsic rate constants that reproduce the experimental data, including details about how the experiment was performed in modeling the data. Rather than discuss how to design and interpret rapid-quench and stopped-flow experiments individually, we now focus on how to fit them simultaneously so that the quench-flow data anchor the interpretation of fluorescence signals. Computer simulation streamlines the fitting of multiple experiments globally to yield a single unifying model to account for all available data.

Topological crossing in the misfolded Tetrahymena ribozyme resolved by cryo-EM

Taking advantage of single-particle cryogenic electron microscopy (cryo-EM) to analyze highly heterogeneous or flexible samples, we obtained long-awaited three-dimensional (3D) structures of the misfolded Tetrahymena ribozyme. These structures provide clear evidence for a previously proposed topological isomer model, in which the stereochemically impossible crossing of two core RNA strands prevents rapid rearrangement of the misfolded state to the native state. Topological isomers may be widespread in misfolding of complex RNA, and these cryo-EM structures set a foundation for dissecting their detailed kinetic mechanisms and functional consequences in a paradigmatic model system.

The Tetrahymena group I intron has been a key system in the understanding of RNA folding and misfolding. The molecule folds into a long-lived misfolded intermediate (M) in vitro, which has been known to form extensive native-like secondary and tertiary structures but is separated by an unknown kinetic barrier from the native state (N). Here, we used cryogenic electron microscopy (cryo-EM) to resolve misfolded structures of the Tetrahymena L-21 ScaI ribozyme. Maps of three M substates (M1, M2, M3) and one N state were achieved from a single specimen with overall resolutions of 3.5 Å, 3.8 Å, 4.0 Å, and 3.0 Å, respectively. Comparisons of the structures reveal that all the M substates are highly similar to N, except for rotation of a core helix P7 that harbors the ribozyme’s guanosine binding site and the crossing of the strands J7/3 and J8/7 that connect P7 to the other elements in the ribozyme core. This topological difference between the M substates and N state explains the failure of 5′-splice site substrate docking in M, supports a topological isomer model for the slow refolding of M to N due to a trapped strand crossing, and suggests pathways for M-to-N refolding.

Molecular chaperones maximize the native state yield on biological times by driving substrates out of equilibrium

Molecular chaperones facilitate the folding of proteins and RNA in vivo. Under physiological conditions, the in vitro folding of Tetrahymena ribozyme by the RNA chaperone CYT-19 behaves paradoxically; increasing the chaperone concentration reduces the yield of native ribozymes. In contrast, the protein chaperone GroEL works as expected; the yield of the native substrate increases with chaperone concentration. The discrepant chaperone-assisted ribozyme folding thus contradicts the expectation that it operates as an efficient annealing machine. To resolve this paradox, we propose a minimal stochastic model based on the Iterative Annealing Mechanism (IAM) that offers a unified description of chaperone-mediated folding of both proteins and RNA. Our theory provides a general relation that quantitatively predicts how the yield of native states depends on chaperone concentration. Although the absolute yield of native states decreases in the Tetrahymena ribozyme, the product of the folding rate and the steady-state native yield increases in both cases. By using energy from ATP hydrolysis, both CYT-19 and GroEL drive their substrate concentrations far out of equilibrium, thus maximizing the native yield in a short time. This also holds when the substrate concentration exceeds that of GroEL. Our findings satisfy the expectation that proteins and RNA be folded by chaperones on biologically relevant time scales, even if the final yield is lower than what equilibrium thermodynamics would dictate. The theory predicts that the quantity of chaperones in vivo has evolved to optimize native state production of the folded states of RNA and proteins in a given time.

Engineering Structurally Interacting RNA (sxRNA)

RNA-based three-way junctions (3WJs) are naturally occurring structures found in many functional RNA molecules including rRNA, tRNA, snRNA and ribozymes. 3WJs are typically characterized as resulting from an RNA molecule folding back on itself in cis but could also form in trans when one RNA, for instance a microRNA binds to a second structured RNA, such as a mRNA. Trans-3WJs can influence the final shape of one or both of the RNA molecules and can thus provide a means for modulating the availability of regulatory motifs including potential protein or microRNA binding sites. Regulatory 3WJs generated in trans represent a newly identified regulatory category that we call structurally interacting RNA or sxRNA for convenience. Here we show that they can be rationally designed using familiar cis-3WJ examples as a guide. We demonstrate that an sxRNA “bait” sequence can be designed to interact with a specific microRNA “trigger” sequence, creating a regulatable RNA-binding protein motif that retains its functional activity. Further, we show that when placed downstream of a coding sequence, sxRNA can be used to switch “ON” translation of that sequence in the presence of the trigger microRNA and the amount of translation corresponded with the amount of microRNA present.

Use of a Fluorescent Aptamer RNA as an Exonic Sequence to Analyze Self-Splicing Ability of aGroup I Intron from Structured RNAs

Group I self-splicing intron constitutes an important class of functional RNA molecules that can promote chemical transformation. Although the fundamental mechanism of the auto-excision from its precursor RNA has been established, convenient assay systems for its splicing activity are still useful for a further understanding of its detailed mechanism and of its application. Because some host RNA sequences, to which group I introns inserted form stable three-dimensional (3D) structures, the effects of the 3D structures of exonic elements on the splicing efficiency of group I introns are important but not a fully investigated issue. We developed an assay system for group I intron self-splicing by employing a fluorescent aptamer RNA (spinach RNA) as a model exonic sequence inserted by the Tetrahymena group I intron. We investigated self-splicing of the intron from spinach RNA, serving as a model exonic sequence with a 3D structure.

RNA Structural Modules Control the Rate and Pathway of RNA Folding and Assembly

Structured RNAs fold through multiple pathways, but we have little understanding of the molecular features that dictate folding pathways and determine rates along a given pathway. Here, we asked whether folding of a complex RNA can be understood from its structural modules. In a two-piece version of the Tetrahymena group I ribozyme, the separated P5abc subdomain folds to local native secondary and tertiary structure in a linked transition and assembles with the ribozyme core via three tertiary contacts: a kissing loop (P14), a metal core-receptor interaction, and a tetraloop-receptor interaction, the first two of which are expected to depend on native P5abc structure from the local transition. Native gel, NMR, and chemical footprinting experiments showed that mutations that destabilize the native P5abc structure slowed assembly up to 100-fold, indicating that P5abc folds first and then assembles with the core by conformational selection. However, rate decreases beyond 100-fold were not observed because an alternative pathway becomes dominant, with nonnative P5abc binding the core and then undergoing an induced-fit rearrangement. P14 is formed in the rate-limiting step along the conformational selection pathway but after the rate-limiting step along the induced-fit pathway. Strikingly, the assembly rate along the conformational selection pathway resembles that of an isolated kissing loop similar to P14, and the rate along the induced-fit pathway resembles that of an isolated tetraloop-receptor interaction. Our results indicate substantial modularity in RNA folding and assembly and suggest that these processes can be understood in terms of underlying structural modules.

Crystal structure of group II intron domain 1 reveals a template for RNA assembly

Although the importance of large noncoding RNAs is increasingly appreciated, our understanding of their structures and architectural dynamics remains limited. In particular, we know little about RNA folding intermediates and how they facilitate the productive assembly of RNA tertiary structures. Here, we report the crystal structure of an obligate intermediate that is required during the earliest stages of group II intron folding. Composed of domain 1 from the Oceanobacillus iheyensis group II intron (266 nucleotides), this intermediate retains native-like features but adopts a compact conformation in which the active site cleft is closed. Transition between this closed and the open (native) conformation is achieved through discrete rotations of hinge motifs in two regions of the molecule. The open state is then stabilized by sequential docking of downstream intron domains, suggesting a 'first come, first folded' strategy that may represent a generalizable pathway for assembly of large RNA and ribonucleoprotein structures.

Suppression of the Arboviruses Dengue and Chikungunya Using a Dual-Acting Group-I Intron Coupled with Conditional Expression of the Bax C-Terminal Domain

In portions of South Asia, vectors and patients co-infected with dengue (DENV) and chikungunya (CHIKV) are on the rise, with the potential for this occurrence in other regions of the world, for example the United States. Therefore, we engineered an antiviral approach that suppresses the replication of both arboviruses in mosquito cells using a single antiviral group I intron. We devised unique configurations of internal, external, and guide sequences that permit homologous recognition and splicing with conserved target sequences in the genomes of both viruses using a single trans-splicing Group I intron, and examined their effectiveness to suppress infections of DENV and CHIKV in mosquito cells when coupled with a proapoptotic 3' exon, ΔN Bax. RT-PCR demonstrated the utility of these introns in trans-splicing the ΔN Bax sequence downstream of either the DENV or CHIKV target site in transformed Aedes albopictus C6/36 cells, independent of the order in which the virus specific targeting sequences were inserted into the construct. This trans-splicing reaction forms DENV or CHIKV ΔN Bax RNA fusions that led to apoptotic cell death as evidenced by annexin V staining, caspase, and DNA fragmentation assays. TCID50-IFA analyses demonstrate effective suppression of DENV and CHIKV infections by our anti-arbovirus group I intron approach. This represents the first report of a dual-acting Group I intron, and demonstrates that we can target DENV and CHIKV RNAs in a sequence specific manner with a single, uniquely configured CHIKV/DENV dual targeting group I intron, leading to replication suppression of both arboviruses, and thus providing a promising single antiviral for the transgenic suppression of multiple arboviruses.

Effective suppression of dengue virus using a novel group-I intron that induces apoptotic cell death upon infection through conditional expression of the Bax C-terminal domain

INTRODUCTION

Approximately 100 million confirmed infections and 20,000 deaths are caused by Dengue virus (DENV) outbreaks annually. Global warming and rapid dispersal have resulted in DENV epidemics in formally non-endemic regions. Currently no consistently effective preventive measures for DENV exist, prompting development of transgenic and paratransgenic vector control approaches. Production of transgenic mosquitoes refractory for virus infection and/or transmission is contingent upon defining antiviral genes that have low probability for allowing escape mutations, and are equally effective against multiple serotypes. Previously we demonstrated the effectiveness of an anti-viral group I intron targeting U143 of the DENV genome in mediating trans-splicing and expression of a marker gene with the capsid coding domain. In this report we examine the effectiveness of coupling expression of ΔN Bax to trans-splicing U143 intron activity as a means of suppressing DENV infection of mosquito cells.

RESULTS

Targeting the conserved DENV circularization sequence (CS) by U143 intron trans-splicing activity appends a 3' exon RNA encoding ΔN Bax to the capsid coding region of the genomic RNA, resulting in a chimeric protein that induces premature cell death upon infection. TCID50-IFA analyses demonstrate an enhancement of DENV suppression for all DENV serotypes tested over the identical group I intron coupled with the non-apoptotic inducing firefly luciferase as the 3' exon. These cumulative results confirm the increased effectiveness of this αDENV-U143-ΔN Bax group I intron as a sequence specific antiviral that should be useful for suppression of DENV in transgenic mosquitoes. Annexin V staining, caspase 3 assays, and DNA ladder observations confirm DCA-ΔN Bax fusion protein expression induces apoptotic cell death.

CONCLUSION

This report confirms the relative effectiveness of an anti-DENV group I intron coupled to an apoptosis-inducing ΔN Bax 3' exon that trans-splices conserved sequences of the 5' CS region of all DENV serotypes and induces apoptotic cell death upon infection. Our results confirm coupling the targeted ribozyme capabilities of the group I intron with the generation of an apoptosis-inducing transcript increases the effectiveness of infection suppression, improving the prospects of this unique approach as a means of inducing transgenic refractoriness in mosquitoes for all serotypes of this important disease.

Structural dynamics of a mitochondrial tRNA possessing weak thermodynamic stability

Folding dynamics are ubiquitously involved in controlling the multivariate functions of RNAs. While the high thermodynamic stabilities of some RNAs favor purely native states at equilibrium, it is unclear whether weakly stable RNAs exist in random, partially folded states or sample well-defined, globally folded conformations. Using a folding assay that precisely tracks the formation of native aminoacylable tRNA, we show that the folding of a weakly stable human mitochondrial (hmt) leucine tRNA is hierarchical with a distinct kinetic folding intermediate. The stabilities of the native and intermediate conformers are separated by only about 1.2 kcal/mol, and the species are readily interconvertible. Comparison of folding dynamics between unmodified and fully modified tRNAs reveals that post-transcriptional modifications produce a more constrained native structure that does not sample intermediate conformations. These structural dynamics may thus be crucial for recognition by some modifying enzymes in vivo, especially those targeting the globular core region, by allowing access to pretransition state conformers. Reduced conformational sampling of the native, modified tRNAs could then permit improved performance in downstream processes of translation. More generally, weak stabilities of small RNAs that fold in the absence of chaperone proteins may facilitate conformational switching that is central to biological function.

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

A two-piece derivative of a group I intron RNA as a platform for designing self-assembling RNA templates to promote Peptide ligation

Multicomponent RNA-peptide complexes are attractive from the viewpoint of artificial design of functional biomacromolecular systems. We have developed self-folding and self-assembling RNAs that serve as templates to assist chemical ligation between two reactive peptides with RNA-binding capabilities. The design principle of previous templates, however, can be applied only to limited classes of RNA-binding peptides. In this study, we employed a two-piece derivative of a group I intron RNA from the Tetrahymena large subunit ribosomal RNA (LSU rRNA) as a platform for new template RNAs. In this group I intron-based self-assembling platform, modules for the recognition of substrate peptides can be installed independently from modules holding the platform structure. The new self-assembling platform allows us to expand the repertoire of substrate peptides in template RNA design.

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.

Kinetics of tRNA folding monitored by aminoacylation

We describe a strategy for tracking Mg²⁺-initiated folding of ³²P-labeled tRNA molecules to their native structures based on the capacity for aminoacylation by the cognate aminoacyl-tRNA synthetase enzyme. The approach directly links folding to function, paralleling a common strategy used to study the folding of catalytic RNAs. Incubation of unfolded tRNA with magnesium ions, followed by the addition of aminoacyl-tRNA synthetase and further incubation, yields a rapid burst of aminoacyl-tRNA formation corresponding to the prefolded tRNA fraction. A subsequent slower increase in product formation monitors continued folding in the presence of the enzyme. Further analysis reveals the presence of a parallel fraction of tRNA that folds more rapidly than the majority of the population. The application of the approach to study the influence of post-transcriptional modifications in folding of Escherichia coli tRNA₁(Gln) reveals that the modified bases increase the folding rate but do not affect either the equilibrium between properly folded and misfolded states or the folding pathway. This assay allows the use of ³²P-labeled tRNA in integrated studies combining folding, post-transcriptional processing, and aminoacylation reactions.

The Azoarcus group I intron ribozyme misfolds and is accelerated for refolding by ATP-dependent RNA chaperone proteins

Structured RNAs traverse complex energy landscapes that include valleys representing misfolded intermediates. In Neurospora crassa and Saccharomyces cerevisiae, efficient splicing of mitochondrial group I and II introns requires the DEAD box proteins CYT-19 and Mss116p, respectively, which promote folding transitions and function as general RNA chaperones. To test the generality of RNA misfolding and the activities of DEAD box proteins in vitro, here we measure native folding of a small group I intron ribozyme from the bacterium Azoarcus by monitoring its catalytic activity. To develop this assay, we first measure cleavage of an oligonucleotide substrate by the prefolded ribozyme. Substrate cleavage is rate-limited by binding and is readily reversible, with an internal equilibrium near unity, such that the amount of product observed is less than the amount of native ribozyme. We use this assay to show that approximately half of the ribozyme folds readily to the native state, whereas the other half forms an intermediate that transitions slowly to the native state. This folding transition is accelerated by urea and increased temperature and slowed by increased Mg(2+) concentration, suggesting that the intermediate is misfolded and must undergo transient unfolding during refolding to the native state. CYT-19 and Mss116p accelerate refolding in an ATP-dependent manner, presumably by disrupting structure in the intermediate. These results highlight the tendency of RNAs to misfold, underscore the roles of CYT-19 and Mss116p as general RNA chaperones, and identify a refolding transition for further dissection of the roles of DEAD box proteins in RNA folding.

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.

Enhanced specificity against misfolding in a thermostable mutant of the Tetrahymena ribozyme

Structured RNAs encode native conformations that are more stable than the vast ensembles of alternative conformations, but how this specificity is evolved is incompletely understood. Here we show that a variant of the Tetrahymena group I intron ribozyme that was generated previously by in vitro selection for enhanced thermostability also displays modestly enhanced specificity against a stable misfolded structure that is globally similar to the native state, despite the absence of selective pressure to increase the energy gap between these structures. The enhanced specificity for native folding arises from mutations in two nucleotides that are close together in space in the native structure, and additional experiments show that these two mutations do not affect the stability of the misfolded conformation relative to the largely unstructured transition state ensemble for interconversion between the native and misfolded conformers. Thus, they selectively stabilize the native state, presumably by strengthening a local tertiary contact network that cannot form in the misfolded conformation. The stabilization is larger in the presence of the peripheral element P5abc, suggesting that cooperative tertiary structure formation plays a key role in the enhanced stability. The increased specificity in the absence of explicit selection suggests that the large energy gap in the wild-type RNA may have arisen analogously, a consequence of selective pressure for stability of the functional structure. More generally, the structural rigidity and intricate networks of contacts in structured RNAs may allow them to evolve substantial structural specificity without explicit negative selection, even against closely related alternative structures.

Trans-acting RNAs as molecular probes for monitoring time-dependent structural change of an RNA complex adapting two structures

As dynamic structural changes are pivotal for the functions of some classes of RNA molecule, it is important to develop methods to monitor structural changes in RNA in a time-dependent manner without chemical modification. Based on previous reports that trans-acting RNAs can be used as probes for analysis and control of 3D structures of target RNAs, we applied this method to monitor time-dependent structural changes in RNA. We designed and performed a proof-of-principle study using a simple model RNA complex that adopts two different structures as a target. The time-dependent structural changes in the target RNA were successfully monitored using two trans-acting RNAs, which stably form a ternary complex with the bimolecular target RNA and act as a catalyst to join two RNA fragments of the target complex, respectively.

A maturase that specifically stabilizes and activates its cognate group I intron at high temperatures

Folding of large structured RNAs into their functional tertiary structures at high temperatures is challenging. Here we show that I-TnaI protein, a small LAGLIDADG homing endonuclease encoded by a group I intron from a hyperthermophilic bacterium, acts as a maturase that is essential for the catalytic activity of this intron at high temperatures and physiological cationic conditions. I-TnaI specifically binds to and induces tertiary packing of the P4-P6 domain of the intron; this RNA-protein complex might serve as a thermostable platform for active folding of the entire intron. Interestingly, the binding affinity of I-TnaI to its cognate intron RNA largely increases with temperature; over 30-fold stronger binding at higher temperatures relative to 37 °C correlates with a switch from an entropy-driven (37 °C) to an enthalpy-driven (55-60 °C) interaction mode. This binding mode may represent a novel strategy how an RNA binding protein can promote the function of its target RNA specifically at high temperatures.

Targeting of highly conserved Dengue virus sequences with anti-Dengue virus trans-splicing group I introns

Background

Dengue viruses (DENV) are one of the most important viral diseases in the world with approximately 100 million infections and 200,000 deaths each year. The current lack of an approved tetravalent vaccine and ineffective insecticide control measures warrant a search for alternatives to effectively combat DENV. The trans-splicing variant of the Tetrahymena thermophila group I intron catalytic RNA, or ribozyme, is a powerful tool for post-transcriptional RNA modification. The nature of the ribozyme and the predictability with which it can be directed makes it a powerful tool for modifying RNA in nearly any cell type without the need for genome-altering gene therapy techniques or dependence on native cofactors.

Results

Several anti-DENV Group I trans-splicing introns (αDENV-GrpIs) were designed and tested for their ability to target DENV-2 NGC genomes in situ. We have successfully targeted two different uracil bases on the positive sense genomic strand within the highly conserved 5'-3' cyclization sequence (CS) region common to all serotypes of DENV with our αDENV-GrpIs. Our ribozymes have demonstrated ability to specifically trans-splice a new RNA sequence downstream of the targeted site in vitro and in transfected insect cells as analyzed by firefly luciferase and RT-PCR assays. The effectiveness of these αDENV-GrpIs to target infecting DENV genomes is also validated in transfected or transformed Aedes mosquito cell lines upon infection with unattenuated DENV-2 NGC.

Conclusions

Analysis shows that our αDENV-GrpIs have the ability to effectively trans-splice the DENV genome in situ. Notably, these results show that the αDENV-GrpI 9v1, designed to be active against all forms of Dengue virus, effectively targeted the DENV-2 NGC genome in a sequence specific manner. These novel αDENV-GrpI introns provide a striking alternative to other RNA based approaches for the transgenic suppression of DENV in transformed mosquito cells and tissues.

Multiple unfolding events during native folding of the Tetrahymena group I ribozyme

Despite the ubiquitous nature of misfolded intermediates in RNA folding, little is known about their physical properties or the folding transitions that allow them to continue folding productively. Folding of the Tetrahymena group I ribozyme includes sequential accumulation of two intermediates, termed I(trap) and misfolded (M). Here, we probe the structure and folding transition of I(trap) and compare them to those of M. Hydroxyl radical and dimethyl sulfate footprinting show that both I(trap) and M are extensively structured and crudely resemble the native RNA. However, regions of the core P3-P8 domain are more exposed to solvent in I(trap) than in M. I(trap) rearranges to continue folding nearly 1000-fold faster than M, and urea accelerates folding of I(trap) much less than M. Thus, the rate-limiting transition from I(trap) requires a smaller increase in exposed surface. Mutations that disrupt peripheral tertiary contacts give large and nearly uniform increases in re-folding of M, whereas the same mutations give at most modest increases in folding from I(trap). Intriguingly, mutations within the peripheral element P5abc give 5- to 10-fold accelerations in escape from I(trap), whereas ablation of P13, which lies on the opposite surface in the native structure, near the P3-P8 domain, has no effect. Thus, the unfolding required from I(trap) appears to be local, whereas the unfolding of M appears to be global. Further, the modest effects from several mutations suggest that there are multiple pathways for escape from I(trap) and that escape is aided by loosening nearby native structural constraints, presumably to facilitate local movements of nucleotides or segments that have not formed native contacts. Overall, these and prior results suggest a model in which the global architecture and peripheral interactions of the RNA are achieved relatively early in folding. Multiple folding and re-folding events occur on the predominant pathway to the native state, with increasing native core interactions and cooperativity as folding progresses.

A rearrangement of the guanosine-binding site establishes an extended network of functional interactions in the Tetrahymena group I ribozyme active site

Protein enzymes appear to use extensive packing and hydrogen bonding interactions to precisely position catalytic groups within active sites. Because of their inherent backbone flexibility and limited side chain repertoire, RNA enzymes face additional challenges relative to proteins in precisely positioning substrates and catalytic groups. Here, we use the group I ribozyme to probe the existence, establishment, and functional consequences of an extended network of interactions in an RNA active site. The group I ribozyme catalyzes a site-specific attack of guanosine on an oligonucleotide substrate. We previously determined that the hydrogen bond between the exocyclic amino group of guanosine and the 2'-hydroxyl group at position A261 of the Tetrahymena group I ribozyme contributes to overall catalysis. We now use functional data, aided by double mutant cycles, to probe this hydrogen bond in the individual reaction steps of the catalytic cycle. Our results indicate that this hydrogen bond is not formed upon guanosine binding to the ribozyme but instead forms at a later stage of the catalytic cycle. Formation of this hydrogen bond is correlated with other structural rearrangements in the ribozyme's active site that are promoted by docking of the oligonucleotide substrate into the ribozyme's active site, and disruption of this interaction has deleterious consequences for the chemical transformation within the ternary complex. These results, combined with earlier results, provide insight into the nature of the multiple conformational steps used by the Tetrahymena group I ribozyme to achieve its active structure and reveal an intricate, extended network of interactions that is used to establish catalytic interactions within this RNA's active site.

Catalytic activity as a probe of native RNA folding

As RNAs fold to functional structures, they traverse complex energy landscapes that include many partially folded and misfolded intermediates. For structured RNAs that possess catalytic activity, this activity can provide a powerful means of monitoring folding that is complementary to biophysical approaches. RNA catalysis can be used to track accumulation of the native RNA specifically and quantitatively, readily distinguishing the native structure from intermediates that resemble it and may not be differentiated by other approaches. Here, we outline how to design and interpret experiments using catalytic activity to monitor RNA folding, and we summarize adaptations of the method that have been used to probe aspects of folding well beyond determination of the folding rates.

Protein roles in group I intron RNA folding: the tyrosyl-tRNA synthetase CYT-18 stabilizes the native state relative to a long-lived misfolded structure without compromising folding kinetics

The Neurospora crassa CYT-18 protein is a mitochondrial tyrosyl-tRNA synthetase that also promotes self-splicing of group I intron RNAs by stabilizing the functional structure in the conserved core. CYT-18 binds the core along the same surface as a common peripheral element, P5abc, suggesting that CYT-18 can replace P5abc functionally. In addition to stabilizing structure generally, P5abc stabilizes the native conformation of the Tetrahymena group I intron relative to a globally similar misfolded conformation that has only local differences within the core and is populated significantly at equilibrium by a ribozyme variant lacking P5abc (E(DeltaP5abc)). Here, we show that CYT-18 specifically promotes formation of the native group I intron core from this misfolded conformation. Catalytic activity assays demonstrate that CYT-18 shifts the equilibrium of E(DeltaP5abc) toward the native state by at least 35-fold, and binding assays suggest an even larger effect. Thus, similar to P5abc, CYT-18 preferentially recognizes the native core, despite the global similarity of the misfolded core and despite forming crudely similar complexes, as revealed by dimethyl sulfate footprinting. Interestingly, the effects of CYT-18 and P5abc on folding kinetics differ. Whereas P5abc inhibits refolding of the misfolded conformation by forming peripheral contacts that must break during refolding, CYT-18 does not display analogous inhibition, most likely because it relies to a greater extent on direct interactions with the core. Although CYT-18 does not encounter this RNA in vivo, our results suggest that it stabilizes its cognate group I introns relative to analogous misfolded intermediates. By specifically recognizing native structural features, CYT-18 may also interact with earlier folding intermediates to avoid RNA misfolding or to trap native contacts as they form. More generally, our results highlight the ability of a protein cofactor to stabilize a functional RNA structure specifically without incurring associated costs in RNA folding kinetics.

Metal ion dependence of cooperative collapse transitions in RNA

Positively charged counterions drive RNA molecules into compact configurations that lead to their biologically active structures. To understand how the valence and size of the cations influences the collapse transition in RNA, small-angle X-ray scattering was used to follow the decrease in the radius of gyration (R(g)) of the Azoarcus and Tetrahymena ribozymes in different cations. Small, multivalent cations induced the collapse of both ribozymes more efficiently than did monovalent ions. Thus, the cooperativity of the collapse transition depends on the counterion charge density. Singular value decomposition of the scattering curves showed that folding of the smaller and more thermostable Azoarcus ribozyme is well described by two components, whereas collapse of the larger Tetrahymena ribozyme involves at least one intermediate. The ion-dependent persistence length, extracted from the distance distribution of the scattering vectors, shows that the Azoarcus ribozyme is less flexible at the midpoint of transition in low-charge-density ions than in high-charge-density ions. We conclude that the formation of sequence-specific tertiary interactions in the Azoarcus ribozyme overlaps with neutralization of the phosphate charge, while tertiary folding of the Tetrahymena ribozyme requires additional counterions. Thus, the stability of the RNA structure determines its sensitivity to the valence and size of the counterions.

Selection of a novel class of RNA-RNA interaction motifs based on the ligase ribozyme with defined modular architecture

To develop molecular tools for the detection and control of RNA molecules whose functions rely on their 3D structures, we have devised a selection system to isolate novel RNA motifs that interact with a target RNA structure within a given structural context. In this system, a GAAA tetraloop and its specific receptor motif (11-ntR) from an artificial RNA ligase ribozyme with modular architecture (the DSL ribozyme) were replaced with a target structure and random sequence, respectively. Motifs recognizing the target structure can be identified by in vitro selection based on ribozyme activity. A model selection targeting GAAA-loop successfully identified motifs previously known as GAAA-loop receptors. In addition, a new selection targeting a C-loop motif also generated novel motifs that interact with this structure. Biochemical analysis of one of the C-loop receptor motifs revealed that it could also function as an independent structural unit.

Kinetic redistribution of native and misfolded RNAs by a DEAD-box chaperone

DExD/H-box proteins are ubiquitously involved in RNA-mediated processes and use ATP to accelerate conformational changes in RNA. However, their mechanisms of action, and what determines which RNA species are targeted, are not well understood. Here we show that the DExD/H-box protein CYT-19, a general RNA chaperone, mediates ATP-dependent unfolding of both the native conformation and a long-lived misfolded conformation of a group I catalytic RNA with efficiencies that depend on the stabilities of the RNA species but not on specific structural features. CYT-19 then allows the RNA to refold, changing the distribution from equilibrium to kinetic control. Because misfolding is favoured kinetically, conditions that allow unfolding of the native RNA yield large increases in the population of misfolded species. Our results suggest that DExD/H-box proteins act with sufficient breadth and efficiency to allow structured RNAs to populate a wider range of conformations than would be present at equilibrium. Thus, RNAs may face selective pressure to stabilize their active conformations relative to inactive ones to avoid significant redistribution by DExD/H-box proteins. Conversely, RNAs whose functions depend on forming multiple conformations may rely on DExD/H-box proteins to increase the populations of less stable conformations, thereby increasing their overall efficiencies.

The molecular evolution and structural organization of group I introns at position 1389 in nuclear small subunit rDNA of myxomycetes

The number of nuclear group I introns from myxomycetes is rapidly increasing in GenBank as more rDNA sequences from these organisms are being sequenced. They represent an interesting and complex group of intervening sequences because several introns are mobile (or inferred to be mobile) and many contain large and unusual insertions in peripheral loops. Here we describe related group I introns at position 1389 in the small subunit rDNA of representatives from the myxomycete family Didymiaceae. Phylogenetic analyses support a common origin and mainly vertical inheritance of the intron. All S1389 introns from the Didymiaceae belong to the IC1 subclass of nuclear group I introns. The central catalytic core region of about 100 nt appears divergent in sequence composition even though the introns reside in closely related species. Furthermore, unlike the majority of group I introns from myxomycetes the S1389 introns do not self-splice as naked RNA in vitro under standard conditions, consistent with a dependence on host factors for folding or activity. Finally, the myxomycete S1389 introns are exclusively found within the family Didymiaceae, which suggests that this group I intron was acquired after the split between the families Didymiaceae and Physaraceae.

Deletion of the P5abc peripheral element accelerates early and late folding steps of the Tetrahymena group I ribozyme

The P5abc peripheral element stabilizes the Tetrahymena group I ribozyme and enhances its catalytic activity. Despite its beneficial effects on the native structure, prior studies have shown that early formation of P5abc structure during folding can slow later folding steps. Here we use a P5abc deletion variant E(deltaP5abc) to systematically probe the role of P5abc throughout tertiary folding. Time-resolved hydroxyl radical footprinting shows that E(deltaP5abc) forms its earliest stable tertiary structure on the millisecond time scale, approximately 5-fold faster than the wild-type ribozyme, and stable structure spreads throughout E(deltaP5abc) in seconds. Nevertheless, activity measurements show that the earliest detectable formation of native E(deltaP5abc) ribozyme is much slower (approximately 0.6 min(-1)), in a manner similar to that of the wild type. Also similar, only a small fraction of E(deltaP5abc) attains the native state on this time scale under standard conditions at 25 degrees C, whereas the remainder misfolds; footprinting experiments show that the misfolded conformer shares structural features with the long-lived misfolded conformer of the wild-type ribozyme. Thus, P5abc does not have a large overall effect on the rate-limiting step(s) along this pathway. However, once misfolded, E(deltaP5abc) refolds to the native state 80-fold faster than the wild-type ribozyme and is less accelerated by urea, indicating that P5abc stabilizes the misfolded structure relative to the less-ordered transition state for refolding. Together, the results suggest that, under these conditions, even the earliest tertiary folding intermediates of the wild-type ribozyme represent misfolded species and that P5abc is principally a liability during the tertiary folding process.

Charge density of divalent metal cations determines RNA stability

RNA molecules are exquisitely sensitive to the properties of counterions. The folding equilibrium of the Tetrahymena ribozyme is measured by nondenaturing gel electrophoresis in the presence of divalent group IIA metal cations. The stability of the folded ribozyme increases with the charge density (zeta) of the cation. Similar scaling is found when the free energy of the RNA folded in small and large metal cations is measured by urea denaturation. Brownian dynamics simulations of a polyelectrolyte show that the experimental observations can be explained by nonspecific ion-RNA interactions in the absence of site-specific metal chelation. The experimental and simulation results establish that RNA stability is largely determined by a combination of counterion charge and the packing efficiency of condensed cations that depends on the excluded volume of the cations.

Nonspecific binding to structured RNA and preferential unwinding of an exposed helix by the CYT-19 protein, a DEAD-box RNA chaperone

We explore the interactions of CYT-19, a DExD/H-box protein that functions in folding of group I RNAs, with a well characterized misfolded species of the Tetrahymena ribozyme. Consistent with its function, CYT-19 accelerates refolding of the misfolded RNA to its native state. Unexpectedly, CYT-19 performs another reaction much more efficiently; it unwinds the 6-bp P1 duplex formed between the ribozyme and its oligonucleotide substrate. Furthermore, CYT-19 performs this reaction 50-fold more efficiently than it unwinds the same duplex free in solution, suggesting that it forms additional interactions with the ribozyme, most likely using a distinct RNA binding site from the one responsible for unwinding. This site can apparently bind double-stranded RNA, as attachment of a simple duplex adjacent to P1 recapitulates much of the activation provided by the ribozyme. Unwinding the native P1 duplex does not accelerate refolding of the misfolded ribozyme, implying that CYT-19 can disrupt multiple contacts on the RNA, consistent with its function in folding of multiple RNAs. Further experiments showed that the P1 duplex unwinding activity is virtually the same whether the ribozyme is misfolded or native but is abrogated by formation of tertiary contacts between the P1 duplex and the body of the ribozyme. Together these results suggest a mechanism for CYT-19 and other general DExD/H-box RNA chaperones in which the proteins bind to structured RNAs and efficiently unwind loosely associated duplexes, which biases the proteins to disrupt nonnative base pairs and gives the liberated strands an opportunity to refold.

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.

Fast, reversible interaction of prion protein with RNA aptamers containing specific sequence patterns

One of the unsolved problems in prion diseases relates to the physiological function of cellular prion protein (PrP), of which a misfolded isoform is the major component of the transmissible spongiform encephalopathies agent. Knowledge of the PrP-binding molecules may help in elucidating its role and understanding the pathological events underlying prion diseases. Because nucleic acids are known to bind PrP, we attempted to identify the preferred RNA sequences that bind to the ovine recombinant PrP. An in vitro selection approach (SELEX) was applied to a pool of 80-nucleotide(nt)-long RNAs containing a randomised 40-nt central region. The most frequently isolated aptamer, RM312, was also the best ligand (20 nM KD value), according to both surface plasmon resonance and filter binding assays. The fast rates of association and dissociation of RM312 with immobilized PrP, which are reminiscent of biologically relevant interactions, could point to a physiological function of PrP towards cellular nucleic acids. The minimal sequence that we found necessary for binding of RM312 to PrP presents a striking similarity with one previously described PrP aptamer of comparable affinity. In addition, we here identify the two lysine clusters contained in the N-terminal part of PrP as its main nucleic-acid binding sites.

Topology of three-way junctions in folded RNAs

The three-way junctions contained in X-ray structures of folded RNAs have been compiled and analyzed. Three-way junctions with two helices approximately coaxially stacked can be divided into three main families depending on the relative lengths of the segments linking the three Watson-Crick helices. Each family has topological characteristics with some conservation in the non-Watson-Crick pairs within the linking segments as well as in the types of contacts between the segments and the helices. The most populated family presents tertiary interactions between two helices as well as extensive shallow/minor groove contacts between a linking segment and the third helix. On the basis of the lengths of the linking segments, some guidelines could be deduced for choosing a topology for a three-way junction on the basis of a secondary structure. Examples and prediction bas'ed on those rules are discussed.

Concordant exploration of the kinetics of RNA folding from global and local perspectives

Time-resolved small-angle X-ray scattering (SAXS) with millisecond time-resolution reveals two discrete phases of global compaction upon Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. Electrostatic relaxation of the RNA occurs rapidly and dominates the first phase of compaction during which the observed radius of gyration (R(g)) decreases from 75 angstroms to 55 angstroms. A further decrease in R(g) to 45 angstroms occurs in a well-defined second phase. An analysis of mutant ribozymes shows that the latter phase depends upon the formation of long-range tertiary contacts within the P4-P6 domain of the ribozyme; disruption of the three remaining long-range contacts linking the peripheral helices has no effect on the 55-45 angstroms compaction transition. A better understanding of the role of specific tertiary contacts in compaction was obtained by concordant time-resolved hydroxyl radical (OH) analyses that report local changes in the solvent accessibility of the RNA backbone. Comparison of the global and local measures of folding shows that formation of a subset of native tertiary contacts (i.e. those defining the ribozyme core) can occur within a highly compact ensemble whose R(g) is close to that of the fully folded ribozyme. Analyses of additional ribozyme mutants and reaction conditions establish the generality of the rapid formation of a partially collapsed state with little to no detectable tertiary structure. These studies directly link global RNA compaction with formation of tertiary structure as the molecule acquires its biologically active structure, and underscore the strong dependence on salt of both local and global measures of folding kinetics.