Structural basis of glmS ribozyme activation by glucosamine-6-phosphate

Some general principles of riboswitch structure and interactions with small-molecule ligands

Riboswitches are RNA elements with a defined structure found in noncoding sections of genes that allow the direct control of gene expression by the binding of small molecules functionally related to the gene product. In most cases, this is a metabolite in the same (typically biosynthetic) pathway as an enzyme (or transporter) encoded by the gene that is controlled. The structures of many riboswitches have been determined and this provides a large database of RNA structure and ligand binding. In this review, we extract general principles of RNA structure and the manner or ligand binding from this resource.

TaRTLEt: Transcriptionally-active Riboswitch Tracer Leveraging Edge deTection

Structured RNAs have emerged as a major component of cellular regulatory systems, but their mechanism of action is often poorly understood. Riboswitches are structured RNAs that allosterically regulate gene expression through any of several different mechanisms. In vitro approaches to characterizing this mechanism are costly, low-throughput, and must be repeated for each individual riboswitch locus of interest. Bioinformatic methods promise higher throughput; despite robust computational identification of riboswitches, however, computational classification of the riboswitch mechanism has so far been both model-bound, relying on identification of sequence motifs known to be required for specific models of riboswitch activity, and empirically untested, with predictions far outpacing biological validation. Here, we introduce TaRTLEt (Transcriptionally-active Riboswitch Tracer Leveraging Edge deTection), a new high-throughput tool that recovers in vivo patterns of riboswitch-mediated transcription termination from paired-end RNA-seq data using edge detection methods. TaRTLEt successfully extracts transcription termination signals despite numerous sources of biological and technical noise. We tested the effectiveness of TaRTLEt on riboswitches identified from a wide range of sequenced bacterial taxa by utilizing publicly available paired-end RNA-seq readsets, finding broad agreement with previously published in vitro characterization results. In addition, we use TaRTLEt to infer the in vivo regulatory mechanism of uncharacterized riboswitch loci from existing public data. TaRTLEt is available on GitHub and can be applied to paired-end RNA-seq datasets from isolates or complex communities.

Transferability of MACE Graph Neural Network for Range Corrected Δ-Machine Learning Potential QM/MM Applications

We previously introduced a "range corrected" Δ-machine learning potential (ΔMLP) that used deep neural networks to improve the accuracy of combined quantum mechanical/molecular mechanical (QM/MM) simulations by correcting both the internal QM and QM/MM interaction energies and forces [J. Chem. Theory Comput. 2021, 17, 6993-7009]. The present work extends this approach to include graph neural networks. Specifically, the approach is applied to the MACE message passing neural network architecture, and a series of AM1/d + MACE models are trained to reproduce PBE0/6-31G* QM/MM energies and forces of model phosphoryl transesterification reactions. Several models are designed to test the transferability of AM1/d + MACE by varying the amount of training data and calculating free energy surfaces of reactions that were not included in the parameter refinement. The transferability is compared to AM1/d + DP models that use the DeepPot-SE (DP) deep neural network architecture. The AM1/d + MACE models are found to reproduce the target free energy surfaces even in instances where the AM1/d + DP models exhibit inaccuracies. We train "end-state" models that include data only from the reactant and product states of the 6 reactions. Unlike the uncorrected AM1/d profiles, the AM1/d + MACE method correctly reproduces a stable pentacoordinated phosphorus intermediate even though the training did not include structures with a similar bonding pattern. Furthermore, the message passing mechanism hyperparameters defining the MACE network are varied to explore their effect on the model's accuracy and performance. The AM1/d + MACE simulations are 28% slower than AM1/d QM/MM when the ΔMLP correction is performed on a graphics processing unit. Our results suggest that the MACE architecture may lead to ΔMLP models with improved transferability.

Sugar-sensing swodkoreceptors and swodkocrine signaling

Sugars are one of the major metabolites and are essential for nucleic acid synthesis and energy production. In addition, sugars can act as signaling molecules. To study sugar signaling at the systemic level, there is an urgent need to systematically identify sugar-sensing proteins and nucleic acids. I propose the terms "swodkoreceptor" and "swodkocrine signaling," derived from the Polish word "słodki" meaning "sweet," to comprise all sugar-sensing proteins and signaling events, respectively, regardless of their cellular location and signaling domains. This proposal is intended to facilitate the inclusion of proteins such as the Escherichia coli LacI repressor as an allolactose receptor, human glucokinase regulatory protein (GCKR) as a fructose receptor, and other sugar-binding based allosterically regulated enzymes and transcription factors as sugar-sensing receptors. In addition, enzyme-interacting proteins whose interaction state is regulated by sugar binding have also been proposed as sugar receptors. The systemic study of protein- and nucleic-acid-based swodkoreceptors may help to identify organelle-specific swodkoreceptors and to also address receptor duality. The study of intra- and inter-organism swodkocrine signaling and its crosstalk with gasocrine signaling may help to understand the etiology of diseases due to dysregulation in sugar homeostasis and signaling.

