General acid–base catalysis is a key mechanistic strategy in protein and RNA enzymes. Ribozymes use hydrated metal ions, nucleobases, and organic cofactors to carry this out. In most small ribozymes, a guanosine is positioned to participate in proton transfer with the nucleophilic 2′-OH. The unshifted pKa values for nucleobases and solvated metal ions are far from neutrality, however, and thus nonideal for general acid–base catalysis. Herein, evidence is provided for cooperative interaction in the hammerhead ribozyme among the guanine that interacts with the nucleophilic 2′-OH, G12, the −1 nucleobase C17, and Mg2+ ions. We introduce global fitting for analyzing ribozyme rate–pH data parametric in Mg2+ concentration and benchmark this method on data from the hepatitis delta virus ribozyme. We then apply global fitting to new rate–pH data for the hammerhead ribozyme using a minimal three-dimensional, four-channel cooperative model. The value for the pKa of G12 that we obtain is channel-dependent and varies from 8.1 to 9.9, shifting closest toward neutrality in the presence of two cationic species: C17H+ and a Mg2+ ion. The value for the pKa of the −1 nucleotide, C17, is increased a remarkable 3.5–5 pKa units toward neutrality. Shifting of the pKa of C17 appears to be driven by an electrostatic sandwich of C17 between carbonyl groups of the 5′-neighboring U and of G12 and involves cation−π interactions. Rate–pH profiles reveal that the major reactive channel under biological Mg2+ and pH involves a cationic C17 rather than a second metal ion. Substitution of a cationic base for a metal underscores the versatility of RNA.

RNA folding has been studied extensively in vitro, typically under dilute solution conditions and abiologically high salt concentrations of 1 M Na+ or 10 mM Mg2+. The cellular environment is very different, with 20–40% crowding and only 10–40 mM Na+, 140 mM K+, and 0.5–2.0 mM Mg2+. As such, RNA structures and functions can be radically altered under cellular conditions. We previously reported that tRNAphe secondary and tertiary structures unfold together in a cooperative two-state fashion under crowded in vivo-like ionic conditions, but in a noncooperative multistate fashion under dilute in vitro ionic conditions unless in nonphysiologically high concentrations of Mg2+. The mechanistic basis behind these effects remains unclear, however. To address the mechanism that drives RNA folding cooperativity, we probe effects of cellular conditions on structures and stabilities of individual secondary structure fragments comprising the full-length RNA. We elucidate effects of a diverse set of crowders on tRNA secondary structural fragments and full-length tRNA at three levels: at the nucleotide level by temperature-dependent in-line probing, at the tertiary structure level by small-angle X-ray scattering, and at the global level by thermal denaturation. We conclude that cooperative RNA folding is induced by two overlapping mechanisms: increased stability and compaction of tertiary structure through effects of Mg2+, and decreased stability of certain secondary structure elements through the effects of molecular crowders. These findings reveal that despite having very different chemical makeups RNA and protein can both have weak secondary structures in vivo leading to cooperative folding.

RNA enzymes, or ribozymes, catalyze internal phosphodiester bond cleavage using diverse catalytic strategies. These include the four classic strategies: in-line nucleophilic attack, deprotonation of the 2′-OH nucleophile, protonation of the 5′-O leaving group, and stabilization of developing charge on the nonbridging oxygen atoms of the scissile phosphate. In addition, we recently identified two additional ribozyme strategies: acidification of the 2′-OH and release of the 2′-OH from inhibitory interactions. Herein, we report inverse thio effects in the presence of glmS ribozyme variants and a 1-deoxyglucosamine 6-phosphate cofactor analogue and demonstrate that activation of the 2′-OH nucleophile is promoted by competitive hydrogen bonding among diverse ribozyme moieties for the pro-RP nonbridging oxygen. We conclude that the glmS ribozyme uses an overdetermined set of competing hydrogen bond donors in its active site to ensure potent activation and regulation by the cofactor. Nucleophile activation through competitive, overdetermined hydrogen bonding could be a general strategy for ribozyme activation and may be applicable for controlling the function of ribozymes and riboswitches in the laboratory.

