Folding of the SAM-I riboswitch: A tale with a twist

Riboswitches are ligand-dependent RNA genetic regulators that control gene expression by altering their structures. The elucidation of riboswitch conformational changes before and after ligand recognition is crucial to understand how riboswitches can achieve high ligand binding affinity and discrimination against cellular analogs. The detailed characterization of riboswitch folding pathways suggest that they may use their intrinsic conformational dynamics to sample a large array of structures, some of which being nearly identical to ligand-bound molecules. Some of these structural conformers can be "captured" upon ligand binding, which is crucial for the outcome of gene regulation. Recent studies about the SAM-I riboswitch have revealed unexpected and previously unknown RNA folding mechanisms. For instance, the observed helical twist of the P1 stem upon ligand binding to the SAM-I aptamer adds a new element in the repertoire of RNA strategies for recognition of small metabolites. From an RNA folding perspective, these findings also strongly indicate that the SAM-I riboswitch could achieve ligand recognition by using an optimized combination of conformational capture and induced-fit approaches, a feature that may be shared by other RNA regulatory sequences.     

Ligand-induced conformational capture of a synthetic tetracycline riboswitch revealed by pulse EPR

RNA aptamers are in vitro-selected binding domains that recognize their respective ligand with high affinity and specificity. They are characterized by complex three-dimensional conformations providing preformed binding pockets that undergo conformational changes upon ligand binding. Small molecule-binding aptamers have been exploited as synthetic riboswitches for conditional gene expression in various organisms. In the present study, double electron-electron resonance (DEER) spectroscopy combined with site-directed spin labeling was used to elucidate the conformational transition of a tetracycline aptamer upon ligand binding. Different sites were selected for post-synthetic introduction of either the (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl) methanethiosulfonate by reaction with a 4-thiouridine modified RNA or of 4-isocyanato-2,6-tetramethylpiperidyl-N-oxid spin label by reaction with 2'-aminouridine modified RNA. The results of the DEER experiments indicate the presence of a thermodynamic equilibrium between two aptamer conformations in the free state and capture of one conformation upon tetracycline binding.     

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA''s secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.     

Structural probing and damage selection of citrulline- and arginine-specific RNA aptamers identify base positions required for binding.

In a recent study, an RNA aptamer for the specific recognition of the amino acid L-arginine was evolved from an in vitro selected L-citrulline binding parent sequence [M. Famulok (1994) J. Am. Chem. Soc. 116, 1698-1706]. We have now carried out a structural analysis of these aptamers by using chemical modification experiments. Footprinting experiments and a damage selection approach were performed to identify those positions protected from modification in the presence of the amino acids and modifications that interfere with the binding of the ligand. It is shown that of the two bulged regions present in both aptamers one can be modified without loss of binding activity whereas in the other bulge nearly every position is shown to be involved in the recognition of the ligands. This might be indicative for non-canonical base pairing to occur within the non-Watson-Crick paired regions which might be stabilized by the complexed amino acid. Binding to the cognate amino acid significantly enhances the conformational stability of the RNA. We also tested the sensitivity of both aptamers towards lead (II) ion induced cleavage and identified a hypersensitive cleavage site within the invariant bulged region. Lead cleavage is inhibited by the complexed amino acid, indicating a conformational change of the aptamer upon ligand binding. NMR titration data obtained with both aptamers and their cognate ligands confirm the proposed conformational changes and indicate the formation of a 1:1 complex of RNA:amino acid.     

On the link between conformational changes,ligand binding and heat capacity

