Adaptive ligand binding by the purine riboswitch in the recognition of guanine and adenine analogs

Advances and Mechanisms of RNA-Ligand Interaction Predictions

The diversity and complexity of RNA include sequence, secondary structure, and tertiary structure characteristics. These elements are crucial for RNA's specific recognition of other molecules. With advancements in biotechnology, RNA-ligand structures allow researchers to utilize experimental data to uncover the mechanisms of complex interactions. However, determining the structures of these complexes experimentally can be technically challenging and often results in low-resolution data. Many machine learning computational approaches have recently emerged to learn multiscale-level RNA features to predict the interactions. Predicting interactions remains an unexplored area. Therefore, studying RNA-ligand interactions is essential for understanding biological processes. In this review, we analyze the interaction characteristics of RNA-ligand complexes by examining RNA's sequence, secondary structure, and tertiary structure. Our goal is to clarify how RNA specifically recognizes ligands. Additionally, we systematically discuss advancements in computational methods for predicting interactions and to guide future research directions. We aim to inspire the creation of more reliable RNA-ligand interaction prediction tools.

Divergent Metabolic Fates of Aromatic Amino Acid-Derived Isomers: Insights from Ex Vivo Metabolomics and HDX-HRMS/MS-Based Resolution of Tautomers

Tautomers are one of the many types of isomers, and differences in tautomeric structures confer altered chemical and biological properties. Using ultrahigh-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS) ex vivo metabolomics, we investigate, in whole blood, the divergent metabolism of enol and keto forms of indole-3-pyruvate (IPyA), a tautomeric product of aromatic amino acid metabolism. Two new compounds resulting from IPyA metabolism were discovered, 3-(1H-indol-3-yl)-2,3-dioxopropanoic acid or "indole-3-oxopyruvic acid" and glutathionyl-indole pyruvate (GSHIPyA), which were characterized via ultraviolet photodissociation (UVPD) and higher-energy collisional dissociation (HCD). Computational calculations support the hypothesis that GSHIPyA forms specifically through the enol form of IPyA. GSHIPyA is also hypothesized to be tautomeric, and a hydrogen-deuterium exchange-high-resolution tandem mass spectrometry (HDX-HRMS/MS) approach is developed to prove the presence of an enol and keto tautomer. HDX of GSHIPyA labels the keto form with an additional deuterium, relative to the enol form. HRMS/MS of the labeled isomers is employed to leverage the relationship of resolving power scaling inversely with the square root of m/z, for Orbitrap mass analyzers. HRMS/MS yields a smaller-molecular-weight deuterated tautomeric product ion, reducing the analyte ion m/z and thus lowering the resolving power necessary to separate the deuterated keto tautomer product ion from the [13]C product ion.

Thermodynamic compensation to temperature extremes in B. subtilis vs T. maritima lysine riboswitches

T. maritima and B. subtilis are bacteria that inhabit significantly different thermal environments, ∼80 vs. ∼40°C, yet employ similar lysine riboswitches to aid in the transcriptional regulation of the genes involved in the synthesis and transport of amino acids. Despite notable differences in G-C basepair frequency and primary sequence, the aptamer moieties of each riboswitch have striking similarities in tertiary structure, with several conserved motifs and long-range interactions. To explore genetic adaptation in extreme thermal environments, we compare the kinetic and thermodynamic behaviors in T. maritima and B. subtilis lysine riboswitches via single-molecule fluorescence resonance energy transfer analysis. Kinetic studies reveal that riboswitch folding rates increase with lysine concentration while the unfolding rates are independent of lysine. This indicates that both riboswitches bind lysine through an induced-fit ("bind-then-fold") mechanism, with lysine binding necessarily preceding conformational changes. Temperature-dependent van't Hoff studies reveal qualitative similarities in the thermodynamic landscapes for both riboswitches in which progression from the open, lysine-unbound state to both transition states (‡) and closed, lysine-bound conformations is enthalpically favored yet entropically penalized, with comparisons of enthalpic and entropic contributions extrapolated to a common [K+] = 100 mM in quantitative agreement. Finally, temperature-dependent Eyring analysis reveals the TMA and BSU riboswitches to have remarkably similar folding/unfolding rate constants when extrapolated to their respective (40 and 80°C) environmental temperatures. Such behavior suggests a shared strategy for ligand binding and aptamer conformational change in the two riboswitches, based on thermodynamic adaptations in number of G-C basepairs and/or modifications in tertiary structure that stabilize the ligand-unbound conformation to achieve biocompetence under both hyperthermophilic and mesothermophilic conditions.

High-content tailoring strategy to improve the multifunctionality of functional nucleic acids

Functional nucleic acids (FNAs) have attracted increasing attention in recent years due to their diverse physiological functions. The understanding of their conformational recognition mechanisms has advanced through nucleic acid tailoring strategies and sequence optimization. With the development of the FNA tailoring techniques, they have become a methodological guide for nucleic acid repurposing. Therefore, it is necessary to systematize the relationship between FNA tailoring strategies and the development of nucleic acid multifunctionality. This review systematically categorizes eight types of FNA multifunctionality, and introduces the traditional FNA tailoring strategy from five aspects, including deletion, substitution, splitting, fusion and elongation. Based on the current state of FNA modification, a new generation of FNA tailoring strategy, called the high-content tailoring strategy, was unprecedentedly proposed to improve FNA multifunctionality. In addition, the multiple applications of rational tailoring-driven FNA performance enhancement in various fields were comprehensively summarized. The limitations and potential of FNA tailoring and repurposing in the future are also explored in this review. In summary, this review introduces a novel tailoring theory, systematically summarizes eight FNA performance enhancements, and provides a systematic overview of tailoring applications across all categories of FNAs. The high-content tailoring strategy is expected to expand the application scenarios of FNAs in biosensing, biomedicine and materials science, thus promoting the synergistic development of various fields.