•Rapid quench is a useful technique for probing fast (<5 s) ribozyme reactions.•Probing the reaction with rapid quench requires a well-behaved RNA system.•Guidelines and considerations are provided for sample preparation and experiment.•Troubleshooting for general rapid-quench procedures and issues are given.Reaction kinetics on the millisecond timescale pervade the protein and RNA fields. To study such reactions, investigators often perturb the system with abiological solution conditions or substrates in order to slow the rate to timescales accessible by hand mixing; however, such perturbations can change the rate-limiting step and obscure key folding and chemical steps that are found under biological conditions. Mechanical methods for collecting data on the millisecond timescale, which allow these perturbations to be avoided, have been developed over the last few decades. These methods are relatively simple and can be conducted on affordable and commercially available instruments. Here, we focus on using the rapid quench-flow technique to study the fast reaction kinetics of RNA enzymes, or ribozymes, which often react on the millisecond timescale under biological conditions. Rapid quench of ribozymes is completely parallel to the familiar hand-mixing approach, including the use of radiolabeled RNAs and fractionation of reactions on polyacrylamide gels. We provide tips on addressing and preventing common problems that can arise with the rapid-quench technique. Guidance is also offered on ensuring the ribozyme is properly folded and fast-reacting. We hope that this article will facilitate the broader use of rapid-quench instrumentation to study fast-reacting ribozymes under biological reaction conditions.

Pathogens are recognized by the innate immune system in part via their unique and complex RNA signatures. A key sensor in human innate immunity is the RNA-activated protein kinase, protein kinase R (PKR), which has two double-stranded RNA (dsRNA) binding motifs (dsRBMs) at its N-terminus. Early studies described PKR as being activated potently by long stretches of perfect dsRNA, a signature typical of viruses. More recently, we and others have found that PKR is also activated by RNAs having structural defects such as bulges and internal loops.This Account describes advances in our understanding of the ability of PKR to detect diverse foreign RNAs and how that recognition plays significant roles in discriminating self from non-self. The experiments discussed employ a wide range of techniques including activation assays, native polyacrylamide gel electrophoresis (PAGE), protein footprinting, and small-angle X-ray scattering (SAXS). We discuss how misfolding and dimerization of RNA lead to activation of PKR. We also present recent findings on the activation of PKR by varied bacterial functional RNAs including ribozymes and riboswitches, which are among the few structured RNAs known to interact with PKR in a site-specific manner. Molecular models for how these structured RNAs activate PKR are provided. Studies by SAXS revealed that PKR straightens bent RNAs. Most external and internal RNA cellular modifications introduced in vitro and found naturally, such as the m7G cap and m6A group, abrogate activation of PKR, but other modifications, such as 5′-ppp and 2′-fluoro groups, are immunostimulatory and potential anticancer agents.Genome-wide studies of RNA folding in vitro and in vivo have provided fresh insights into general differences in RNA structure among bacteria, viruses, and human. These studies suggest that in vivo, cellular human RNAs are less folded than once thought, unwound by helicases, destabilized by m6A modifications, and often bound up with proteins, all conditions known to abrogate activation of PKR. It thus appears that non-self RNAs are detected as unmodified, naked RNAs with appreciable secondary and tertiary structure. Observation that PKR is activated by structured but otherwise diverse RNAs is consistent both with the broad-spectrum nature of innate immunity and the nonspecific recognition of RNA by the dsRBM family. These findings provide a possible explanation for the apparent absence of protein-free structured human RNAs, such as ribozymes and riboswitches.

The innate immune system provides the first line of defense against pathogens through the recognition of nonspecific patterns in RNA to protect the cell in a generalized way. The human RNA-activated protein kinase, PKR, is a dsRNA binding protein and an essential sensor in the innate immune response, which recognizes viral and bacterial pathogens through their RNAs. Upon activation via RNA-dependent autophosphorylation, PKR phosphorylates the eukaryotic initiation factor eIF2α, leading to termination of translation. PKR has a well-characterized role in recognizing viral RNA, where it binds long stretches of double-stranded RNA nonsequence specifically to promote activation; however, the mechanism by which bacterial RNA activates PKR and the mode by which self RNA avoids activating PKR are unknown. We characterized activation of PKR by three functional bacterial RNAs with pseudoknots and extensive tertiary structure: the cyclic di-GMP riboswitch, the glmS riboswitch-ribozyme, and the twister ribozyme, two of which are ligand-activated. These RNAs were found to activate PKR with comparable potency to long dsRNA. Enzymatic structure mapping in the absence and presence of PKR reveals a clear PKR footprint and provides a structural basis for how these bacterial RNAs activate PKR. In the case of the cyclic di-GMP riboswitch and the glmS riboswitch-ribozyme, PKR appears to dimerize on the peripheral double-stranded regions of the native RNA tertiary structure. Overall, these results provide new insights into how PKR acts as an innate immune signaling protein for the presence of bacteria and suggest a reason for the apparent absence of protein-free riboswitches and ribozymes in the human genome.

