Time-resolved RNA SHAPE chemistry

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry yields quantitative RNA secondary and tertiary structure information at single nucleotide resolution. SHAPE takes advantage of the discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl group is modulated by local nucleotide flexibility in the RNA backbone. Flexible nucleotides are reactive toward hydroxyl-selective electrophiles, whereas constrained nucleotides are unreactive. Initial versions of SHAPE chemistry, which employ isatoic anhydride derivatives that react on the minute time scale, are emerging as the ideal technology for monitoring equilibrium structures of RNA in a wide variety of biological environments. Here, we extend SHAPE chemistry to a benzoyl cyanide scaffold to make possible facile time-resolved kinetic studies of RNA in approximately 1 s snapshots. We then use SHAPE chemistry to follow the time-dependent folding of an RNase P specificity domain RNA. Tertiary interactions form in two distinct steps with local tertiary contacts forming an order of magnitude faster than long-range interactions. Rate-determining tertiary folding requires minutes despite that no non-native interactions must be disrupted to form the native structure. Instead, overall folding is limited by simultaneous formation of interactions approximately 55 A distant in the RNA. Time-resolved SHAPE holds broad potential for understanding structural biogenesis and the conformational interconversions essential to the functions of complex RNA molecules at single nucleotide resolution.     

Structural insight into RNA hairpin folding intermediates

Hairpins are a ubiquitous secondary structure motif in RNA molecules. Despite their simple structure, there is some debate over whether they fold in a two-state or multi-state manner. We have studied the folding of a small tetraloop hairpin using a serial version of replica exchange molecular dynamics on a distributed computing environment. On the basis of these simulations, we have identified a number of intermediates that are consistent with experimental results. We also find that folding is not simply the reverse of high-temperature unfolding and suggest that this may be a general feature of biomolecular folding.     

2'-mercaptonucleotide interference reveals regions of close packing within folded RNA molecules

The 2'-hydroxyl group makes essential contributions to RNA structure and function. As an approach to assess the ability of a mercapto group to serve as a functional analogue for the 2'-hydroxyl group, we synthesized 2'-mercaptonucleotides for use in nucleotide analogue interference mapping. To correlate the observed interference effects with tertiary structure, we used the independently folding DeltaC209 P4-P6 domain from the Tetrahymena group I intron. We generated populations of DeltaC209 P4-P6 molecules containing 2'-mercaptonucleotides located randomly throughout the domain and separated the folded molecules from the unfolded molecules by nondenaturing gel electrophoresis. Iodine-induced cleavage of the RNA molecules revealed the sites at which 2'-mercaptonucleotides interfere with folding. These interferences cluster in the most densely packed regions of the tertiary structure, occurring only at sites that lack the space and flexibility to accommodate a sulfur atom. Interference mapping with 2'-mercaptonucleotides therefore provides a method by which to identify structurally rigid and densely packed regions within folded RNA molecules.     

A powerful approach for the selection of 2-aminopurine substitution sites to investigate RNA folding

A precise tertiary structure must be adopted to allow the function of many RNAs in cells. Accordingly, increasing resources have been devoted to the elucidation of RNA structures and the folding of RNAs. 2-Aminopurine (2AP), a fluorescent nucleobase analogue, can be substituted in strategic positions of DNA or RNA molecules to act as site-specific probe to monitor folding and folding dynamics of nucleic acids. Recent studies further demonstrated the potential of 2AP modifications in the assessment of folding kinetics during ligand-induced secondary and tertiary RNA structure rearrangements. However, an efficient way to unambiguously identify reliable positions for 2AP sensors is as yet unavailable and would represent a major asset, especially in the absence of crystallographic or NMR structural data for a target molecule. We report evidence of a novel and direct correlation between the 2'-OH flexibility of nucleotides, observed by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probing and the fluorescence response following nucleotide substitutions by 2AP. This correlation leads to a straightforward method, using SHAPE probing with benzoyl cyanide, to select appropriate nucleotide sites for 2AP substitution. This clear correlation is presented for three model RNAs of biological significance: the SAM-II, adenine (addA), and preQ(1) class II (preQ(1)cII) riboswitches.     

The earliest events in protein folding: a structural requirement for ultrafast folding in cytochrome C

The folding dynamics of reduced cytochrome c (redcyt c) obtained from tuna heart, which contains a tryptophan residue at the site occupied by His33 in horse heart cytochrome c, was studied using nanosecond time-resolved optical rotatory dispersion spectroscopy. As observed previously for horse heart redcyt c, two time regimes were observed for secondary structure formation in tuna redcyt c: a fast (microseconds) and a slow (milliseconds) phase. However, the fast phase of tuna redcyt c folding was much slower and smaller in amplitude than the same phase in horse. The differences in the fast folding phases suggest that for horse heart redcyt c, the conformers that undergo the fastest observed folding have the His18-Fe-His33 heme configuration, which appears to be necessary, but not sufficient, to poise an unfolded chain conformation for fastest folding in redcyt c.     

