Effects of polyvalent cations on the folding of an rRNA three-way junction and binding of ribosomal protein S15

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) predicted over 10 inner-shell Mg2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminishes the selectivity between SAM and SAH. One Mg2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg2+ thus plays vital roles in SAM/SAH riboswitch function.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg 2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) identified over 10 inner-shell Mg 2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminish the selectivity between SAM and SAH. One Mg 2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg 2+ thus plays vital roles in SAM/SAH riboswitch function.

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl2), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg2+ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl2) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals ΔG≃ 10 kcal/mol. Besides stabilizing the 3WJ, Mg2+ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.

Shape changes and cooperativity in the folding of the central domain of the 16S ribosomal RNA

Both the small and large subunits of the ribosome, the molecular machine that synthesizes proteins, are complexes of ribosomal RNAs (rRNAs) and a number of proteins. In bacteria, the small subunit has a single 16S rRNA whose folding is the first step in its assembly. The central domain of the 16S rRNA folds independently, driven either by Mg²⁺ ions or by interaction with ribosomal proteins. To provide a quantitative description of ion-induced folding of the ∼350-nucleotide rRNA, we carried out extensive coarse-grained molecular simulations spanning Mg²⁺ concentration between 0 and 30 mM. The Mg²⁺ dependence of the radius of gyration shows that globally the rRNA folds cooperatively. Surprisingly, various structural elements order at different Mg²⁺ concentrations, indicative of the heterogeneous assembly even within a single domain of the rRNA. Binding of Mg²⁺ ions is highly specific, with successive ion condensation resulting in nucleation of tertiary structures. We also predict the Mg²⁺-dependent protection factors, measurable in hydroxyl radical footprinting experiments, which corroborate the specificity of Mg²⁺-induced folding. The simulations, which agree quantitatively with several experiments on the folding of a three-way junction, show that its folding is preceded by formation of other tertiary contacts in the central junction. Our work provides a starting point in simulating the early events in the assembly of the small subunit of the ribosome.

Real-time assembly of ribonucleoprotein complexes on nascent RNA transcripts

Cellular protein-RNA complexes assemble on nascent transcripts, but methods to observe transcription and protein binding in real time and at physiological concentrations are not available. Here, we report a single-molecule approach based on zero-mode waveguides that simultaneously tracks transcription progress and the binding of ribosomal protein S15 to nascent RNA transcripts during early ribosome biogenesis. We observe stable binding of S15 to single RNAs immediately after transcription for the majority of the transcripts at 35 °C but for less than half at 20 °C. The remaining transcripts exhibit either rapid and transient binding or are unable to bind S15, likely due to RNA misfolding. Our work establishes the foundation for studying transcription and its coupled co-transcriptional processes, including RNA folding, ligand binding, and enzymatic activity such as in coupling of transcription to splicing, ribosome assembly or translation.

The early steps of ribosome assembly occur co-transcriptionally on the nascent ribosomal RNA. Here the authors demonstrate an approach that allows simultaneous monitoring of transcription and ribosomal protein assembly at the single-molecule level in real time.

Evolution of protein-coupled RNA dynamics during hierarchical assembly of ribosomal complexes

Assembly of 30S ribosomes involves the hierarchical addition of ribosomal proteins that progressively stabilize the folded 16S rRNA. Here, we use three-color single molecule FRET to show how combinations of ribosomal proteins uS4, uS17 and bS20 in the 16S 5′ domain enable the recruitment of protein bS16, the next protein to join the complex. Analysis of real-time bS16 binding events shows that bS16 binds both native and non-native forms of the rRNA. The native rRNA conformation is increasingly favored after bS16 binds, explaining how bS16 drives later steps of 30S assembly. Chemical footprinting and molecular dynamics simulations show that each ribosomal protein switches the 16S conformation and dampens fluctuations at the interface between rRNA subdomains where bS16 binds. The results suggest that specific protein-induced changes in the rRNA dynamics underlie the hierarchy of 30S assembly and simplify the search for the native ribosome structure.

Ribosomes assemble through the hierarchical addition of proteins to a ribosomal RNA scaffold. Here the authors use three-color single-molecule FRET to show how the dynamics of the rRNA dictate the order in which multiple proteins assemble on the 5′ domain of the E. coli 16S rRNA.

Protein-RNA Dynamics in the Central Junction Control 30S Ribosome Assembly

Interactions between ribosomal proteins (rproteins) and ribosomal RNA (rRNA) facilitate the formation of functional ribosomes. S15 is a central domain primary binding protein that has been shown to trigger a cascade of conformational changes in 16S rRNA, forming the functional structure of the central domain. Previous biochemical and structural studies in vitro have revealed that S15 binds a three-way junction of helices 20, 21, and 22, including nucleotides 652-654 and 752-754. All junction nucleotides except 653 are highly conserved among the Bacteria. To identify functionally important motifs within the junction, we subjected nucleotides 652-654 and 752-754 to saturation mutagenesis and selected and analyzed functional mutants. Only 64 mutants with greater than 10% ribosome function in vivo were isolated. S15 overexpression complemented mutations in the junction loop in each of the partially active mutants, although mutations that produced inactive ribosomes were not complemented by overexpression of S15. Single-molecule Förster or fluorescence resonance energy transfer (smFRET) was used to study the Mg(2+)- and S15-induced conformational dynamics of selected junction mutants. Comparison of the structural dynamics of these mutants with the wild type in the presence and absence of S15 revealed specific sequence and structural motifs in the central junction that are important in ribosome function.

Structural Analysis of Multi-Helical RNAs by NMR-SAXS/WAXS: Application to the U4/U6 di-snRNA