Multidimensional Applications and Challenges of Riboswitches in Biosensing and Biotherapy

Riboswitches have received significant attention over the last two decades for their multiple functionalities and great potential for applications in various fields. This article highlights and reviews the recent advances in biosensing and biotherapy. These fields involve a wide range of applications, such as food safety detection, environmental monitoring, metabolic engineering, live cell imaging, wearable biosensors, antibacterial drug targets, and gene therapy. The discovery, origin, and optimization of riboswitches are summarized to help readers better understand their multidimensional applications. Finally, this review discusses the multidimensional challenges and development of riboswitches in order to further expand their potential for novel applications.

Mechanism of Ligand Discrimination by the NMT1 Riboswitch

Riboswitches are conserved functional domains in mRNA that almost exclusively exist in bacteria. They regulate the biosynthesis and transport of amino acids and essential metabolites such as coenzymes, nucleobases, and their derivatives by specifically binding small molecules. Due to their ability to precisely discriminate between different cognate molecules as well as their common existence in bacteria, riboswitches have become potential antibacterial drug targets that could deliver urgently needed antibiotics with novel mechanisms of action. In this work, we report the recognition mechanisms of four oxidization products (XAN, AZA, UAC, and HPA) generated during purine degradation by an RNA motif termed the NMT1 riboswitch. Specifically, we investigated the physical interactions between the riboswitch and the oxidized metabolites by computing the changes in the free energy on mutating key nucleobases in the ligand binding pocket of the riboswitch. We discovered that the electrostatic interactions are central to ligand discrimination by this riboswitch. The relative binding free energies of the mutations further indicated that some of the mutations can also strengthen the binding affinities of the ligands (AZA, UAC, and HPA). These mechanistic details are also potentially relevant in the design of novel compounds targeting riboswitches.

Translational and rotational diffusion of short ribonucleic acids

Inevitable precondition for ribonucleic acids to regulate gene expression and to perform gene editing is diffusion. Free three-dimensional translational diffusion velocity of RNA of up to 200 nucleotides could be predicted with high accuracy by the empirical model D = 4.58 10^-10 N^-0.39 m²s^-1. Furthermore, the biological function of ribonucleic acids is determined by rotational diffusion. In the presented work, an empirical model is derived applying atom-level shell-modeling of electron density maps, D_r = 1.62 10⁹ N^-1.20 s^-1, to predict the rotational diffusion coefficient of short ribonucleic acids based on the polymer size.

Discovery of RNA-targeted small molecules through the merging of experimental and computational technologies

INTRODUCTION

The field of RNA-targeted small molecules is rapidly evolving, owing to the advances in experimental and computational technologies. With the identification of several bioactive small molecules that target RNA, including the FDA-approved risdiplam, the biopharmaceutical industry is gaining confidence in the field. This review, based on the literature obtained from PubMed, aims to disseminate information about the various technologies developed for targeting RNA with small molecules and propose areas for improvement to develop drugs more efficiently, particularly those linked to diseases with unmet medical needs.

AREAS COVERED

The technologies for the identification of RNA targets, screening of chemical libraries against RNA, assessing the bioactivity and target engagement of the hit compounds, structure determination, and hit-to-lead optimization are reviewed. Along with the description of the technologies, their strengths, limitations, and examples of how they can impact drug discovery are provided.

EXPERT OPINION

Many existing technologies employed for protein targets have been repurposed for use in the discovery of RNA-targeted small molecules. In addition, technologies tailored for RNA targets have been developed. Nevertheless, more improvements are necessary, such as artificial intelligence to dissect important RNA structures and RNA-small-molecule interactions and more powerful chemical probing and structure prediction techniques.

Brief considerations on targeting RNA with small molecules

For more than three decades, RNA has been known to be a relevant and attractive macromolecule to target but figuring out which RNA should be targeted and how remains challenging. Recent years have seen the confluence of approaches for screening, drug optimization, and target validation that have led to the approval of a few RNA-targeting therapeutics for clinical applications. This focused perspective aims to highlight - but not exhaustively review - key factors accounting for these successes while pointing at crucial aspects worth considering for further breakthroughs.

Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch

Riboswitches normally regulate gene expression through structural changes in response to the specific binding of cellular metabolites or metal ions. Taking add adenine riboswitch as an example, we explore the influences of metal ions (especially for K⁺ and Mg²⁺ ions) on the structure and dynamics of riboswitch aptamer (with and without ligand) by using molecular dynamic (MD) simulations. Our results show that a two-state transition marked by the structural deformation at the connection of J12 and P1 (C_J12-P1) is not only related to the binding of cognate ligands, but also strongly coupled with the change of metal ion environments. Moreover, the deformation of the structure at C_J12-P1 can be transmitted to P1 directly connected to the expression platform in multiple ways, which will affect the structure and stability of P1 to varying degrees, and finally change the regulation state of this riboswitch.

Molecular dynamic simulations are employed to assess the influence of metal ions on riboswitch structure and dynamics, suggesting a conformational control of riboswitch aptamers by metal ions before ligand binding.

SILCS-RNA: Toward a Structure-Based Drug Design Approach for Targeting RNAs with Small Molecules

RNA molecules can act as potential drug targets in different diseases, as their dysregulated expression or misfolding can alter various cellular processes. Noncoding RNAs account for ∼70% of the human genome, and these molecules can have complex tertiary structures that present a great opportunity for targeting by small molecules. In the present study, the site identification by ligand competitive saturation (SILCS) computational approach is extended to target RNA, termed SILCS-RNA. Extensions to the method include an enhanced oscillating excess chemical potential protocol for the grand canonical Monte Carlo calculations and individual simulations of the neutral and charged solutes from which the SILCS functional group affinity maps (FragMaps) are calculated for subsequent binding site identification and docking calculations. The method is developed and evaluated against seven RNA targets and their reported small molecule ligands. SILCS-RNA provides a detailed characterization of the functional group affinity pattern in the small molecule binding sites, recapitulating the types of functional groups present in the ligands. The developed method is also shown to be useful for identification of new potential binding sites and identifying ligand moieties that contribute to binding, granular information that can facilitate ligand design. However, limitations in the method are evident including the ability to map the regions of binding sites occupied by ligand phosphate moieties and to fully account for the wide range of conformational heterogeneity in RNA associated with binding of different small molecules, emphasizing inherent challenges associated with applying computer-aided drug design methods to RNA. While limitations are present, the current study indicates how the SILCS-RNA approach may enhance drug discovery efforts targeting RNAs with small molecules.