CHiTA: A scarless high-throughput pipeline for characterization of ribozymes

Small self-cleaving ribozymes are catalytic RNAs that cleave their phosphodiester backbone rapidly and site-specifically, without the assistance of proteins. Their catalytic properties make them ideal targets for applications in RNA pharmaceuticals and bioengineering. Consequently, computational pipelines that predict or design thousands of self-cleaving ribozyme candidates have been developed. Traditional experimental techniques for verifying the activity of these putative ribozymes, however, are low-throughput and time intensive. High-throughput (HT) pipelines that employ next-generation sequencing (NGS) analyze the activity of these thousands of ribozymes simultaneously. Until recently, the application of these HT pipelines has been limited to studying all single and double mutants of a select representative ribozyme. Unfortunately, this prevents the exploration of candidates having different lengths, circular permutations, and auxiliary stem-loops. Moreover, pipelines that analyze ribozymes en masse often include transcription of non-native flanking sequences that preclude accurate assessment of the intrinsic rate of ribozyme self-cleavage. To overcome these limitations, we developed a HT pipeline, "Cleavage High-Throughput Assay (CHiTA)", which employs NGS and massively parallel oligonucleotide synthesis (MPOS) to characterize ribozyme activity for thousands of candidates in a scarless fashion. Herein, we describe detailed strategies and protocols to implement CHiTA to measure the activity of putative ribozymes from a wide range of ribozyme classes.

Structure and catalytic activity of the SAM-utilizing ribozyme SAMURI

Ribozymes that catalyze site-specific RNA modification have recently gained increasing interest for their ability to mimic methyltransferase enzymes and for their application to install molecular tags. Recently, we reported SAMURI as a site-specific alkyltransferase ribozyme using S-adenosylmethionine (SAM) or a stabilized analog to transfer a methyl or propargyl group to N3 of an adenosine. Here, we report the crystal structures of SAMURI in the postcatalytic state. The structures reveal a three-helix junction with the catalytic core folded into four stacked layers, harboring the cofactor and the modified nucleotide. Detailed structure-activity analyses explain the cofactor scope and the structural basis for site selectivity. A structural comparison of SAMURI with SAM riboswitches sheds light on how the synthetic ribozyme overcomes the strategies of natural riboswitches to avoid self-methylation. Our results suggest that SAM and its analogs may serve as substrates for various RNA-catalyzed reactions, for which the corresponding ribozymes remain to be identified.

GlmS plays a key role in the virulence factor expression and biofilm formation ability of Staphylococcus aureus promoted by advanced glycation end products

Staphylococcus aureus (S. aureus) is well known for its biofilm formation ability and is responsible for serious, chronic refractory infections worldwide. We previously demonstrated that advanced glycation end products (AGEs), a hallmark of chronic hyperglycaemia in diabetic tissues, enhanced biofilm formation by promoting eDNA release via sigB upregulation in S. aureus, contributing to the high morbidity and mortality of patients presenting a diabetic foot ulcer infection. However, the exact regulatory network has not been completely described. Here, we used pull-down assay and LC-MS/MS to identify the GlmS as a candidate regulator of sigB in S. aureus stimulated by AGEs. Dual-luciferase assays and electrophoretic mobility shift assays (EMSAs) revealed that GlmS directly upregulated the transcriptional activity of sigB. We constructed NCTC 8325 ∆glmS for further validation. qRT-PCR analysis revealed that AGEs promoted both glmS and sigB expression in the NCTC 8325 strain but had no effect on NCTC 8325 ∆glmS. NCTC 8325 ∆glmS showed a significant attenuation in biofilm formation and virulence factor expression, accompanied by a decrease in sigB expression, even under AGE stimulation. All of the changes, including pigment deficiency, decreased haemolysis ability, downregulation of hla and hld expression, and less and sparser biofilms, indicated that sigB and biofilm formation ability no longer responded to AGEs in NCTC 8325 ∆glmS. Our data extend the understanding of GlmS in the global regulatory network of S. aureus and demonstrate a new mechanism by which AGEs can upregulate GlmS, which directly regulates sigB and plays a significant role in mediating biofilm formation and virulence factor expression.

Optimized periphery-core interface increases fitness of the Bacillus subtilis glmS ribozyme

Like other functional RNAs, ribozymes encode a conserved catalytic center supported by peripheral domains that vary among ribozyme sub-families. To understand how core-periphery interactions contribute to ribozyme fitness, we compared the cleavage kinetics of all single base substitutions at 152 sites across the Bacillus subtilis glmS ribozyme by high-throughput sequencing (k-seq). The in vitro activity map mirrored phylogenetic sequence conservation in glmS ribozymes, indicating that biological fitness reports all biochemically important positions. The k-seq results and folding assays showed that most deleterious mutations lower activity by impairing ribozyme self-assembly. All-atom molecular dynamics simulations of the complete ribozyme revealed how individual mutations in the core or the IL4 peripheral loop introduce a non-native tertiary interface that rewires the catalytic center, eliminating activity. We conclude that the need to avoid non-native helix packing powerfully constrains the evolution of tertiary structure motifs in RNA.

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

