Sequential folding of transfer RNA. A nuclear magnetic resonance study of successively longer tRNA fragments with a common 5' end

Kinetic pathway of HIV-1 TAR cotranscriptional folding

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

Regulation of Gene Expression Through Effector-dependent Conformational Switching by Cobalamin Riboswitches

Riboswitches are an outstanding example of genetic regulation mediated by RNA conformational switching. In these non-coding RNA elements, the occupancy status of a ligand-binding domain governs the mRNA's decision to form one of two mutually exclusive structures in the downstream expression platform. Temporal constraints upon the function of many riboswitches, requiring folding of complex architectures and conformational switching in a limited co-transcriptional timeframe, make them ideal model systems for studying these processes. In this review, we focus on the mechanism of ligand-directed conformational changes in one of the most widely distributed riboswitches in bacteria: the cobalamin family. We describe the architectural features of cobalamin riboswitches whose structures have been determined by x-ray crystallography, which suggest a direct physical role of cobalamin in effecting the regulatory switch. Next, we discuss a series of experimental approaches applied to several model cobalamin riboswitches that interrogate these structural models. As folding is central to riboswitch function, we consider the differences in folding landscapes experienced by RNAs that are produced in vitro and those that are allowed to fold co-transcriptionally. Finally, we highlight a set of studies that reveal the difficulties of studying cobalamin riboswitches outside the context of transcription and that co-transcriptional approaches are essential for developing a more accurate picture of their structure-function relationships in these switches. This understanding will be essential for future advancements in the use of small-molecule guided RNA switches in a range of applications such as biosensors, RNA imaging tools, and nucleic acid-based therapies.

Landscape Zooming toward the Prediction of RNA Cotranscriptional Folding

RNA molecules fold as they are transcribed. Cotranscriptional folding of RNA plays a critical role in RNA functions in vivo. Present computational strategies focus on simulations where large structural changes may not be completely sampled. Here, we describe an alternative approach to predicting cotranscriptional RNA folding by zooming in and out of the RNA folding energy landscape. By classifying the RNA structural ensemble into "partitions" based on long, stable helices, we zoom out of the landscape and predict the overall slow folding kinetics from the interpartition kinetic network, and for each interpartition transition, we zoom in on the landscape to simulate the kinetics. Applications of the model to the 117-nucleotide E. coli SRP RNA and the 59-nucleotide HIV-1 TAR RNA show agreements with the experimental data and new structural and kinetic insights into biologically significant conformational switches and pathways for these important systems. This approach, by zooming in/out of an RNA folding landscape at different resolutions, might allow us to treat large RNAs in vivo with transcriptional pause, transcription speed, and other in vivo effects.

Transcription complexes as RNA chaperones

To exert their functions, RNAs adopt diverse structures, ranging from simple secondary to complex tertiary and quaternary folds. In vivo, RNA folding starts with RNA transcription, and a wide variety of processes are coupled to co-transcriptional RNA folding events, including the regulation of fundamental transcription dynamics, gene regulation by mechanisms like attenuation, RNA processing or ribonucleoprotein particle formation. While co-transcriptional RNA folding and associated co-transcriptional processes are by now well accepted as pervasive regulatory principles in all organisms, investigations into the role of the transcription machinery in co-transcriptional folding processes have so far largely focused on effects of the order in which RNA regions are produced and of transcription kinetics. Recent structural and structure-guided functional analyses of bacterial transcription complexes increasingly point to an additional role of RNA polymerase and associated transcription factors in supporting co-transcriptional RNA folding by fostering or preventing strategic contacts to the nascent transcripts. In general, the results support the view that transcription complexes can act as RNA chaperones, a function that has been suggested over 30 years ago. Here, we discuss transcription complexes as RNA chaperones based on recent examples from bacterial transcription.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

The Role of RNA Secondary Structure in Regulation of Gene Expression in Bacteria

Due to the high exposition to changing environmental conditions, bacteria have developed many mechanisms enabling immediate adjustments of gene expression. In many cases, the required speed and plasticity of the response are provided by RNA-dependent regulatory mechanisms. This is possible due to the very high dynamics and flexibility of an RNA structure, which provide the necessary sensitivity and specificity for efficient sensing and transduction of environmental signals. In this review, we will discuss the current knowledge about known bacterial regulatory mechanisms which rely on RNA structure. To better understand the structure-driven modulation of gene expression, we describe the basic theory on RNA structure folding and dynamics. Next, we present examples of multiple mechanisms employed by RNA regulators in the control of bacterial transcription and translation.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.

Transient Protein-RNA Interactions Guide Nascent Ribosomal RNA Folding

Ribosome assembly is an efficient but complex and heterogeneous process during which ribosomal proteins assemble on the nascent rRNA during transcription. Understanding how the interplay between nascent RNA folding and protein binding determines the fate of transcripts remains a major challenge. Here, using single-molecule fluorescence microscopy, we follow assembly of the entire 3' domain of the bacterial small ribosomal subunit in real time. We find that co-transcriptional rRNA folding is complicated by the formation of long-range RNA interactions and that r-proteins self-chaperone the rRNA folding process prior to stable incorporation into a ribonucleoprotein (RNP) complex. Assembly is initiated by transient rather than stable protein binding, and the protein-RNA binding dynamics gradually decrease during assembly. This work questions the paradigm of strictly sequential and cooperative ribosome assembly and suggests that transient binding of RNA binding proteins to cellular RNAs could provide a general mechanism to shape nascent RNA folding during RNP assembly.