Stretches of guanines in DNA and RNA can fold into guanine quadruplex structures (GQSs). These structures protect telomeres in DNA and regulate gene expression in RNA. GQSs have an intrinsic fluorescence that is sensitive to different parameters, including loop sequence and length. However, the dependence of GQS fluorescence on solution and sequence parameters and the origin of this fluorescence are poorly understood. Herein we examine effects of dangling nucleotides and cosolute conditions on GQS fluorescence using both steady-state and time-resolved fluorescence spectroscopy. The quantum yield of dGGGTGGGTGGGTGGG, termed “dG3T”, is found to be modest at ∼2 × 10–3. Nevertheless, dG3T and its variants are significantly brighter than the common nucleic acid fluorophore 2-aminopurine (2AP) largely due to their sizable extinction coefficients. Dangling 5′-end nucleotides generally reduce emission and blue-shift the resultant spectrum, whereas dangling 3′-end nucleotides slightly enhance fluorescence, particularly on the red side of the emission band. Time-resolved fluorescence decays are broadly distributed in time and require three exponential components for accurate fits. Time-resolved emission spectra suggest the presence of two emitting populations centered at ∼330 and ∼390 nm, with the redder component being a well-defined long-lived (∼1 ns) entity. Insights into GQS fluorescence obtained here should be useful in designing brighter intrinsic RNA and DNA quadruplexes for use in label-free biotechnological applications.

Phase separation of aqueous solutions containing polyelectrolytes can lead to formation of dense, solute-rich liquid droplets referred to as coacervates, surrounded by a dilute continuous phase of much larger volume. This type of liquid–liquid phase separation is thought to help explain the appearance of polyelectrolyte-rich intracellular droplets in the cytoplasm and nucleoplasm of extant biological cells and may be relevant to protocellular compartmentalization of nucleic acids on the early Earth. Here we describe complex coacervates formed upon mixing the polycation poly(allylamine) (PAH, 15 kDa) with the anionic nucleotides adenosine 5′-mono-, di-, and triphosphate (AMP, ADP, and ATP). Droplet formation was observed over a wide range of pH and MgCl2 concentrations. The nucleotides themselves as well as Mg2+ and RNA oligonucleotides were all extremely concentrated within the coacervates. Nucleotides present at just 2.5 mM in bulk solution had concentrations greater than 1 M inside the coacervate droplets. A solution with a total Mg2+ concentration of 10 mM had 1–5 M Mg2+ in the coacervates, and RNA random sequence (N54) partitioned ∼10 000-fold into the coacervates. Coacervate droplets are thus rich in nucleotides, Mg2+, and RNA, providing a medium favorable for generating functional RNAs. Compartmentalization of nucleotides at high concentrations could have facilitated their polymerization to form oligonucleotides, which preferentially accumulate in the droplets. Locally high Mg2+ concentrations could have aided folding and catalysis in an RNA world, making coacervate droplets an appealing platform for exploring protocellular environments.

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates. The RP substrate gives large normal thio effects in the presence of all monovalent ions. The SP substrate also gives normal or no thio effects, but only for smaller monovalent and divalent cations, such as Li+, Mg2+, Ca2+, and Sr2+; in contrast, sizable inverse thio effects are found for larger monovalent and divalent cations, including Na+, K+, NH4+, and Ba2+. Proton inventories are found to be unity in the presence of the larger monovalent and divalent ions, but two in the presence of Mg2+. Additionally, rate–pH profiles are inverted for the low charge density ions, and only imidazole plus ammonium ions rescue an inactive C75Δ variant in the absence of Mg2+. Results from the thio effect experiments, rate–pH profiles, proton inventories, and ammonium/imidazole rescue experiments, combined with QM/MM free energy simulations, support a change in the mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized to stepwise and proton transfer-stabilized as the charge density of the metal ion decreases.