Formation of a Ligand‐Assisted Complex of Two RNA Hairpin Loops

The hairpin structure is one of the most common secondary structures in RNA and holds a central position in the stream of RNA folding from a non‐structured RNA to structurally complex and functional ribonucleoproteins. Since the RNA secondary structure is strongly correlated to the function and can be modulated by the binding of small molecules, we have investigated the modulation of RNA folding by a ligand‐assisted formation of loop–loop complexes of two RNA hairpin loops. With a ligand (NCT6), designed based on the ligand binding to the G–G mismatches in double‐stranded DNA, we successfully demonstrated the formation of both inter‐ and intra‐molecular NCT6‐assisted complex of two RNA hairpin loops. NCT6 selectively bound to the two hairpin loops containing (CGG)₃ in the loop region. Native polyacrylamide gel electrophoresis analysis of two doubly‐labeled RNA hairpin loops clearly showed the formation of intermolecular NCT6‐assisted loop–loop complex. Förster resonance energy‐transfer studies of RNA constructs containing two hairpin loops, in which each hairpin was labeled with Alexa488 and Cy3 fluorophores, showed the conformational change of the RNA constructs upon binding of NCT6. These experimental data showed that NCT6 simultaneously bound to two hairpin RNAs at the loop region, and can induce the conformational change of the RNA molecule. These data strongly support that NCT6 functions as molecular glue for two hairpin RNAs.     

Determining the Mg2+ stoichiometry for folding an RNA metal ion core

The folding and catalytic function of RNA molecules depend on their interactions with divalent metal ions, such as magnesium. As with every molecular process, the most basic knowledge required for understanding the close relationship of an RNA with its metal ions is the stoichiometry of the interaction. Unfortunately, inventories of the numbers of divalent ions associated with unfolded and folded RNA states have been unattainable. A common approach has been to interpret Hill coefficients fit to folding equilibria as the number of metal ions bound upon folding. However, this approach is vitiated by the presence of diffusely associated divalent ions in a dynamic ion atmosphere and by the likelihood of multiple transitions along a folding pathway. We demonstrate that the use of molar concentrations of background monovalent salt can alleviate these complications. These simplifying solution conditions allow a precise determination of the stoichiometry of the magnesium ions involved in folding the metal ion core of the P4-P6 domain of the Tetrahymena group I ribozyme. Hill analysis of hydroxyl radical footprinting data suggests that the P4-P6 RNA core folds cooperatively upon the association of two metal ions. This unexpectedly small stoichiometry is strongly supported by counting magnesium ions associated with the P4-P6 RNA via fluorescence titration and atomic emission spectroscopy. By pinpointing the metal ion stoichiometry, these measurements provide a critical but previously missing step in the thermodynamic dissection of the coupling between metal ion binding and RNA folding.     

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.     

Length-Dependent Deep Learning Model for RNA Secondary Structure Prediction

Deep learning methods for RNA secondary structure prediction have shown higher performance than traditional methods, but there is still much room to improve. It is known that the lengths of RNAs are very different, as are their secondary structures. However, the current deep learning methods all use length-independent models, so it is difficult for these models to learn very different secondary structures. Here, we propose a length-dependent model that is obtained by further training the length-independent model for different length ranges of RNAs through transfer learning. 2dRNA, a coupled deep learning neural network for RNA secondary structure prediction, is used to do this. Benchmarking shows that the length-dependent model performs better than the usual length-independent model.     

Direct measurement of tertiary contact cooperativity in RNA folding

All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.     

Selective stabilization of natively folded RNA structure by DNA constraints

Learning how native RNA conformations can be stabilized relative to unfolded states is an important objective, for both understanding natural RNAs and improving the design of artificial functional RNAs. Here we show that covalently attached double-stranded DNA constraints (ca. 14 base pairs in length) can significantly stabilize the native conformation of an RNA molecule. Using the P4-P6 domain of the Tetrahymena group I intron as the test system, we identified pairs of RNA sites where attaching a DNA duplex is predicted to be structurally compatible with only the folded state of the RNA. The DNA-constrained RNAs were synthesized and shown by nondenaturing polyacrylamide gel electrophoresis (native PAGE) to have substantial decreases in their Mg2+ midpoints ([Mg2+]1/2 values). These changes are equivalent to free energy stabilizations as large as DeltaDeltaGdegrees = -2.5 kcal/mol, which is approximately 14% of the total tertiary folding energy. For comparison, the sole modification of P4-P6 previously reported to stabilize this RNA is a single-nucleotide deletion (DeltaC209) that provides only 1.1 kcal/mol of stabilization. Our findings indicate that nature has not completely optimized P4-P6 RNA folding. Furthermore, the DNA constraints are designed not to interact directly and extensively with the RNA, but rather more indirectly to modulate the relative stabilities of folded and unfolded RNA states. The successful implementation of this strategy to further stabilize a natively folded RNA conformation suggests an important element of modularity in stabilization of RNA structure, with implications for how nature might use other molecules such as proteins to stabilize specific RNA conformations.     