NMR and SAXS (small-angle X-ray scattering)/WAXS (wide-angle X-ray scattering) are highly complementary approaches for the analysis of RNA structure in solution. Here we describe an efficient NMR-SAXS/WAXS approach for structural investigation of multi-helical RNAs. We illustrate this approach by determining the overall fold of a 92-nt 3-helix junction from the U4/U6 di-snRNA. The U4/U6 di-snRNA is conserved in eukaryotes and is part of the U4/U6.U5 tri-snRNP, a large ribonucleoprotein complex that comprises a major subunit of the assembled spliceosome. Helical orientations can be determined by X-ray scattering data alone, but the addition of NMR RDC (residual dipolar coupling) restraints improves the structure models. RDCs were measured in two different external alignment media and also by magnetic susceptibility anisotropy. The resulting alignment tensors are collinear, which is a previously noted problem for nucleic acids. Including WAXS data in the calculations produces models with significantly better fits to the scattering data. In solution, the U4/U6 di-snRNA forms a 3-helix junction with a planar Y-shaped structure and has no detectable tertiary interactions. Single-molecule Förster resonance energy transfer data support the observed topology. A comparison with the recently determined cryo-electron microscopy structure of the U4/U6.U5 tri-snRNP illustrates how proteins scaffold the RNA and dramatically alter the geometry of the U4/U6 3-helix junction.

Assembly constraints drive co-evolution among ribosomal constituents

Ribosome biogenesis, a central and essential cellular process, occurs through sequential association and mutual co-folding of protein-RNA constituents in a well-defined assembly pathway. Here, we construct a network of co-evolving nucleotide/amino acid residues within the ribosome and demonstrate that assembly constraints are strong predictors of co-evolutionary patterns. Predictors of co-evolution include a wide spectrum of structural reconstitution events, such as cooperativity phenomenon, protein-induced rRNA reconstitutions, molecular packing of different rRNA domains, protein-rRNA recognition, etc. A correlation between folding rate of small globular proteins and their topological features is known. We have introduced an analogous topological characteristic for co-evolutionary network of ribosome, which allows us to differentiate between rRNA regions subjected to rapid reconstitutions from those hindered by kinetic traps. Furthermore, co-evolutionary patterns provide a biological basis for deleterious mutation sites and further allow prediction of potential antibiotic targeting sites. Understanding assembly pathways of multicomponent macromolecules remains a key challenge in biophysics. Our study provides a 'proof of concept' that directly relates co-evolution to biophysical interactions during multicomponent assembly and suggests predictive power to identify candidates for critical functional interactions as well as for assembly-blocking antibiotic target sites.

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.

RNA and its ionic cloud: solution scattering experiments and atomically detailed simulations

RNA molecules play critical roles in many cellular processes. Traditionally viewed as genetic messengers, RNA molecules were recently discovered to have diverse functions related to gene regulation and expression. RNA also has great potential as a therapeutic and a tool for further investigation of gene regulation. Metal ions are an integral part of RNA structure and should be considered in any experimental or theoretical study of RNA. Here, we report a multidisciplinary approach that combines anomalous small-angle x-ray scattering and molecular-dynamics (MD) simulations with explicit solvent and ions around RNA. From experiment and simulation results, we find excellent agreement in the number and distribution of excess monovalent and divalent ions around a short RNA duplex. Although similar agreement can be obtained from a continuum description of the solvent and mobile ions (by solving the Poisson-Boltzmann equation and accounting for finite ion size), the use of MD is easily extended to flexible RNA systems with thermal fluctuations. Therefore, we also model a short RNA pseudoknot and find good agreement between the MD results and the experimentally derived solution structures. Surprisingly, both deviate from crystal structure predictions. These favorable comparisons of experiment and simulations encourage work on RNA in all-atom dynamic models.

Molecular sensing by the aptamer domain of the FMN riboswitch: a general model for ligand binding by conformational selection

Understanding the nature of the free state of riboswitch aptamers is important for illuminating common themes in gene regulation by riboswitches. Prior evidence indicated the flavin mononucleotide (FMN)-binding riboswitch aptamer adopted a 'bound-like' structure in absence of FMN, suggesting only local conformational changes upon ligand binding. In the scope of pinpointing the general nature of such changes at the nucleotide level, we performed SHAPE mapping experiments using the aptamer domain of two phylogenetic variants, both in absence and in presence of FMN. We also solved the crystal structures of one of these domains both free (3.3 Å resolution) and bound to FMN (2.95 Å resolution). Our comparative study reveals that structural rearrangements occurring upon binding are restricted to a few of the joining regions that form the binding pocket in both RNAs. This type of binding event with minimal structural perturbations is reminiscent of binding events by conformational selection encountered in other riboswitches and various RNAs.

Mitochondrial enzyme rhodanese is essential for 5 S ribosomal RNA import into human mitochondria

5 S rRNA is an essential component of ribosomes. In eukaryotic cells, it is distinguished by particularly complex intracellular traffic, including nuclear export and re-import. The finding that in mammalian cells 5 S rRNA can eventually escape its usual circuit toward nascent ribosomes to get imported into mitochondria has made the scheme more complex, and it has raised questions about both the mechanism of 5 S rRNA mitochondrial targeting and its function inside the organelle. Previously, we showed that import of 5 S rRNA into mitochondria requires unknown cytosolic proteins. Here, one of them was identified as mitochondrial thiosulfate sulfurtransferase, rhodanese. Rhodanese in its misfolded form was found to possess a strong and specific 5 S rRNA binding activity, exploiting sites found earlier to function as signals of 5 S rRNA mitochondrial localization. The interaction with 5 S rRNA occurs cotranslationally and results in formation of a stable complex in which rhodanese is preserved in a compact enzymatically inactive conformation. Human 5 S rRNA in a branched Mg(2+)-free form, upon its interaction with misfolded rhodanese, demonstrates characteristic functional traits of Hsp40 cochaperones implicated in mitochondrial precursor protein targeting, suggesting that it may use this mechanism to ensure its own mitochondrial localization. Finally, silencing of the rhodanese gene caused not only a proportional decrease of 5 S rRNA import but also a general inhibition of mitochondrial translation, indicating the functional importance of the imported 5 S rRNA inside the organelle.