The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid–base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2′), which would be followed by proton transfer from G40(N1) to A-1(O2′). After this initial deprotonation, A-1(O2′) starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2′) to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2′) to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg2+ ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa’s of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

The hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage reaction using a combination of nucleobase and metal ion catalysis. Both divalent and monovalent ions can catalyze this reaction, although the rate is slower with monovalent ions alone. Herein, we use quantum mechanical/molecular mechanical (QM/MM) free energy simulations to investigate the mechanism of this ribozyme and to elucidate the roles of the catalytic metal ion. With Mg2+ at the catalytic site, the self-cleavage mechanism is observed to be concerted with a phosphorane-like transition state and a free energy barrier of ∼13 kcal/mol, consistent with free energy barrier values extrapolated from experimental studies. With Na+ at the catalytic site, the mechanism is observed to be sequential, passing through a phosphorane intermediate, with free energy barriers of 2–4 kcal/mol for both steps; moreover, proton transfer from the exocyclic amine of protonated C75 to the nonbridging oxygen of the scissile phosphate occurs to stabilize the phosphorane intermediate in the sequential mechanism. To explain the slower rate observed experimentally with monovalent ions, we hypothesize that the activation of the O2′ nucleophile by deprotonation and orientation is less favorable with Na+ ions than with Mg2+ ions. To explore this hypothesis, we experimentally measure the pKa of O2′ by kinetic and NMR methods and find it to be lower in the presence of divalent ions rather than only monovalent ions. The combined theoretical and experimental results indicate that the catalytic Mg2+ ion may play three key roles: assisting in the activation of the O2′ nucleophile, acidifying the general acid C75, and stabilizing the nonbridging oxygen to prevent proton transfer to it.

Charged nucleobases exist in RNA and DNA at neutral pH owing to pKa shifting. These bases can affect polymerase fidelity and participate in ribozyme general acid–base catalysis. Protonated RNA bases further influence miRNA processing and viral frameshifting. It is therefore important to have a simple and rapid method for determining the pKa of nucleobases in RNA and DNA. Here we describe the application of 2-aminopurine (2AP), a fluorescent isomer of adenine, to report on the pKa of a nearby ionizing base both in DNA secondary structure and RNA tertiary structure. We observe large, up to 5-fold quenching in fluorescence upon protonation of a nearby base. Using this method, we identify highly shifted pKa’s of 7.6 for adenine in a DNA oligonucleotide and 8.15 for cytidine in a tertiary structure element from beet western yellows virus (BWYV) RNA. These pKa values, which were corroborated by 31P NMR measurements and comparison to literature, are shifted over 4 units from their standard values. This fluorescence method can be used to determine pKa’s for ionization of both A and C and reveals that shifted pKa’s are prevalent in DNA and RNA secondary and tertiary structures.

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pKa shifted toward neutrality. Less is known about the role of metal ions in the mechanism. A recent crystal structure of the precleavage ribozyme identified a Mg2+ ion that interacts through its partial hydration sphere with the G25·U20 reverse wobble. In addition, this Mg2+ ion is in position to directly coordinate the nucleophile, the 2′-hydroxyl of U(−1), suggesting it can serve as a Lewis acid to facilitate deprotonation of the 2′-hydroxyl. To test the role of the active site Mg2+ ion, we replaced the G25·U20 reverse wobble with an isosteric A25·C20 reverse wobble. This change was found to significantly reduce the negative potential at the active site, as supported by electrostatics calculations, suggesting that active site Mg2+ binding could be adversely affected by the mutation. The kinetic analysis and molecular dynamics of the A25·C20 double mutant suggest that this variant stably folds into an active structure. However, pH–rate profiles of the double mutant in the presence of Mg2+ are inverted relative to the profiles for the wild-type ribozyme, suggesting that the A25·C20 double mutant has lost the active site metal ion. Overall, these studies support a model in which the partially hydrated Mg2+ positioned at the G25·U20 reverse wobble is catalytic and could serve as a Lewis acid, a Brønsted base, or both to facilitate deprotonation of the nucleophile.

