Enhanced group II intron retrohoming in magnesium-deficient Escherichia coli via selection of mutations in the ribozyme core

Hydrolytic endonucleolytic ribozyme (HYER) is programmable for sequence-specific DNA cleavage

Ribozymes are catalytic RNAs with diverse functions including self-splicing and polymerization. This work aims to discover natural ribozymes that behave as hydrolytic and sequence-specific DNA endonucleases, which could be repurposed as DNA manipulation tools. Focused on bacterial group II-C introns, we found that many systems without intron-encoded protein propagate multiple copies in their resident genomes. These introns, named HYdrolytic Endonucleolytic Ribozymes (HYERs), cleaved RNA, single-stranded DNA, bubbled double-stranded DNA (dsDNA), and plasmids in vitro. HYER1 generated dsDNA breaks in the mammalian genome. Cryo-electron microscopy analysis revealed a homodimer structure for HYER1, where each monomer contains a Mg2+-dependent hydrolysis pocket and captures DNA complementary to the target recognition site (TRS). Rational designs including TRS extension, recruiting sequence insertion, and heterodimerization yielded engineered HYERs showing improved specificity and flexibility for DNA manipulation.

Thermodynamic determination of RNA duplex stability in magnesium solutions

The prediction of RNA secondary structure and thermodynamics from sequence relies on free energy minimization and nearest neighbor parameters. Currently, algorithms used to make these predictions are based on parameters from optical melting studies performed in 1 M NaCl. However, many physiological and biochemical buffers containing RNA include much lower concentrations of monovalent cations and the presence of divalent cations. In order to improve these algorithms, thermodynamic data was previously collected for RNA duplexes in solutions containing 71, 121, 221, and 621 mM Na+. From this data, correction factors for free energy (ΔG°37) and melting temperature (Tm) were derived. Despite these newly derived correction factors for sodium, the stabilizing effects of magnesium have been ignored. Here, the same RNA duplexes were melted in solutions containing 0.5, 1.5, 3.0, and 10.0 mM Mg2+ in the absence of monovalent cations. Correction factors for Tm and ΔG°37 were derived to scale the current parameters to a range of magnesium concentrations. The Tm correction factor predicts the melting temperature within 1.2°C, and the ΔG°37 correction factor predicts the free energy within 0.30 kcalmol. These newly derived magnesium correction factors can be incorporated into algorithms that predict RNA secondary structure and stability from sequence.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Functional Roles of Chelated Magnesium Ions in RNA Folding and Function

RNA regulates myriad cellular events such as transcription, translation, and splicing. To perform these essential functions, RNA often folds into complex tertiary structures in which its negatively charged ribose-phosphate backbone interacts with metal ions. Magnesium, the most abundant divalent metal ion in cells, neutralizes the backbone, thereby playing essential roles in RNA folding and function. This has been known for more than 50 years, and there are now thousands of in vitro studies, most of which have used ≥10 mM free Mg²⁺ ions to achieve optimal RNA folding and function. In the cell, however, concentrations of free Mg²⁺ ions are much lower, with most Mg²⁺ ions chelated by metabolites. In this Perspective, we curate data from a number of sources to provide extensive summaries of cellular concentrations of metabolites that bind Mg²⁺ and to estimate cellular concentrations of metabolite-chelated Mg²⁺ species, in the representative prokaryotic and eukaryotic systems Escherichia coli, Saccharomyces cerevisiae, and iBMK cells. Recent research from our lab and others has uncovered the fact that such weakly chelated Mg²⁺ ions can enhance RNA function, including its thermodynamic stability, chemical stability, and catalysis. We also discuss how metabolite-chelated Mg²⁺ complexes may have played roles in the origins of life. It is clear from this analysis that bound Mg²⁺ should not be simply considered non-RNA-interacting and that future RNA research, as well as protein research, could benefit from considering chelated magnesium.

Unraveling the role of the enigmatic MatK maturase in chloroplast group IIA intron excision