Water-Controlled Keto-Enol Tautomerization of a Prebiotic Nucleobase

Barbituric acid is believed to be a proto-RNA nucleobase that was used for biological information transfer on prebiotic earth before DNA and RNA in their present forms evolved. Nucleobases have various tautomeric forms and the relative stability of these forms is critical to their biological function. It has been shown that barbituric acid has a tri-keto form in the gas phase and an enol form in the solid state. However, its dominant tautomeric form in aqueous medium that is most relevant for biology has been investigated only to a limited extent and the findings are inconclusive. We have used multiple approaches, namely, molecular dynamics, quantum chemistry, NMR, and IR spectroscopy to determine the most stable tautomer of barbituric acid in solution. We find a delicate balance in the stability of the two tautomers, tri-keto and enol, which is tipped toward the enol as the extent of solvation by water increases.

Synthesis, DFT Molecular Geometry and Anticancer Activity of Symmetrical 2,2'-(2-Oxo-1H-benzo[d]imidazole-1,3(2H)-diyl) Diacetate and Its Arylideneacetohydrazide Derivatives

To identify new candidate anticancer compounds, we here report the synthesis of benzimidazole derivatives: diethyl 2,2'-(2-oxo-1H-benzo[d]imidazole-1,3(2H)-diyl) diacetate and its arylideneacetohydrazide derivatives, using ultrasonic irradiation and conventional heating. The compounds were confirmed by Nuclear magnetic resonance (NMR) (JEOL, Tokyo, Japan) and Fourier transform infrared spectroscopy (FTIR) spectroscopy (Thermoscientific, Waltham, MA, USA). The molecular structure and electronic properties of the studied compounds were predicted for the acetohydrazide hydrazones. These compounds exist as a mixture of configurational and conformational isomerism as well as amido-amidic acid tautomerism. The NMR spectral data proved the predominance of syn-E amido isomers. In addition, density functional theory (DFT) predicted stability in the gas phase and showed that syn-E amido isomers are the most stable in the presence of an electron donating group, while the anti-isomer is the most stable in the presence of electron-attracting substituents. The anticancer activity of these synthetic compounds 6a, 6b and 6c towards both colon cancer (HCT-116) and cervical cancer (HeLa) cells was examined by MTT assay and DAPI staining. The MTT assay revealed a strong antiproliferative effect against the cancer cells at low concentrations, and interestingly, no significant inhibitory action against the non-cancerous cell line, HEK-293. The IC50 values for HCT-116 were 29.5 + 4.53 µM, 57.9 + 7.01 µM and 40.6 + 5.42 µM for 6a, 6b, and 6c, respectively. The IC50 values for HeLa cells were 57.1 + 6.7 µM, 65.6 + 6.63 µM and 33.8 + 3.54 µM for 6a, 6b, and 6c, respectively. DAPI staining revealed that these synthesized benzimidazole derivatives caused apoptotic cell death in both the colon and cervical cancer cells. Thus, these synthetic compounds demonstrate encouraging anticancer activity as well as being safe for normal human cells, making them attractive candidates as anticancer agents.

Structure and mechanism of the methyltransferase ribozyme MTR1

RNA-catalysed RNA methylation was recently shown to be part of the catalytic repertoire of ribozymes. The methyltransferase ribozyme MTR1 catalyses the site-specific synthesis of 1-methyladenosine (m¹A) in RNA, using O⁶-methylguanine (m⁶G) as methyl group donor. Here we report the crystal structure of MTR1 at a resolution of 2.8 Å, which reveals a guanine binding site reminiscent of natural guanine riboswitches. The structure represents the postcatalytic state of a split ribozyme in complex with the m¹A-containing RNA product and the demethylated cofactor guanine. The structural data suggest the mechanistic involvement of a protonated cytidine in the methyl transfer reaction. A synergistic effect of two 2'-O-methylated ribose residues in the active site results in accelerated methyl group transfer. Supported by these results, it seems plausible that modified nucleotides may have enhanced early RNA catalysis and that metabolite-binding riboswitches may resemble inactivated ribozymes that have lost their catalytic activity during evolution.

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.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

Insights into xanthine riboswitch structure and metal ion-mediated ligand recognition

Riboswitches are conserved functional domains in mRNA that mostly exist in bacteria. They regulate gene expression in response to varying concentrations of metabolites or metal ions. Recently, the NMT1 RNA motif has been identified to selectively bind xanthine and uric acid, respectively, both are involved in the metabolic pathway of purine degradation. Here, we report a crystal structure of this RNA bound to xanthine. Overall, the riboswitch exhibits a rod-like, continuously stacked fold composed of three stems and two internal junctions. The binding-pocket is determined by the highly conserved junctional sequence J1 between stem P1 and P2a, and engages a long-distance Watson-Crick base pair to junction J2. Xanthine inserts between a G-U pair from the major groove side and is sandwiched between base triples. Strikingly, a Mg2+ ion is inner-sphere coordinated to O6 of xanthine and a non-bridging oxygen of a backbone phosphate. Two further hydrated Mg2+ ions participate in extensive interactions between xanthine and the pocket. Our structure model is verified by ligand binding analysis to selected riboswitch mutants using isothermal titration calorimetry, and by fluorescence spectroscopic analysis of RNA folding using 2-aminopurine-modified variants. Together, our study highlights the principles of metal ion-mediated ligand recognition by the xanthine riboswitch.

