Structure and ligand binding of the SAM-V riboswitch

SAM-III aptamer enables detection of enzymatic SAM analogue generation

S-Adenosyl-L-methionine (SAM) is crucial for methylation and tightly controlled in cells. We examined SAM-III riboswitch response to 17 SAM analogues and used a Spinach/SAM aptasensor to monitor their enzymatic formation in situ. Most SAM analogues were recognized, unless they featured an ortho-substituted benzyl ring, indicating potential regulatory effects via SAM riboswitches.

Linker-Mediated Inactivation of the SAM-II Domain in the Tandem SAM-II/SAM-V Riboswitch

Tandem SAM-II/SAM-V riboswitch belongs to a class of riboswitches found in the marine bacterium 'Candidatus Pelagibacter ubique'. Previous studies have demonstrated that these riboswitches have the potential for digital modulation of gene expression at both the transcriptional and translational levels. In this study, we investigate the conformational changes in the tandem SAM-II/SAM-V riboswitch binding to S-adenosylmethionine (SAM) using selective 2'-hydroxyl acylation analyzed by the primer extension (SHAPE) assay, small-angle X-ray scattering (SAXS), and oligos depressing probing. Our findings reveal that the linker between SAM-II/SAM-V aptamers blocks the SAM response of the SAM-II domain. This result proposes a new mechanism for gene expression regulation, where the ligand-binding functions of tandem riboswitches can be selectively masked or released through a linker.

Bifidobacterium bifidum SAM-VI Riboswitch Conformation Change Requires Peripheral Helix Formation

The Bifidobacterium bifidum SAM-VI riboswitch undergoes dynamic conformational changes that modulate downstream gene expression. Traditional structural methods such as crystallography capture the bound conformation at high resolution, and additional efforts would reveal details from the dynamic transition. Here, we revealed a transcription-dependent conformation model for Bifidobacterium bifidum SAM-VI riboswitch. In this study, we combine small-angle X-ray scattering, chemical probing, and isothermal titration calorimetry to unveil the ligand-binding properties and conformational changes of the Bifidobacterium bifidum SAM-VI riboswitch and its variants. Our results suggest that the SAM-VI riboswitch contains a pre-organized ligand-binding pocket and stabilizes into the bound conformation upon binding to SAM. Whether the P1 stem formed and variations in length critically influence the conformational dynamics of the SAM-VI riboswitch. Our study provides the basis for artificially engineering the riboswitch by manipulating its peripheral sequences without modifying the SAM-binding core.

Biophysical Investigation of RNA ⋅ DNA : DNA Triple Helix and RNA : DNA Heteroduplex Formation by the lncRNAs MEG3 and Fendrr

Long non-coding RNAs (lncRNAs) are important regulators of gene expression and can associate with DNA as RNA : DNA heteroduplexes or RNA ⋅ DNA : DNA triple helix structures. Here, we review in vitro biochemical and biophysical experiments including electromobility shift assays (EMSA), circular dichroism (CD) spectroscopy, thermal melting analysis, microscale thermophoresis (MST), single-molecule Förster resonance energy transfer (smFRET) and nuclear magnetic resonance (NMR) spectroscopy to investigate RNA ⋅ DNA : DNA triple helix and RNA : DNA heteroduplex formation. We present the investigations of the antiparallel triplex-forming lncRNA MEG3 targeting the gene TGFB2 and the parallel triplex-forming lncRNA Fendrr with its target gene Emp2. The thermodynamic properties of these oligonucleotides lead to concentration-dependent heterogeneous mixtures, where a DNA duplex, an RNA : DNA heteroduplex and an RNA ⋅ DNA : DNA triplex coexist and their relative populations are modulated in a temperature-dependent manner. The in vitro data provide a reliable readout of triplex structures, as RNA ⋅ DNA : DNA triplexes show distinct features compared to DNA duplexes and RNA : DNA heteroduplexes. Our experimental results can be used to validate computationally predicted triple helix formation between novel disease-relevant lncRNAs and their DNA target genes.

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.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Structure and stability of different triplets involving artificial nucleobases: clues for the formation of semisynthetic triple helical DNA