Kinetic cooperativity in Escherichia coli 30S ribosomal subunit reconstitution reveals additional complexity in the assembly landscape

The Escherichia coli 30S ribosomal subunit self-assembles in vitro in a hierarchical manner, with the RNA binding by proteins enabled by the prior binding of others under equilibrium conditions. Early 16S rRNA binding proteins also bind faster than late-binding proteins, but the specific causes for the slow binding of late proteins remain unclear. Previously, a pulse-chase monitored by quantitative mass spectrometry method was developed for monitoring 30S subunit assembly kinetics, and here a modified experimental scheme was used to probe kinetic cooperativity by including a step where subsets of ribosomal proteins bind and initiate assembly prior to the pulse-chase kinetics. In this work, 30S ribosomal subunit kinetic reconstitution experiments revealed that thermodynamic dependency does not always correlate with kinetic cooperativity. Some folding transitions that cause subsequent protein binding to be more energetically favorable do not result in faster protein binding. Although 3(') domain primary protein S7 is required for RNA binding by both proteins S9 and S19, prior binding of S7 accelerates the binding of S9, but not S19, indicating there is an additional mechanistic step required for S19 to bind. Such data on kinetic cooperativity and the presence of multiphasic assembly kinetics reveal complexity in the assembly landscape that was previously hidden.

Salt-dependent folding energy landscape of RNA three-way junction

RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg(2+) ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na(+) and Mg(2+) ion concentrations. For Mg(2+) solutions, the TBI model, which accounts for the Mg(2+) ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg(2+) in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.

RNA-protein binding kinetics in an automated microfluidic reactor

Microfluidic chips can automate biochemical assays on the nanoliter scale, which is of considerable utility for RNA-protein binding reactions that would otherwise require large quantities of proteins. Unfortunately, complex reactions involving multiple reactants cannot be prepared in current microfluidic mixer designs, nor is investigation of long-time scale reactions possible. Here, a microfluidic 'Riboreactor' has been designed and constructed to facilitate the study of kinetics of RNA-protein complex formation over long time scales. With computer automation, the reactor can prepare binding reactions from any combination of eight reagents, and is optimized to monitor long reaction times. By integrating a two-photon microscope into the microfluidic platform, 5-nl reactions can be observed for longer than 1000 s with single-molecule sensitivity and negligible photobleaching. Using the Riboreactor, RNA-protein binding reactions with a fragment of the bacterial 30S ribosome were prepared in a fully automated fashion and binding rates were consistent with rates obtained from conventional assays. The microfluidic chip successfully combines automation, low sample consumption, ultra-sensitive fluorescence detection and a high degree of reproducibility. The chip should be able to probe complex reaction networks describing the assembly of large multicomponent RNPs such as the ribosome.

Comparative gel electrophoresis analysis of helical junctions in RNA

Comparative gel electrophoresis provides information on the relative angles subtended between helical arms at a branchpoint in RNA. It is based upon the comparison of electrophoretic mobility in polyacrylamide gels of species containing two long arms, with the remaining one(s) being significantly shorter. Although the method currently lacks a really well-established basis of physical theory, it is very powerful, yet simple to apply. It has had a number of significant successes in RNA, DNA and DNA-protein complexes, and in all cases to date the results have stood the test of time and eventual comparison with crystallographic analysis.

Single-molecule mechanical unfolding and folding of a pseudoknot in human telomerase RNA

RNA unfolding and folding reactions in physiological conditions can be facilitated by mechanical force one molecule at a time. By using force-measuring optical tweezers, we studied the mechanical unfolding and folding of a hairpin-type pseudoknot in human telomerase RNA in a near-physiological solution, and at room temperature. Discrete two-state folding transitions of the pseudoknot are seen at approximately 10 and approximately 5 piconewtons (pN), with ensemble rate constants of approximately 0.1 sec(-1), by stepwise force-drop experiments. Folding studies of the isolated 5'-hairpin construct suggested that the 5'-hairpin within the pseudoknot forms first, followed by formation of the 3'-stem. Stepwise formation of the pseudoknot structure at low forces are in contrast with the one-step unfolding at high forces of approximately 46 pN, at an average rate of approximately 0.05 sec(-1). In the constant-force folding trajectories at approximately 10 pN and approximately 5 pN, transient formation of nonnative structures were observed, which is direct experimental evidence that folding of both the hairpin and pseudoknot takes complex pathways. Possible nonnative structures and folding pathways are discussed.

Recognition elements in rRNA for the tylosin resistance methyltransferase RlmA(II)

The methyltransferase RlmA(II) (formerly TlrB) is found in many Gram-positive bacteria, and methylates the N-1 position of nucleotide G748 within the loop of hairpin 35 in 23S rRNA. Methylation of the rRNA by RlmA(II) confers resistance to tylosin and other mycinosylated 16-membered ring macrolide antibiotics. We have previously solved the solution structure of hairpin 35 in the conformation that is recognized by the RlmA(II) methyltransferase from Streptococcus pneumoniae. It was shown that while essential recognition elements are located in hairpin 35, the interactions between RlmA(II) and hairpin 35 are insufficient on their own to support the methylation reaction. Here we use biochemical techniques in conjunction with heteronuclear/homonuclear nuclear magnetic resonance spectroscopy to define the RNA structures that are required for efficient methylation by RlmA(II). Progressive truncation of the rRNA substrate indicated that multiple contacts occur between RlmA(II) and nucleotides in stem-loops 33, 34 and 35. RlmA(II) appears to recognize its rRNA target through specific surface shape complementarity at the junction formed by these three helices. This means of recognition is highly similar to that of the orthologous Gram-negative methyltransferase, RlmA(I) (formerly RrmA), which also interacts with hairpin 35, but methylates at the adjacent nucleotide G745.

Ion atmosphere around nucleic acid

The poyion-ion preferential interaction coefficient Gamma describes the exclusion of coions and accumulations of counterions in the vicinity of a polyion in an aqueous solution. We give tight upper and lower bounds for Gamma when the polyion can be modeled by a cylinder of infinite length but of arbitrary charge density. This case can be used as a model for long strands of DNA or RNA in an aqueous solution containing univalent cations. The salt dependence of Gamma is predicted from low to intermediate and high salt concentrations. We also indicate how the bounds for the infinite polyion can be exploited to place bounds for polyions of length greater than a constant on the order of the inverse Debye screening length.

Ion atmosphere of three-way junction nucleic acid

The ion atmosphere of three-armed symmetric Y-shaped and asymmetric y-shaped A-RNA junctions in aqueous solution containing multivalent ions is described within the framework of a polyelectrolyte model. The fraction of "screening counterions" per polyion charge that shield the residual unneutralized charges from interacting with one another and the condensed counterions per polyion charge as a function of sodium and magnesium ion concentrations are determined. The predictions for the slope of log(k(o)/k(f)) as a function of Na+ and Mg2+ concentration, where k(o) and k(f) are the opening and folding rates of the three-helix junction molecule, respectively, are compared with experimental data (Kim et al. Proc. Nat. Acad. Sci. U.S.A. 2002, 96, 9077-9082).

H/ACA small nucleolar RNA pseudouridylation pockets bind substrate RNA to form three-way junctions that position the target U for modification

During the biogenesis of eukaryotic ribosomal RNA (rRNA) and spliceosomal small nuclear RNA (snRNA), uridines at specific sites are converted to pseudouridines by H/ACA ribonucleoprotein particles (RNPs). Each H/ACA RNP contains a substrate-specific H/ACA RNA and four common proteins, the pseudouridine synthase Cbf5, Nop10, Gar1, and Nhp2. The H/ACA RNA contains at least one pseudouridylation (psi) pocket, which is complementary to the sequences flanking the target uridine. In this article, we show structural evidence that the psi pocket can form the predicted base pairs with substrate RNA in the absence of protein components. We report the solution structure of the complex between an RNA hairpin derived from the 3' psi pocket of human U65 H/ACA small nucleolar RNA (snoRNA) and the substrate rRNA. The snoRNA-rRNA substrate complex has a unique structure with two offset parallel pairs of stacked helices and two unusual intermolecular three-way junctions, which together organize the substrate for docking into the active site of Cbf5. The substrate RNA interacts on one face of the snoRNA in the complex, forming a structure that easily could be accommodated in the H/ACA RNP, and explains how successive substrate RNAs could be loaded onto and unloaded from the H/ACA RNA in the RNP.

An allosteric-feedback mechanism for protein-assisted group I intron splicing

The I-AniI maturase facilitates self-splicing of a mitochondrial group I intron in Aspergillus nidulans. Binding occurs in at least two steps: first, a specific but labile encounter complex rapidly forms and then this intermediate is slowly resolved into a native, catalytically active RNA/protein complex. Here we probe the structure of the RNA throughout the assembly pathway. Although inherently unstable, the intron core, when bound by I-AniI, undergoes rapid folding to a near-native state in the encounter complex. The next transition includes the slow destabilization and docking into the core of the peripheral stacked helix that contains the 5' splice site. Mutational analyses confirm that both transitions are important for native complex formation. We propose that protein-driven destabilization and docking of the peripheral stacked helix lead to subtle changes in the I-AniI binding site that facilitate native complex formation. These results support an allosteric-feedback mechanism of RNA-protein recognition in which proteins engaged in an intermediate complex can influence RNA structure far from their binding sites. The linkage of these changes to stable binding ensures that the protein and RNA do not get sequestered in nonfunctional complexes.

The conformational landscape of the ribosomal protein S15 and its influence on the protein interaction with 16S RNA

The interaction between the ribosomal protein S15 and its binding sites in the 16S RNA was examined from two points of view. First, the isolated protein S15 was studied by comparing NMR conformer sets, available in the PDB and recalculated using the CNS-ARIA protocol. Molecular dynamics (MD) trajectories were then recorded starting from a conformer of each set. The recalculation of the S15 NMR structure, as well as the recording of MD trajectories, reveals that several orientations of the N-terminal alpha-helix alpha1 with respect to the structure core are populated. MD trajectories of the complex between the ribosomal protein S15 and RNA were also recorded in the presence and absence of Mg(2+) ions. The Mg(2+) ions are hexacoordinated by water and RNA oxygens. The coordination spheres mainly interact with the RNA phosphodiester backbone, reducing the RNA mobility and inducing electrostatic screening. When the Mg(2+) ions are removed, the internal mobility of the RNA and of the protein increases at the interaction interface close to the RNA G-U/G-C motif as a result of a gap between the phosphate groups in the UUCG capping tetraloop and of the disruption of S15-RNA hydrogen bonds in that region. On the other hand, several S15-RNA hydrogen bonds are reinforced, and water bridges appear between the three-way junction region and S15. The network of hydrogen bonds observed in the loop between alpha1 and alpha2 is consequently reorganized. In the absence of Mg(2+), this network has the same pattern as the network observed in the isolated protein, where the helix alpha1 is mobile with respect to the protein core. The presence of Mg(2+) ions may thus play a role in stabilizing the orientation of the helix alpha1 of S15.

Riboswitches: fold and function

Riboswitches in the 5'-untranslated regions of mRNAs cotranscriptionally couple ligand binding and gene regulation. In this issue of Chemistry & Biology, Lemay et al. describe folding of a key tertiary interaction in the adenine riboswitch and its mechanistic consequences.

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.

The binding interface between Bacillus stearothermophilus ribosomal protein S15 and its 5'-translational operator mRNA

The Bacillus stearothermophilus ribosomal protein S15 (BS15) binds a purine-rich three-helix junction motif in the central domain of 16S ribosomal RNA (rRNA) as well as a translational operator located in the 5'-untranslated region (5'-UTR) of its cognate messenger RNA (mRNA). An in-frame fusion between the 5'-UTR of the BS15 gene and beta-galactosidase (lacZ) was prepared, and tested for BS15-dependent translational repression of lacZ activity in Escherichia coli. The presence of BS15 in trans represses lacZ activity 24-fold. A series of detailed point mutations in BS15 were tested for their effects upon translational repression of lacZ activity. These point mutations demonstrated that the 5'-UTR-BS15 binding interface utilizes many of the same conserved amino acid residues implicated in the binding of BS15 to 16S rRNA. The data demonstrate that the S15 protein can bind to an RNA target motif based primarily upon appropriate minor groove and sugar-phosphate backbone contacts, irrespective of the specific RNA sequence.