Monitoring co-transcriptional folding of riboswitches through helicase unwinding

In the cell, RNAs fold and begin to function as they are being transcribed. In contrast, in the laboratory, RNAs are typically studied after transcription is completed. Co-transcriptional folding can regulate the function of riboswitches and ribozymes and dictate the order of ribonucleoprotein assembly. Methods to observe and investigate RNA folding and activity during transcription are therefore desirable, yet synchronizing RNA polymerases and incorporating labels at specific sites for biophysical studies can be challenging. A recent methodological advance has been to harness highly processive, engineered "super-helicases" to unwind hybrid RNA-DNA duplexes, thereby releasing the RNA 5'-3'. When combined with single-molecule fluorescence detection, RNA folding and concomitant activity can be studied in vitro in a manner that mimics vectorial folding during transcription. Herein, we describe methods for designing and preparing fluorescently labeled RNA-DNA duplex substrates for sequential helicase-dependent RNA folding experiments.

Oritatami: A Computational Model for Molecular Co-Transcriptional Folding

We introduce and study the computational power of Oritatami, a theoretical model that explores greedy molecular folding, whereby a molecular strand begins to fold before its production is complete. This model is inspired by our recent experimental work demonstrating the construction of shapes at the nanoscale from RNA, where strands of RNA fold into programmable shapes during their transcription from an engineered sequence of synthetic DNA. In the model of Oritatami, we explore the process of folding a single-strand bit by bit in such a way that the final fold emerges as a space-time diagram of computation. One major requirement in order to compute within this model is the ability to program a single sequence to fold into different shapes dependent on the state of the surrounding inputs. Another challenge is to embed all of the computing components within a contiguous strand, and in such a way that different fold patterns of the same strand perform different functions of computation. Here, we introduce general design techniques to solve these challenges in the Oritatami model. Our main result in this direction is the demonstration of a periodic Oritatami system that folds upon itself algorithmically into a prescribed set of shapes, depending on its current local environment, and whose final folding displays the sequence of binary integers from 0 to N=2k−1 with a seed of size O(k). We prove that designing Oritatami is NP-hard in the number of possible local environments for the folding. Nevertheless, we provide an efficient algorithm, linear in the length of the sequence, that solves the Oritatami design problem when the number of local environments is a small fixed constant. This shows that this problem is in fact fixed parameter tractable (FPT) and can thus be solved in practice efficiently. We hope that the numerous structural strategies employed in Oritatami enabling computation will inspire new architectures for computing in RNA that take advantage of the rapid kinetic-folding of RNA.

Predicting Cotranscriptional Folding Kinetics For Riboswitch

On the basis of a helix-based transition rate model, we developed a new method for sampling cotranscriptional RNA conformational ensemble and the prediction of cotranscriptional folding kinetics. Applications to E. coli. SRP RNA and pbuE riboswitch indicate that the model may provide reliable predictions for the cotranscriptional folding pathways and population kinetics. For E. coli. SRP RNA, the predicted population kinetics and the folding pathway are consistent with the SHAPE profiles in the recent cotranscriptional SHAPE-seq experiments. For the pbuE riboswitch, the model predicts the transcriptional termination efficiency as a function of the force. The theoretical results show (a) a force-induced transition from the aptamer (antiterminator) to the terminator structure and (b) the different folding pathways for the riboswitch with and without the ligand (adenine). More specifically, without adenine, the aptamer structure emerges as a short-lived kinetic transient state instead of a thermodynamically stable intermediate state. Furthermore, from the predicted extension-time curves, the model identifies a series of conformational switches in the pulling process, where the predicted relative residence times for the different structures are in accordance with the experimental data. The model may provide a new tool for quantitative predictions of cotranscriptional folding kinetics, and results can offer useful insights into cotranscriptional folding-related RNA functions such as regulation of gene expression with riboswitches.

Efficient computation of co-transcriptional RNA-ligand interaction dynamics

Riboswitches form an abundant class of cis-regulatory RNA elements that mediate gene expression by binding a small metabolite. For synthetic biology applications, they are becoming cheap and accessible systems for selectively triggering transcription or translation of downstream genes. Many riboswitches are kinetically controlled, hence knowledge of their co-transcriptional mechanisms is essential. We present here an efficient implementation for analyzing co-transcriptional RNA-ligand interaction dynamics. This approach allows for the first time to model concentration-dependent metabolite binding/unbinding kinetics. We exemplify this novel approach by means of the recently studied I-A 2^'-deoxyguanosine (2^'dG)-sensing riboswitch from Mesoplasma florum.

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

RNA folds into intricate structures that are crucial for its functions and regulations. To date, a multitude of approaches for probing structures of the whole transcriptome, i.e., RNA structuromes, have been developed. Applications of these approaches to different cell lines and tissues have generated a rich resource for the study of RNA structure–function relationships at a systems biology level. In this review, we first introduce the designs of these methods and their applications to study different RNA structuromes. We emphasize their technological differences especially their unique advantages and caveats. We then summarize the structural insights in RNA functions and regulations obtained from the studies of RNA structuromes. And finally, we propose potential directions for future improvements and studies.

In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding

Defining the in vivo folding pathway of cellular RNAs is essential to understand how they reach their final native conformation. We here introduce a novel method, named Structural Probing of Elongating Transcripts (SPET-seq), that permits single-base resolution analysis of transcription intermediates' secondary structures on a transcriptome-wide scale, enabling base-resolution analysis of the RNA folding events. Our results suggest that cotranscriptional RNA folding in vivo is a mixture of cooperative folding events, in which local RNA secondary structure elements are formed as they get transcribed, and non-cooperative events, in which 5'-halves of long-range helices get sequestered into transient non-native interactions until their 3' counterparts have been transcribed. Together our work provides the first transcriptome-scale overview of RNA cotranscriptional folding in a living organism.

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

