Single-molecule force spectroscopy of the add adenine riboswitch relates folding to regulatory mechanism

Riboswitches regulate gene expression via ligand binding to an aptamer domain which induces conformational changes in a regulatory expression platform. By unfolding and refolding single add adenine riboswitch molecules in an optical trap, an integrated picture of the folding was developed and related to the regulatory mechanism. Force-extension curves (FECs) and constant-force folding trajectories measured on the aptamer alone revealed multiple partially-folded states, including several misfolded states not on the native folding pathway. All states were correlated to key structural components and interactions within hierarchical folding pathways. FECs of the full-length riboswitch revealed that the thermodynamically stable conformation switches upon ligand binding from a structure repressing translation to one permitting it. Along with rapid equilibration of the two structures in the absence of adenine, these results support a thermodynamically-controlled regulatory mechanism, in contrast with the kinetic control of the closely-related pbuE adenine riboswitch. Comparison of the folding of these riboswitches revealed many similarities arising from shared structural features but also essential differences related to their different regulatory mechanisms.     

Mg2+-induced conformational changes in the btuB riboswitch from E. coli

Mg²⁺ has been shown to modulate the function of riboswitches by facilitating the ligand-riboswitch interactions. The btuB riboswitch from Escherichia coli undergoes a conformational change upon binding to its ligand, coenzyme B₁₂ (adenosyl-cobalamine, AdoCbl), and down-regulates the expression of the B₁₂ transporter protein BtuB in order to control the cellular levels of AdoCbl. Here, we discuss the structural folding attained by the btuB riboswitch from E. coli in response to Mg²⁺ and how it affects the ligand binding competent conformation of the RNA. The btuB riboswitch notably adopts different conformational states depending upon the concentration of Mg²⁺. With the help of in-line probing, we show the existence of at least two specific conformations, one being achieved in the complete absence of Mg²⁺ (or low Mg²⁺ concentration) and the other appearing above ∼0.5 mM Mg²⁺. Distinct regions of the riboswitch exhibit different dissociation constants toward Mg²⁺, indicating a stepwise folding of the btuB RNA. Increasing the Mg²⁺ concentration drives the transition from one conformation toward the other. The conformational state existing above 0.5 mM Mg²⁺ defines the binding competent conformation of the btuB riboswitch which can productively interact with the ligand, coenzyme B₁₂, and switch the RNA conformation. Moreover, raising the Mg²⁺ concentration enhances the ratio of switched RNA in the presence of AdoCbl. The lack of a AdoCbl-induced conformational switch experienced by the btuB riboswitch in the absence of Mg²⁺ indicates a crucial role played by Mg²⁺ for defining an active conformation of the riboswitch.     

Riboswitch structure: an internal residue mimicking the purine ligand

The adenine and guanine riboswitches regulate gene expression in response to their purine ligand. X-ray structures of the aptamer moiety of these riboswitches are characterized by a compact fold in which the ligand forms a Watson–Crick base pair with residue 65. Phylogenetic analyses revealed a strict restriction at position 39 of the aptamer that prevents the G39–C65 and A39–U65 combinations, and mutational studies indicate that aptamers with these sequence combinations are impaired for ligand binding. In order to investigate the rationale for sequence conservation at residue 39, structural characterization of the U65C mutant from Bacillus subtilis pbuE adenine riboswitch aptamer was undertaken. NMR spectroscopy and X-ray crystallography studies demonstrate that the U65C mutant adopts a compact ligand-free structure, in which G39 occupies the ligand-binding site of purine riboswitch aptamers. These studies present a remarkable example of a mutant RNA aptamer that adopts a native-like fold by means of ligand mimicking and explain why this mutant is impaired for ligand binding. Furthermore, this work provides a specific insight into how the natural sequence has evolved through selection of nucleotide identities that contribute to formation of the ligand-bound state, but ensures that the ligand-free state remains in an active conformation.     

Molecular mechanism for preQ1-II riboswitch function revealed by molecular dynamics