Determination of the intrinsic affinities of multiple site-specific Mg(2+) ions coordinated to domain 6 of a group II intron ribozyme

Group II introns are large metallo-ribozymes that use divalent metal ions in folding and catalysis. The 3'-terminal domain 6 (D6) contains a conserved adenosine whose 2'-OH group acts as the nucleophile in the first splicing step. In the hierarchy of folding, D6 binds last into the active site. In order to investigate and understand the folding process to the catalytically active intron structure, it is important to know the individual binding affinities of Mg2+ ions to D6. We recently studied the solution structure of a 27 nucleotide long D6 (D6-27) from the mitochondrial yeast group II intron Sc.ai5gamma, also identifying five Mg2+ binding sites including the one at the 5'-terminal phosphate residues. Mg2+ coordination to the 5'-terminal di- and triphosphate groups is strongest (e.g., log KA,TP = 4.55 +/- 0.10) and is evaluated here in detail for the first time. The other four binding sites within D6-27 are filled simultaneously (e.g., log KA,BR = 2.38 +/- 0.06) and thus compete for the free Mg2+ ions in solution, having a distinct influence on the individual affinities of the various sites. For the first time, we take this competition into account to obtain the intrinsic binding constants, describing a method that is generally applicable. Our data illustrates that any RNA molecule undergoing tertiary contacts to a second RNA molecule first needs to be loaded evenly and specifically with metal ions to compensate for the repulsion between the negatively charged RNA molecules.     

Metal-induced specific and nonspecific oligonucleotide folding studied by FRET and related biophysical and bioanalytical implications

Metal induced nucleic acid folding has been extensively studied with ribozymes, DNAzymes, tRNA and riboswitches. These RNA/DNA molecules usually have a high content of double-stranded regions to support a rigid scaffold. On the other hand, such rigid structural features are not available for many in vitro selected or rationally designed DNA aptamers; they adopt flexible random coil structures in the absence of target molecules. Upon target binding, these aptamers adaptively fold into a compact structure with a reduced end-to-end distance, making fluorescence resonance energy transfer (FRET) a popular signaling mechanism. However, nonspecific folding induced by mono- or divalent metal ions can also reduce the end-to-end distance and thus lead to false positive results. In this study we used a FRET pair labeled Hg(II) binding DNA and monitored metal-induced folding in the presence of various cations. While nonspecific electrostatically mediated folding can be very significant, at each tested salt condition, Hg(II) induced folding was still observed with a similar sensitivity. We also studied the biophysical meaning of the acceptor/donor fluorescence ratio that allowed us to explain the experimental observations. Potential solutions for this ionic strength problem have been discussed. For example, probes designed to signal the formation of double-stranded DNA showed a lower dependency on ionic strength.     

Slow conformational dynamics at C2'-endo nucleotides in RNA

RNA molecules undergo local conformational dynamics on timescales spanning picoseconds to minutes. Slower local motions have the greater potential to govern RNA folding, ligand recognition, and ribonucleoprotein assembly reactions but are difficult to detect in large RNAs with complex structures. RNA SHAPE chemistry employs acylation of the ribose 2'-hydroxyl position to measure local nucleotide flexibility in RNA and is well-characterized by a mechanism in which each nucleotide samples unreactive (closed) and reactive (open) states. We monitor RNA conformational dynamics over distinct time domains by varying the electrophilicity of the acylating reagent. Select C2'-endo nucleotides are nonreactive toward fast reagents but reactive toward slower SHAPE reagents in both model RNAs and in a large RNA with a tertiary fold. We conclude, first, that the C2'-endo conformation by itself does not govern SHAPE reactivity. However, some C2'-endo nucleotides undergo extraordinarily slow conformational changes, on the order of 10(-4) s(-1). Due to their distinctive local dynamics, C2'-endo nucleotides have the potential to function as rate-determining molecular switches and are likely to play central, currently unexplored, roles in RNA folding and function.     

Optimum folding pathways of proteins: their determination and properties