The elusive quest for RNA knots

Physical entanglement, and particularly knots arise spontaneously in equilibrated polymers that are sufficiently long and densely packed. Biopolymers are no exceptions: knots have long been known to occur in proteins as well as in encapsidated viral DNA. The rapidly growing number of RNA structures has recently made it possible to investigate the incidence of physical knots in this type of biomolecule, too. Strikingly, no knots have been found to date in the known RNA structures. In this Point of View Article we discuss the absence of knots in currently available RNAs and consider the reasons why knots in RNA have not yet been found, despite the expectation that they should exist in Nature. We conclude by singling out a number of RNA sequences that, based on the properties of their predicted secondary structures, are good candidates for knotted RNAs.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.

Helix capping in RNA structure

Helices are an essential element in defining the three-dimensional architecture of structured RNAs. While internal basepairs in a canonical helix stack on both sides, the ends of the helix stack on only one side and are exposed to the loop side, thus susceptible to fraying unless they are protected. While coaxial stacking has long been known to stabilize helix ends by directly stacking two canonical helices coaxially, based on analysis of helix-loop junctions in RNA crystal structures, herein we describe helix capping, topological stacking of a helix end with a basepair or an unpaired nucleotide from the loop side, which in turn protects helix ends. Beyond the topological protection of helix ends against fraying, helix capping should confer greater stability onto the resulting composite helices. Our analysis also reveals that this general motif is associated with the formation of tertiary structure interactions. Greater knowledge about the dynamics at the helix-junctions in the secondary structure should enhance the prediction of RNA secondary structure with a richer set of energetic rules and help better understand the folding of a secondary structure into its three-dimensional structure. These together suggest that helix capping likely play a fundamental role in driving RNA folding.

On the importance of cotranscriptional RNA structure formation

The expression of genes, both coding and noncoding, can be significantly influenced by RNA structural features of their corresponding transcripts. There is by now mounting experimental and some theoretical evidence that structure formation in vivo starts during transcription and that this cotranscriptional folding determines the functional RNA structural features that are being formed. Several decades of research in bioinformatics have resulted in a wide range of computational methods for predicting RNA secondary structures. Almost all state-of-the-art methods in terms of prediction accuracy, however, completely ignore the process of structure formation and focus exclusively on the final RNA structure. This review hopes to bridge this gap. We summarize the existing evidence for cotranscriptional folding and then review the different, currently used strategies for RNA secondary-structure prediction. Finally, we propose a range of ideas on how state-of-the-art methods could be potentially improved by explicitly capturing the process of cotranscriptional structure formation.

MD simulation studies to investigate iso-energetic conformational behaviour of modified nucleosides m(2)G and m(2) 2G present in tRNA

Modified nucleic acid bases are most commonly found in tRNA. These may contain modifications from simple methylation to addition of bulky groups. Methylation of the four canonical nucleotide bases at a wide variety of positions is particularly prominent among the known modification. Methylation of N2 group of guanine is a relatively common modification in tRNA and rRNA. N2-methylguanosine (m²G) is the second most often encountered nucleoside in E. coli tRNAs. N2, N2- dimethylguanosine (m²₂G) is found in the majority of eukaryotic tRNAs and involved in forming base pair interactions with adjacent bases. Hence, in order to understand the structural significance of these methylated nucleic acid bases we have carried out molecular dynamics simulation to see the salvation effect. The results obtained shows iso-energetic conformational behaviors for m²G and m²₂G. The simulation trajectory of m²G shows regular periodical fluctuations suggesting that m²G is equally stable as either s-cis or s-trans rotamers. The two rotamers of m²G may interact canonically or non-canonically with opposite base as s-trans m²G26:C/A/U44 and s-cis m²G26:A/U44. The free rotations around the C-N bond could be the possible reason for these iso-energetic conformations. Dimethylation of G has almost no influence on base pairing with either A or U. Thus, these results reveal that modified nucleosides m²G and m²₂G may play an important role to prevent tRNA from adopting the unusual mitochondrial like conformation.

Transient RNA structure features are evolutionarily conserved and can be computationally predicted

Functional RNA structures tend to be conserved during evolution. This finding is, for example, exploited by comparative methods for RNA secondary structure prediction that currently provide the state-of-art in terms of prediction accuracy. We here provide strong evidence that homologous RNA genes not only fold into similar final RNA structures, but that their folding pathways also share common transient structural features that have been evolutionarily conserved. For this, we compile and investigate a non-redundant data set of 32 sequences with known transient and final RNA secondary structures and devise a dedicated computational analysis pipeline.

COFOLD: an RNA secondary structure prediction method that takes co-transcriptional folding into account

Existing state-of-the-art methods that take a single RNA sequence and predict the corresponding RNA secondary structure are thermodynamic methods. These aim to predict the most stable RNA structure. There exists by now ample experimental and theoretical evidence that the process of structure formation matters and that sequences in vivo fold while they are being transcribed. None of the thermodynamic methods, however, consider the process of structure formation. Here, we present a conceptually new method for predicting RNA secondary structure, called CoFold, that takes effects of co-transcriptional folding explicitly into account. Our method significantly improves the state-of-art in terms of prediction accuracy, especially for long sequences of >1000 nt in length.