Antisense RNA C9orf72 hexanucleotide repeat associated with amyotrophic lateral sclerosis and frontotemporal dementia forms a triplex-like structure and binds small synthetic ligand

The abnormal expansion of GGGGCC/GGCCCC hexanucleotide repeats (HR) in C9orf72 is associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Structural polymorphisms of HR result in the multifactorial pathomechanism of ALS/FTD. Consequently, many ongoing studies are focused at developing therapies targeting pathogenic HR RNA. One of them involves small molecules blocking sequestration of important proteins, preventing formation of toxic nuclear foci. However, rational design of potential therapeutics is hindered by limited number of structural studies of RNA-ligand complexes. We determined the crystal structure of antisense HR RNA in complex with ANP77 ligand (1.1 Å resolution) and in the free form (0.92 and 1.5 Å resolution). HR RNA folds into a triplex structure composed of four RNA chains. ANP77 interacted with two neighboring single-stranded cytosines to form pseudo-canonical base pairs by adopting sandwich-like conformation and adjusting the position of its naphthyridine units to the helical twist of the RNA. In the unliganded structure, the cytosines formed a peculiar triplex i-motif, assembled by trans C•C+ pair and a third cytosine located at the Hoogsteen edge of the C•C+ pair. These results extend our knowledge of the structural polymorphisms of HR and can be used for rational design of small molecules targeting disease-related RNAs.

Flipping the script: Understanding riboswitches from an alternative perspective

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

RNA-based controllers for engineering gene and cell therapies

Engineered RNA-based genetic controllers provide compact, tunable, post-transcriptional gene regulation. As RNA devices are generally small, these devices are portable to DNA and RNA viral vectors. RNA tools have recently expanded to allow reading and editing of endogenous RNAs for profiling and programming of transcriptional states. With their expanded capabilities and highly compact, modular, and programmable nature, RNA-based controllers will support greater safety, efficacy, and performance in gene and cell-based therapies. In this review, we highlight RNA-based controllers and their potential as user-guided and autonomous systems for control of gene and cell-based therapies.

Co-transcriptional folding of the glmS ribozyme enables a rapid response to metabolite

The glmS ribozyme riboswitch, located in the 5' untranslated region of the Bacillus subtilis glmS messenger RNA (mRNA), regulates cell wall biosynthesis through ligand-induced self-cleavage and decay of the glmS mRNA. Although self-cleavage of the refolded glmS ribozyme has been studied extensively, it is not known how early the ribozyme folds and self-cleaves during transcription. Here, we combine single-molecule fluorescence with kinetic modeling to show that self-cleavage can occur during transcription before the ribozyme is fully synthesized. Moreover, co-transcriptional folding of the RNA at a physiological elongation rate allows the ribozyme catalytic core to react without the downstream peripheral stability domain. Dimethyl sulfate footprinting further revealed how slow sequential folding favors formation of the native core structure through fraying of misfolded helices and nucleation of a native pseudoknot. Ribozyme self-cleavage at an early stage of transcription may benefit glmS regulation in B. subtilis, as it exposes the mRNA to exoribonuclease before translation of the open reading frame can begin. Our results emphasize the importance of co-transcriptional folding of RNA tertiary structure for cis-regulation of mRNA stability.

Gas-sensing riboceptors

Understanding how cells sense gases or gaseous solutes is a fundamental question in biology and is pivotal for the evolution of molecular and organismal life. In numerous organisms, gases can diffuse into cells, be transported, generated, and sensed. Controlling gases in the cellular environment is essential to prevent cellular and molecular damage due to interactions with gas-dependent free radicals. Consequently, the mechanisms governing acute gas sensing are evolutionarily conserved and have been experimentally elucidated in various organisms. However, the scientific literature on direct gas sensing is largely based on hemoprotein-based gasoreceptors (or sensors). As RNA-based G-quadruplex (G4) structures can also bind to heme, I propose that some ribozymes can act as gas-sensing riboceptors (ribonucleic acid receptors). Additionally, I present a few other ideas for non-heme metal ion- or metal cluster-based gas-sensing riboceptors. Studying riboceptors can help understand the evolutionary origins of cellular and gasocrine signaling.

Computational design and experimental verification of pseudoknotted ribozymes