A triple helical DNA can control gene expression, help in homologous recombination, induce mutations to facilitate DNA repair mechanisms, suppress oncogene formations, etc. However, the structure and function of semisynthetic triple helical DNA are not known. To understand this, various triplets formed between eight artificial nucleobases (P, Z, J, V, B, S, X, and K) and four natural DNA bases (G, C, A, and T) are studied herein by employing a reliable density functional theoretic (DFT) method. Initially, the triple helix-forming artificial nucleobases interacted with the duplex DNA containing GC and AT base pairs, and subsequently, triple helix-forming natural bases (G and C) interacted with artificial duplex DNA containing PZ, JV, BS, and XK base pairs. Among the different triplets formed in the first category, the C-JV triplet is found to be the most stable with a binding energy of about - 31 kcal/mol. Similarly, among the second category of triplets, the Z-GC and V-GC triplets are the most stable. Interestingly, Z-GC and V-GC are found to be isoenergetic with a binding energy of about - 30 kcal/mol. The C-JV, and Z-GC or V-GC triplets are about 12-14 kcal/mol more stable than the JV and GC base pairs respectively. Microsolvation of these triplets in 5 explicit water molecules further enhanced their stability by 16-21 kcal/mol. These results along with the consecutive stacking of the C-JV triplet (C-JV/C-JV) data indicate that the synthetic nucleobases can form stable semisynthetic triple helical DNA. However, consideration of a full-length DNA containing one or more semisynthetic bases or base pairs is necessary to understand the formation of semisynthetic DNA in living cells.

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

Structure and function analysis of a type III preQ1-I riboswitch from Escherichia coli reveals direct metabolite sensing by the Shine-Dalgarno sequence

Riboswitches are small noncoding RNAs found primarily in the 5' leader regions of bacterial messenger RNAs where they regulate expression of downstream genes in response to binding one or more cellular metabolites. Such noncoding RNAs are often regulated at the translation level, which is thought to be mediated by the accessibility of the Shine-Dalgarno sequence (SDS) ribosome-binding site. Three classes (I-III) of prequeuosine1 (preQ1)-sensing riboswitches are known that control translation. Class I is divided into three subtypes (types I-III) that have diverse mechanisms of sensing preQ1, which is involved in queuosine biosynthesis. To provide insight into translation control, we determined a 2.30 Å-resolution cocrystal structure of a class I type III preQ1-sensing riboswitch identified in Escherichia coli (Eco) by bioinformatic searches. The Eco riboswitch structure differs from previous preQ1 riboswitch structures because it has the smallest naturally occurring aptamer and the SDS directly contacts the preQ1 metabolite. We validated structural observations using surface plasmon resonance and in vivo gene-expression assays, which showed strong switching in live E. coli. Our results demonstrate that the Eco riboswitch is relatively sensitive to mutations that disrupt noncanonical interactions that form the pseudoknot. In contrast to type II preQ1 riboswitches, a kinetic analysis showed that the type III Eco riboswitch strongly prefers preQ1 over the chemically similar metabolic precursor preQ0. Our results reveal the importance of noncanonical interactions in riboswitch-driven gene regulation and the versatility of the class I preQ1 riboswitch pseudoknot as a metabolite-sensing platform that supports SDS sequestration.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibility and the role of co-transcriptional folding in these processes.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

Observation of structural switch in nascent SAM-VI riboswitch during transcription at single-nucleotide and single-molecule resolution

Growing RNAs fold differently as they are transcribed, which modulates their finally adopted structures. Riboswitches regulate gene expression by structural change, which are sensitive to co-transcriptionally structural biology. Here we develop a strategy to track the structural change of RNAs during transcription at single-nucleotide and single-molecule resolution and use it to monitor individual transcripts of the SAM-VI riboswitch (riboSAM) as transcription proceeds, observing co-existence of five states in riboSAM. We report a bifurcated helix in one newly identified state from NMR and single-molecule FRET (smFRET) results, and its presence directs the translation inhibition in our cellular translation experiments. A model is proposed to illustrate the distinct switch patterns and gene-regulatory outcome of riboSAM when SAM is present or absent. Our strategy enables the precise mapping of RNAs' conformational landscape during transcription, and may combine with detection methods other than smFRET for structural studies of RNAs in general.