Cotranscriptional folding kinetics of ribonucleic acid secondary structures

We develop a systematic helix-based computational method to predict RNA folding kinetics during transcription. In our method, the transcription is modeled as stepwise process, where each step is the transcription of a nucleotide. For each step, the kinetics algorithm predicts the population kinetics, transition pathways, folding intermediates, and the transcriptional folding products. The folding pathways, rate constants, and the conformational populations for cotranscription folding show contrastingly different features than the refolding kinetics for a fully transcribed chain. The competition between the transcription speed and rate constants for the transitions between the different nascent structures determines the RNA folding pathway and the end product of folding. For example, fast transcription favors the formation of branch-like structures than rod-like structures and chain elongation in the folding process may reduce the probability of the formation of misfolded structures. Furthermore, good theory-experiment agreements suggest that our method may provide a reliable tool for quantitative prediction for cotranscriptional RNA folding, including the kinetics for the population distribution for the whole conformational ensemble.

Conformational preferences of modified nucleoside N(2)-methylguanosine (m(2)G) and its derivative N(2), N(2)-dimethylguanosine (m(2)(2)G) occur at 26th position (hinge region) in tRNA

Conformational preferences of the modified nucleosides N(2)-methylguanosine (m(2)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G) have been studied theoretically by using quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. Automated complete geometry optimization using semiempirical quantum chemical RM1, along with ab initio molecular orbital Hartree-Fock (HF-SCF), and density functional theory (DFT) calculations has also been made to compare the salient features. Single-point energy calculation studies have been made on various models of m(2)G26:C/A/U44 and m(2)(2)G26:C/A/U44. The glycosyl torsion angle prefers "syn" (χ = 286°) conformation for m(2)G and m(2)(2)G molecules. These conformations are stabilized by N(3)-HC2' and N(3)-HC3' by replacing weak interaction between O5'-HC(8). The N(2)-methyl substituent of (m(2)G26) prefers "proximal" or s-trans conformation. It may also prefer "distal" or s-cis conformation that allows base pairing with A/U44 instead of C at the hinge region. Thus, N(2)-methyl group of m(2)G may have energetically two stable s-trans m(2)G:C/A/U or s-cis m(2)G:A/U rotamers. This could be because of free rotations around C-N bond. Similarly, N(2), N(2)-dimethyl substituent of (m(2)(2)G) prefers "distal" conformation that may allow base pairing with A/U instead of C at 44th position. Such orientations of m(2)G and m(2)(2)G could play an important role in base-stacking interactions at the hinge region of tRNA during protein biosynthesis process.

TRANSAT-- method for detecting the conserved helices of functional RNA structures, including transient, pseudo-knotted and alternative structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment.

Many non-coding genes exert their function via an RNA structure which starts emerging while the RNA sequence is being transcribed from the genome. The resulting folding pathway is known to depend on a variety of features such as the transcription speed, the concentration of various ions and the binding of proteins and other molecules. Not all of these influences can be adequately captured by the existing computational methods which try to replicate what happens in vivo. So far, it has been challenging to experimentally investigate co-transcriptional folding pathways in vivo and only little data from in vitro experiments exists. In order to investigate if functionally similar RNA sequences from different organisms fold in a similar way, we have developed a new computational method, called Transat, which does not require the detailed computational modeling of the cellular environment. We show in a comprehensive analysis that our method is capable of detecting known structural features and provide evidence that structural features of the in vivo folding pathways have been conserved for several biologically interesting classes of RNA sequences.

RNA structure prediction

We give an overview of RNA structure predictions in this chapter. We discuss here the main approaches to RNA structure prediction: combinatorial approaches, comparative approaches, and kinetic approaches. The main algorithms and mathematical concepts such as transformational grammars will be briefly introduced.

Dynamic extended folding: modeling the RNA secondary structures during co-transcriptional folding

For RNA secondary structure prediction, it is an important issue that how to deal with co-transcriptional folding during the RNA synthesis in the cell. On one hand, co-transcriptional folding, leads to the correct final structure of the whole RNA molecule. On the other hand, it may form the recognition sites for the progress of the transcription. Considering the hurdles in the experimental determination of RNA folding structures, we proposed a so-called "dynamic extended folding simulation" approach. We used two human pre-mRNA samples, the first functional alpha-gene HBZ and the fifth beta-gene HBB, to "display" the co-transcriptional folding images in detail. The modeling process starts from the prediction of a 30-nucleotide (nt) sequence, then in each update 30 nts was extended, say, 1-30, 1-60, 1-90, 1-120,..., 1-1651 nts (for HBB, 1-1606 nts). We selected the RNAstructure program to predict the folding secondary structures of all the segments. We defined "hairpin" as the unit of the secondary structure and analyzed the states of such unit during the sequential dynamic extended folding processes. We found that some hairpins are "conserved", i.e., after its appearance, it always is there in the followed foldings. Some hairpins present partially in the folding segments, and some hairpins appear for only once or twice. This phenomenon vividly depicts the generation and adjusting of the temporal structural units during the co-transcriptional folding process. It is these "hairpins" that support the thermodynamically stable structure at the end of the RNA synthesis. They may also play a role in RNA splicing process and even in the folding structure of the synthesized protein.

The jerky and knotty dynamics of RNA