The design of new RNA sequences that retain the function of a model RNA structure is a challenge in bioinformatics because of the structural complexity of these molecules. RNA can fold into its secondary and tertiary structures by forming stem-loops and pseudoknots. A pseudoknot is a set of base pairs between a region within a stem-loop and nucleotides outside of this stem-loop; this motif is very important for numerous functional structures. It is important for any computational design algorithm to take into account these interactions to give a reliable result for any structures that include pseudoknots. In our study, we experimentally validated synthetic ribozymes designed by Enzymer, which implements algorithms allowing for the design of pseudoknots. Enzymer is a program that uses an inverse folding approach to design pseudoknotted RNAs; we used it in this study to design two types of ribozymes. The ribozymes tested were the hammerhead and the glmS, which have a self-cleaving activity that allows them to liberate the new RNA genome copy during rolling-circle replication or to control the expression of the downstream genes, respectively. We demonstrated the efficiency of Enzymer by showing that the pseudoknotted hammerhead and glmS ribozymes sequences it designed were extensively modified compared to wild-type sequences and were still active.

Antagonistic transcriptome profile reveals potential mechanisms of action on Xanthomonas oryzae pv. oryzicola by the cell-free supernatants of Bacillus velezensis 504, a versatile plant probiotic bacterium

Bacterial leaf streak (BLS) of rice is a severe disease caused by the bacterial pathogen Xanthomonas oryzae pv. oryzicola (Xoc) that has gradually become the fourth major disease on rice in some rice-growing regions in southern China. Previously, we isolated a Bacillus velezensis strain 504 that exhibited apparent antagonistic activity against the Xoc wild-type strain RS105, and found that B. velezensis 504 was a potential biocontrol agent for BLS. However, the underlying mechanisms of antagonism and biocontrol are not completely understood. Here we mine the genomic data of B. velezensis 504, and the comparative transcriptomic data of Xoc RS105 treated by the cell-free supernatants (CFSs) of B. velezensis 504 to define differentially expressed genes (DEGs). We show that B. velezensis 504 shares over 89% conserved genes with FZB42 and SQR9, two representative model strains of B. velezensis, but 504 is more closely related to FZB42 than SQR9, as well as B. velezensis 504 possesses the secondary metabolite gene clusters encoding the essential anti-Xoc agents difficidin and bacilysin. We conclude that approximately 77% of Xoc RS105 coding sequences are differentially expressed by the CFSs of B. velezensis 504, which significantly downregulates genes involved in signal transduction, oxidative phosphorylation, transmembrane transport, cell motility, cell division, DNA translation, and five physiological metabolisms, as well as depresses an additional set of virulence-associated genes encoding the type III secretion, type II secretion system, type VI secretion system, type IV pilus, lipopolysaccharides and exopolysaccharides. We also show that B. velezensis 504 is a potential biocontrol agent for bacterial blight of rice exhibiting relative control efficiencies over 70% on two susceptible cultivars, and can efficiently antagonize against some important plant pathogenic fungi including Colletotrichum siamense and C. australisinense that are thought to be the two dominant pathogenic species causing leaf anthracnose of rubber tree in Hainan province of China. B. velezensis 504 also harbors some characteristics of plant growth-promoting rhizobacterium such as secreting protease and siderophore, and stimulating plant growth. This study reveals the potential biocontrol mechanisms of B. velezensis against BLS, and also suggests that B. velezensis 504 is a versatile plant probiotic bacterium.

Real-Time Assessment of Intracellular Metabolites in Single Cells through RNA-Based Sensors

Quantification of the concentration of particular cellular metabolites reports on the actual utilization of metabolic pathways in physiological and pathological conditions. Metabolite concentration also constitutes the readout for screening cell factories in metabolic engineering. However, there are no direct approaches that allow for real-time assessment of the levels of intracellular metabolites in single cells. In recent years, the modular architecture of natural bacterial RNA riboswitches has inspired the design of genetically encoded synthetic RNA devices that convert the intracellular concentration of a metabolite into a quantitative fluorescent signal. These so-called RNA-based sensors are composed of a metabolite-binding RNA aptamer as the sensor domain, connected through an actuator segment to a signal-generating reporter domain. However, at present, the variety of available RNA-based sensors for intracellular metabolites is still very limited. Here, we go through natural mechanisms for metabolite sensing and regulation in cells across all kingdoms, focusing on those mediated by riboswitches. We review the design principles underlying currently developed RNA-based sensors and discuss the challenges that hindered the development of novel sensors and recent strategies to address them. We finish by introducing the current and potential applicability of synthetic RNA-based sensors for intracellular metabolites.

CRISPR-dCas12a-mediated genetic circuit cascades for multiplexed pathway optimization