A natural riboswitch scaffold with self-methylation activity

Methylation is a prevalent post-transcriptional modification encountered in coding and non-coding RNA. For RNA methylation, cells use methyltransferases and small organic substances as methyl-group donors, such as S-adenosylmethionine (SAM). SAM and other nucleotide-derived cofactors are viewed as evolutionary leftovers from an RNA world, in which riboswitches have regulated, and ribozymes have catalyzed essential metabolic reactions. Here, we disclose the thus far unrecognized direct link between a present-day riboswitch and its inherent reactivity for site-specific methylation. The key is O6-methyl pre-queuosine (m6preQ1), a potentially prebiotic nucleobase which is recognized by the native aptamer of a preQ1 class I riboswitch. Upon binding, the transfer of the ligand's methyl group to a specific cytidine occurs, installing 3-methylcytidine (m3C) in the RNA pocket under release of pre-queuosine (preQ1). Our finding suggests that nucleic acid-mediated methylation is an ancient mechanism that has offered an early path for RNA epigenetics prior to the evolution of protein methyltransferases. Furthermore, our findings may pave the way for the development of riboswitch-descending methylation tools based on rational design as a powerful alternative to in vitro selection approaches.

Multiscale Modeling of Wobble to Watson-Crick-Like Guanine-Uracil Tautomerization Pathways in RNA

Energetically unfavorable Watson–Crick (WC)-like tautomeric forms of nucleobases are known to introduce spontaneous mutations, and contribute to replication, transcription, and translation errors. Recent NMR relaxation dispersion techniques were able to show that wobble (w) G•U mispair exists in equilibrium with the short-lived, low-population WC-like enolic tautomers. Presently, we have investigated the wG•U → WC-like enolic reaction pathway using various theoretical methods: quantum mechanics (QM), molecular dynamics (MD), and combined quantum mechanics/molecular mechanics (QM/MM). The previous studies on QM gas phase calculations were inconsistent with experimental data. We have also explored the environmental effects on the reaction energies by adding explicit water. While the QM-profile clearly becomes endoergic in the presence of water, the QM/MM-profile remains consistently endoergic in the presence and absence of water. Hence, by including microsolvation and QM/MM calculations, the experimental data can be explained. For the G•Uenol→ Genol•U pathway, the latter appears to be energetically more favorable throughout all computational models. This study can be considered as a benchmark of various computational models of wG•U to WC-like tautomerization pathways with and without the environmental effects, and may contribute on further studies of other mispairs as well.

Deaminated purine bypass by DNA polymerase η

A recent work by Jung and colleagues (Biochem J.477, 4797–4810) provides an explanation of how DNA polymerase η replicates through deaminated purine bases such as xanthine and hypoxanthine. This commentary discusses the crystal structures of the polymerase η complexes that implicate the role of tautomerism in the bypass of these DNA lesions.

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.

ITScore-NL: An Iterative Knowledge-Based Scoring Function for Nucleic Acid-Ligand Interactions

Nucleic acid-ligand complexes underlie numerous cellular processes, such as gene function expression and regulation, in which their three-dimensional structures are important to understand their functions and thus to develop therapeutic interventions. Given the high cost and technical difficulties in experimental methods, computational methods such as molecular docking have been actively used to investigate nucleic acid-ligand interactions in which an accurate scoring function is crucial. However, because of the limited number of experimental nucleic acid-ligand binding data and structures, the scoring function development for nucleic acid-ligand interactions falls far behind that for protein-protein and protein-ligand interactions. Here, based on our statistical mechanics-based iterative approach, we have developed an iterative knowledge-based scoring function for nucleic acid-ligand interactions, named as ITScore-NL, by explicitly including stacking and electrostatic potentials. Our ITScore-NL scoring function was extensively evaluated for its ability in the binding mode and binding affinity predictions on three diverse test sets and compared with state-of-the-art scoring functions. Overall, ITScore-NL obtained significantly better performance than the other 12 scoring functions and predicted near-native poses with rmsd ≤ 1.5 Å for 71.43% of the cases when the top three binding modes were considered and a good correlation of R = 0.64 in binding affinity prediction on the large test set of 77 nucleic acid-ligand complexes. These results suggested the accuracy of ITScore-NL and the necessity of explicitly including stacking and electrostatic potentials.

Site-specific RNA methylation by a methyltransferase ribozyme

Nearly all classes of coding and non-coding RNA undergo post-transcriptional modification, including RNA methylation. Methylated nucleotides are among the evolutionarily most-conserved features of transfer (t)RNA and ribosomal (r)RNA1,2. Many contemporary methyltransferases use the universal cofactor S-adenosylmethionine (SAM) as a methyl-group donor. SAM and other nucleotide-derived cofactors are considered to be evolutionary leftovers from an RNA world, in which ribozymes may have catalysed essential metabolic reactions beyond self-replication3. Chemically diverse ribozymes seem to have been lost in nature, but may be reconstructed in the laboratory by in vitro selection. Here we report a methyltransferase ribozyme that catalyses the site-specific installation of 1-methyladenosine in a substrate RNA, using O6-methylguanine as a small-molecule cofactor. The ribozyme shows a broad RNA-sequence scope, as exemplified by site-specific adenosine methylation in various RNAs. This finding provides fundamental insights into the catalytic abilities of RNA, serves a synthetic tool to install 1-methyladenosine in RNA and may pave the way to in vitro evolution of other methyltransferase and demethylase ribozymes.

SPRINT: a Cas13a-based platform for detection of small molecules