BackgroundConformational changes coupled to ligand binding constitute the structural and energetics basis underlying cooperativity, allostery and, in general, protein regulation. These conformational rearrangements are associated with heat capacity changes. ITC is a unique technique for studying binding interactions because of the simultaneous determination of the binding affinity and enthalpy, and for providing the best estimates of binding heat capacity changes.Scope of reviewStill controversial issues in ligand binding are the discrimination between the “conformational selection model” and the “induced fit model”, and whether or not conformational changes lead to temperature dependent apparent binding heat capacities. The assessment of conformational changes associated with ligand binding by ITC is discussed. In addition, the “conformational selection” and “induced fit” models are reconciled, and discussed within the context of intrinsically (partially) unstructured proteins.Major conclusionsConformational equilibrium is a major contribution to binding heat capacity changes. A simple model may explain both conformational selection and induced fit scenarios. A temperature-independent binding heat capacity does not necessarily indicate absence of conformational changes upon ligand binding. ITC provides information on the energetics of conformational changes associated with ligand binding (and other possible additional coupled equilibria).General significancePreferential ligand binding to certain protein states leads to an equilibrium shift that is reflected in the coupling between ligand binding and additional equilibria. This represents the structural/energetic basis of the widespread dependence of ligand binding parameters on temperature, as well as pH, ionic strength and the concentration of other chemical species. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Elucidating the Role of Microprocessor Protein DGCR8 in Bending RNA Structures

Although conformational dynamics of RNA molecules are potentially important in microRNA (miRNA) processing, the role of the protein binding partners in facilitating the requisite structural changes is not well understood. In previous work, we and others have demonstrated that nonduplex structural elements and the conformational flexibility they support are necessary for efficient RNA binding and cleavage by the proteins associated with the two major stages of miRNA processing. However, recent studies showed that the protein DGCR8 binds primary miRNA and duplex RNA with similar affinities. Here, we study RNA binding by a small recombinant construct of the DGCR8 protein and the RNA conformation changes that result. This construct, the DGCR8 core, contains two double-stranded RNA-binding domains (dsRBDs) and a C-terminal tail. To assess conformational changes resulting from binding, we applied small-angle x-ray scattering with contrast variation to detect conformational changes of primary-miR-16-1 in complex with the DGCR8 core. This method reports only on the RNA conformation within the complex and suggests that the protein bends the RNA upon binding. Supporting work using smFRET to study the conformation of RNA duplexes bound to the core also shows bending. Together, these studies elucidate the role of DGCR8 in interacting with RNA during the early stages of miRNA processing.     

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.     

Riboswitch Detection Using Profile Hidden Markov Models

Background  

Riboswitches are a type of noncoding RNA that regulate gene expression by switching from one structural conformation to another on ligand binding. The various classes of riboswitches discovered so far are differentiated by the ligand, which on binding induces a conformational switch. Every class of riboswitch is characterized by an aptamer domain, which provides the site for ligand binding, and an expression platform that undergoes conformational change on ligand binding. The sequence and structure of the aptamer domain is highly conserved in riboswitches belonging to the same class. We propose a method for fast and accurate identification of riboswitches using profile Hidden Markov Models (pHMM). Our method exploits the high degree of sequence conservation that characterizes the aptamer domain.     

A theophylline responsive riboswitch based on helix slipping controls gene expression in vivo

Riboswitches are newly discovered regulatory elements which control a wide set of basic metabolic pathways. They consist solely of RNA, sense their ligand in a preformed binding pocket and perform a conformational switch in response to ligand binding resulting in altered gene expression. We have utilized the enormous potential of RNA for molecular sensing and conformational changes to develop novel molecular switches with predetermined structural transitions in response to the binding of a small molecule. To validate these in vivo, we exploit the distance-dependent inhibitory potential of secondary structure elements placed close to the bacterial ribosome binding site. We created a translational control element by combining the theophylline aptamer with a helical communication module for which a ligand-dependent one-nucleotide slipping mechanism had been proposed. This structural element was inserted at a position just interfering with translation in the non ligand-bound form. Addition of the ligand then shifts the inhibitory element to a distance which permits efficient translation. We present here a novel regulatory mechanism in the first rationally designed, in vivo active RNA switch. Its use of a slippage mechanism to control gene expression makes it different from natural riboswitches which are based on sequestration or antitermination.     

Purposely engineered drug-target mismatches for entropy-based drug optimization

Proteins are dynamic objects that often undergo significant structural change and reduce their conformational possibilities upon ligand binding. Thus, unless dynamic information is incorporated, structure-based drug design becomes of limited applicability. Even within a dynamic approach, a rarely visited scenario arises as proteins increase their entropy content upon ligand binding by locally enhancing conformational exploration in the complex. In this opinion piece, we argue that this binding mode is of primary importance in drug development because it allows for drugs that are not optimized in the conventional way but feature mismatches with the target. Thus, we advocate entropy optimization that exploits dynamic information for drug design.     