Riboswitches are RNA molecules that regulate gene expression using conformational change, affected by binding of small molecule ligands. A crystal structure of a ligand-bound class II preQ₁ riboswitch has been determined in a previous structural study. To gain insight into the dynamics of this riboswitch in solution, eight total molecular dynamic simulations, four with and four without ligand, were performed using the Amber force field. In the presence of ligand, all four of the simulations demonstrated rearranged base pairs at the 3′ end, consistent with expected base-pairing from comparative sequence analysis in a prior bioinformatic analysis; this suggests the pairing in this region was altered by crystallization. Additionally, in the absence of ligand, three of the simulations demonstrated similar changes in base-pairing at the ligand binding site. Significantly, although most of the riboswitch architecture remained intact in the respective trajectories, the P3 stem was destabilized in the ligand-free simulations in a way that exposed the Shine–Dalgarno sequence. This work illustrates how destabilization of two major groove base triples can influence a nearby H-type pseudoknot and provides a mechanism for control of gene expression by a fold that is frequently found in bacterial riboswitches.     

Direct evidence for RNA–RNA interactions at the 3′ end of the Hepatitis C virus genome using surface plasmon resonance

Surface plasmon resonance was used to investigate two previously described interactions analyzed by reverse genetics and complementation mutation experiments, involving 5BSL3.2, a stem–loop located in the NS5B coding region of HCV. 5BSL3.2 was immobilized on a sensor chip by streptavidin-biotin coupling, and its interaction either with the SL2 stem–loop of the 3′ end or with an upstream sequence centered on nucleotide 9110 (referred to as Seq9110) was monitored in real-time. In contrast with previous results obtained by NMR assays with the same short RNA sequences that we used or SHAPE analysis with longer RNAs, we demonstrate that recognition between 5BSL3.2 and SL2 can occur in solution through a kissing-loop interaction. We show that recognition between Seq9110 and the internal loop of 5BSL3.2 does not prevent binding of SL2 on the apical loop of 5BSL3.2 and does not influence the rate constants of the SL2-5BSL3.2 complex. Therefore, the two binding sites of 5BSL3.2, the apical and internal loops, are structurally independent and both interactions can coexist. We finally show that the stem–loop SL2 is a highly dynamic RNA motif that fluctuates between at least two conformations: One is able to hybridize with 5BSL3.2 through loop–loop interaction, and the other one is capable of self-associating in the absence of protein, reinforcing the hypothesis of SL2 being a dimerization sequence. This result suggests also that the conformational dynamics of SL2 could play a crucial role for controlling the destiny of the genomic RNA.     

The structure of NoRC-associated RNA is crucial for targeting the chromatin remodelling complex NoRC to the nucleolus

Silencing of ribosomal RNA genes (rDNA) requires binding of the chromatin remodelling complex NoRC to RNA that is complementary to the rDNA promoter. NoRC-associated RNA (pRNA) folds into a conserved stem–loop structure that is required for nucleolar localization and rDNA silencing. Mutations that disrupt the stem–loop structure impair binding of TIP5, the large subunit of NoRC, to pRNA and abolish targeting of NoRC to nucleoli. Binding to pRNA results in a conformational change of TIP5, as shown by enhanced sensitivity of TIP5 towards trypsin digestion. Our results indicate an RNA-dependent mechanism that targets NoRC to chromatin and facilitates the interaction with co-repressors that promote heterochromatin formation and rDNA silencing.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Ligand binding by the tandem glycine riboswitch depends on aptamer dimerization but not double ligand occupancy