Crystal structures of the NAD+-II riboswitch reveal two distinct ligand-binding pockets

We present crystal structures of a new NAD+-binding riboswitch termed NAD+-II, bound to nicotinamide mononucleotide (NMN), nicotinamide adenine dinucleotide (NAD+) and nicotinamide riboside (NR). The RNA structure comprises a number of structural features including three helices, one of which forms a triple helix by interacting with an A5 strand in its minor-groove, and another formed from a long-range pseudoknot. The core of the structure (centrally located and coaxial with the triplex and the pseudoknot) includes two consecutive quadruple base interactions. Unusually the riboswitch binds two molecules of ligand, bound at distinct, non-overlapping sites in the RNA. Binding occurs primarily through the nicotinamide moiety of each ligand, held by specific hydrogen bonding and stacking interactions with the pyridyl ring. The mode of binding is the same for NMN, NR and the nicotinamide moiety of NAD+. In addition, when NAD+ is bound into one site it adopts an elongated conformation such that its diphosphate linker occupies a groove on the surface of the RNA, following which the adenine portion inserts into a pocket and makes specific hydrogen bonding interactions. Thus the NAD+-II riboswitch is distinct from the NAD+-I riboswitch in that it binds two molecules of ligand at separate sites, and that binding occurs principally through the nicotinamide moiety.

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.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

Insights into the structural stability of major groove RNA triplexes by WAXS-guided MD simulations

RNA triple helices are commonly observed tertiary motifs that are associated with critical biological functions, including signal transduction. Because the recognition of their biological importance is relatively recent, their full range of structural properties has not yet been elucidated. The integration of solution wide-angle X-ray scattering (WAXS) with molecular dynamics (MD) simulations, described here, provides a new way to capture the structures of major-groove RNA triplexes that evade crystallographic characterization. This method yields excellent agreement between measured and computed WAXS profiles and allows for an atomically detailed visualization of these motifs. Using correlation maps, the relationship between well-defined features in the scattering profiles and real space characteristics of RNA molecules is defined, including the subtle conformational variations in the double-stranded RNA upon the incorporation of a third strand by base triples. This readily applicable approach has the potential to provide insight into interactions that stabilize RNA tertiary structure that enables function.

Subsite Ligand Recognition and Cooperativity in the TPP Riboswitch: Implications for Fragment-Linking in RNA Ligand Discovery

RNA molecules can show high levels of cooperativity in their global folding and interactions with divalent ions. However, cooperativity at individual ligand-RNA interaction sites remains poorly understood. Here, we investigated the binding of thiamine and methylene diphosphonic acid (MDP, a soluble structural analogue of pyrophosphate) to the thiamine pyrophosphate riboswitch. These ligands each bind weakly at proximal subsites, with 10 μM and 1 mM affinities, respectively. The affinity of MDP moderately improves when thiamine or thiamine-like fragments are pre-bound to the RNA. Covalent linking of thiamine and MDP substantially increases riboswitch binding to a notable high affinity of 20 nM. Crystal structures and single-molecule correlated chemical probing revealed favorable induced fit effects upon binding of individual ligands and, unexpectedly, a substantial thermodynamically unfavorable RNA structural rearrangement upon binding of the linked thiamine-MDP ligand. Thus, linking of two ligands of modest affinity, accompanied by an unfavorable structural rearrangement, still yields a potent linked RNA-binding compound. Since complex ligands often bind riboswitches and other RNAs at proximal subsites, principles derived from this work inform and support fragment-linking strategies for identifying small molecules that interact with RNA specifically and with high affinity.

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.

Charge Assisted S/Se Chalcogen Bonds in SAM Riboswitches: A Combined PDB and ab Initio Study