Recent efforts in biological engineering have made detection of nucleic acids in samples more rapid, inexpensive and sensitive using CRISPR-based approaches. We expand one of these Cas13a-based methods to detect small molecules in a one-batch assay. Using SHERLOCK-based profiling of in vitrotranscription (SPRINT), in vitro transcribed RNA sequence-specifically triggers the RNase activity of Cas13a. This event activates its non-specific RNase activity, which enables cleavage of an RNA oligonucleotide labeled with a quencher/fluorophore pair and thereby de-quenches the fluorophore. This fluorogenic output can be measured to assess transcriptional output. The use of riboswitches or proteins to regulate transcription via specific effector molecules is leveraged as a coupled assay that transforms effector concentration into fluorescence intensity. In this way, we quantified eight different compounds, including cofactors, nucleotides, metabolites of amino acids, tetracycline and monatomic ions in samples. In this manner, hundreds of reactions can be easily quantified in a few hours. This increased throughput also enables detailed characterization of transcriptional regulators, synthetic compounds that inhibit transcription, or other coupled enzymatic reactions. These SPRINT reactions are easily adaptable to portable formats and could therefore be used for the detection of analytes in the field or at point-of-care situations.

Rising trend on the halogen and non-halogen derivatives (Br, Cl, CF3, F, CH3 and NH2) in ruminal β-d-Xylopyranose - a quantum chemical perspective

The utmost aim of the current study is to find significance of the binding affinity in the halogen and non-halogen derivatives: Br, Cl, CF₃, F, CH₃ and NH₂ of β-d-Xylopyranose with the hinge region amino acids of ruminant-β-glycosidase. The interaction energy analysis was carried out in detail through various density functional studies as M062X/def2-QZVP, M062X/LANL2DZ, B3LYP/LANL2DZ and M06HF/LANL2DZ level of theories. The total interaction energy of halogen derivatives: Br, Cl, F and CF₃ are -618.21, -599.00, -720.45 and -553.08 kcal/mol respectively, and non-halogen derivative: amine group (NH₂) is -87.96 kcal/mol at M062X/def2-QZVP level of theory, which exist with strong binding affinity. Ligand properties: dipole moment, polarizability, volume, molecular mass, electrostatic potential map was evaluated to understand its electrostatic and structural behavior. The nature of the bonds was inferred from the electrostatic potential map for all the halogen and non-halogen derivatives ligand. The stabilization energy from NBO analysis reveals the stability of single hydrogen and halogen bonds (N-H…Br, C-Br…O, N-H…Cl, C-Cl…O, O-H…F, C-H…F, N-H…F, C-F…O, N-H…O, O-H…O, N-H…N, O-H…N) in β-d-Xylopyranose and its derivatives. Overall, this study paves way for scientist and medicinal chemist in modelling new drugs. Further, it suggests mutations that increase the binding and may enhance the catalytic action and strengthen the complex diet in animals and hence recommended for experimental synthesis.Communicated by Ramaswamy H. Sarma.

Parallel Discovery Strategies Provide a Basis for Riboswitch Ligand Design

Riboswitches are mRNA domains that make gene-regulatory decisions upon binding their cognate ligands. Bacterial riboswitches that specifically recognize 5-aminoimidazole-4-carboxamide riboside 5'-monophosphate (ZMP) and 5'-triphosphate (ZTP) regulate genes involved in folate and purine metabolism. Now, we have developed synthetic ligands targeting ZTP riboswitches by replacing the sugar-phosphate moiety of ZMP with various functional groups, including simple heterocycles. Despite losing hydrogen bonds from ZMP, these analogs bind ZTP riboswitches with similar affinities as the natural ligand, and activate transcription more strongly than ZMP in vitro. The most active ligand stimulates gene expression ∼3 times more than ZMP in a live Escherichia coli reporter. Co-crystal structures of the Fusobacterium ulcerans ZTP riboswitch bound to synthetic ligands suggest stacking of their pyridine moieties on a conserved RNA nucleobase primarily determines their higher activity. Altogether, these findings guide future design of improved riboswitch activators and yield insights into how RNA-targeted ligand discovery may proceed.

Crystal structure of an atypical cobalamin riboswitch reveals RNA structural adaptability as basis for promiscuous ligand binding

Cobalamin riboswitches encompass a structurally diverse group of cis-acting, gene regulatory elements found mostly in bacterial messenger RNA and are classified into subtypes based on secondary and tertiary characteristics. An unusual variant of the cobalamin riboswitch with predicted structural features was identified in Bacillus subtilis over a decade ago, but its structure and mechanisms of cobalamin selectivity and translational control have remained unsolved. We present the crystal structure of the aptamer domain of this atypical cobalamin riboswitch and a model for the complete riboswitch, including its expression platform domain. We demonstrate that this riboswitch binds to multiple cobalamin derivatives and correlate its promiscuous behavior to its structure and unique arrangement of peripheral elements. Comparative structural analyses between conventional cobalamin riboswitches and the B. subtilis cobalamin riboswitch reveal that the likely basis for this promiscuous ligand binding is intrinsic structural adaptability encoded in the RNA structure. It suggests that cobalamin selectivity might ultimately be viewed as existing on a spectrum of affinity for each derivative rather than as belonging to distinct types based on ligand specificities. Our work provides an interesting and notable example of functional coupling of ligand-sensing and adaptive folding by a structured RNA molecule.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Systematic analysis of the interactions driving small molecule-RNA recognition

RNA molecules are becoming an important target class in drug discovery. However, the principles for designing RNA-binding small molecules are yet to be fully uncovered. In this study, we examined the Protein Data Bank (PDB) to highlight privileged interactions underlying small molecule-RNA recognition. By comparing this analysis with previously determined small molecule-protein interactions, we find that RNA recognition is driven mostly by stacking and hydrogen bonding interactions, while protein recognition is instead driven by hydrophobic effects. Furthermore, we analyze patterns of interactions to highlight potential strategies to tune RNA recognition, such as stacking and cation-π interactions that favor purine and guanine recognition, and note an unexpected paucity of backbone interactions, even for cationic ligands. Collectively, this work provides further understanding of RNA-small molecule interactions that may inform the design of small molecules targeting RNA.