Maturases are prokaryotic enzymes that aid self-excision of introns in precursor RNAs and have evolutionary ties to the nuclear spliceosome. Both the mitochondria and chloroplast, due to their prokaryotic origin, encode a single intron maturase, MatR for the mitochondria and MatK for the chloroplast. MatK is proposed to aid excision of seven different chloroplast group IIA introns that reside within precursor RNAs for essential elements of chloroplast function. We have developed an in vitro activity assay to test chloroplast group IIA intron excision. Using this assay, we demonstrate self-excision of the group IIA intron of the second intron of rps12 and the group IIA intron of rpl2. We further show that the addition of heterologously expressed MatK protein increases efficiency of group IIA intron self-splicing for the second intron of rps12 but not the group IIA intron of rpl2. These data, to our knowledge, provide the first direct evidence of MatK's maturase activity.

Molecular Mechanism for Folding Cooperativity of Functional RNAs in Living Organisms

A diverse set of organisms has adapted to live under extreme conditions. The molecular origin of the stability is unclear, however. It is not known whether the adaptation of functional RNAs, which have intricate tertiary structures, arises from strengthening of tertiary or secondary structure. Herein we evaluate effects of sequence changes on the thermostability of tRNAphe using experimental and computational approaches. To separate out effects of secondary and tertiary structure on thermostability, we modify base pairing strength in the acceptor stem, which does not participate in tertiary structure. In dilute solution conditions, strengthening secondary structure leads to non-two-state thermal denaturation curves and has small effects on thermostability, or the temperature at which tertiary structure and function are lost. In contrast, under cellular conditions with crowding and Mg2+-chelated amino acids, where two-state cooperative unfolding is maintained, strengthening secondary structure enhances thermostability. Investigation of stabilities of each tRNA stem across 44 organisms with a range of optimal growing temperatures revealed that organisms that grow in warmer environments have more stable stems. We also used Shannon entropies to identify positions of higher and lower information content, or sequence conservation, in tRNAphe and found that secondary structures have modest information content allowing them to drive thermal adaptation, while tertiary structures have maximal information content hindering them from participating in thermal adaptation. Base-paired regions with no tertiary structure and modest information content thus offer a facile evolutionary route to enhancing the thermostability of functional RNA by the simple molecular rules of base pairing.

Cooperative RNA Folding under Cellular Conditions Arises From Both Tertiary Structure Stabilization and Secondary Structure Destabilization

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 Mg²⁺. 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 Mg²⁺. As such, RNA structures and functions can be radically altered under cellular conditions. We previously reported that tRNA^phe 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 Mg²⁺. 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 Mg²⁺, 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.

Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

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 pK_a 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 Mg²⁺ ions. We introduce global fitting for analyzing ribozyme rate-pH data parametric in Mg²⁺ 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 pK_a 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 Mg²⁺ ion. The value for the pK_a of the -1 nucleotide, C17, is increased a remarkable 3.5-5 pK_a units toward neutrality. Shifting of the pK_a 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 Mg²⁺ 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.

Mobile self-splicing introns and inteins as environmental sensors

Self-splicing introns and inteins are often mobile at the level of the genome. Although these RNA and protein elements, respectively, are generally considered to be selfish parasites, group I and group II introns and inteins can be triggered by environmental cues to splice and/or to mobilize. These cues include stressors such as oxidizing agents, reactive oxygen and nitrogen species, starvation, temperature, osmolarity and DNA damage. Their sensitivity to these stimuli leads to a carefully choreographed dance between the mobile element and its host that is in tune with the cellular environment. This responsiveness to a changing milieu provides strong evidence that these diverse, self-splicing mobile elements have adapted to react to prevailing conditions, to the potential advantage of both the element and its host.

Programmable Genome Editing Tools and their Regulation for Efficient Genome Engineering