The production efficiency of microbial cell factories is sometimes limited by the lack of effective methods to regulate multiple targets in a coordinated manner. Here taking the biosynthesis of glucosamine-6-phosphate (GlcN6P) in Bacillus subtilis as an example, a 'design-build-test-learn' framework was proposed to achieve efficient multiplexed optimization of metabolic pathways. A platform strain was built to carry biosensor signal-amplifying circuits and two genetic regulation circuits. Then, a synthetic CRISPR RNA array blend for boosting and leading (ScrABBLE) device was integrated into the platform strain, which generated 5,184 combinatorial assemblies targeting three genes. The best GlcN6P producer was screened and engineered for the synthesis of valuable pharmaceuticals N-acetylglucosamine and N-acetylmannosamine. The N-acetylglucosamine titer reached 183.9 g liter-1 in a 15-liter bioreactor. In addition, the potential generic application of the ScrABBLE device was also verified using three fluorescent proteins as a case study.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

Probing RNA Structures and Interactions Using Fluorescence Lifetime Analyses

Structural analyses of large, complex noncoding RNAs continue to lag behind their rapid discovery and functional descriptions. Site-specifically incorporated, minimally invasive fluorescent probes such as 2-aminopurine (2AP) and pyrrolo-cytosine (PyC) have provided essential complementary information about local RNA structure, conformational dynamics, and interactions. Here I describe a protocol that benchmarks and correlates local RNA conformations with their respective fluorescence lifetimes, as a general technique that confers key advantages over fluorescence intensity-based methods. The observation that fluorescence lifetimes are more sensitive to local structures than sequence contexts suggests broad utility across diverse RNA and ribonucleoprotein systems.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

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

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

A Riboswitch-Driven Era of New Antibacterials

Riboswitches are structured non-coding RNAs found in the 5′ UTR of important genes for bacterial metabolism, virulence and survival. Upon the binding of specific ligands that can vary from simple ions to complex molecules such as nucleotides and tRNAs, riboswitches change their local and global mRNA conformations to affect downstream transcription or translation. Due to their dynamic nature and central regulatory role in bacterial metabolism, riboswitches have been exploited as novel RNA-based targets for the development of new generation antibacterials that can overcome drug-resistance problems. During recent years, several important riboswitch structures from many bacterial representatives, including several prominent human pathogens, have shown that riboswitches are ideal RNA targets for new compounds that can interfere with their structure and function, exhibiting much reduced resistance over time. Most interestingly, mainstream antibiotics that target the ribosome have been shown to effectively modulate the regulatory behavior and capacity of several riboswitches, both in vivo and in vitro, emphasizing the need for more in-depth studies and biological evaluation of new antibiotics. Herein, we summarize the currently known compounds that target several main riboswitches and discuss the role of mainstream antibiotics as modulators of T-box riboswitches, in the dawn of an era of novel inhibitors that target important bacterial regulatory RNAs.

Key players in regulatory RNA realm of bacteria

Precise regulation of gene expression is crucial for living cells to adapt for survival in diverse environmental conditions. Among the common cellular regulatory mechanisms, RNA-based regulators play a key role in all domains of life. Discovery of regulatory RNAs have made a paradigm shift in molecular biology as many regulatory functions of RNA have been identified beyond its canonical roles as messenger, ribosomal and transfer RNA. In the complex regulatory RNA network, riboswitches, small RNAs, and RNA thermometers can be identified as some of the key players. Herein, we review the discovery, mechanism, and potential therapeutic use of these classes of regulatory RNAs mainly found in bacteria. Being highly adaptive organisms that inhabit a broad range of ecological niches, bacteria have adopted tight and rapid-responding gene regulation mechanisms. This review aims to highlight how bacteria utilize versatile RNA structures and sequences to build a sophisticated gene regulation network.

•

The three major classes of prokaryotic ncRNAs and their characterized mechanisms of operation in gene regulation.

•

sRNAs emerging as major players in global gene regulatory networks.

•

Riboswitch mediated gene control mechanisms through on/off switches in response to ligand binding.

•

RNA thermo sensors for temperature-dependent gene expression.

•

Therapeutic importance of ncRNAs and computational approaches involved in the discovery of ncRNAs.

Overview of Methods for Large-Scale RNA Synthesis

In recent years, it has become clear that RNA molecules are involved in almost all vital cellular processes and pathogenesis of human disorders. The functional diversity of RNA comes from its structural richness. Although composed of only four nucleotides, RNA molecules present a plethora of secondary and tertiary structures critical for intra and intermolecular contacts with other RNAs and ligands (proteins, small metabolites, etc.). In order to fully understand RNA function it is necessary to define its spatial structure. Crystallography, nuclear magnetic resonance and cryogenic electron microscopy have demonstrated considerable success in determining the structures of biologically important RNA molecules. However, these powerful methods require large amounts of sample. Despite their limitations, chemical synthesis and in vitro transcription are usually employed to obtain milligram quantities of RNA for structural studies, delivering simple and effective methods for large-scale production of homogenous samples. The aim of this paper is to provide an overview of methods for large-scale RNA synthesis with emphasis on chemical synthesis and in vitro transcription. We also present our own results of testing the efficiency of these approaches in order to adapt the material acquisition strategy depending on the desired RNA construct.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

Architectures and complex functions of tandem riboswitches