High Affinity Binding of N2-Modified Guanine Derivatives Significantly Disrupts the Ligand Binding Pocket of the Guanine Riboswitch

Riboswitches are important model systems for the development of approaches to search for RNA-targeting therapeutics. A principal challenge in finding compounds that target riboswitches is that the effector ligand is typically almost completely encapsulated by the RNA, which severely limits the chemical space that can be explored. Efforts to find compounds that bind the guanine/adenine class of riboswitches with a high affinity have in part focused on purines modified at the C6 and C2 positions. These studies have revealed compounds that have low to sub-micromolar affinity and, in a few cases, have antimicrobial activity. To further understand how these compounds interact with the guanine riboswitch, we have performed an integrated structural and functional analysis of representative guanine derivatives with modifications at the C8, C6 and C2 positions. Our data indicate that while modifications of guanine at the C6 position are generally unfavorable, modifications at the C8 and C2 positions yield compounds that rival guanine with respect to binding affinity. Surprisingly, C2-modified guanines such as N2-acetylguanine completely disrupt a key Watson–Crick pairing interaction between the ligand and RNA. These compounds, which also modulate transcriptional termination as efficiently as guanine, open up a significant new chemical space of guanine modifications in the search for antimicrobial agents that target purine riboswitches.

Structural basis for 2'-deoxyguanosine recognition by the 2'-dG-II class of riboswitches

A recent bioinformatic analysis of well-characterized classes of riboswitches uncovered subgroups unable to bind to the regulatory molecule of the parental class. Within the guanine/adenine class, seven groups of RNAs were identified that deviate from the consensus sequence at one or more of three positions directly involved purine nucleobase recognition, one of which was validated as a second class of 2'-deoxyguanosine riboswitch (called 2'-dG-II). To understand how 2'-dG-II riboswitches recognize their cognate ligand and how they differ from a previously identified class of 2'-deoxyguanosine binding riboswitches, we have solved the crystal structure of a 2'-dG-II aptamer domain bound to 2'-deoxyguanosine. This structure reveals a global architecture similar to other members of the purine riboswitch family, but contains key differences within the ligand binding core. Defining the 2'-dG-II riboswitches is a two-nucleotide insertion in the three-way junction that promotes novel base-base interactions. Unlike 2'-dG-I riboswitches, the 2'-dG-II class only requires local changes to the ligand binding pocket of the guanine/adenine class to achieve a change in ligand preference. Notably, members of the 2'-dG-II family have variable ability to discriminate between 2'-deoxyguanosine and riboguanosine, suggesting that a subset of 2'-dG-II riboswitches may bind either molecule to regulate gene expression.

Ligand Binding Mechanism and Its Relationship with Conformational Changes in Adenine Riboswitch

Riboswitches are naturally occurring RNA aptamers that control the expression of essential bacterial genes by binding to specific small molecules. The binding with both high affinity and specificity induces conformational changes. Thus, riboswitches were proposed as a possible molecular target for developing antibiotics and chemical tools. The adenine riboswitch can bind not only to purine analogues but also to pyrimidine analogues. Here, long molecular dynamics (MD) simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) computational methodologies were carried out to show the differences in the binding model and the conformational changes upon five ligands binding. The binding free energies of the guanine riboswitch aptamer with C74U mutation complexes were compared to the binding free energies of the adenine riboswitch (AR) aptamer complexes. The calculated results are in agreement with the experimental data. The differences for the same ligand binding to two different aptamers are related to the electrostatic contribution. Binding dynamical analysis suggests a flexible binding pocket for the pyrimidine ligand in comparison with the purine ligand. The 18 μs of MD simulations in total indicate that both ligand-unbound and ligand-bound aptamers transfer their conformation between open and closed states. The ligand binding obviously affects the conformational change. The conformational states of the aptamer are associated with the distance between the mass center of two key nucleotides (U51 and A52) and the mass center of the other two key nucleotides (C74 and C75). The results suggest that the dynamical character of the binding pocket would affect its biofunction. To design new ligands of the adenine riboswitch, it is recommended to consider the binding affinities of the ligand and the conformational change of the ligand binding pocket.

Engineering Cellular Signal Sensors based on CRISPR-sgRNA Reconstruction Approaches

The discovery of the CRISPR systems has enriched the application of gene therapy and biotechnology. As a type of robust and simple toolbox, the CRISPR system has greatly promoted the development of cellular signal sensors at the genomic level. Although CRISPR systems have demonstrated that they can be used in eukaryotic and even mammalian cells after extraction from prokaryotic cells, controlling their gene-editing activity remains a challenge. Here we summarize the advantages and disadvantages of building a CRIRPR-based signal sensor through sgRNA reconstruction, as well as possible ways to reprogram the signal network of cells. We also propose how to further improve the design of the current signal sensors based on sgRNA-riboswitch. We believe that the development of these technologies and the construction of platforms can further promote the development of environment detection, disease diagnosis, and gene therapy by means of synthetic biology.

Molecular dynamics simulation of the binding process of ligands to the add adenine riboswitch aptamer

Riboswitches are RNA-structured elements that modulate gene expression through changing their conformation in response to specific metabolite binding. However, the regulation mechanisms of riboswitches by ligand binding are still not well understood. At present two possible ligand-binding mechanisms have been proposed: conformational selection and induced fit. Based on explicit-solvent molecular dynamics (MD) simulations, we have studied the process of the binding of ligands (adenines) to add adenine riboswitch aptamer (AARA) in detail. Our results show that the relative high flexibility of the junction J23 of AARA allows the ligands to be captured by the binding pocket of AARA in a near-native state, which may be driven by hydrophobic and base-stacking interactions. In addition, the binding of a ligand makes the stem P1 and the junction J23 of AARA more stable than in the absence of the ligand. Generally, our analyses show that the ligand-binding process of the add adenine riboswitch may be partially embodied by a conformational selection mechanism.

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations

Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