Targeted genome editing has become a powerful genetic tool for studying gene function or for modifying genomes by correcting defective genes or introducing genes. A variety of reagents have been developed in recent years that can generate targeted double-stranded DNA cuts which can be repaired by the error-prone, non-homologous end joining repair system or via the homologous recombination-based double-strand break repair pathway provided a suitable template is available. These genome editing reagents require components for recognizing a specific DNA target site and for DNA-cleavage that generates the double-stranded break. In order to reduce potential toxic effects of genome editing reagents, it might be desirable to control the in vitro or in vivo activity of these reagents by incorporating regulatory switches that can reduce off-target activities and/or allow for these reagents to be turned on or off. This review will outline the various genome editing tools that are currently available and describe the strategies that have so far been employed for regulating these editing reagents. In addition, this review will examine potential regulatory switches/strategies that can be employed in the future in order to provide temporal control for these reagents.

Genome-Wide Analysis of RNA Secondary Structure

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

Host Factors Influencing the Retrohoming Pathway of Group II Intron RmInt1, Which Has an Intron-Encoded Protein Naturally Devoid of Endonuclease Activity

Bacterial group II introns are self-splicing catalytic RNAs and mobile retroelements that have an open reading frame encoding an intron-encoded protein (IEP) with reverse transcriptase (RT) and RNA splicing or maturase activity. Some IEPs carry a DNA endonuclease (En) domain, which is required to cleave the bottom strand downstream from the intron-insertion site for target DNA-primed reverse transcription (TPRT) of the inserted intron RNA. Host factors complete the insertion of the intron. By contrast, the major retrohoming pathway of introns with IEPs naturally lacking endonuclease activity, like the Sinorhizobium meliloti intron RmInt1, is thought to involve insertion of the intron RNA into the template for lagging strand DNA synthesis ahead of the replication fork, with possible use of the nascent strand to prime reverse transcription of the intron RNA. The host factors influencing the retrohoming pathway of such introns have not yet been described. Here, we identify key candidates likely to be involved in early and late steps of RmInt1 retrohoming. Some of these host factors are common to En⁺ group II intron retrohoming, but some have different functions. Our results also suggest that the retrohoming process of RmInt1 may be less dependent on the intracellular free Mg²⁺ concentration than those of other group II introns.

Bridging the gap between in vitro and in vivo RNA folding

Deciphering the folding pathways and predicting the structures of complex three-dimensional biomolecules is central to elucidating biological function. RNA is single-stranded, which gives it the freedom to fold into complex secondary and tertiary structures. These structures endow RNA with the ability to perform complex chemistries and functions ranging from enzymatic activity to gene regulation. Given that RNA is involved in many essential cellular processes, it is critical to understand how it folds and functions in vivo. Within the last few years, methods have been developed to probe RNA structures in vivo and genome-wide. These studies reveal that RNA often adopts very different structures in vivo and in vitro, and provide profound insights into RNA biology. Nonetheless, both in vitro and in vivo approaches have limitations: studies in the complex and uncontrolled cellular environment make it difficult to obtain insight into RNA folding pathways and thermodynamics, and studies in vitro often lack direct cellular relevance, leaving a gap in our knowledge of RNA folding in vivo. This gap is being bridged by biophysical and mechanistic studies of RNA structure and function under conditions that mimic the cellular environment. To date, most artificial cytoplasms have used various polymers as molecular crowding agents and a series of small molecules as cosolutes. Studies under such in vivo-like conditions are yielding fresh insights, such as cooperative folding of functional RNAs and increased activity of ribozymes. These observations are accounted for in part by molecular crowding effects and interactions with other molecules. In this review, we report milestones in RNA folding in vitro and in vivo and discuss ongoing experimental and computational efforts to bridge the gap between these two conditions in order to understand how RNA folds in the cell.

Bacterial Riboswitches and Ribozymes Potently Activate the Human Innate Immune Sensor PKR

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.

Mobile Bacterial Group II Introns at the Crux of Eukaryotic Evolution