RNA tertiary structure and cooperative assembly of a large ribonucleoprotein complex

The mechanisms that govern the ordered assembly of multiprotein ribonucleoprotein complexes are not well understood. The in vitro reconstitution of the small subunit of the bacterial ribosome provides a tractable system for the detailed study of ordered assembly. We present a quantitative thermodynamic description of the hierarchical binding of ribosomal proteins to 16S rRNA during assembly of the platform of the 30S ribosomal subunit. The binding of S8, S11, S15, and the S6:S18 heterodimer to the central domain of 16S rRNA has been measured both individually and in combination using isothermal titration calorimetry and gel mobility shift assays. Both enthalpy and free energy measurements demonstrate the cooperative binding of S15 and the S6:S18 heterodimer, but no cooperativity is observed for either S8 or S11. The results define a thermodynamic framework that describes cooperative platform assembly.

The N-glycosidase activity of the ribosome-inactivating protein ME1 targets single-stranded regions of nucleic acids independent of sequence or structural motifs

ME(1), a type I ribosome-inactivating protein (RIP), belongs to a family of enzymes long believed to possess rRNA N-glycosidase activity directed solely at the universally conserved residue A4324 in the sarcin/ricin loop of large eukaryotic and prokaryotic rRNAs. We have investigated the effect of modifying the structure of nonribosomal RNA substrates on their interaction with ME(1) and other RIPs. ME(1) was shown to depurinate a variety of partially denatured nucleic acids, randomly removing adenine residues from single-stranded regions and, to a lesser extent, guanine residues from wobble base-pairs in hairpin stems. A defined sequence motif was not required for recognition of non-paired adenosines and cleavage of the N-glycosidic bond. Substrate recognition and ME(1) activity appeared to depend on the physical availability of nucleotides, and denaturation of nucleic acid substrates increased their interaction with ME(1). Pretreatment of mRNA at 75 degrees C rather than 60 degrees C, for example, lowered the apparent K(D) from 87.1 to 73.9 nm, making it more vulnerable to depurination by RIPs. Exposure to ME(1) in vitro completely abolished the infectivity of partially denatured RNA transcripts of the potato spindle tuber viroid, suggesting that RIPs may target invading nucleic acids before they reach host ribosomes in vivo. Our data suggest that the extensive folding of many potential substrates interferes with their ability to interact with RIPs, thereby blocking their inactivation by ME(1) (or other RIPs).

Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations

The goal of the current study is to utilize molecular dynamic (MD) simulations to investigate the dynamic behavior of 16S rRNA in the presence and absence of S15 and to identify the binding interactions between these two molecules. The simulations show that: (i) 16S rRNA remains in a highly folded structure when it is bound to S15; (ii) in the absence of S15, 16S rRNA significantly alters its conformation and transiently forms conformations that are similar to the bound structure that make it available for binding with S15; (iii) the unbound rRNA spends the majority of its time in extended conformations. The formation of the extended conformations is a result of the molecule reaching a lower electrostatic energy and the formation of the highly folded, crystal-like conformation is a result of achieving a lower solvation energy. In addition, our MD simulations show that 16S rRNA and S15 bind across the major groove of helix 22 (H22) via electrostatic interactions. The negatively charged phosphate groups of G658, U740, G741 and G742 bind to the positively charged S15 residues Lys7, Arg34 and Arg37. The current study provides a dynamic view of the binding of 16S rRNA with S15.

The identification of spermine binding sites in 16S rRNA allows interpretation of the spermine effect on ribosomal 30S subunit functions

A photoreactive analogue of spermine, N1-azidobenzamidino (ABA)-spermine, was covalently attached after irradiation to Escherichia coli 30S ribosomal subunits or naked 16S rRNA. By means of RNase H digestion and primer extension, the cross-linking sites of ABA-spermine in naked 16S rRNA were characterised and compared with those identified in 30S subunits. The 5' domain, the internal and terminal loops of helix H24, as well as the upper part of helix H44 in naked 16S rRNA, were found to be preferable binding sites for polyamines. Association of 16S rRNA with ribosomal proteins facilitated its interaction with photoprobe, except for 530 stem-loop nt, whose modification by ABA-spermine was abolished. Association of 30S with 50S subunits, poly(U) and AcPhe-tRNA (complex C) further altered the susceptibility of ABA-spermine cross-linking to 16S rRNA. Complex C, modified in its 30S subunit by ABA-spermine, reacted with puromycin similarly to non-photolabelled complex. On the contrary, poly(U)-programmed 70S ribosomes reconstituted from photolabelled 30S subunits and untreated 50S subunits bound AcPhe-tRNA more efficiently than untreated ribosomes, but were less able to recognise and reject near cognate aminoacyl-tRNA. The above can be interpreted in terms of conformational changes in 16S rRNA, induced by the incorporation of ABA-spermine.

The folding pathway of the genomic hepatitis delta virus ribozyme is dominated by slow folding of the pseudoknots