RNA is known to exhibit a jerky dynamics, as intramolecular thermal motion, on <0.1 micros time scales, is punctuated by infrequent structural rearrangements on much longer time scales, i.e. from >10 micros up to a few minutes or even hours. These rare stochastic events correspond to the formation or dissociation of entire stems through cooperative base pairing/unpairing transitions. Such a clear separation of time scales in RNA dynamics has made it possible to implement coarse grained RNA simulations, which predict RNA folding and unfolding pathways including kinetically trapped structures on biologically relevant time scales of seconds to minutes. RNA folding simulations also enable to predict the formation of pseudoknots, that is, helices interior to loops, which mechanically restrain the relative orientations of other non-nested helices. But beyond static structural constraints, pseudoknots can also strongly affect the folding and unfolding dynamics of RNA, as the order by which successive helices are formed and dissociated can lead to topologically blocked transition intermediates. The resulting knotty dynamics can enhance the stability of RNA switches, improve the efficacy of co-transcriptional folding pathways and lead to unusual self-assembly properties of RNA.

Dynamics of co-transcriptional pre-mRNA folding influences the induction of dystrophin exon skipping by antisense oligonucleotides

Antisense oligonucleotides (AONs) mediated exon skipping offers potential therapy for Duchenne muscular dystrophy. However, the identification of effective AON target sites remains unsatisfactory for lack of a precise method to predict their binding accessibility. This study demonstrates the importance of co-transcriptional pre-mRNA folding in determining the accessibility of AON target sites for AON induction of selective exon skipping in DMD. Because transcription and splicing occur in tandem, AONs must bind to their target sites before splicing factors. Furthermore, co-transcriptional pre-mRNA folding forms transient secondary structures, which redistributes accessible binding sites. In our analysis, to approximate transcription elongation, a “window of analysis” that included the entire targeted exon was shifted one nucleotide at a time along the pre-mRNA. Possible co-transcriptional secondary structures were predicted using the sequence in each step of transcriptional analysis. A nucleotide was considered “engaged” if it formed a complementary base pairing in all predicted secondary structures of a particular step. Correlation of frequency and localisation of engaged nucleotides in AON target sites accounted for the performance (efficacy and efficiency) of 94% of 176 previously reported AONs. Four novel insights are inferred: (1) the lowest frequencies of engaged nucleotides are associated with the most efficient AONs; (2) engaged nucleotides at 3′ or 5′ ends of the target site attenuate AON performance more than at other sites; (3) the performance of longer AONs is less attenuated by engaged nucleotides at 3′ or 5′ ends of the target site compared to shorter AONs; (4) engaged nucleotides at 3′ end of a short target site attenuates AON efficiency more than at 5′ end.

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

SimulFold: simultaneously inferring RNA structures including pseudoknots, alignments, and trees using a Bayesian MCMC framework

Computational methods for predicting evolutionarily conserved rather than thermodynamic RNA structures have recently attracted increased interest. These methods are indispensable not only for elucidating the regulatory roles of known RNA transcripts, but also for predicting RNA genes. It has been notoriously difficult to devise them to make the best use of the available data and to predict high-quality RNA structures that may also contain pseudoknots. We introduce a novel theoretical framework for co-estimating an RNA secondary structure including pseudoknots, a multiple sequence alignment, and an evolutionary tree, given several RNA input sequences. We also present an implementation of the framework in a new computer program, called SimulFold, which employs a Bayesian Markov chain Monte Carlo method to sample from the joint posterior distribution of RNA structures, alignments, and trees. We use the new framework to predict RNA structures, and comprehensively evaluate the quality of our predictions by comparing our results to those of several other programs. We also present preliminary data that show SimulFold's potential as an alignment and phylogeny prediction method. SimulFold overcomes many conceptual limitations that current RNA structure prediction methods face, introduces several new theoretical techniques, and generates high-quality predictions of conserved RNA structures that may include pseudoknots. It is thus likely to have a strong impact, both on the field of RNA structure prediction and on a wide range of data analyses.

Encoding folding paths of RNA switches

RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions (UTRs). Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and inter-molecular base pair interactions can effectively regulate the folding of nascent RNA molecules into different native structures, provided limited coding requirements, as discussed from an information theory perspective. This constitutive coupling between RNA synthesis and RNA folding regulation may have enabled the early emergence of autonomous RNA-based regulation networks.

A counterintuitive Mg2+-dependent and modification-assisted functional folding of mitochondrial tRNAs

Mitochondrial tRNAs (mtRNAs) often lack domains and posttranscriptional modifications that are found in cytoplasmic tRNAs. These structural and chemical elements normally stabilize the folding of cytoplasmic tRNAs into canonical structures that are competent for aminoacylation and translation. For example, the dihydrouridine (D) stem and loop domain is involved in the tertiary structure of cytoplasmic tRNAs through hydrogen bonds and a Mg2+ bridge to the ribothymidine (T) stem and loop domain. These interactions are often absent in mtRNA because the D-domain is truncated or missing. Using gel mobility shift analyses, UV, circular dichroism and NMR spectroscopies and aminoacylation assays, we have investigated the functional folding interactions of chemically synthesized and site-specifically modified mitochondrial and cytoplasmic tRNAs. We found that Mg2+ is critical for folding of the truncated D-domain of bovine mtRNAMet with the tRNA's T-domain. Contrary to the expectation that Mg2+ stabilizes RNA folding, the mtRNAMet D-domain structure was unfolded and relaxed, rather than stabilized in the presence of Mg2+. Because the D-domain is transcribed prior to the T-domain, we conclude that Mg2+ prevents misfolding of the 5'-half of bovine mtRNAMet facilitating its correct interaction with the T-domain. The interaction of the mtRNAMet D-domain with the T-domain was enhanced by a pseudouridine located in either the D or T-domains compared to that of the unmodified RNAs (Kd=25.3, 24.6 and 44.4 microM, respectively). Mg2+ also affected the folding interaction of a yeast mtRNALeu1, but had minimal effect on the folding of an Escherichia coli cytoplasmic tRNALeu. The D-domain modification, dihydrouridine, facilitated mtRNALeu folding. These data indicate that conserved modifications assist and stabilize the formation of the functional mtRNA tertiary structure.