This review focuses on recent developments in our understanding of group II intron function, the relationships of these introns to retrotransposons and spliceosomes, and how their common features have informed thinking about bacterial group II introns as key elements in eukaryotic evolution. Reverse transcriptase-mediated and host factor-aided intron retrohoming pathways are considered along with retrotransposition mechanisms to novel sites in bacteria, where group II introns are thought to have originated. DNA target recognition and movement by target-primed reverse transcription infer an evolutionary relationship among group II introns, non-LTR retrotransposons, such as LINE elements, and telomerase. Additionally, group II introns are almost certainly the progenitors of spliceosomal introns. Their profound similarities include splicing chemistry extending to RNA catalysis, reaction stereochemistry, and the position of two divalent metals that perform catalysis at the RNA active site. There are also sequence and structural similarities between group II introns and the spliceosome's small nuclear RNAs (snRNAs) and between a highly conserved core spliceosomal protein Prp8 and a group II intron-like reverse transcriptase. It has been proposed that group II introns entered eukaryotes during bacterial endosymbiosis or bacterial-archaeal fusion, proliferated within the nuclear genome, necessitating evolution of the nuclear envelope, and fragmented giving rise to spliceosomal introns. Thus, these bacterial self-splicing mobile elements have fundamentally impacted the composition of extant eukaryotic genomes, including the human genome, most of which is derived from close relatives of mobile group II introns.

Using Group II Introns for Attenuating the In Vitro and In Vivo Expression of a Homing Endonuclease

In Chaetomium thermophilum (DSM 1495) within the mitochondrial DNA (mtDNA) small ribosomal subunit (rns) gene a group IIA1 intron interrupts an open reading frame (ORF) encoded within a group I intron (mS1247). This arrangement offers the opportunity to examine if the nested group II intron could be utilized as a regulatory element for the expression of the homing endonuclease (HEase). Constructs were generated where the codon-optimized ORF was interrupted with either the native group IIA1 intron or a group IIB type intron. This study showed that the expression of the HEase (in vivo) in Escherichia coli can be regulated by manipulating the splicing efficiency of the HEase ORF-embedded group II introns. Exogenous magnesium chloride (MgCl₂) stimulated the expression of a functional HEase but the addition of cobalt chloride (CoCl₂) to growth media antagonized the expression of HEase activity. Ultimately the ability to attenuate HEase activity might be useful in precision genome engineering, minimizing off target activities, or where pathways have to be altered during a specific growth phase.

Retrohoming of a Mobile Group II Intron in Human Cells Suggests How Eukaryotes Limit Group II Intron Proliferation

Mobile bacterial group II introns are evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of an autocatalytic intron RNA (a “ribozyme”) and an intron-encoded reverse transcriptase, which function together to promote intron integration into new DNA sites by a mechanism termed “retrohoming”. Although mobile group II introns splice and retrohome efficiently in bacteria, all examined thus far function inefficiently in eukaryotes, where their ribozyme activity is limited by low Mg²⁺ concentrations, and intron-containing transcripts are subject to nonsense-mediated decay (NMD) and translational repression. Here, by using RNA polymerase II to express a humanized group II intron reverse transcriptase and T7 RNA polymerase to express intron transcripts resistant to NMD, we find that simply supplementing culture medium with Mg²⁺ induces the Lactococcus lactis Ll.LtrB intron to retrohome into plasmid and chromosomal sites, the latter at frequencies up to ~0.1%, in viable HEK-293 cells. Surprisingly, under these conditions, the Ll.LtrB intron reverse transcriptase is required for retrohoming but not for RNA splicing as in bacteria. By using a genetic assay for in vivo selections combined with deep sequencing, we identified intron RNA mutations that enhance retrohoming in human cells, but <4-fold and not without added Mg²⁺. Further, the selected mutations lie outside the ribozyme catalytic core, which appears not readily modified to function efficiently at low Mg²⁺ concentrations. Our results reveal differences between group II intron retrohoming in human cells and bacteria and suggest constraints on critical nucleotide residues of the ribozyme core that limit how much group II intron retrohoming in eukaryotes can be enhanced. These findings have implications for group II intron use for gene targeting in eukaryotes and suggest how differences in intracellular Mg²⁺ concentrations between bacteria and eukarya may have impacted the evolution of introns and gene expression mechanisms.