Hepatitis delta virus (HDV) replicates by a double rolling-circle mechanism that requires self-cleavage by closely related genomic and antigenomic versions of a ribozyme. We have previously shown that the uncleaved genomic ribozyme is subject to a variety of alternative (Alt) pairings. Sequence upstream of the ribozyme can regulate self-cleavage activity by formation of an Alt 1 ribozyme-containing structure that severely inhibits self-cleavage, or a P(-1) self-structure that permits rapid self-cleavage. Here, we test three other alternative pairings: Alt P1, Alt 2, and Alt 3. Alt P1 and Alt 3 contain primarily ribozyme-ribozyme interactions, while Alt 2 involves ribozyme-flanking sequence interaction. A number of single and double mutant ribozymes were prepared to increase or decrease the stability of the alternative pairings, and rates of self-cleavage were determined. Results of these experiments were consistent with the existence of the proposed alternative pairings and their ability to inhibit the overall rate of native ribozyme folding. Local misfolds are treated as internal equilibrium constants in a binding polynomial that modulates the intrinsic rate of cleavage. This model of equilibrium effects of misfolds should be general and apply to other ribozymes. All of the alternative pairings sequester a pseudoknot-forming strand. Folding of ribozymes containing Alt P1 and Alt 2 was accelerated by urea as long as the native ribozyme fold was sufficiently stable, while folding of mutants in which both of these alternative pairings had been removed were not stimulated by urea. This behavior suggests that the pseudoknots form by capture of an unfolded or appropriately rearranged alternative pairing by its complementary native strand. Fast-folding mutants were prepared by either weakening alternative pairings or by strengthening native pairings. A kinetic model was developed that accommodates these features and explains the position of the rate-limiting step for the G11C mutant. Implications of these results for structural and dynamic studies of the uncleaved HDV ribozyme are discussed.

Mg2+-dependent conformational change of RNA studied by fluorescence correlation and FRET on immobilized single molecules

Fluorescence correlation spectroscopy (FCS) of fluorescence resonant energy transfer (FRET) on immobilized individual fluorophores was used to study the Mg2+-facilitated conformational change of an RNA three-helix junction, a structural element that initiates the folding of the 30S ribosomal subunit. Transitions of the RNA junction between open and folded conformations resulted in fluctuations in fluorescence by FRET. Fluorescence fluctuations occurring between two FRET states on the millisecond time scale were found to be dependent on Mg2+ and Na+ concentrations. Correlation functions of the fluctuations were used to determine transition rates between the two conformations as a function of Mg2+ or Na+ concentration. Both the opening and folding rates were found to vary with changing salt conditions. Assuming specific binding of divalent ions to RNA, the Mg2+ dependence of the observed rates cannot be explained by conformational change induced by Mg2+ binding/unbinding, but is consistent with a model in which the intrinsic conformational change of the RNA junction is altered by uptake of Mg2+ ion(s). This version of FCS/FRET on immobilized single molecules is demonstrated to be a powerful technique in the study of conformational dynamics of biomolecules over time scales ranging from microseconds to seconds.

Design, construction, and analysis of a novel class of self-folding RNA

RNA can play multiple biological roles through use of its three-dimensional (3-D) structures. Recent advances in RNA structural biology have revealed that complex RNA 3D structures are assemblages of double-stranded helices with a variety of tertiary structural motifs. By employing RNA tertiary structural motifs together with the helices, we designed a novel class of self-folding RNA. In RNA composed of three helices (P1, P2, and P3), P1 interacts with P3 via a tetraloop-receptor interaction and P2 forms consecutive base-triples. Two designed RNAs of this class were prepared and their folding properties indicate that they form defined tertiary structures as designed. These RNAs may be used as modular units for constructing artificial ribozymes or nanometer-scale materials.

Molecular dynamics simulations of the denaturation and refolding of an RNA tetraloop

Tetraloops are very abundant structural elements of RNA that are formed by four nucleotides in a hairpin loop which is closed by a double stranded helical stem with some Watson-Crick base pairs. A tetraloop r(GCGAAGGC) was identified from the crystal structure of the central domain of 16S rRNA (727-730) in the Thermus thermophilus 30S ribosomal complex. The crystal structure of the 30S complex includes a total of 104 nucleotides from the central domain of the 16S rRNA and three ribosomal proteins S6, S15 and S18. Independent biochemical experiments have demonstrated that protein S15 plays the role in initiating the formation of the central domain of this complex. In the crystal, the tetraloop interacts with the protein S15 at two sites: one of them is associated with hydrogen bond interactions between residue His50 and nucleotide G730, and the other is associated with the occurrence of residue Arg53 beside A728. This paper uses molecular dynamics (MD) simulations to investigate the protein-dependent structural stability of the tetraloop and demonstrates the folding pathway of this tetraloop via melting denaturation and its subsequent refolding. Three important results are derived from these simulations: (i) The stability of nucleotide A728 appears to be protein dependent. Without the interaction with S15, A728 flips away from stacking with A729. (ii) The melting temperature demonstrated by the simulations is analogous to the results of thermodynamic experiments. In addition, the simulated folding of the tetraloop is stepwise: the native shape of the backbone is formed first; this is then followed by the formation of the Watson- Crick base pairs in the stem; and finally the hydrogen bonds and base stacking in the loop are formed. (iii) The tetraloop structure is similar to the crystal structure at salt concentrations of 0.1 M and 1.0 M used for the simulations, but the refolded structure at 0.1 M salt is more comparable to the crystal structure than at 1.0 M. The results from the simulations using both the Generalized Born continuum model and explicit solvent model (Particle Mesh Ewald) generate a similar pathway for unfolding/refolding of the tetraloop.

Interaction of the Bacillus stearothermophilus ribosomal protein S15 with its 5'-translational operator mRNA