Co-transcriptional folding is encoded within RNA genes

Background

Most of the existing RNA structure prediction programs fold a completely synthesized RNA molecule. However, within the cell, RNA molecules emerge sequentially during the directed process of transcription. Dedicated experiments with individual RNA molecules have shown that RNA folds while it is being transcribed and that its correct folding can also depend on the proper speed of transcription.

Methods

The main aim of this work is to study if and how co-transcriptional folding is encoded within the primary and secondary structure of RNA genes. In order to achieve this, we study the known primary and secondary structures of a comprehensive data set of 361 RNA genes as well as a set of 48 RNA sequences that are known to differ from the originally transcribed sequence units. We detect co-transcriptional folding by defining two measures of directedness which quantify the extend of asymmetry between alternative helices that lie 5' and those that lie 3' of the known helices with which they compete.

Results

We show with statistical significance that co-transcriptional folding strongly influences RNA sequences in two ways: (1) alternative helices that would compete with the formation of the functional structure during co-transcriptional folding are suppressed and (2) the formation of transient structures which may serve as guidelines for the co-transcriptional folding pathway is encouraged.

Conclusions

These findings have a number of implications for RNA secondary structure prediction methods and the detection of RNA genes.

Effect of transcription on folding of the Tetrahymena ribozyme

Sequential formation of RNA interactions during transcription can bias the folding pathway and ultimately determine the functional state of a transcript. The kinetics of cotranscriptional folding of the Tetrahymena L-21 ribozyme was compared with refolding of full-length transcripts under the same conditions. Sequential folding after transcription by phage T7 or Escherichia coli polymerase is only twice as fast as refolding, and the yield of native RNA is the same. By contrast, a greater fraction of circularly permuted variants folded correctly at early times during transcription than during refolding. Hybridization of complementary oligonucleotides suggests that cotranscriptional folding enables a permuted RNA beginning at G303 to escape non-native interactions in P3 and P9. We propose that base pairing of upstream sequences during transcription elongation favors branched secondary structures that increase the probability of forming the native ribozyme structure.

Detection and analysis of hairpin II, an essential metastable structural element in viroid replication intermediates

In (-)-stranded replication intermediates of the potato spindle tuber viroid (PSTVd) a thermodynamically metastable structure containing a specific hairpin structure (HP II) has been proposed to be essential for viroid replication. In the present work a method was devised allowing the direct detection of the HP II structure in vitro and in vivo using a biophysical approach. An RNA oligonucleotide was constructed which specifically binds to the HP II loop region in transient (-)-strand intermediates. Analysis of the resulting oligonucleotide/HP II complexes on temperature-gradient gels enabled us to follow the formation of HP II during in vitro transcription by T7 RNA polymerase. Moreover, we were able to demonstrate the formation of HP II during viroid replication in potato (Solanum tuberosum) cells.

A physical origin for functional domain structure in nucleic acids as evidenced by cross-linking entropy: II

In Part I, cross-linking entropy (CLE) was proposed as a mechanism that limits the size of functional domains of RNA. To test this hypothesis, the theory is developed into an RNA secondary structure prediction filter which is applied to nearest-neighbor secondary structure (NNSS) algorithms that utilize a free energy (FE) minimization strategy. (The NNSS strategies are also referred to as the dynamic programming algorithm in the literature.) The cross-linking entropy for RNA is derived from a generalized Gaussian polymer chain model where the entropic contributions caused by the formation of base pairs (stacking) in RNA are analysed globally. Local entropic contributions are associated with the freezing out of degrees of freedom in the links. Both global and local entropic effects are strongly influenced by the persistence length. The cross-linking entropy provides a physical origin for the size of functional domains in long nucleic acid sequences and may go further to explain as to why the majority of the domain regions in typical sequences tend to be less than 600 nucleotides in length. In addition, improvements were observed in the "best guess" predictive capacity over NNSS prediction strategies. The thermodynamic distribution is more representative of the expected structures and is strongly governed by such physical parameters as the persistence length and the excluded volume. The CLE appears to generalize the tabulated penalties used in NNSS algorithms. The principal parameter influencing this entropy is the persistence length. The model is shown to accomodate a variable persistence length and is capable of describing the folding dynamics of RNA. A two-state kinetic model based on the CLE principle is used to help elucidate the folding kinetics of functional domains in the group I introns.

A role for upstream RNA structure in facilitating the catalytic fold of the genomic hepatitis delta virus ribozyme