Mobile group II introns are bacterial retrotransposons that are evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of an autocatalytic intron RNA (a ribozyme) and an intron-encoded reverse transcriptase, which together promote intron mobility to new DNA sites by a mechanism called retrohoming. Although found in bacteria, archaea and eukaryotic organelles, group II introns are absent from eukaryotic nuclear genomes, where host defenses impede their expression and lower intracellular Mg²⁺ concentrations limit their ribozyme activity. Here, we developed a mobile group II intron expression system that bypasses expression barriers and show that simply adding Mg²⁺ to culture medium enables group II intron retrohoming into plasmid and chromosomal target sites in human cells at appreciable frequencies. Genetic selections and deep sequencing identified intron RNA mutations that moderately enhance retrohoming in human cells, but not without added Mg²⁺. Thus, low Mg²⁺ concentrations in human cells are a natural barrier to efficient retrohoming that is not readily overcome by mutational variation and selection. Our results have implications for group II intron use for gene targeting in higher organisms and highlight the impact of different intracellular environments on intron evolution and gene expression mechanisms in bacteria and eukarya.

Inactivation of group II intron RmInt1 in the Sinorhizobium meliloti genome

Group II introns are self-splicing catalytic RNAs that probably originated in bacteria and act as mobile retroelements. The dispersal and dynamics of group II intron spread within a bacterial genome are thought to follow a selection-driven extinction model. Likewise, various studies on the evolution of group II introns have suggested that they are evolving toward an inactive form by fragmentation, with the loss of the intron 3′-terminus, but with some intron fragments remaining and continuing to evolve in the genome. RmInt1 is a mobile group II intron that is widespread in natural populations of Sinorhizobium meliloti, but some strains of this species have no RmInt1 introns. We studied the splicing ability and mobility of the three full-length RmInt1 copies harbored by S. meliloti 1021, and obtained evidence suggesting that specific mutations may lead to the impairment of intron splicing and retrohoming. Our data suggest that the RmInt1 copies in this strain are undergoing a process of inactivation.

Mechanistic Analysis of Activation of the Innate Immune Sensor PKR by Bacterial RNA

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.5mM, 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.

The role of Mg(II) in DNA cleavage site recognition in group II intron ribozymes: solution structure and metal ion binding sites of the RNA-DNA complex

Group II intron ribozymes catalyze the cleavage of (and their reinsertion into) DNA and RNA targets using a Mg2(+)-dependent reaction. The target is cleaved 3' to the last nucleotide of intron binding site 1 (IBS1), one of three regions that form base pairs with the intron's exon binding sites (EBS1 to -3).We solved the NMR solution structure of the d3' hairpin of the Sc.ai5γ intron containing EBS1 in its 11-nucleotide loop in complex with the dIBS1 DNA 7-mer and compare it with the analogous RNA-RNA contact. The EBS1-dIBS1 helix is slightly flexible and non-symmetric. NMR data reveal two major groove binding sites for divalent metal ions at the EBS1-dIBS1 helix, and surface plasmon resonance experiments show that low concentrations of Mg2(+) considerably enhance the affinity of dIBS1 for EBS1. Our results indicate that identification of both RNA and DNA IBS1 targets, presentation of the scissile bond, and stabilization of the structure by metal ions are governed by the overall structure of EBS1-dIBS1 and the surrounding loop nucleotides but are irrespective of different EBS1-(d)IBS1 geometries and interstrand affinities.

Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: quantum mechanical/molecular mechanical free energy simulations

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 Mg(2+) 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.

IRE mRNA riboregulators use metabolic iron (Fe(2+)) to control mRNA activity and iron chemistry in animals