The Bacillus stearothermophilus ribosomal protein S15 (BS15) binds both a three-helix junction in the central domain of 16 S ribosomal RNA and its cognate mRNA. Native gel mobility-shift assays show that BS15 interacts specifically and with high affinity to the 5'-untranslated region (5'-UTR) of this cognate mRNA with an apparent dissociation constant of 3(+/-0.3) nM. In order to localize the structural elements that are essential for BS15 recognition, a series of deletion mutants of the full cognate mRNA were prepared and tested in the same gel-shift assay. The minimal binding site for BS15 is a 50 nucleotide RNA showing a close secondary structure resemblance to the BS15 binding region from 16 S rRNA. There are two major structural motifs that must be maintained for high-affinity binding. The first being a purine-rich three-helix junction, and the second being an internal loop. The sequence identity of the internal loops differs greatly between the BS15 mRNA and rRNA sites, and this difference is correlated to discrimination between wild-type BS15 and a BS15(H45R) mutant. The association and dissociation kinetics measured for the 5'-UTR-BS15 interaction are quite slow, but are typical for a ribosomal protein-RNA interaction. The BS15 mRNA and 16 S rRNA binding sites share a common secondary structure yet have little sequence identity. The mRNA and rRNA may in fact present similar if not identical structural elements that confer BS15 recognition.

Central domain assembly: thermodynamics and kinetics of S6 and S18 binding to an S15-RNA complex

The 30 S ribosomal subunit assembles in vitro through the hierarchical binding of 21 ribosomal proteins to 16 S rRNA. The central domain of 16 S rRNA becomes the platform of the 30 S subunit upon binding of ribosomal proteins S6, S8, S11, S15, S18 and S21. The assembly of the platform is nucleated by binding of S15 to 16 S rRNA, followed by the cooperative binding of S6 and S18. The prior binding of S6 and S18 is required for binding of S11 and S21. We have studied the mechanism of the cooperative binding of S6 and S18 to the S15-rRNA complex by isothermal titration calorimetry and gel mobility shift assays with rRNA and proteins from the hyperthermophilic bacterium Aquifex aeolicus. S6 and S18 form a stable heterodimer in solution with an apparent dissociation constant of 8.7 nM at 40 degrees C. The S6:S18 heterodimer binds to the S15-rRNA complex with an equilibrium dissociation constant of 2.7 nM at 40 degrees C. Consistent with previous studies using rRNA and proteins from Escherichia coli, we observed no binding of S6 or S18 in the absence of the other protein or S15. The presence of S15 increases the affinity of S6:S18 for the RNA by at least four orders of magnitude. The kinetics of S6:S18 binding to the S15-rRNA complex are slow, with an apparent bimolecular rate constant of 8.0 x 10(4) M(-1) s(-1) and an apparent unimolecular dissociation rate of 1.6 x 10(-4) s(-1). These results, which are consistent with a model in which S6 and S18 bind as a heterodimer to the S15-rRNA complex, provide a mechanistic framework to describe the previously observed S15-mediated cooperative binding of S6 and S18 in the ordered assembly of a multi-protein ribonucleoprotein complex.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Beyond kinetic traps in RNA folding

Large RNAs often have rugged folding energy landscapes that result in severe misfolding and slow folding kinetics. Several interdependent parameters that contribute to misfolding are now well understood and examples of large RNAs and ribonucleoproteins that avoid kinetic traps have been reported. These advances have facilitated the exploration of fundamental RNA folding processes that were previously inaccessible.

Effects of cobalt hexammine on folding and self-cleavage of the Neurospora VS ribozyme

We have investigated the effects of Co(NH3)6(3+), an analog of hexahydrated Mg2+, on folding and catalysis of the Neurospora VS ribozyme. Most of the metal ion-induced changes detected by chemical modification structure probing in either metal ion are similar, but occur at approximately 33-fold lower concentrations of Co(NH3)6(3+) than Mg2+. However, Co(NH3)6(3+) is not as effective at inducing two functionally important structural changes: stabilizing the pseudoknot interaction between loops I and V, and rearranging the secondary structure of helix Ib. Comparison of the folding of the precursor and the downstream cleavage product, which lacks helix Ia, shows that helix Ia inhibits stable pseudoknot formation and rearrangement of helix Ib. The VS ribozyme does not self-cleave with Co(NH3)6(3+) as the sole polyvalent cation; however, mixed-metal kinetic experiments show that Co(NH3)6(3+) does not inhibit Mg2+-induced self-cleavage. In contrast, at sub-saturating concentrations of Mg2+, Co(NH3)6(3+) increases the rate of Mg2+-induced self-cleavage, indicating that Co(NH3)6(3+) contributes to the functionally relevant folding of the VS ribozyme.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Folding of branched RNA species

The global structures of branched RNA species are important to their function. Branched RNA species are defined as molecules in which double-helical segments are interrupted by abrupt discontinuities. These include helical junctions of different orders, and base bulges and loops. Common helical junctions are three- and four-way junctions, often interrupted by mispairs or additional nucleotides. There are many interesting examples of functional RNA junctions, including the hammerhead and hairpin ribozymes, and junctions that serve as binding sites for proteins. The junctions display some common structural properties. These include a tendency to undergo pairwise helical stacking and ion-induced conformational transitions. Helical branchpoints can act as key architectural components and as important sites for interactions with proteins. Copyright 1999 John Wiley & Sons, Inc.

Role of conserved nucleotides in building the 16 S rRNA binding site for ribosomal protein S15

Ribosomal protein S15 recognizes a highly conserved target on 16 S rRNA, which consists of two distinct binding regions. Here, we used extensive site-directed mutagenesis on a Escherichia coli 16 S rRNA fragment containing the S15 binding site, to investigate the role of conserved nucleotides in protein recognition and to evaluate the relative contribution of the two sites. The effect of mutations on S15 recognition was studied by measuring the relative binding affinity, RNA probing and footprinting. The crystallographic structure of the Thermus thermophilus complex allowed molecular modelling of the E. coli complex and facilitated interpretation of biochemical data. Binding is essentially driven by site 1, which includes a three-way junction constrained by a conserved base triple and cross-strand stacking. Recognition is based mainly on shape complementarity, and the role of conserved nucleotides is to maintain a unique backbone geometry. The wild-type base triple is absolutely required for protein interaction, while changes in the conserved surrounding nucleotides are partially tolerated. Site 2, which provides functional groups in a conserved G-U/G-C motif, contributes only modestly to the stability of the interaction. Binding to this motif is dependent on binding at site 1 and is allowed only if the two sites are in the correct relative orientation. Non-conserved bulged nucleotides as well as a conserved purine interior loop, although not directly involved in recognition, are used to provide an appropriate flexibility between the two sites. In addition, correct binding at the two sites triggers conformational adjustments in the purine interior loop and in a distal region, which are known to be involved for subsequent binding of proteins S6 and S18. Thus, the role of site 1 is to anchor S15 to the rRNA, while binding at site 2 is aimed to induce a cascade of events required for subunit assembly.