Riboswitch architectures that involve the binding of a single ligand to a single RNA aptamer domain result in ordinary dose-response curves that require approximately a 100-fold change in ligand concentration to cover nearly the full dynamic range for gene regulation. However, by using multiple riboswitches or aptamer domains in tandem, these ligand-sensing structures can produce additional, complex gene control outcomes. In the current study, we have computationally searched for tandem riboswitch architectures in bacteria to provide a more complete understanding of the diverse biological and biochemical functions of gene control elements that are made exclusively of RNA. Numerous different arrangements of tandem homologous riboswitch architectures are exploited by bacteria to create more 'digital' gene control devices, which operate over a narrower ligand concentration range. Also, two heterologous riboswitch aptamers are sometimes employed to create two-input Boolean logic gates with various types of genetic outputs. These findings illustrate the sophisticated genetic decisions that can be made by using molecular sensors and switches based only on RNA.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

Riboswitch Mechanisms: New Tricks for an Old Dog

Discovered almost twenty years ago, riboswitches turned out to be one of the most common regulatory systems in bacteria, with representatives found in eukaryotes and archaea. Unlike many other regulatory elements, riboswitches are entirely composed of RNA and capable of modulating expression of genes by direct binding of small cellular molecules. While bacterial riboswitches had been initially thought to control production of enzymes and transporters associated with small organic molecules via feedback regulatory circuits, later findings identified riboswitches directing expression of a wide range of genes and responding to various classes of molecules, including ions, signaling molecules, and others. The 5'-untranslated mRNA regions host a vast majority of riboswitches, which modulate transcription or translation of downstream genes through conformational rearrangements in the ligand-sensing domains and adjacent expression-controlling platforms. Over years, the repertoire of regulatory mechanisms employed by riboswitches has greatly expanded; most recent studies have highlighted the importance of alternative mechanisms, such as RNA degradation, for the riboswitch-mediated genetic circuits. This review discusses the plethora of bacterial riboswitch mechanisms and illustrates how riboswitches utilize different features and approaches to elicit various regulatory responses.

The potential versatility of RNA catalysis

It is commonly thought that in the early development of life on this planet RNA would have acted both as a store of genetic information and as a catalyst. While a number of RNA enzymes are known in contemporary cells, they are largely confined to phosphoryl transfer reactions, whereas an RNA based metabolism would have required a much greater chemical diversity of catalysis. Here we discuss how RNA might catalyze a wider variety of chemistries, and particularly how information gleaned from riboswitches could suggest how ribozymes might recruit coenzymes to expand their chemical range. We ask how we might seek such activities in modern biology. This article is categorized under: RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions Regulatory RNAs/RNAi/Riboswitches > Riboswitches RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Peptidoglycan Salvage Enables the Periodontal Pathogen Tannerella forsythia to Survive within the Oral Microbial Community

Tannerella forsythia is an anaerobic, fusiform Gram-negative oral pathogen strongly associated with periodontitis, a multibacterial inflammatory disease that leads to the destruction of the teeth-supporting tissue, ultimately causing tooth loss. To survive in the oral habitat, T. forsythia depends on cohabiting bacteria for the provision of nutrients. For axenic growth under laboratory conditions, it specifically relies on the external supply of N-acetylmuramic acid (MurNAc), which is an essential constituent of the peptidoglycan (PGN) of bacterial cell walls. T. forsythia comprises a typical Gram-negative PGN; however, as evidenced by genome sequence analysis, the organism lacks common enzymes required for the de novo synthesis of precursors of PGN, which rationalizes its MurNAc auxotrophy. Only recently insights were obtained into how T. forsythia gains access to MurNAc in its oral habitat, enabling synthesis of the own PGN cell wall. This report summarizes T. forsythia's strategies to survive in the oral habitat by means of PGN salvage pathways, including recovery of exogenous MurNAc and PGN-derived fragments but also polymeric PGN, which are all derived from cohabiting bacteria either via cell wall turnover or decay of cells. Salvage of polymeric PGN presumably requires the removal of peptides from PGN by an unknown amidase, concomitantly with the translocation of the polymer across the outer membrane. Two recently identified exo-lytic N-acetylmuramidases (Tf_NamZ1 and Tf_NamZ2) specifically cleave the peptide-free, exogenous (nutrition source) PGN in the periplasm and release the MurNAc and disaccharide substrates for the transporters Tf_MurT and Tf_AmpG, respectively, whereas the peptide-containing, endogenous (the self-cell wall) PGN stays unattached. This review also outlines how T. forsythia synthesises the PGN precursors UDP-MurNAc and UDP-N-acetylglucosamine (UDP-GlcNAc), involving homologs of the Pseudomonas sp. recycling enzymes AmgK/MurU and a monofunctional uridylyl transferase (named Tf_GlmU*), respectively.

The case of the missing allosteric ribozymes