We develop a dynamic optimization technique for determining optimum folding pathways of proteins starting from different initial configurations. A coarse-grained Go model is used. Forces acting on each bead are (i) the friction force, (ii) forces from bond length constraints, (iii) excluded volume constraints, and (iv) attractive forces between residue pairs that are in contact in the native state. An objective function is defined as the total attractive energy between nonbonded residues, which are neighbors in the native state. The objective function is minimized over all feasible paths, satisfying bond length and excluded volume constraints. The optimization problem is nonconvex and contains a large number of constraints. An augmented Lagrangian method with a penalty barrier function was used to solve the problem. The method is applied to a 36-residue protein, chicken villin headpiece. Sequences of events during folding of the protein are determined for various pathways and analyzed. The relative time scales are compared and scaled according to experimentally measured events. Formation times of the helices, turn, and the loop agree with experimental data. We obtain the overall folding time of the protein in the range of 600 ns-1.2 micros that is smaller than the experimental result of 4-5 micros, showing that the optimal folding times that we obtain may be possible lower bounds. Time dependent variables during folding and energies associated with short- and long-range interactions between secondary structures are analyzed in modal space using Karhunen-Loeve expansion.     

A Periodic System of Supramolecular Elements

Chemistry “beyond the molecule” is based on weak, noncovalent, and reversible interactions. As a consequence of these bonds being weak, structural organization by folding and self‐assembly can only be fully exploited with larger molecules that can provide multiple binding sites. Such “supramolecules” can now be synthesized and their folding into desired conformations predicted. A new level of chemistry can now be realized through the creation of non‐natural entities composed of molecular building blocks with defined secondary structures. Herein we define these building blocks as “supramolecular elements”. We anticipate that further research on such large molecules will reveal construction principles dictated by recurring motifs that govern structure formation through folding and self‐assembly. These principles are comparable to the organization of atoms in the Periodic Table of Chemical Elements and may lead to the establishment of a Periodic System of Supramolecular Elements.     

Biologically important reactions catalyzed by RNA molecules

The last few years have seen a considerable increase in our understanding of catalysis by naturally occurring RNA molecules called ribozymes. The biological functions of RNA molecules depend upon their adoption of appropriate three-dimensional structures. The structure of RNA has a very important electrostatic component, which results from the presence of charged phosphodiester bonds. Metal ions are usually required to stabilize the folded structures and/or catalysis. Some ribozymes utilize metal ions as catalysts, whereas others use the ions to maintain appropriate three-dimensional structures. In the latter case, the correct folding of the RNA structures can perturb the pKa values of the nucleotide(s) within a catalytic pocket such that they act as general acid/bases catalysts.     

Folding of 8-17 deoxyribozyme studied by three-color alternating-laser excitation of single molecules

The folding of 8-17 deoxyribozyme was investigated by three-color alternating-laser excitation (3c-ALEX), a new single-molecule fluorescence resonance energy transfer (FRET) method we recently developed. Since 3c-ALEX has the capability of simultaneously sorting fluorescent molecules based on their labeling status and monitoring three interprobe distances of a biomolecule by employing three-color FRET, it is an ideal tool to study folding of multibranched molecules. The 8-17 deoxyribozyme, a DNA enzyme that cleaves a specific RNA substrate, is a good model system for a multibranched molecule, since it has the structure of a three-way DNA junction with a bulge. Labeling all three branches of the 8-17 with different fluorescent probes, we studied its [Mg2+]-dependent folding in a Na+ buffer solution. With the stoichiometric sorting capability of 3c-ALEX, we first selected only the triply labeled 8-17 in a solution of all heterogeneous mixtures and then simultaneously measured all three interprobe distances of the selected species. Our results show that the 8-17 folds into a pyramidal form upon increasing [Mg2+], in a similar way with [Zn2+] as found in an earlier study conducted at the ensemble level. The apparent dissociation constant of Mg2+ was more than 100 times larger than that of Zn2+ and showed considerable variance with buffer concentration. No clear sign of two-step folding was observed for Mg2+, in contrast to the case of Zn2+. Compared with the hammerhead ribozyme, the 8-17 was found to require 10 times higher [Mg2+] to undergo folding. By comparison with the folding of several inactive 8-17 analogues, we found that the two conserved sequences (A and G) of the triad loop of the shortest branch are critical elements for folding, especially for the folding at low [Mg2+]. Our results suggest that the role of the stem loop is to provide a scaffold for the two bases to be properly positioned for the necessary interaction and that the two bases are directly involved in the interaction that plays a critical role in folding. This work demonstrates that 3c-ALEX is a powerful single-molecule method to study the structure and folding of complex and multibranched biomolecules.