The glycine riboswitch predominantly exists as a tandem structure, with two adjacent, homologous ligand-binding domains (aptamers), followed by a single expression platform. The recent identification of a leader helix, the inclusion of which eliminates cooperativity between the aptamers, has reopened the debate over the purpose of the tandem structure of the glycine riboswitch. An equilibrium dialysis-based assay was combined with binding-site mutations to monitor glycine binding in each ligand-binding site independently to understand the role of each aptamer in glycine binding and riboswitch tertiary interactions. A series of mutations disrupting the dimer interface was used to probe how dimerization impacts ligand binding by the tandem glycine riboswitch. While the wild-type tandem riboswitch binds two glycine equivalents, one for each aptamer, both individual aptamers are capable of binding glycine when the other aptamer is unoccupied. Intriguingly, glycine binding by aptamer-1 is more sensitive to dimerization than glycine binding by aptamer-2 in the context of the tandem riboswitch. However, monomeric aptamer-2 shows dramatically weakened glycine-binding affinity. In addition, dimerization of the two aptamers in trans is dependent on glycine binding in at least one aptamer. We propose a revised model for tandem riboswitch function that is consistent with these results, wherein ligand binding in aptamer-1 is linked to aptamer dimerization and stabilizes the P1 stem of aptamer-2, which controls the expression platform.     

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.     

RNA stem–loop enhanced expression of previously non-expressible genes

The key step in bacterial translation is formation of the pre-initiation complex. This requires initial contacts between mRNA, fMet-tRNA and the 30S subunit of the ribosome, steps that limit the initiation of translation. Here we report a method for improving translational initiation, which allows expression of several previously non-expressible genes. This method has potential applications in heterologous protein synthesis and high-throughput expression systems. We introduced a synthetic RNA stem–loop (stem length, 7 bp; ΔG₀ = –9.9 kcal/mol) in front of various gene sequences. In each case, the stem–loop was inserted 15 nt downstream from the start codon. Insertion of the stem–loop allowed in vitro expression of five previously non-expressible genes and enhanced the expression of all other genes investigated. Analysis of the RNA structure proved that the stem–loop was formed in vitro, and demonstrated that stabilization of the ribosome binding site is due to stem–loop introduction. By theoretical RNA structure analysis we showed that the inserted RNA stem–loop suppresses long-range interactions between the translation initiation domain and gene-specific mRNA sequences. Thus the inserted RNA stem–loop supports the formation of a separate translational initiation domain, which is more accessible to ribosome binding.     

Conformational capture of the SAM-II riboswitch

Riboswitches are gene regulation elements in mRNA that function by specifically responding to metabolites. Although the metabolite-bound states of riboswitches have proven amenable to structure determination efforts, knowledge of the structural features of riboswitches in their ligand-free forms and their ligand-response mechanisms giving rise to regulatory control is lacking. Here we explore the ligand-induced folding process of the S-adenosylmethionine type II (SAM-II) riboswitch using chemical and biophysical methods, including NMR and fluorescence spectroscopy, and single-molecule fluorescence imaging. The data reveal that the unliganded SAM-II riboswitch is dynamic in nature, in that its stem-loop element becomes engaged in a pseudoknot fold through base-pairing with nucleosides in the 3' overhang containing the Shine-Dalgarno sequence. Although the pseudoknot structure is highly transient in the absence of its ligand, S-adenosylmethionine (SAM), it becomes conformationally restrained upon ligand recognition, through a conformational capture mechanism. These insights provide a molecular understanding of riboswitch dynamics that shed new light on the mechanism of riboswitch-mediated translational regulation.     

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.     

Role of the Adenine Ligand on the Stabilization of the Secondary and Tertiary Interactions in the Adenine Riboswitch

Riboswitches are RNA-based genetic control elements that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket. To facilitate an atomic detail interpretation of experimental investigations on the role of the adenine ligand on the conformational properties and kinetics of folding of the add adenine riboswitch, we performed molecular dynamics simulations in both the presence and the absence of the ligand. In the absence of ligand, structural deviations were observed in the J23 junction and the P1 stem. Destabilization of the P1 stem in the absence of ligand involves the loss of direct stabilizing interactions with the ligand, with additional contributions from the J23 junction region. The J23 junction of the riboswitch is found to be more flexible, and the tertiary contacts among the junction regions are altered in the absence of the adenine ligand; results suggest that the adenine ligand associates and dissociates from the riboswitch in the vicinity of J23. Good agreement was obtained with the experimental data with the results indicating dynamic behavior of the adenine ligand on the nanosecond time scale to be associated with the dynamic behavior of hydrogen bonding with the riboswitch. Results also predict that direct interactions of the adenine ligand with U74 of the riboswitch are not essential for stable binding although it is crucial for its recognition. The possibility of methodological artifacts and force-field inaccuracies impacting the present observations was checked by additional molecular dynamics simulations in the presence of 2,6-diaminopurine and in the crystal environment.     

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.     