Guanine quadruplex structures (GQSs) exhibit unique spectroscopic features, including an inverse melting profile at 295 nm, distinctive circular dichroism features, and intrinsic fluorescence. Herein, we investigate effects of loop sequence and loop length on the intrinsic fluorescence of 13 DNA GQSs. We report label-free fluorescence enhancements upon intramolecular GQS formation of up to 16-fold and a shift in the emission maximum to the visible portion of the spectrum. Effects can be understood in the context of available nuclear magnetic resonance GQS structures. The intrinsic fluorescence of GQSs may be useful for nucleic acid studies and for the development of label-free detection methods.

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg2+ ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn2+. Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd2+. Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg2+ or Cd2+, supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ∼1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd2+ of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd2+ inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg2+ binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2′-hydroxyl, supporting the experimental findings.

In an effort to relate RNA folding to function under cellular-like conditions, we monitored the self-cleavage reaction of the human hepatitis delta virus-like CPEB3 ribozyme in the background of physiological ionic concentrations and various crowding and cosolute agents. We found that at physiological free Mg2+ concentrations (∼0.1–0.5 mM), both crowders and cosolutes stimulate the rate of self-cleavage, up to ∼6-fold, but that in 10 mM Mg2+ (conditions widely used for in vitro ribozyme studies) these same additives have virtually no effect on the self-cleavage rate. We further observe a dependence of the self-cleavage rate on crowder size, wherein the level of rate stimulation is diminished for crowders larger than the size of the unfolded RNA. Monitoring effects of crowding and cosolute agents on rates in biological amounts of urea revealed additive-promoted increases at both low and high Mg2+ concentrations, with a maximal stimulation of more than 10-fold and a rescue of the rate to its urea-free values. Small-angle X-ray scattering experiments reveal a structural basis for this stimulation in that higher-molecular weight crowding agents favor a more compact form of the ribozyme in 0.5 mM Mg2+ that is essentially equivalent to the form under standard ribozyme conditions of 10 mM Mg2+ without a crowder. This finding suggests that at least a portion of the rate enhancement arises from favoring the native RNA tertiary structure. We conclude that cellular-like crowding supports ribozyme reactivity by favoring a compact form of the ribozyme, but only under physiological ionic and cosolute conditions.

Shifting of pKa’s in RNA is important for many biological processes; however, the driving forces responsible for shifting are not well understood. Herein, we determine how structural environments surrounding protonated bases affect pKa shifting in double-stranded RNA (dsRNA). Using 31P NMR, we determined the pKa of the adenine in an A+·C base pair in various sequence and structural environments. We found a significant dependence of pKa on the base pairing strength of nearest neighbors and the location of a nearby bulge. Increasing nearest neighbor base pairing strength shifted the pKa of the adenine in an A+·C base pair higher by an additional 1.6 pKa units, from 6.5 to 8.1, which is well above neutrality. The addition of a bulge two base pairs away from a protonated A+·C base pair shifted the pKa by only ∼0.5 units less than a perfectly base paired hairpin; however, positioning the bulge just one base pair away from the A+·C base pair prohibited formation of the protonated base pair as well as several flanking base pairs. Comparison of data collected at 25 °C and 100 mM KCl to biological temperature and Mg2+ concentration revealed only slight pKa changes, suggesting that similar sequence contexts in biological systems have the potential to be protonated at biological pH. We present a general model to aid in the determination of the roles protonated bases may play in various dsRNA-mediated processes including ADAR editing, miRNA processing, programmed ribosomal frameshifting, and general acid–base catalysis in ribozymes.

Single-stranded DNA (ssDNA) ligation is a crucial step in many biochemical assays. Efficient ways of carrying out this reaction, however, are lacking. We show here that existing ssDNA ligation methods suffer from slow kinetics, poor yield, and severe nucleotide preference. To resolve these issues, we introduce a hybridization-based strategy that provides efficient and low-bias ligation of ssDNA. Our method uses a hairpin DNA to hybridize to any incoming acceptor ssDNA with low bias, with ligation of these strands mediated by T4 DNA ligase. This technique potentially can be applied in protocols that require ligation of ssDNA, including ligation-mediated polymerase chain reaction (LMPCR) and complementary DNA (cDNA) library construction.

Aptamers are RNA molecules that bind small molecules. They were originally isolated from random libraries and then found in bacteria. In this issue of Chemistry & Biology, Vu et al. demonstrate that a motif for an adenosine aptamer occurs in human and bacterial genomes and binds AMP and ATP in vitro.