Structure of the S15,S6,S18-rRNA complex: assembly of the 30S ribosome central domain

The crystal structure of a 70-kilodalton ribonucleoprotein complex from the central domain of the Thermus thermophilus 30S ribosomal subunit was solved at 2.6 angstrom resolution. The complex consists of a 104-nucleotide RNA fragment composed of two three-helix junctions that lie at the end of a central helix, and the ribosomal proteins S15, S6, and S18. S15 binds the ribosomal RNA early in the assembly of the 30S ribosomal subunit, stabilizing a conformational reorganization of the two three-helix junctions that creates the RNA fold necessary for subsequent binding of S6 and S18. The structure of the complex demonstrates the central role of S15-induced reorganization of central domain RNA for the subsequent steps of ribosome assembly.

The identification of the A-type RNA helices in a 55mer RNA by selective incorporation of deuterium-labelled nucleotide residues (Uppsala NMR-window concept)

The 55-nt long RNA, modelling a three-way junction, with non-uniformly incorporated deuterated nucleotides has been synthesised in a pure form. The NMR-window part in this partially deuterated 55mer RNA consists of natural non-enriched nucleotide blocks at the three-way junction (shown in a square box in Fig. 2), whereas all other nucleotides of the rest of the molecule are partially deuterated (> 97 atom% 2H at C2', C3', C5', C5, and approximately 50 atom% 2H at C4'). The secondary structure of this 55mer RNA was determined by 2D 1H NOESY spectroscopy in D2O or in 10% D2O-H2O mixture. The use of deuterated building blocks in the specific region of the 55mer RNA allowed us to identify two distinct A-type RNA helices in a straightforward manner by observing connectivities of H1' with the basepaired imino and the aromatic H2 of all adenosine nucleotides as the first step for the determination of its tertiary structure in a cost- and time-effective manner without employing any 13C/15N labelling. These two decameric helices involve 40 nucleotides, for which all non-exchangeable H1', H6, H2, H8 and H5 protons (all 40 H1', all 40 H6 or H8 aromatics, all seven H2 of adenine nucleotide and all four non-deuterated H5 of cytosines) as well as all 16 exchangeable imino protons (with the exception of four terminal basepairs) and 16 amino protons of cytosines have been assigned. Since all aromatic-H2', H3' as well as H5'/5'' crosspeaks from partially deuterated residues have been eliminated from the NMR spectra, the observation of natural nucleotide residues in the NMR window part has essentially been simplified. It has been found that the crosspeaks from the natural nucleotides located at the three-way junction in the NMR-window part show different degrees of line-broadening, thereby indicating that the various nucleotide residues have very different mobilities with respect to themselves as well as compared to other nucleotides in the helices. The assignment of H2' and H3' in the NMR-window part has been made based on NOESY and DQF-COSY crosspeaks. It is noteworthy that, even in this preliminary study, it has been possible to identify 10 H2' out of total 14 and 9 H3' out of 14. The data show that expanded AU containing a tract of 55mer RNA does not self-organise into a tight third helix, as the two decameric A-type helices, across the three-way junction which is evident from the absence of any additional imino protons, except those that already have been assigned for the two decameric helices.

Thermodynamics of stabilization of RNA pseudoknots by cobalt(III) hexaammine

Equilibrium unfolding (folding) studies reveal that the autoregulatory RNA pseudoknots derived from the bacteriophage T2 and T4 gene 32 mRNAs exhibit significant stabilization by increasing concentrations of divalent metal ions in solution. In this report, the apparent affinities of exchange inert trivalent Co(NH(3))(3+)(6) have been determined, relative to divalent Mg(2+), for the folded, partially folded (K(f)), and fully unfolded (K(u)) conformations of these molecules. A general nonspecific, delocalized ion binding model was developed and applied to the analysis of the metal ion concentration dependence of individual two-state unfolding transitions. Trivalent Co(NH(3))(3+)(6) was found to associate with the fully folded and partially unfolded pseudoknotted forms of these RNAs with a K(f) of 5-8 x 10(4) M(-1) in a background of 0.10 M K(+), or 3- to 5-fold larger than the K(f) obtained for two model RNA hairpins and hairpin unfolding intermediates, and approximately 40-50-fold larger than K(f) for Mg(2+). The magnitude of K(f) was found to be strongly dependent on the monovalent salt concentration in a manner qualitatively consistent with polyelectrolyte theory, with K(f) reaching 1.2 x 10(5) M(-1) in 50 mM K(+). Two RNA hairpins were found to have affinities for Co(NH(3))(3+)(6) and Ru(NH(3))(3+)(6) of 1-2 x10(4) M(-1), or approximately 15-fold larger than the K(f) of approximately 1000 M(-1) observed for Mg(2+). Additionally, the K(u) of 4,800 M(-1) for the trivalent ligands is approximately 8-fold larger than the K(u) of 600 M(-1) observed for Mg(2+). These findings suggest that the T2 and T4 gene 32 mRNA pseudoknots possess a site(s) for Mg(2+) and Co(NH(3))(3+)(6) binding of significantly higher affinity than a "duplexlike" delocalized ion binding site that is strongly linked to the thermodynamic stability of these molecules. Imino proton perturbation nmr spectroscopy suggests that this site(s) lies near the base of the pseudoknot stem S2, near a patch of high negative electrostatic potential associated with the region where the single loop L1 adenosine crosses the major groove of stem S2.