Hepatitis delta virus (HDV) has a circular RNA genome that replicates by a double rolling-circle mechanism. The genomic and antigenomic versions of HDV contain a ribozyme that undergoes cis-cleavage, thereby processing the transcript into unit-length monomers. A genomic HDV transcript containing 30 nucleotides immediately upstream of the cleavage site was found to have attenuated self-cleavage. Structure mapping and site-directed mutagenesis revealed an inhibitory stretch consisting of upstream nucleotides -24 to -15 that forms a long-range pairing, termed Alt 1, with the 3' strand of P2 (P2(3')) located at the very 3'-end of the ribozyme. Two other alternative pairings were found, Alt 2, which involves upstream nucleotide-ribozyme interactions, and Alt 3, which involves ribozyme-ribozyme interactions. Self-cleavage was rescued 2700 to 20,000-fold by adding DNA oligomers, which sequester the -24/-15 inhibitory stretch in trans. Surprisingly, co-transcriptional self-cleavage occurs when the number of upstream nucleotides is increased to 54. Computer prediction and structure mapping support the existence of an unusually stable upstream hairpin involving nucleotides -54 to -18, termed P(-1)/L(-1), which sequesters the majority of the -24/-15 inhibitory stretch in cis. This hairpin is followed by a stretch of single-stranded pyrimidine-rich nucleotides, termed J(-1/1). Sequence comparison suggests that the P(-1)/L(-1)/J(-1/1) motif is conserved among known genomic HDV isolates, and that the J(-1/1) stretch is conserved among antigenomic HDV isolates. Lastly, the secondary structure of the Alt 1-containing ribozyme provides insight into possible folding intermediates of the ribozyme.

Characterization of the initiation sites of both polarity strands of a viroid RNA reveals a motif conserved in sequence and structure

Viroids replicate through a rolling-circle mechanism in which the infecting circular RNA and its complementary (-) strand are transcribed. The precise site at which transcription starts was investigated for the avocado sunblotch viroid (ASBVd), the type species of the family of viroids with hammerhead ribozymes. Linear ASBVd (+) and (-) RNAs begin with a UAAAA sequence that maps to similar A+U-rich terminal loops in their predicted quasi-rod-like secondary structures. The sequences around the initiation sites of ASBVd, which replicates and accumulates in the chloroplast, are similar to the promoters of a nuclear-encoded chloroplastic RNA polymerase (NEP), supporting the involvement of an NEP-like activity in ASBVd replication. Since RNA folding appears to be kinetically determined, the specific location of both ASBVd initiation sites provides a mechanistic insight into how the nascent ASBVd strands may fold in vivo. The approach used here, in vitro capping and RNase protection assays, may be useful for investigating the initiation sites of other small circular RNA replicons.

Sequential folding of the genomic ribozyme of the hepatitis delta virus: structural analysis of RNA transcription intermediates

The structures of the model oligoribonucleotides that mimic the consecutive stages in the transcription of genomic HDV ribozyme have been analyzed by the Pb(2+)-induced cleavage method, partial digestion with specific nucleases and chemical probing. In the transcription intermediates, the P1 and P4 helical segments are found to be present in the final folded forms in which they exist in the full-length transcript. However, the region corresponding to the central hairpin forms another thermodynamically stable hairpin structure. Its correct folding requires the presence of a ribozyme 3'-terminal sequence and the formation of helix P2. This confirms the ribozyme structure of the pseudoknot type and points to the crucial role of helix P2 in its overall folding. Moreover, we show that the J4/2 region can be specifically cleaved in the presence of selected divalent metal ions in the full-length transcript, but not in a shorter one lacking six 3'-terminal nucleotides, which cannot form the pseudoknotted structure. Thus, a particular RNA conformation around that cleavage site is required for specific hydrolysis, and the J4/2 region seems to be involved in the formation of a general metal ion binding site. Recently, it has been proposed that, in the antigenomic ribozyme, a four nucleotide sequence within the J1/2 region may contribute to the folding pathway, being part of a mechanism responsible for controlling ribozyme cleavage activity. Our study shows that in the genomic ribozyme the central hairpin region may contribute to a similar mechanism, providing a barrier to the formation of an active structure in the ribozyme folding pathway.

Prediction of sequentially optimal RNA secondary structures

A rigorous mathematical modeling of the RNA sequential folding process during transcription is proposed. It is based, at each transcription step, on a homogeneous markovian jump process, the state space of which is the set of structures constructible on the part of the RNA already transcribed. A theoretical formula permitting the computation of the structures probabilities at the end of the RNA transcription is derived. Successive approximations, aimed at reducing the size of the state space, permit the design of a prediction algorithm. The algorithm is tested on some structural RNAs (tRNA, 5S, 16S, hammerhead, ...), results are discussed and possible improvements are proposed.

Hierarchy and dynamics of RNA folding

The evidence showing that the self-assembly of complex RNAs occurs in discrete transitions, each relating to the folding of sub-systems of increasing size and complexity starting from a state with most of the secondary structure, is reviewed. The reciprocal influence of the concentration of magnesium ions and nucleotide mutations on tertiary structure is analyzed. Several observations demonstrate that detrimental mutations can be rescued by high magnesium concentrations, while stabilizing mutations lead to a lesser dependence on magnesium ion concentration. Recent data point to the central controlling and monitoring roles of RNA-binding proteins that can bind to the different folding stages, either before full establishment of the secondary structure or at the molten globule state before the cooperative transition to the final three-dimensional structure.

Enzymatic formation of N2,N2-dimethylguanosine in eukaryotic tRNA: importance of the tRNA architecture