Riboswitch distribution, structure, and function in bacteria

Riboswitches are gene control elements that directly bind to specific ligands to regulate gene expression without the need for proteins. They are found in all three domains of life, including Bacteria, Archaea, and Eukaryota. Riboswitches are mostly spread in bacteria and archaea. In this paper, we discuss the general distribution, structure, and function of 28 different riboswitch classes as we focus our attention on riboswitches in bacteria. Bacterial riboswitches regulate gene expression by four distinct mechanisms. They regulate the expression of a limited number of genes. However, most of these genes are responsible for the synthesis of essential metabolites without which the cell cannot function. Therefore, riboswitch distribution is also important for antibacterial drug development.

Synthetic ligands for PreQ1 riboswitches provide structural and mechanistic insights into targeting RNA tertiary structure

Riboswitches are naturally occurring RNA aptamers that regulate gene expression by binding to specific small molecules. Riboswitches control the expression of essential bacterial genes and are important models for RNA-small molecule recognition. Here, we report the discovery of a class of synthetic small molecules that bind to PreQ₁ riboswitch aptamers. These molecules bind specifically and reversibly to the aptamers with high affinity and induce a conformational change. Furthermore, the ligands modulate riboswitch activity through transcriptional termination despite no obvious chemical similarity to the cognate ligand. X-ray crystallographic studies reveal that the ligands share a binding site with the cognate ligand but make different contacts. Finally, alteration of the chemical structure of the ligand causes changes in the mode of RNA binding and affects regulatory function. Thus, target- and structure-based approaches can be used to identify and understand the mechanism of synthetic ligands that bind to and regulate complex, folded RNAs.

RNA sensors—Riboswitches—respond to the binding of small molecules ligands through structure modification. Here the authors identify synthetic small molecules that bind and regulate the activity of PreQ₁ Riboswitches despite having no obvious chemical similarity to the cognate ligand.

Synthesis, characterization, photoluminescence, anti-tumor activity, DFT calculations and molecular docking with proteins of zinc(ii) halogen substituted terpyridine compounds

Study on the synthesis, characterization, photoluminescence and anti-tumor activity of a series of zinc(ii) halogen substituted terpyridine complexes.

Blueprints for Biosensors: Design, Limitations, and Applications

Biosensors are enabling major advances in the field of analytics that are both facilitating and being facilitated by advances in synthetic biology. The ability of biosensors to rapidly and specifically detect a wide range of molecules makes them highly relevant to a range of industrial, medical, ecological, and scientific applications. Approaches to biosensor design are as diverse as their applications, with major biosensor classes including nucleic acids, proteins, and transcription factors. Each of these biosensor types has advantages and limitations based on the intended application, and the parameters that are required for optimal performance. Specifically, the choice of biosensor design must consider factors such as the ligand specificity, sensitivity, dynamic range, functional range, mode of output, time of activation, ease of use, and ease of engineering. This review discusses the rationale for designing the major classes of biosensor in the context of their limitations and assesses their suitability to different areas of biotechnological application.

Principles for targeting RNA with drug-like small molecules

Recent studies have indicated the potential to develop small-molecule drugs that act on RNA targets, leading to burgeoning interest in the field. This article discusses general principles for discovering small-molecule drugs that target RNA and argues that the overarching challenge is to identify appropriate target structures in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets.

RNA molecules are essential for cellular information transfer and gene regulation, and RNAs have been implicated in many human diseases. Messenger and non-coding RNAs contain highly structured elements, and evidence suggests that many of these structures are important for function. Targeting these RNAs with small molecules offers opportunities to therapeutically modulate numerous cellular processes, including those linked to 'undruggable' protein targets. Despite this promise, there is currently only a single class of human-designed small molecules that target RNA used clinically — the linezolid antibiotics. However, a growing number of small-molecule RNA ligands are being identified, leading to burgeoning interest in the field. Here, we discuss principles for discovering small-molecule drugs that target RNA and argue that the overarching challenge is to identify appropriate target structures — namely, in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets. If focus is placed on such druggable binding sites in RNA, extensive knowledge of the typical physicochemical properties of drug-like small molecules could then enable small-molecule drug discovery for RNA targets to become (only) roughly as difficult as for protein targets.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Sensitive and specific detection of ligands using engineered riboswitches

Riboswitches are RNA elements found in non-coding regions of messenger RNAs that regulate gene expression through a ligand-triggered conformational change. Riboswitches typically bind tightly and specifically to their ligands, so they have the potential to serve as highly effective sensors in vitro. In B. subtilis and other gram-positive bacteria, purine nucleotide synthesis is regulated by riboswitches that bind to guanine. We modified the xpt-pbuX guanine riboswitch for use in a fluorescence quenching assay that allowed us to specifically detect and quantify guanine in vitro. Using this assay, we reproducibly detected as little as 5 nM guanine. We then produced sensors for 2'-deoxyguanosine and cyclic diguanylate (c-diGMP) by appending the P1 stem of the guanine riboswitch to the ligand-binding domains of a 2'-deoxyguanosine riboswitch and a c-diGMP riboswitch. These hybrid sensors could detect 15 nM 2'-deoxyguanosine and 3 nM c-diGMP, respectively. Each sensor retained the ligand specificity of its corresponding natural riboswitch. In order to extend the utility of our approach, we developed a strategy for the in vitro selection of sensors with novel ligand specificity. Here we report a proof-of-principle experiment that demonstrated the feasibility of our selection strategy.

Allosteric mechanism of the V. vulnificus adenine riboswitch resolved by four-dimensional chemical mapping