Conformational dynamics and thermodynamics of protein-ligand binding studied by NMR relaxation

Protein conformational dynamics can be critical for ligand binding in two ways that relate to kinetics and thermodynamics respectively. First, conformational transitions between different substates can control access to the binding site (kinetics). Secondly, differences between free and ligand-bound states in their conformational fluctuations contribute to the entropy of ligand binding (thermodynamics). In the present paper, I focus on the second topic, summarizing our recent results on the role of conformational entropy in ligand binding to Gal3C (the carbohydrate-recognition domain of galectin-3). NMR relaxation experiments provide a unique probe of conformational entropy by characterizing bond-vector fluctuations at atomic resolution. By monitoring differences between the free and ligand-bound states in their backbone and side chain order parameters, we have estimated the contributions from conformational entropy to the free energy of binding. Overall, the conformational entropy of Gal3C increases upon ligand binding, thereby contributing favourably to the binding affinity. Comparisons with the results from isothermal titration calorimetry indicate that the conformational entropy is comparable in magnitude to the enthalpy of binding. Furthermore, there are significant differences in the dynamic response to binding of different ligands, despite the fact that the protein structure is virtually identical in the different protein-ligand complexes. Thus both affinity and specificity of ligand binding to Gal3C appear to depend in part on subtle differences in the conformational fluctuations that reflect the complex interplay between structure, dynamics and ligand interactions.     

π–π Interactions in Structural Stability: Role in RNA Binding Proteins

RNA binding proteins play significant roles in many bio-macromolecular systems. Aromatic amino acid residues are vital for several biological functions. In the present work, the influences of π–π interactions in RNA binding proteins are analyzed. There are a total of 3,396 π-residues in RNA binding proteins out of which 1,547, 1,241, and 608 are phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp), respectively. Among these 945, 634, and 356 Phe, Tyr, and Trp residues, respectively, are involved in π–π interactions. The observations indicate that majority of the aromatic residues in RNA binding proteins are involved in π–π interactions. Side chain–side chain π–π interactions are the predominant type of interactions in RNA binding proteins. These π–π interactions stabilize the core regions within RNA binding proteins. π–π interacting residues are evolutionary conserved. Residue-wise analysis indicates that π–π interacting residues have higher long-range contacts and hence they are important in the global conformational stability of these proteins.     

Trans-acting glmS catalytic riboswitch: locked and loaded

A recently discovered class of gene regulatory RNAs, coined riboswitches, are commonly found in noncoding segments of bacterial and some eukaryotic mRNAs. Gene up- or down-regulation is triggered by binding of a small organic metabolite, which typically induces an RNA conformational change. Unique among these noncoding RNAs is the glmS catalytic riboswitch, or ribozyme, found in the 5'-untranslated region of the glmS gene in Gram-positive bacteria. It is activated by glucosamine-6-phosphate (GlcN6P), leading to site-specific backbone cleavage of the mRNA and subsequent repression of the glmS gene, responsible for cellular GlcN6P production. Recent biochemical and structural evidence suggests that the GlcN6P ligand acts as a coenzyme and participates in the cleavage reaction without inducing a conformational change. To better understand the role of GlcN6P in solution structural dynamics and function, we have separated the glmS riboswitch core from Bacillus subtilis into a trans-cleaving ribozyme and an externally cleaved substrate. We find that trans cleavage is rapidly activated by nearly 5000-fold to a rate of 4.4 min(-1) upon addition of 10 mM GlcN6P, comparable to the cis-acting ribozyme. Fluorescence resonance energy transfer suggests that this ribozyme-substrate complex does not undergo a global conformational change upon ligand binding in solution. In addition, footprinting at nucleotide resolution using terbium(III) and RNase V1 indicates no significant changes in secondary and tertiary structure upon ligand binding. These findings suggest that the glmS ribozyme is fully folded in solution prior to binding its activating ligand, supporting recent observations in the crystalline state.     

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.