The malachite green aptamer binds two closely related ligands, malachite green (MG) and tetramethylrosamine (TMR), with nearly equal affinity. The MG ligand consists of three phenyl rings emanating from a central carbon, while TMR has two of the three rings connected by an ether linkage. The binding pockets for MG and TMR in the aptamer, known from high-resolution structures, differ only in the conformation of a few nucleotides. Herein, we applied isothermal titration calorimetry (ITC) to compare the thermodynamics of binding of MG and TMR to the aptamer. Binding heat capacities were obtained from ITC titrations over the temperature range of 15–60 °C. Two temperature regimes were found for MG binding: one from 15 to 45 °C where MG bound with a large negative heat capacity and an apparent stoichiometry (n) of ∼0.4 and another from 50 to 60 °C where MG bound with a positive heat capacity and an n of ∼1.1. The binding of TMR, on the other hand, revealed only one temperature regime for binding, with a more modest negative heat capacity and an n of ∼1.2. The large difference in heat capacity between the two ligands suggests that significantly more conformational rearrangement occurs upon the binding of MG than that of TMR, which is consistent with differences in solvent accessible surface area calculated for available ligand-bound structures. Lastly, we note that the binding stoichiometry of MG was improved not only by an increase in the temperature but also by a decrease in the concentration of Mg2+ or an increase in the time between ITC injections. These studies suggest that binding of a dynamical ligand to a functional RNA requires the RNA itself to have significant dynamics.

The interferon-inducible, double-stranded (ds) RNA-activated protein kinase (PKR) contains a dsRNA-binding domain (dsRBD) and plays key roles in viral pathogenesis and innate immunity. Activation of PKR is typically mediated by long dsRNA, and regulation of PKR is disfavored by most RNA imperfections, including bulges and internal loops. Herein, we combine isothermal titration calorimetry (ITC), electrophoretic mobility shift assays, and small-angle X-ray scattering (SAXS) to dissect the thermodynamic basis for the specificity of the dsRBD termed “p20” for various RNAs and to detect any RNA conformational changes induced upon protein binding. We monitor binding of p20 to chimeric duplexes containing terminal RNA–DNA hybrid segments and a central dsRNA segment, which was either unbulged (“perfect”) or bulged. The ITC data reveal strong binding of p20 to the perfect duplex (Kd ∼ 30 nM) and weaker binding to the bulged duplex (Kd ∼ 2–5 μM). SAXS reconstructions and p(r) distance distribution functions further uncover that p20 induces no significant conformational change in perfect dsRNA but largely straightens bulged dsRNA. Together, these observations support the dsRBD’s ability to tightly bind to only A-form RNA and suggest that in a noninfected cell, PKR may be buffered via weak interactions with various bulged and looped RNAs, which it may straighten. This work suggests that PKR-regulating RNAs with complex secondary and tertiary structures likely mimic dsRNA and/or engage portions of PKR outside of the dsRBD.

Self-cleaving RNAs have recently been identified in mammalian genomes. A small ribozyme related in structure to the hepatitis delta virus (HDV) ribozyme occurs in a number of mammals, including chimpanzees and humans, within an intron of the CPEB3 gene. The catalytic mechanisms for the CPEB3 and HDV ribozymes appear to be similar, generating cleavage products with 5′-hydroxyl and 2′,3′-cyclic phosphate termini; nonetheless, the cleavage rate reported for the CPEB3 ribozyme is more than 6000-fold slower than for the fastest HDV ribozyme. Herein, we use full-length RNA and cotranscriptional self-cleavage assays to compare reaction rates among human CPEB3, chimp CPEB3, and HDV ribozymes. Our data reveal that a single base change of the upstream flanking sequence, which sequesters an intrinsically weak P1.1 pairing in a misfold, increases the rate of the wild-type human CPEB3 ribozyme by ∼250-fold; thus, the human ribozyme is intrinsically fast-reacting. Secondary structure determination and native gel analyses reveal that the cleaved population of the CPEB3 ribozyme has a single, secondary structure that closely resembles the HDV ribozyme. In contrast, the precleavage population of the CPEB3 ribozyme appears to have a more diverse secondary structure, possibly reflecting misfolding with the upstream sequence and dynamics intrinsic to the ribozyme. Prior identification of expressed sequence tags (ESTs) in human cells indicated that cleavage activity of the human ribozyme is tissue-specific. It is therefore possible that cellular factors interact with regions upstream of the CPEB3 ribozyme to unmask its high intrinsic reactivity.