The RNA World theory encompasses the hypothesis that sophisticated ribozymes and riboswitches were the primary drivers of metabolic processes in ancient organisms. Several types of catalytic RNAs and many classes of ligand-sensing RNA switches still exist in modern cells. Curiously, allosteric ribozymes formed by the merger of RNA enzyme and RNA switch components are largely absent in today's biological systems. This is true despite the striking abundances of various classes of both self-cleaving ribozymes and riboswitch aptamers. Here we present the known types of ligand-controlled ribozymes and riboswitches and discuss the possible reasons why fused ribozyme-aptamer constructs have been disfavored through evolution.

Synthetic biology-driven microbial production of folates: Advances and perspectives

With the development and application of synthetic biology, significant progress has been made in the production of folate by microbial fermentation using cell factories, especially for using generally regarded as safe (GRAS) microorganism as production host. In this review, the physiological functions and applications of folates were firstly discussed. Second, the current advances of folate-producing GRAS strains development were summarized. Third, the applications of synthetic biology-based metabolic regulatory tools in GRAS strains were introduced, and the progress in the application of these tools for folate production were summarized. Finally, the challenges to folates efficient production and corresponding emerging strategies to overcome them by synthetic biology were discussed, including the construction of biosensors using tetrahydrofolate riboswitches to regulate metabolic pathways, adaptive evolution to overcome the flux limitations of the folate pathway. The combination of new strategies and tools of synthetic biology is expected to further improve the efficiency of microbial folate synthesis.

Riboswitches as Drug Targets for Antibiotics

Riboswitches reside in the untranslated region of RNA and regulate genes involved in the biosynthesis of essential metabolites through binding of small molecules. Since their discovery at the beginning of this century, riboswitches have been regarded as potential antibacterial targets. Using fragment screening, high-throughput screening and rational ligand design guided by X-ray crystallography, lead compounds against various riboswitches have been identified. Here, we review the current status and suitability of the thiamine pyrophosphate (TPP), flavin mononucleotide (FMN), glmS, guanine, and other riboswitches as antibacterial targets and discuss them in a biological context. Further, we highlight challenges in riboswitch drug discovery and emphasis the need to develop riboswitch specific high-throughput screening methods.

Improving RNA Crystal Diffraction Quality by Postcrystallization Treatment

The crystallization and structural determination of large RNAs and their complexes remain major bottlenecks in the mechanistic analysis of cellular and viral RNAs. Here, we describe a protocol that combines postcrystallization dehydration and ion replacement that dramatically improved the diffraction quality of crystals of a large gene-regulatory tRNA-mRNA complex. Through this method, the resolution limit of X-ray data extended from 8.5 to 3.2 Å, enabling structure determination. Although this protocol was developed for a particular RNA complex, the general importance of solvent and counterions in nucleic acid structure may render it generally useful for crystallographic analysis of other RNAs.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

Design of a programmable biosensor-CRISPRi genetic circuits for dynamic and autonomous dual-control of metabolic flux in Bacillus subtilis

Dynamic regulation is an effective strategy for fine-tuning metabolic pathways in order to maximize target product synthesis. However, achieving dynamic and autonomous up- and down-regulation of the metabolic modules of interest simultaneously, still remains a great challenge. In this work, we created an autonomous dual-control (ADC) system, by combining CRISPRi-based NOT gates with novel biosensors of a key metabolite in the pathway of interest. By sensing the levels of the intermediate glucosamine-6-phosphate (GlcN6P) and self-adjusting the expression levels of the target genes accordingly with the GlcN6P biosensor and ADC system enabled feedback circuits, the metabolic flux towards the production of the high value nutraceutical N-acetylglucosamine (GlcNAc) could be balanced and optimized in Bacillus subtilis. As a result, the GlcNAc titer in a 15-l fed-batch bioreactor increased from 59.9 g/l to 97.1 g/l with acetoin production and 81.7 g/l to 131.6 g/l without acetoin production, indicating the robustness and stability of the synthetic circuits in a large bioreactor system. Remarkably, this self-regulatory methodology does not require any external level of control such as the use of inducer molecules or switching fermentation/environmental conditions. Moreover, the proposed programmable genetic circuits may be expanded to engineer other microbial cells and metabolic pathways.

Computational study of switching mechanism in add A-riboswitch

Riboswitch can bind small molecules to regulate gene expression. Unlike other RNAs, riboswitch relies on its conformational switching for regulation. However, the understanding of the switching mechanism is still limited. Here, we focussed on the add A-riboswitch to illustrate the dynamical switching mechanism as an example. We performed molecular dynamics simulation, conservation and co-evolution calculations to infer the dynamical motions and evolutionary base pairings. The results suggest that the binding domain is stable for molecule recognition and binding, whereas the switching base pairings are co-evolutionary for translation. The understanding of the add A-riboswitch switching mechanism provides a potential solution for riboswitch drug design.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Classification of the nucleolytic ribozymes based upon catalytic mechanism