In this study, we provide experimental (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the involvement of charge assisted chalcogen bonding (ChB) interactions in the recognition and folding mechanisms of S-adenosylmethionine (SAM) riboswitches. Concretely, an initial PDB search revealed several examples where ChBs between S-adenosyl methionine (SAM)/adenosyl selenomethionine (EEM) molecules and uracil (U) bases belonging to RNA take place. While these interactions are usually described as a merely Coulombic attraction between the positively charged S/Se group and RNA, theoretical calculations indicated that the σ holes of S and Se are involved. Moreover, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader’s “Atoms in Molecules” (AIM) theory, Non-Covalent Interaction plot (NCIplot) visual index, and Natural Bonding Orbital (NBO) analyses. As far as our knowledge extends, this is the first time that ChBs in SAM–RNA complexes have been systematically analyzed, and we believe the results might be useful for scientists working in the field of RNA engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

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.

Siblings or doppelgängers? Deciphering the evolution of structured cis-regulatory RNAs beyond homology

Structured cis-regulatory RNAs have evolved across all domains of life, highlighting the utility and plasticity of RNA as a regulatory molecule. Homologous RNA sequences and structures often have similar functions, but homology may also be deceiving. The challenges that derive from trying to assign function to structure and vice versa are not trivial. Bacterial riboswitches, viral and eukaryotic IRESes, CITEs, and 3′ UTR elements employ an array of mechanisms to exert their effects. Bioinformatic searches coupled with biochemical and functional validation have elucidated some shared and many unique ways cis-regulators are employed in mRNA transcripts. As cis-regulatory RNAs are resolved in greater detail, it is increasingly apparent that shared homology can mask the full spectrum of mRNA cis-regulator functional diversity. Furthermore, similar functions may be obscured by lack of obvious sequence similarity. Thus looking beyond homology is crucial for furthering our understanding of RNA-based regulation.

Recent Insights into the Structure and Function of Mycobacterial Membrane Proteins Facilitated by Cryo-EM

Mycobacterium tuberculosis (Mtb) is one of the deadliest pathogens encountered by humanity. Over the decades, its characteristic membrane organization and composition have been understood. However, there is still limited structural information and mechanistic understanding of the constituent membrane proteins critical for drug discovery pipelines. Recent advances in single-particle cryo-electron microscopy and cryo-electron tomography have provided the much-needed impetus towards structure determination of several vital Mtb membrane proteins whose structures were inaccessible via X-ray crystallography and NMR. Important insights into membrane composition and organization have been gained via a combination of electron tomography and biochemical and biophysical assays. In addition, till the time of writing this review, 75 new structures of various Mtb proteins have been reported via single-particle cryo-EM. The information obtained from these structures has improved our understanding of the mechanisms of action of these proteins and the physiological pathways they are associated with. These structures have opened avenues for structure-based drug design and vaccine discovery programs that might help achieve global-TB control. This review describes the structural features of selected membrane proteins (type VII secretion systems, Rv1819c, Arabinosyltransferase, Fatty Acid Synthase, F-type ATP synthase, respiratory supercomplex, ClpP1P2 protease, ClpB disaggregase and SAM riboswitch), their involvement in physiological pathways, and possible use as a drug target. Tuberculosis is a deadly disease caused by Mycobacterium tuberculosis. The Cryo-EM and tomography have simplified the understanding of the mycobacterial membrane organization. Some proteins are located in the plasma membrane; some span the entire envelope, while some, like MspA, are located in the mycomembrane. Cryo-EM has made the study of such membrane proteins feasible.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

Imaginary Ribozymes

Biocatalysis is dominated by protein enzymes, and only a few classes of ribozymes are known to contribute to the task of promoting biochemical transformations. The RNA World theory encompasses the notion that earlier forms of life made use of a much greater diversity of ribozymes and other functional RNAs to guide complex metabolic states long before proteins had emerged in evolution. In recent years, the discoveries of various classes of ribozymes, riboswitches, and other noncoding RNAs in bacteria have provided additional support for the hypothesis that RNA molecules indeed have the catalytic competence to promote diverse chemical reactions without the aid of protein enzymes. Herein, some of the most striking observations made from examinations of natural riboswitches that bind small ligands are highlighted and used as a basis to imagine the characteristics and functions of long-extinct ribozymes from the RNA World.

Crystal structure and ligand-induced folding of the SAM/SAH riboswitch

While most SAM riboswitches strongly discriminate between SAM and SAH, the SAM/SAH riboswitch responds to both ligands with similar apparent affinities. We have determined crystal structures of the SAM/SAH riboswitch bound to SAH, SAM and other variant ligands at high resolution. The riboswitch forms an H-type pseudoknot structure with coaxial alignment of the stem–loop helix (P1) and the pseudoknot helix (PK). An additional three base pairs form at the non-open end of P1, and the ligand is bound at the interface between the P1 extension and the PK helix. The adenine nucleobase is stacked into the helix and forms a trans Hoogsteen–Watson–Crick base pair with a uridine, thus becoming an integral part of the helical structure. The majority of the specific interactions are formed with the adenosine. The methionine or homocysteine chain lies in the groove making a single hydrogen bond, and there is no discrimination between the sulfonium of SAM or the thioether of SAH. Single-molecule FRET analysis reveals that the riboswitch exists in two distinct conformations, and that addition of SAM or SAH shifts the population into a stable state that likely corresponds to the form observed in the crystal. A model for translational regulation is presented whereby in the absence of ligand the riboswitch is largely unfolded, lacking the PK helix so that translation can be initiated at the ribosome binding site. But the presence of ligand stabilizes the folded conformation that includes the PK helix, so occluding the ribosome binding site and thus preventing the initiation of translation.

Structure and ligand binding of the ADP-binding domain of the NAD+ riboswitch

The nadA motif is the first known NAD+-dependent riboswitch, comprising two similar tandem bulged stem-loop structures. We have determined the structure of the 5' domain 1 of the riboswitch. It has three coaxial helical segments, separated by an ACANCCCC bulge and by an internal loop, with a tertiary contact between them that includes two C:G base pairs. We have determined the structure with a number of ligands related to NADH, but in each case only the ADP moiety is observed. The adenosine adopts an anti conformation, forms multiple hydrogen bonds across the width of the sugar edge of the penultimate C:G base pair of the helix preceding the bulge, and the observed contacts have been confirmed by mutagenesis and calorimetry. Two divalent metal ions play a key structural role at the narrow neck of the bulge. One makes direct bonding contacts to the diphosphate moiety, locking it into position. Thus the nucleobase, ribose, and phosphate groups of the ADP moiety are all specifically recognized by the RNA. The NAD+ riboswitch is modular. Domain 1 is an ADP binding domain that may be ancient and could potentially be used in combination with other ligand binding motifs such as CoA.

Molecular structure of a U•A-U-rich RNA triple helix with 11 consecutive base triples

Three-dimensional structures have been solved for several naturally occurring RNA triple helices, although all are limited to six or fewer consecutive base triples, hindering accurate estimation of global and local structural parameters. We present an X-ray crystal structure of a right-handed, U•A-U-rich RNA triple helix with 11 continuous base triples. Due to helical unwinding, the RNA triple helix spans an average of 12 base triples per turn. The double helix portion of the RNA triple helix is more similar to both the helical and base step structural parameters of A′-RNA rather than A-RNA. Its most striking features are its wide and deep major groove, a smaller inclination angle and all three strands favoring a C3′-endo sugar pucker. Despite the presence of a third strand, the diameter of an RNA triple helix remains nearly identical to those of DNA and RNA double helices. Contrary to our previous modeling predictions, this structure demonstrates that an RNA triple helix is not limited in length to six consecutive base triples and that longer RNA triple helices may exist in nature. Our structure provides a starting point to establish structural parameters of the so-called ‘ideal’ RNA triple helix, analogous to A-RNA and B-DNA double helices.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

SAM-VI riboswitch structure and signature for ligand discrimination

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

Structure and ligand binding of the glutamine-II riboswitch

We have determined the structure of the glutamine-II riboswitch ligand binding domain using X-ray crystallography. The structure was solved using a novel combination of homology modeling and molecular replacement. The structure comprises three coaxial helical domains, the central one of which is a pseudoknot with partial triplex character. The major groove of this helix provides the binding site for L-glutamine, which is extensively hydrogen bonded to the RNA. Atomic mutation of the RNA at the ligand binding site leads to loss of binding shown by isothermal titration calorimetry, explaining the specificity of the riboswitch. A metal ion also plays an important role in ligand binding. This is directly bonded to a glutamine carboxylate oxygen atom, and its remaining inner-sphere water molecules make hydrogen bonding interactions with the RNA.

Cryo-EM structure of a 40 kDa SAM-IV riboswitch RNA at 3.7 Å resolution

Specimens below 50 kDa have generally been considered too small to be analyzed by single-particle cryo-electron microscopy (cryo-EM). The high flexibility of pure RNAs makes it difficult to obtain high-resolution structures by cryo-EM. In bacteria, riboswitches regulate sulfur metabolism through binding to the S-adenosylmethionine (SAM) ligand and offer compelling targets for new antibiotics. SAM-I, SAM-I/IV, and SAM-IV are the three most commonly found SAM riboswitches, but the structure of SAM-IV is still unknown. Here, we report the structures of apo and SAM-bound SAM-IV riboswitches (119-nt, ~40 kDa) to 3.7 Å and 4.1 Å resolution, respectively, using cryo-EM. The structures illustrate homologies in the ligand-binding core but distinct peripheral tertiary contacts in SAM-IV compared to SAM-I and SAM-I/IV. Our results demonstrate the feasibility of resolving small RNAs with enough detail to enable detection of their ligand-binding pockets and suggest that cryo-EM could play a role in structure-assisted drug design for RNA.

Stability of an RNA•DNA-DNA triple helix depends on base triplet composition and length of the RNA third strand

Recent studies suggest noncoding RNAs interact with genomic DNA, forming an RNA•DNA–DNA triple helix that regulates gene expression. However, base triplet composition of pyrimidine motif RNA•DNA–DNA triple helices is not well understood beyond the canonical U•A–T and C•G–C base triplets. Using native gel-shift assays, the relative stability of 16 different base triplets at a single position, Z•X–Y (where Z = C, U, A, G and X–Y = A–T, G–C, T–A, C–G), in an RNA•DNA–DNA triple helix was determined. The canonical U•A–T and C•G–C base triplets were the most stable, while three non-canonical base triplets completely disrupted triple-helix formation. We further show that our RNA•DNA–DNA triple helix can tolerate up to two consecutive non-canonical A•G–C base triplets. Additionally, the RNA third strand must be at least 19 nucleotides to form an RNA•DNA–DNA triple helix but increasing the length to 27 nucleotides does not increase stability. The relative stability of 16 different base triplets in DNA•DNA–DNA and RNA•RNA–RNA triple helices was distinctly different from those in RNA•DNA–DNA triple helices, showing that base triplet stability depends on strand composition being DNA and/or RNA. Multiple factors influence the stability of triple helices, emphasizing the importance of experimentally validating formation of computationally predicted triple helices.

Limits in accuracy and a strategy of RNA structure prediction using experimental information

RNA structural complexity and flexibility present a challenge for computational modeling efforts. Experimental information and bioinformatics data can be used as restraints to improve the accuracy of RNA tertiary structure prediction. Regarding utilization of restraints, the fundamental questions are: (i) What is the limit in prediction accuracy that one can achieve with arbitrary number of restraints? (ii) Is there a strategy for selection of the minimal number of restraints that would result in the best structural model? We address the first question by testing the limits in prediction accuracy using native contacts as restraints. To address the second question, we develop an algorithm based on the distance variation allowed by secondary structure (DVASS), which ranks restraints according to their importance to RNA tertiary structure prediction. We find that due to kinetic traps, the greatest improvement in the structure prediction accuracy is achieved when we utilize only 40-60% of the total number of native contacts as restraints. When the restraints are sorted by DVASS algorithm, using only the first 20% ranked restraints can greatly improve the prediction accuracy. Our findings suggest that only a limited number of strategically selected distance restraints can significantly assist in RNA structure modeling.