•Bacterial mRNA 5′-UTR activates the innate immune sensor PKR.•Multiple secondary structural features in the 5′-UTR of trp mRNA from B. subtilis activate PKR.•Various bacterial RNA segments display a 5′-triphosphate dependence for PKR activation.•Tertiary structure in the 5′-UTR contributes to PKR activation.•Bacterial RNAs activate PKR under physiological magnesium concentrations.The protein kinase PKR (protein kinase R) is a sensor in innate immunity. PKR autophosphorylates in the presence of double-stranded RNA enabling it to phosphorylate its substrate, eIF2α (eukaryotic initiation factor 2α), halting cellular translation. Classical activators of PKR are long viral double-stranded RNAs, but recently, PKR has been found to be activated by bacterial RNA. However, the features of bacterial RNA that activate PKR are unknown. We studied the Bacillus subtilis trp 5′-UTR (untranslated region), which is an indirect riboswitch with secondary and tertiary RNA structures that regulate gene function. Additionally, the trp 5′-UTR binds a protein, TRAP (tryptophan RNA-binding attenuation protein), which recognizes l-tryptophan. We present the first evidence that multiple structural features in this RNA, which are typical of bacterial RNAs, activate PKR in TRAP-free and TRAP/l-Trp-bound forms. Segments from the 5′-UTR, including the terminator 5′-stem–loop and Shine–Dalgarno blocking hairpins, demonstrated 5′-triphosphate and flanking RNA tail dependence on PKR activation. Disruption of long-distance tertiary interactions in the 5′-UTR led to partial loss in activation, consistent with highly base-paired regions in bacterial RNA activating PKR. One physiological change a bacterial RNA would face in a human cell is a decrease in the concentration of free magnesium. Upon lowering the magnesium concentration to human physiological conditions of 0.5 mM, the trp 5′-UTR continued to activate PKR potently. Moreover, total RNA from Escherichia coli, depleted of rRNA, also activated PKR under these ionic conditions. This study demonstrates that PKR can signal the presence of bacterial RNAs under physiological ionic conditions and offers a potential explanation for the apparent absence of riboswitches in the human genome.Download high-res image (64KB)Download full-size image

•Single-stranded A and C nucleotides are reactive to DMS (dimethyl sulfate) in vivo.•mRNA regions encoding protein domains have high in vivo DMS reactivity.•This relationship is especially prominent for enzymes.•mRNA regions encoding ordered protein regions in intrinsically disordered proteins have high in vivo DMS reactivity.•Structural relationships between mRNAs and encoded proteins may facilitate protein folding.We assessed whether in vivo mRNA structural reactivity and the structure of the encoded protein are related. This is the first investigation of such a relationship that utilizes information on RNA structure obtained in living cells. Based on our recent genome-wide Structure-seq analysis in Arabidopsis thaliana, we report that, as a meta property, regions of individual mRNAs that code for protein domains generally have higher reactivity to DMS (dimethyl sulfate), a chemical that covalently modifies accessible As and Cs, than regions that encode protein domain junctions. This relationship is prominent for proteins annotated for catalytic activity and reversed in proteins annotated for binding and transcription regulatory activity. Upon analyzing intrinsically disordered proteins, we found a similar pattern for disordered regions as compared to ordered regions: regions of individual mRNAs that code for ordered regions have significantly higher DMS reactivity than regions that code for intrinsically disordered regions. Based on these effects, we hypothesize that the decreased DMS reactivity of RNA regions that encode protein domain junctions or intrinsically disordered regions may reflect increased RNA structure that may slow translation, allowing time for the nascent protein domain or ordered region of the protein to fold, thereby reducing protein misfolding. In addition, a drop in DMS reactivity was observed on portions of mRNA sequences that correspond to the C-termini of protein domains, suggesting ribosome protection at these mRNA regions. Structural relationships between mRNAs and their encoded proteins may have evolved to allow efficient and accurate protein folding.Download high-res image (83KB)Download full-size image