The structural interconversions that mediate the gene regulatory functions of RNA molecules may be different from classic models of allostery, but the relevant structural correlations have remained elusive in even intensively studied systems. Here, we present a four-dimensional expansion of chemical mapping called lock-mutate-map-rescue (LM²R), which integrates multiple layers of mutation with nucleotide-resolution chemical mapping. This technique resolves the core mechanism of the adenine-responsive V. vulnificus add riboswitch, a paradigmatic system for which both Monod-Wyman-Changeux (MWC) conformational selection models and non-MWC alternatives have been proposed. To discriminate amongst these models, we locked each functionally important helix through designed mutations and assessed formation or depletion of other helices via compensatory rescue evaluated by chemical mapping. These LM²R measurements give strong support to the pre-existing correlations predicted by MWC models, disfavor alternative models, and suggest additional structural heterogeneities that may be general across ligand-free riboswitches.

Using Genome Sequence to Enable the Design of Medicines and Chemical Probes

Rapid progress in genome sequencing technology has put us firmly into a postgenomic era. A key challenge in biomedical research is harnessing genome sequence to fulfill the promise of personalized medicine. This Review describes how genome sequencing has enabled the identification of disease-causing biomolecules and how these data have been converted into chemical probes of function, preclinical lead modalities, and ultimately U.S. Food and Drug Administration (FDA)-approved drugs. In particular, we focus on the use of oligonucleotide-based modalities to target disease-causing RNAs; small molecules that target DNA, RNA, or protein; the rational repurposing of known therapeutic modalities; and the advantages of pharmacogenetics. Lastly, we discuss the remaining challenges and opportunities in the direct utilization of genome sequence to enable design of medicines.

Purine analogs targeting the guanine riboswitch as potential antibiotics against Clostridioides difficile

Riboswitches recently emerged as possible targets for the development of alternative antimicrobial approaches. Guanine-sensing riboswitches in the bacterial pathogen Clostridioides difficile (formerly known as Clostridium difficile) constitute potential targets based on their involvement in the regulation of basal metabolic control of purine compounds. In this study, we deciphered the structure-activity relationship of several guanine derivatives on the guanine riboswitch and determined their antimicrobial activity. We describe the synthesis of purine analogs modified in ring B as well as positions 2 and 6. Their biological activity was determined by measuring their affinity for the C. difficile guanine riboswitch and their inhibitory effect on bacterial growth, including a counter-screen to discriminate against riboswitch-independent antibacterial effects. Altogether, our results suggest that improvements in riboswitch binding affinity in vitro do not necessarily translate into improved antibacterial activity in bacteria, despite the fact that some structure-activity relationship was observed at least with respect to binding affinity.

Cobalamin riboswitches exhibit a broad range of ability to discriminate between methylcobalamin and adenosylcobalamin

Riboswitches are a widely distributed class of regulatory RNAs in bacteria that modulate gene expression via small-molecule-induced conformational changes. Generally, these RNA elements are grouped into classes based upon conserved primary and secondary structure and their cognate effector molecule. Although this approach has been very successful in identifying new riboswitch families and defining their distributions, small sequence differences between structurally related RNAs can alter their ligand selectivity and regulatory behavior. Herein, we use a structure-based mutagenic approach to demonstrate that cobalamin riboswitches have a broad spectrum of preference for the two biological forms of cobalamin in vitro using isothermal titration calorimetry. This selectivity is primarily mediated by the interaction between a peripheral element of the RNA that forms a T-loop module and a subset of nucleotides in the cobalamin-binding pocket. Cell-based fluorescence reporter assays in Escherichia coli revealed that mutations that switch effector preference in vitro lead to differential regulatory responses in a biological context. These data demonstrate that a more comprehensive analysis of representative sequences of both previously and newly discovered classes of riboswitches might reveal subgroups of RNAs that respond to different effectors. Furthermore, this study demonstrates a second distinct means by which tertiary structural interactions in cobalamin riboswitches dictate ligand selectivity.

Ligand Selectivity Mechanism and Conformational Changes in Guanine Riboswitch by Molecular Dynamics Simulations and Free Energy Calculations

Riboswitches regulate gene expression through direct and specific interactions with small metabolite molecules. Binding of a ligand to its RNA target is high selectivity and affinity and induces conformational changes of the RNA's secondary and tertiary structure. The structural difference of two purine riboswitches aptamers is caused by only one single mutation, where cytosine 74 in the guanine riboswitch is corresponding to a uracil 74 in adenine riboswitch. Here we employed molecular dynamics (MD) simulation, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and thermodynamic integration computational methodologies to evaluate the energetic and conformational changes of ligands binding to purine riboswitches. The snapshots used in MM-PBSA calculation were extracted from ten 50 ns MD simulation trajectories for each complex. These free energy results are in consistent with the experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that the loss in binding free energy upon mutation is mainly electrostatic in guanine (GUA) and riboswitch complex. Furthermore, new hydrogen bonds are found in mutated complexes. To reveal the conformational properties of guanine riboswitch, we performed a total of 6 μs MD simulations in both the presence and the absence of the ligand GUA. The MD simulations suggest that the conformation of guanine riboswitch depends on the distance of two groups in the binding pocket of ligand. The conformation is in a close conformation when U51-A52 is close to C74-U75.

Engineering and In Vivo Applications of Riboswitches

Riboswitches are common gene regulatory units mostly found in bacteria that are capable of altering gene expression in response to a small molecule. These structured RNA elements consist of two modular subunits: an aptamer domain that binds with high specificity and affinity to a target ligand and an expression platform that transduces ligand binding to a gene expression output. Significant progress has been made in engineering novel aptamer domains for new small molecule inducers of gene expression. Modified expression platforms have also been optimized to function when fused with both natural and synthetic aptamer domains. As this field expands, the use of these privileged scaffolds has permitted the development of tools such as RNA-based fluorescent biosensors. In this review, we summarize the methods that have been developed to engineer new riboswitches and highlight applications of natural and synthetic riboswitches in enzyme and strain engineering, in controlling gene expression and cellular physiology, and in real-time imaging of cellular metabolites and signals.