In eukaryotic tRNA, guanosine at position 26 in the junction between the D-stem and the anticodon stem is mostly modified to N2,N2-dimethylguanosine (m2(2)G26). Here we review the available information on the enzyme catalyzing the formation of this modified nucleoside, the SAM-dependent tRNA (m2(2)G26)-methyltransferase, and our attemps to identify the parameters in tRNA needed for efficient enzymatic dimethylation of guanosine-26. The required identity elements in yeast tRNA for dimethylation under in vitro conditions by the yeast tRNA(m2(2)G26)-methyltransferase (the TRM1 gene product) are comprised of two G-C base pairs at positions G10-C25 and C11-G24 in the D-stem together with a variable loop of at least five nucleotides. These positive determinants do not seem to act via base specific interactions with the methyltransferase; they instead ensure that G26 is presented to the enzyme in a favorable orientation, within the central 3D-core of the tRNA molecule. The anticodon stem and loop is not involved in such an interaction with the enzyme. In a heterologous in vivo system, consisting of yeast tRNAs microinjected into Xenopus laevis oocytes, the requirements for modification of G26 are less stringent than in the yeast homologous in vitro system. Indeed, G26 in several microinjected tRNAs becomes monomethylated, while in yeast extracts it stays unmethylated, even after extensive incubation. Thus either the X laevis tRNA(m2(2)G26)-methyltransferase has a more relaxed specificity than its yeast homolog, or there exist two distinct G26-methylating activities, one for G26-monomethylation, and one for dimethylation of G26.(ABSTRACT TRUNCATED AT 250 WORDS)

Modeling RNA secondary structures. I. Mathematical structural model for predicting RNA secondary structures

A mathematical model for analyzing the secondary structures of RNA is developed that is based on the connection matrix associated with the planar p-h graph. The classification of the elementary structures allows the introduction of the basis of structural space from which to build the global secondary structure. All admissible solutions belong to the configuration space and can be obtained directly from its basis.

Visualizing the higher order folding of a catalytic RNA molecule

The higher order folding process of the catalytic RNA derived from the self-splicing intron of Tetrahymena thermophila was monitored with the use of Fe(II)-EDTA-induced free radical chemistry. The overall tertiary structure of the RNA molecule forms cooperatively with the uptake of at least three magnesium ions. Local folding transitions display different metal ion dependencies, suggesting that the RNA tertiary structure assembles through a specific folding intermediate before the catalytic core is formed. Enzymatic activity, assayed with an RNA substrate that is complementary to the catalytic RNA active site, coincides with the cooperative structural transition. The higher order RNA foldings produced by Mg(II), Ca(II), and Sr(II) are similar; however, only the Mg(II)-stabilized RNA is catalytically active. Thus, these results directly demonstrate that divalent metal ions participate in general folding of the ribozyme tertiary structure, and further indicate a more specific involvement of Mg(II) in catalysis.

Characterization of the lead(II)-induced cleavages in tRNAs in solution and effect of the Y-base removal in yeast tRNAPhe

The specificity of lead(II)-induced hydrolysis of yeast tRNA(Phe) was studied as a function of concentration of Pb2+ ions. The major cut was localized in the D-loop and minor cleavages were detected in the anticodon and T-loops at high metal ion concentration. The effects of pH, temperature, and urea were also analyzed, revealing a basically unchanged specificity of hydrolysis. In the isolated 5'-half-molecule of yeast tRNAPhe not cut was found in the D-loop, indicating its stringent dependence on T-D-loop interaction. Comparison of hydrolysis patterns and efficiencies observed in yeast tRNA(Phe) with those found in other tRNAs suggests that the presence of a U59-C60 sequence in the T-loop is responsible for the highly efficient and specific hydrolysis in the spatially close region of the D-loop. The efficiencies of D-loop cleavage in intact yeast tRNA(Phe) and in tRNA(Phe) deprived of the Y base next to the anticodon were also compared at various Pb2+ ion concentrations. Kinetics of the D-loop hydrolysis analyzed at 0, 25, and 37 degrees C showed a 6 times higher susceptibility of tRNA(Phe) minus Y base (tRNA(Phe)-Y) to lead(II)-induced hydrolysis than in tRNA(Phe). The observed effect is discussed in terms of a long-distance conformational transition in the region of the interacting D- and T-loops triggered by the Y-base excision.

Morpholino spin-labeling for base-pair sequencing of a 3'-terminal RNA stem by proton homonuclear Overhauser enhancements: yeast ribosomal 5S RNA

Base-pair sequences for 5S and 5.8S RNAs are not readily extracted from proton homonuclear nuclear Overhauser enhancement (NOE) connectivity experiments alone, due to extensive peak overlap in the downfield (11-15 ppm) proton NMR spectrum. In this paper, we introduce a new method for base-pair proton peak assignment for ribosomal RNAs, based upon the distance-dependent broadening of the resonances of base-pair protons spatially proximal to a paramagnetic group. Introduction of a nitroxide spin-label covalently attached to the 3'-terminal ribose provides an unequivocal starting point for base-pair hydrogen-bond proton NMR assignment. Subsequent NOE connectivities then establish the base-pair sequence for the terminal stem of a 5S RNA. Periodate oxidation of yeast 5S RNA, followed by reaction with 4-amino-2,2,6,6-tetramethylpiperidinyl-1-oxy (TEMPO-NH2) and sodium borohydride reduction, produces yeast 5S RNA specifically labeled with a paramagnetic nitroxide group at the 3'-terminal ribose. Comparison of the 500-MHz 1H NMR spectra of native and 3'-terminal spin-labeled yeast 5S RNA serves to identify the terminal base pair (G1 . C120) and its adjacent base pair (G2 . U119) on the basis of their proximity to the 3'-terminal spin-label. From that starting point, we have then identified (G . C, A . U, or G . U) and sequenced eight of the nine base pairs in the terminal helix via primary and secondary NOE's.