The nucleolytic ribozymes carry out site-specific RNA cleavage reactions by nucleophilic attack of the 2'-oxygen atom on the adjacent phosphorus with an acceleration of a million-fold or greater. A major part of this arises from concerted general acid-base catalysis. Recent identification of new ribozymes has expanded the group to a total of nine and this provides a new opportunity to identify sub-groupings according to the nature of the general base and acid. These include nucleobases, hydrated metal ions, and 2'-hydroxyl groups. Evolution has selected a number of different combinations of these elements that lead to efficient catalysis. These differences provide a new mechanistic basis for classifying these ribozymes.

Cleaning Up Mechanistic Debris Generated by Twister Ribozymes Using Computational RNA Enzymology

The catalytic properties of RNA have been a subject of fascination and intense research since their discovery over 30 years ago. Very recently, several classes of nucleolytic ribozymes have emerged and been characterized structurally. Among these, the twister ribozyme has been center-stage, and a topic of debate about its architecture and mechanism owing to conflicting interpretations of different crystal structures, and in some cases conflicting interpretations of the same functional data. In the present work, we attempt to clean up the mechanistic "debris" generated by twister ribozymes using a comprehensive computational RNA enzymology approach aimed to provide a unified interpretation of existing structural and functional data. Simulations in the crystalline environment and in solution provide insight into the origins of observed differences in crystal structures, and coalesce on a common active site architecture, and dynamical ensemble in solution. We use GPU-accelerated free energy methods with enhanced sampling to ascertain microscopic nucleobase pK a values of the implicated general acid and base, from which predicted activity-pH profiles can be compared directly with experiments. Next, ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with full dynamic solvation under periodic boundary conditions are used to determine mechanistic pathways through multi-dimensional free energy landscapes for the reaction. We then characterize the rate-controlling transition state, and make predictions about kinetic isotope effects and linear free energy relations. Computational mutagenesis is performed to explain the origin of rate effects caused by chemical modifications and make experimentally testable predictions. Finally, we provide evidence that helps to resolve conflicting issues related to the role of metal ions in catalysis. Throughout each stage, we highlight how a conserved L-platform structural motif, to- gether with a key L-anchor residue, forms the characteristic active site scaffold enabling each of the catalytic strategies to come together not only for the twister ribozyme, but the majority of the known small nucleolytic ribozyme classes.

Hatchet ribozyme structure and implications for cleavage mechanism

Self-cleaving ribozymes are RNAs that catalyze position-specific cleavage of their phosphodiester backbone. The cleavage site of the newly discovered hatchet ribozyme is located at the very 5′ end of its consensus secondary structure motif. Here we report on the 2.1-Å crystal structure of the hatchet ribozyme in the product state, which defines its intricate tertiary fold and identifies key residues lining the catalytic pocket. This in turn has allowed us to propose a model of the precatalytic state structure and a role in catalysis for a conserved guanine. This study therefore provides a structure-based platform toward an improved understanding of the catalytic mechanism of hatchet ribozymes.

Small self-cleaving ribozymes catalyze site-specific cleavage of their own phosphodiester backbone with implications for viral genome replication, pre-mRNA processing, and alternative splicing. We report on the 2.1-Å crystal structure of the hatchet ribozyme product, which adopts a compact pseudosymmetric dimeric scaffold, with each monomer stabilized by long-range interactions involving highly conserved nucleotides brought into close proximity of the scissile phosphate. Strikingly, the catalytic pocket contains a cavity capable of accommodating both the modeled scissile phosphate and its flanking 5′ nucleoside. The resulting modeled precatalytic conformation incorporates a splayed-apart alignment at the scissile phosphate, thereby providing structure-based insights into the in-line cleavage mechanism. We identify a guanine lining the catalytic pocket positioned to contribute to cleavage chemistry. The functional relevance of structure-based insights into hatchet ribozyme catalysis is strongly supported by cleavage assays monitoring the impact of selected nucleobase and atom-specific mutations on ribozyme activity.

Detection of Low-Abundance Metabolites in Live Cells Using an RNA Integrator

Genetically encoded biosensors are useful tools for detecting the presence and levels of diverse biomolecules in living cells. However, low-abundance targets are difficult to detect because they are often unable to bind and activate enough biosensors to detect using standard microscopic imaging approaches. Here we describe a type of RNA-based biosensor, an RNA integrator, which enables detection of low-abundance targets in vitro and in living cells. The RNA integrator is an RNA sequence comprising a ribozyme and an unfolded form of the fluorogenic aptamer Broccoli. Upon binding its target, the ribozyme undergoes cleavage and releases Broccoli, which subsequently folds and becomes fluorescent. Importantly, each target molecule can bind and induce cleavage of multiple copies of the integrator sensor, resulting in an amplified signal. We show that this approach can be generalized to numerous different ribozyme types for the detection of various small molecules.