A family of noncoding RNAs bind Fe(2+) to increase protein synthesis. The structures occur in messenger RNAs encoding animal proteins for iron metabolism. Each mRNA regulatory sequence, ∼30 ribonucleotides long, is called an IRE (Iron Responsive Element), and folds into a bent, A-RNA helix with a terminal loop. Riboregulatory RNAs, like t-RNAs, r-RNAs micro-RNAs, etc. contrast with DNA, since single-stranded RNA can fold into a variety of complex, three-dimensional structures. IRE-RNAs bind two types of proteins: (1) IRPs which are protein repressors, sequence-related to mitochondrial aconitases. (2) eIF-4F, which bind ribosomes and enhances general protein biosynthesis. The competition between IRP and eIF-4F binding to IRE-RNA is controlled by Fe(2+)-induced changes in the IRE-RNA conformation. Mn(2+), which also binds to IRE-RNA in solution, is a convenient experimental proxy for air-sensitive Fe(2+) studies of in vitro protein biosynthesis and protein binding. However, only Fe(2+) has physiological effects on protein biosynthesis directed by IRE-mRNAs. The structures of the IRE-RNA riboregulators is known indirectly from effects of base substitutions on function, from solution NMR of the free RNA, and of X-ray crystallography of the IRE-RNA-IRP repressor complex. However, the inability to date, to crystallize the free IRE-RNA, and the dissociation of the IRE-RNA-IRP complex when metal binds, have hampered direct identification and characterization of the RNA-metal binding sites. The high conservation of the primary sequence in IRE-mRNA control elements has facilitated their identification and analysis of metal-assisted riboregulator function. Expansion of RNA search analyses beyond primary will likely reveal other, metal-dependent families of mRNA riboregulators.

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg(2+) ion concentrations are low, K(+) concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo-like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg(2+) (0.5-2 mM) and K(+) (140 mM) if the solution is supplemented with physiological amounts (∼ 20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.

Biotechnological applications of mobile group II introns and their reverse transcriptases: gene targeting, RNA-seq, and non-coding RNA analysis

Mobile group II introns are bacterial retrotransposons that combine the activities of an autocatalytic intron RNA (a ribozyme) and an intron-encoded reverse transcriptase to insert site-specifically into DNA. They recognize DNA target sites largely by base pairing of sequences within the intron RNA and achieve high DNA target specificity by using the ribozyme active site to couple correct base pairing to RNA-catalyzed intron integration. Algorithms have been developed to program the DNA target site specificity of several mobile group II introns, allowing them to be made into ‘targetrons.’ Targetrons function for gene targeting in a wide variety of bacteria and typically integrate at efficiencies high enough to be screened easily by colony PCR, without the need for selectable markers. Targetrons have found wide application in microbiological research, enabling gene targeting and genetic engineering of bacteria that had been intractable to other methods. Recently, a thermostable targetron has been developed for use in bacterial thermophiles, and new methods have been developed for using targetrons to position recombinase recognition sites, enabling large-scale genome-editing operations, such as deletions, inversions, insertions, and ‘cut-and-pastes’ (that is, translocation of large DNA segments), in a wide range of bacteria at high efficiency. Using targetrons in eukaryotes presents challenges due to the difficulties of nuclear localization and sub-optimal magnesium concentrations, although supplementation with magnesium can increase integration efficiency, and directed evolution is being employed to overcome these barriers. Finally, spurred by new methods for expressing group II intron reverse transcriptases that yield large amounts of highly active protein, thermostable group II intron reverse transcriptases from bacterial thermophiles are being used as research tools for a variety of applications, including qRT-PCR and next-generation RNA sequencing (RNA-seq). The high processivity and fidelity of group II intron reverse transcriptases along with their novel template-switching activity, which can directly link RNA-seq adaptor sequences to cDNAs during reverse transcription, open new approaches for RNA-seq and the identification and profiling of non-coding RNAs, with potentially wide applications in research and biotechnology.