Structural and biophysical characterization of the tandem substrate-binding domains of the ABC importer GlnPQ

The ATP-binding cassette transporter GlnPQ is an essential uptake system that transports glutamine, glutamic acid and asparagine in Gram-positive bacteria. It features two extra-cytoplasmic substrate-binding domains (SBDs) that are linked in tandem to the transmembrane domain of the transporter. The two SBDs differ in their ligand specificities, binding affinities and their distance to the transmembrane domain. Here, we elucidate the effects of the tandem arrangement of the domains on the biochemical, biophysical and structural properties of the protein. For this, we determined the crystal structure of the ligand-free tandem SBD1-2 protein from Lactococcus lactis in the absence of the transporter and compared the tandem to the isolated SBDs. We also used isothermal titration calorimetry to determine the ligand-binding affinity of the SBDs and single-molecule Förster resonance energy transfer (smFRET) to relate ligand binding to conformational changes in each of the domains of the tandem. We show that substrate binding and conformational changes are not notably affected by the presence of the adjoining domain in the wild-type protein, and changes only occur when the linker between the domains is shortened. In a proof-of-concept experiment, we combine smFRET with protein-induced fluorescence enhancement (PIFE–FRET) and show that a decrease in SBD linker length is observed as a linear increase in donor-brightness for SBD2 while we can still monitor the conformational states (open/closed) of SBD1. These results demonstrate the feasibility of PIFE–FRET to monitor protein–protein interactions and conformational states simultaneously.     

Ligand-dependent folding of the three-way junction in the purine riboswitch

Riboswitches are highly structured cis-acting elements located in the 5'-untranslated region of messenger RNAs that directly bind small molecule metabolites to regulate gene expression. Structural and biochemical studies have revealed riboswitches experience significant ligand-dependent conformational changes that are coupled to regulation. To monitor the coupling of ligand binding and RNA folding within the aptamer domain of the purine riboswitch, we have chemically probed the RNA with N-methylisatoic anhydride (NMIA) over a broad temperature range. Analysis of the temperature-dependent reactivity of the RNA in the presence and absence of hypoxanthine reveals that a limited set of nucleotides within the binding pocket change their conformation in response to ligand binding. Our data demonstrate that a distal loop-loop interaction serves to restrict the conformational freedom of a significant portion of the three-way junction, thereby promoting ligand binding under physiological conditions.     

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.     

Conformational stability of the epidermal growth factor (EGF) receptor as influenced by glycosylation,dimerization and EGF hormone binding

The epidermal growth factor receptor (EGFR) is an important transmembrane glycoprotein kinase involved the initiation or perpetuation of signal transduction cascades within cells. These processes occur after EGFR binds to a ligand [epidermal growth factor (EGF)], thus inducing its dimerization and tyrosine autophosphorylation. Previous publications have highlighted the importance of glycosylation and dimerization for promoting proper function of the receptor and conformation in membranes; however, the effects of these associations on the protein conformational stability have not yet been described. Molecular dynamics simulations were performed to characterize the conformational preferences of the monomeric and dimeric forms of the EGFR extracellular domain upon binding to EGF in the presence and absence of N‐glycan moieties. Structural stability analyses revealed that EGF provides the most conformational stability to EGFR, followed by glycosylation and dimerization, respectively. The findings also support that EGF–EGFR binding takes place through a large‐scale induced‐fitting mechanism. Proteins 2017; 85:561–570. © 2016 Wiley Periodicals, Inc.