MD simulations of ligand-bound and ligand-free aptamer: molecular level insights into the binding and switching mechanism of the add A-riboswitch

Computer aided aptamer selection for fabrication of electrochemical sensor to detect Aflatoxin B1

Aflatoxin B1 (AFB1) is a naturally occurring toxin produced by Aspergillus flavus and Aspergillus parasiticus. The AFB1 is classified as a potent carcinogen and poses significant health risks both to humans and animals. Early detection of the toxin in post-harvest agricultural products will save lives and promote healthy food production. In this study, stratified docking approach was utilized to screen and identify potential aptamers that can bind to AFB1. ssDNA sequences were acquired from the Mendeley dataset, secondary and tertiary structures were predicted through a series of bioinformatics pipelines. Further, the final DNA tertiary structures were minimized and SiteMap algorithm was used to probe and locate binding cavities. According to the final XP docking result, a 34 nt sequence (5'-ATCCTGTGAGGAATGCTCATGCATAGCAAGGGCT-3') aptamer with a docking score of -5.959 kcal/mol was considered for 200 ns MD Simulation. Finally, the screened DNA-aptamer was immobilized over the gold surface based on Au-S chemistry and utilized for the detection of AFB1. The fabricated DNA-aptamer electrode demonstrated a good analytical performance including wide linear range (1.0 to 1000 ng L-1), detection limit (1.0 ng L-1), high stability, and reproducibility.

Exploring the Binding Process of Cognate Ligand to Add Adenine Riboswitch Aptamer by Using Explicit Solvent Molecular Dynamics (MD) Simulation

Riboswitches are RNA-structured elements that modulate gene expression by changing their conformation in response to specific metabolite ligand binding. Therefore, the biological functions of riboswitches mainly depend on the switching of secondary and three-dimensional structures in the presence and absence of the metabolite ligands. However, the binding mechanisms of cognate ligands to riboswitches are still not well understood. Here, we have introduced how to use explicit solvent molecular dynamics (MD) simulation to observe the binding process of cognate ligand to add adenine riboswitch aptamer at the atomic level. In addition, we have analyzed the driving factors of the binding process and calculated the binding free energy based on the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method.

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.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.

Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands

Guanidinium ion is a toxic cellular metabolite. The ykkC-III riboswitch, an mRNA stretch, regulates the gene expression by undergoing a conformational change in response to the binding of a free guanidinium ion and thereby plays a potentially important role in alleviating guanidinium toxicity in cells. An experimental crystal structure of the guanidinium-bound aptamer domain of the riboswitch from Thermobifida Fusca revealed the overall RNA architecture and mapped the specific noncovalent interactions that stabilize the ligand within the binding pocket aptamer. However, details of how the aptamer domain discriminates the cognate ligand from its closest structurally analogous physiological metabolites (arginine and urea), and how the binding of cognate ligand arrays information from the aptamer domain to the expression platform for regulating the gene expression, are not well understood. To fill this void, we perform a cumulative of 2 μs all-atom explicit-solvent molecular dynamics (MD) simulations on the full aptamer domain, augmented with quantum-chemical calculations on the ligand-binding pocket, to compare the structural and dynamical details of the guanidinium-bound state with the arginine or urea bound states, as well as the unbound (open) state. Analysis of the ligand-binding pocket reveals that due to unfavorable interactions with the binding-pocket residues, urea cannot bind the aptamer domain and thereby cannot alter the gene expression. Although interaction of the guanidyl moiety of arginine within the binding pocket is either comparable or stronger than the guanidinium ion, additional non-native hydrogen-bonding networks, as well as differences in the dynamical details of the arginine-bound state, explain why arginine cannot transmit the information from the aptamer domain to the expression platform. Based on our simulations, we propose a mechanism of how the aptamer domain communicates with the expression platform. Overall, our work provides interesting insights into the ligand recognition by a specific class of riboswitches and may hopefully inspire future studies to further understand the gene regulation by riboswitches.

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.

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.

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.

Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: crystal structure analysis, classification, optimal structures, and accurate interaction energies

In the present work, 67 crystal structures of the aptamer domains of RNA riboswitches are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of 80 unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson-Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate, and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that 67 out of 80 pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal a rich hydrogen bonding landscape, ranging from weak interactions (base:other, -3 kcal/mol) to strong (base:phosphate, -48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in the context of RNA:ligand recognition. This can, in turn, contribute toward current efforts to develop antimicrobials that target RNAs.

Structural insights into the interactions of flavin mononucleotide (FMN) and riboflavin with FMN riboswitch: a molecular dynamics simulation study

Antibiotics resistance is becoming a serious problem associated with fatalities and suffering patients. New antibiotics that can target the broader spectrum of cellular processes are warranted. One of the recent approaches in this regard is to target the special type of RNA riboswitches in bacteria. In this report, we have explored the mechanistic pathways of ligand-dependent conformational changes of flavin mononucleotide (FMN) riboswitch using molecular dynamics (MD) simulation studies. Cognate ligands FMN and riboflavin (RBF) have shown very different behavior with FMN riboswitch in terms of their role in the gene regulation process. These two ligands have similar scaffold, except the terminal phosphate group in FMN ligand. The MD simulations reveal that the binding of FMN ligand with the riboswitch does not lead to global folding of structure, rather lead to local changes in riboswitch structure. The binding free energy calculated with molecular mechanics Poisson-Boltzmann surface area method suggests the stronger binding of FMN than RBF to the riboswitch and electrostatic energy contributes chiefly to stabilize the complex. Further, the hydrogen bonding analysis identified the key binding site residues G11, G32, G62 of the riboswitch with FMN and RBF. The critical role of the phosphate group in the FMN ligand for binding with the active site of a riboswitch is also borne out in this study. These results unravel the importance of functional groups in natural ligands on designing newer ligands for FMN riboswitch as new antibiotics in the future.Communicated by Ramaswamy H. Sarma.

Conformational flexibility of adenine riboswitch aptamer in apo and bound states using NMR and an X-ray free electron laser

Riboswitches are structured cis-regulators mainly found in the untranslated regions of messenger RNA. The aptamer domain of a riboswitch serves as a sensor for its ligand, the binding of which triggers conformational changes that regulate the behavior of its expression platform. As a model system for understanding riboswitch structures and functions, the add adenine riboswitch has been studied extensively. However, there is a need for further investigation of the conformational dynamics of the aptamer in light of the recent real-time crystallographic study at room temperature (RT) using an X-ray free electron laser (XFEL) and femtosecond X-ray crystallography (SFX). Herein, we investigate the conformational motions of the add adenine riboswitch aptamer domain, in the presence or absence of adenine, using nuclear magnetic resonance relaxation measurements and analysis of RT atomic displacement factors (B-factors). In the absence of ligand, the P1 duplex undergoes a fast exchange where the overall molecule exhibits a motion at k_ex ~ 319 s^-1, based on imino signals. In the presence of ligand, the P1 duplex adopts a highly ordered conformation, with k_ex~ 83 s^-1, similar to the global motion of the molecule, excluding the loops and binding pocket, at 84 s^-1. The µs-ms motions in both the apo and bound states are consistent with RT B-factors. Reduced spatial atomic fluctuation, ~ 50%, in P1 upon ligand binding coincides with significantly attenuated temporal dynamic exchanges. The binding pocket is structured in the absence or presence of ligand, as evidenced by relatively low and similar RT B-factors. Therefore, despite the dramatic rearrangement of the binding pocket, those residues exhibit similar spatial thermal fluctuation before and after binding.

Going beyond base-pairs: topology-based characterization of base-multiplets in RNA

Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.

TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg²⁺ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

Conformational Dynamics of thiM Riboswitch To Understand the Gene Regulation Mechanism Using Markov State Modeling and the Residual Fluctuation Network Approach

Thiamine pyrophosphate (TPP) riboswitch is a cis-regulatory element in the noncoding region of mRNA. The aptamer domain of TPP riboswitch detects the high abundance of coenzyme thiamine pyrophosphate (TPP) and modulates the gene expression for thiamine synthetic gene. The mechanistic understanding in recognition of TPP in aptamer domain and ligand-induced compactness for folding of expression platform are most important to designing novel modulators. To understand the dynamic behavior of TPP riboswitch upon TPP binding, molecular dynamics simulations were performed for 400 ns in both apo and TPP bound forms of thiM riboswitch from E. coli and analyzed in terms of eRMSD-based Markov state modeling and residual fluctuation network. Markov state models show good correlations in transition probability among metastable states from simulated trajectory and generated models. Structural compactness in TPP bound form is observed which is correlated with SAXS experiment. The importance of junction of P4 and P5 is evident during dynamics, which correlates with FRET analysis. The dynamic nature of two sensor forearms is due to the flexible P1 helix, which is its intrinsic property. The transient state in TPP-bound form was observed in the Markov state model, along with stable states. We believe that this transient state is responsible to assist the influx and outflux of ligand molecule by creating a solvent channel around the junction region of P4 and P5 and such a structure was anticipated in FRET analysis. The dynamic nature of riboswitch is dependent on the interaction between residues on distal loops L3 and L5/P3 and junction P4 and P5, J3/2 which stabilize the J2/4. It helps in the transfer of allosteric information between J2/4 and P3/L5 tertiary docking region through the active site residues. Understanding such information flow will benefit in highlighting crucial residues in highly dynamic and kinetic systems. Here, we report the residues and segments in riboswitch that play vital roles in providing stability and this can be exploited in designing inhibitors to regulate the functioning of riboswitches.

Effect of PDGF-B aptamer on PDGFRβ/PDGF-B interaction: Molecular dynamics study

PDGFRβ/PDGF-B interaction plays a role in angiogenesis, and is mandatory in wound healing and cancer treatment. It has been reported that the PDGF-B aptamer was able to bind to PDGF-B, thus regulating the angiogenesis. However, the binding interaction between the aptamer and the growth factor, including the binding sites, has not been well investigated. This study applied a molecular dynamics (MD) simulation to investigate the aptamer-growth factor interaction in the presence or absence of a receptor (PDGFRβ). Characterization of the structure of an aptamer-growth factor complex revealed binding sites from each section in the complex. Upon the complex formation, PDGF-B and its aptamer exhibited less flexibility in their molecular movement, as indicated by the minimum values of RMSD, RMSF, loop-to-loop distance, and the summation of PCA eigenvalues. Our study of residue pairwise interaction demonstrated that the binding interaction was mainly contributed by electrostatic interaction between the positively-charged amino acid and the negatively-charged phosphate backbone. The role of the PDGF-B aptamer in PDGFRβ/PDGF-B interaction was also investigated. We demonstrated that the stability of the Apt-PDGF-B complex could prevent the presence of a competitor, of PDGFRβ, interrupting the binding process. Because the aptamer was capable of binding with PDGF-B, and blocking the growth factor from the PDGFRβ, it could down regulate the consequent signaling pathway. We provide evidence that the PDGF-BB aptamer is a promising molecule for regulation of angiogenesis. The MD study provides a molecular understanding to modification of the aptamer binding interaction, which could be used in a number of medical applications.

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.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches

Riboswitches are genetic control elements within non-coding regions of mRNA. These self-regulatory elements have been found to sense a range of small metabolites, ions, and other physical signals to exert regulatory control of transcription, translation, and splicing. To date, more than a dozen riboswitch classes have been characterized that vary widely in size and secondary structure. Extensive experiments and theoretical studies have made great strides in understanding the general structures, genetic mechanisms, and regulatory activities of individual riboswitches. As the ligand-dependent co-transcriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolites under the transcription. This review will provide a brief summary of the studies about the regulation mechanisms of the pbuE, S_MK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the riboswitches.

Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch

Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.

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.

Approach to the unfolding and folding dynamics of add A-riboswitch upon adenine dissociation using a coarse-grained elastic network model

Riboswitches are noncoding mRNA segments that can regulate the gene expression via altering their structures in response to specific metabolite binding. We proposed a coarse-grained Gaussian network model (GNM) to examine the unfolding and folding dynamics of adenosine deaminase (add) A-riboswitch upon the adenine dissociation, in which the RNA is modeled by a nucleotide chain with interaction networks formed by connecting adjoining atomic contacts. It was shown that the adenine binding is critical to the folding of the add A-riboswitch while the removal of the ligand can result in drastic increase of the thermodynamic fluctuations especially in the junction regions between helix domains. Under the assumption that the native contacts with the highest thermodynamic fluctuations break first, the iterative GNM simulations showed that the unfolding process of the adenine-free add A-riboswitch starts with the denature of the terminal helix stem, followed by the loops and junctions involving ligand binding pocket, and then the central helix domains. Despite the simplified coarse-grained modeling, the unfolding dynamics and pathways are shown in close agreement with the results from atomic-level MD simulations and the NMR and single-molecule force spectroscopy experiments. Overall, the study demonstrates a new avenue to investigate the binding and folding dynamics of add A-riboswitch molecule which can be readily extended for other RNA molecules.

Ligand-Induced Stabilization of a Duplex-like Architecture Is Crucial for the Switching Mechanism of the SAM-III Riboswitch

Riboswitches are structured RNA motifs that control gene expression by sensing the concentrations of specific metabolites and make up a promising new class of antibiotic targets. S-Adenosylmethionine (SAM)-III riboswitch, mainly found in lactic acid bacteria, is involved in regulating methionine and SAM biosynthetic pathways. SAM-III riboswitch regulates the gene expression by switching the translation process on and off with respect to the absence and presence of the SAM ligand, respectively. In this study, an attempt is made to understand the key conformational transitions involved in ligand binding using atomistic molecular dynamics (MD) simulations performed in an explicit solvent environment. G26 is found to recognize the SAM ligand by forming hydrogen bonds, whereas the absence of the ligand leads to opening of the binding pocket. Consistent with experimental results, the absence of the SAM ligand weakens the base pairing interactions between the nucleobases that are part of the Shine-Dalgarno (SD) and anti-Shine-Dalgarno (aSD) sequences, which in turn facilitates recognition of the SD sequence by ribosomes. Detailed analysis reveals that a duplex-like structure formed by nucleotides from different parts of the RNA and the adenine base of the ligand is crucial for the stability of the completely folded state in the presence of the ligand. Previous experimental studies have shown that the SAM-III riboswitch exists in equilibrium between the unfolded and partially folded states in the absence of the ligand, which completely folds upon binding of the ligand. Comparison of the results presented here to the available experimental data indicates the structures obtained using the MD simulations resemble the partially folded state. Thus, this study provides a detailed understanding of the fully and partially folded structures of the SAM-III riboswitch in the presence and absence of the ligand, respectively. This study hypothesizes a dual role for the SAM ligand, which facilitates conformational switching between partially and fully folded states by forming a stable duplex-like structure and strengthening the interactions between SD and aSD nucleotides.

Biophysical Approaches to Bacterial Gene Regulation by Riboswitches

The last decade has witnessed the discovery of a variety of non-coding RNA sequences that perform a broad range of crucial biological functions. Among these, the ability of certain RNA sequences, so-called riboswitches, has attracted considerable interest. Riboswitches control gene expression in response to the concentration of particular metabolites to which they bind without the need for any protein. These RNA switches not only need to adopt a very specific tridimensional structure to perform their function, but also their sequence has been evolutionary optimized to recognize a particular metabolite with high affinity and selectivity. Thus, riboswitches offer a unique opportunity to get fundamental insights into RNA plasticity and how folding dynamics and ligand recognition mechanisms have been efficiently merged to control gene regulation. Because riboswitch sequences have been mostly found in bacterial organisms controlling the expression of genes associated to the synthesis, degradation or transport of crucial metabolites for bacterial survival, they offer exciting new routes for antibiotic development in an era where bacterial resistance is more than ever challenging conventional drug discovery strategies. Here, we give an overview of the architecture, diversity and regulatory mechanisms employed by riboswitches with particular emphasis on the biophysical methods currently available to characterise their structure and functional dynamics.

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 preQ1 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.

Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations

INTRODUCTION

Riboswitches are cis-acting regulatory RNA elements prevalently located in the leader sequences of bacterial mRNA. An adenine sensing riboswitch cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The structural mechanism regulating its conformational changes upon ligand binding mostly remains to be elucidated. In this open framework it has been suggested that the ligand stabilizes the interaction of the distal "kissing loop" complex. Using accurate full-atom molecular dynamics with explicit solvent in combination with enhanced sampling techniques and advanced analysis methods it could be possible to provide a more detailed perspective on the formation of these tertiary contacts.

METHODS

In this work, we used umbrella sampling simulations to study the thermodynamics of the kissing loop complex in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction restraining the distance between the two loops. We also assessed the convergence of the results by using two alternative initialization protocols. A structural analysis was performed using a novel approach to analyze base contacts.

RESULTS

Contacts between the two loops were progressively lost when larger inter-loop distances were enforced. Inter-loop Watson-Crick contacts survived at larger separation when compared with non-canonical pairing and stacking interactions. Intra-loop stacking contacts remained formed upon loop undocking. Our simulations qualitatively indicated that the ligand could stabilize the kissing loop complex. We also compared with previously published simulation studies.

DISCUSSION AND CONCLUSIONS

Kissing complex stabilization given by the ligand was compatible with available experimental data. However, the dependence of its value on the initialization protocol of the umbrella sampling simulations posed some questions on the quantitative interpretation of the results and called for better converged enhanced sampling simulations.

RNABP COGEST: a resource for investigating functional RNAs

Structural bioinformatics of RNA has evolved mainly in response to the rapidly accumulating evidence that non-(protein)-coding RNAs (ncRNAs) play critical roles in gene regulation and development. The structures and functions of most ncRNAs are however still unknown. Most of the available RNA structural databases rely heavily on known 3D structures, and contextually correlate base pairing geometry with actual 3D RNA structures. None of the databases provide any direct information about stabilization energies. However, the intrinsic interaction energies of constituent base pairs can provide significant insights into their roles in the overall dynamics of RNA motifs and structures. Quantum mechanical (QM) computations provide the only approach toward their accurate quantification and characterization. ‘RNA Base Pair Count, Geometry and Stability’ (http://bioinf.iiit.ac.in/RNABPCOGEST) brings together information, extracted from literature data, regarding occurrence frequency, experimental and quantum chemically optimized geometries, and computed interaction energies, for non-canonical base pairs observed in a non-redundant dataset of functional RNA structures. The database is designed to enable the QM community, on the one hand, to identify appropriate biologically relevant model systems and also enable the biology community to easily sift through diverse computational results to gain theoretical insights which could promote hypothesis driven biological research.

Database URL:http://bioinf.iiit.ac.in/RNABPCOGEST

Force field dependence of riboswitch dynamics

Riboswitches are noncoding regulatory elements that control gene expression in response to the presence of metabolites, which bind to the aptamer domain. Metabolite binding appears to occur through a combination of conformational selection and induced fit mechanism. This demands to characterize the structural dynamics of the apo state of aptamer domains. In principle, molecular dynamics (MD) simulations can give insights at the atomistic level into the dynamics of the aptamer domain. However, it is unclear to what extent contemporary force fields can bias such insights. Here, we show that the Amber force field ff99 yields the best agreement with detailed experimental observations on differences in the structural dynamics of wild type and mutant aptamer domains of the guanine-sensing riboswitch (Gsw), including a pronounced influence of Mg2+. In contrast, applying ff99 with parmbsc0 and parmχOL modifications (denoted ff10) results in strongly damped motions and overly stable tertiary loop-loop interactions. These results are based on 58 MD simulations with an aggregate simulation time>11 μs, careful modeling of Mg2+ ions, and thorough statistical testing. Our results suggest that the moderate stabilization of the χ-anti region in ff10 can have an unwanted damping effect on functionally relevant structural dynamics of marginally stable RNA systems. This suggestion is supported by crystal structure analyses of Gsw aptamer domains that reveal χ torsions with high-anti values in the most mobile regions. We expect that future RNA force field development will benefit from considering marginally stable RNA systems and optimization toward good representations of dynamics in addition to structural characteristics.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

Binding site preorganization and ligand discrimination in the purine riboswitch

The progress of RNA research has suggested a wide variety of RNA molecules as possible targets for pharmaceutical drug molecules. Structure-based computational methods for predicting binding modes and affinities are now important tools in drug discovery, but these methods have mainly been focused on protein targets. Here we employ molecular dynamics free-energy perturbation calculations and the linear interaction energy method to analyze the energetics of ligand binding to purine riboswitches. Calculations are carried out for 14 different purine complexes with the guanine and adenine riboswitches in order to examine their ligand recognition principles. The simulations yield binding affinities in good agreement with experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that these receptors have an unusually high degree of electrostatic preorganization for their cognate ligands, and this effect is further quantified by explicit free-energy calculations, which show that the standard electrostatic linear interaction energy parametrization is suboptimal in this case. The adenine riboswitch specifically uses the electrostatic preorganization to discriminate against guanine by preventing the formation of a G-U wobble base pair.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

The importance of helix P1 stability for structural pre-organization and ligand binding affinity of the adenine riboswitch aptamer domain

We report here an in-depth characterization of the aptamer domain of the transcriptional adenine-sensing riboswitch (pbuE) by NMR and fluorescence spectroscopy. By NMR studies, the structure of two aptamer sequences with different lengths of the helix P1, the central element involved in riboswitch conformational switching, was characterized. Hydrogen-bond interactions could be mapped at nucleotide resolution providing information about secondary and tertiary structure, structure homogeneity and dynamics. Our study reveals that the elongation of helix P1 has pronounced effects not only on the local but on the global structure of the apo aptamer domain. The structural differences induced by stabilizing helix P1 were found to be linked to changes of the ligand binding affinity as revealed from analysis of kinetic and thermodynamic data obtained from stopped-flow fluorescence studies. The results provide new insight into the sequence-dependent fine tuning of the structure and function of purine-sensing riboswitches.

DFT and MD study of the divalent-cation-mediated interaction of ochratoxin A with DNA nucleosides

Aptamers are ligand-binding nucleic acids with affinities and selectivities that make them useful for the detection of a variety of compounds, including ochratoxin A. Theoretical methods can be applied to study the recognition interaction between aptamers and the ochratoxin A molecule. In this work, molecular dynamics simulations and quantum chemical calculations performed at the DFT level of theory were used to study the structures and energies of aptamers and aptamer-ochratoxin A complexes. The optimal structures as well as the interaction energies of these structures were elucidated. Divalent cations in the water solvent were shown to be an important influence on the structures and stabilities of the complexes.

Fluorescence tools to investigate riboswitch structural dynamics

Riboswitches are novel regulatory elements that respond to cellular metabolites to control gene expression. They are constituted of highly conserved domains that have evolved to recognize specific metabolites. Such domains, so-called aptamers, are folded into intricate structures to enable metabolite recognition. Over the years, the development of ensemble and single-molecule fluorescence techniques has allowed to probe most of the mechanistic aspects of aptamer folding and ligand binding. In this review, we summarize the current fluorescence toolkit available to study riboswitch structural dynamics. We fist describe those methods based on fluorescent nucleotide analogues, mostly 2-aminopurine (2AP), to investigate short-range conformational changes, including some key steady-state and time-resolved examples that exemplify the versatility of fluorescent analogues as structural probes. The study of long-range structural changes by Förster resonance energy transfer (FRET) is mostly discussed in the context of single-molecule studies, including some recent developments based on the combination of single-molecule FRET techniques with controlled chemical denaturation methods. This article is part of a Special Issue entitled: Riboswitches.

Insights into ligand binding to PreQ1 Riboswitch Aptamer from molecular dynamics simulations

Riboswitches play roles in transcriptional or translational regulation through specific ligand binding of their aptamer domains. Although a number of ligand-bound aptamer complex structures have been solved, it is important to know ligand-free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, preQ₁ riboswitch aptamer domain from Bacillus subtilis is studied by overall 1.5 μs all-atom molecular dynamics simulations We found that the ligand-free aptamer has a stable state with a folded P1-L3 and open binding pocket. The latter forms a cytosine-rich pool in which the nucleotide C19 oscillates between close and open positions, making it a potential conformation for preQ₁ entrance. The dynamic picture further suggests that the specific recognition of preQ₁ by the aptamer domain is not only facilitated by the key nucleotide C19 but also aided and enhanced by other cytosines around the binding pocket. These results should help to understand the details of preQ₁ binding.

Methods to enable the design of bioactive small molecules targeting RNA

RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.

Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch

Riboswitches are structured mRNA elements that modulate gene expression. They undergo conformational changes triggered by highly specific interactions with sensed metabolites. Among the structural rearrangements engaged by riboswitches, the forming and melting of the aptamer terminal helix, the so-called P1 stem, is essential for genetic control. The structural mechanisms by which this conformational change is modulated upon ligand binding mostly remain to be elucidated. Here, we used pulling molecular dynamics simulations to study the thermodynamics of the P1 stem in the add adenine riboswitch. The P1 ligand-dependent stabilization was quantified in terms of free energy and compared with thermodynamic data. This comparison suggests a model for the aptamer folding in which direct P1-ligand interactions play a minor role on the conformational switch when compared with those related to the ligand-induced aptamer preorganization.

RNA structure and dynamics: a base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.

The impact of a ligand binding on strand migration in the SAM-I riboswitch

Riboswitches sense cellular concentrations of small molecules and use this information to adjust synthesis rates of related metabolites. Riboswitches include an aptamer domain to detect the ligand and an expression platform to control gene expression. Previous structural studies of riboswitches largely focused on aptamers, truncating the expression domain to suppress conformational switching. To link ligand/aptamer binding to conformational switching, we constructed models of an S-adenosyl methionine (SAM)-I riboswitch RNA segment incorporating elements of the expression platform, allowing formation of an antiterminator (AT) helix. Using Anton, a computer specially developed for long timescale Molecular Dynamics (MD), we simulated an extended (three microseconds) MD trajectory with SAM bound to a modeled riboswitch RNA segment. Remarkably, we observed a strand migration, converting three base pairs from an antiterminator (AT) helix, characteristic of the transcription ON state, to a P1 helix, characteristic of the OFF state. This conformational switching towards the OFF state is observed only in the presence of SAM. Among seven extended trajectories with three starting structures, the presence of SAM enhances the trend towards the OFF state for two out of three starting structures tested. Our simulation provides a visual demonstration of how a small molecule (<500 MW) binding to a limited surface can trigger a large scale conformational rearrangement in a 40 kDa RNA by perturbing the Free Energy Landscape. Such a mechanism can explain minimal requirements for SAM binding and transcription termination for SAM-I riboswitches previously reported experimentally.

Single-molecule chemical denaturation of riboswitches

To date, single-molecule RNA science has been developed almost exclusively around the effect of metal ions as folding promoters and stabilizers of the RNA structure. Here, we introduce a novel strategy that combines single-molecule Förster resonance energy transfer (FRET) and chemical denaturation to observe and manipulate RNA dynamics. We demonstrate that the competing interplay between metal ions and denaturant agents provides a platform to extract information that otherwise will remain hidden with current methods. Using the adenine-sensing riboswitch aptamer as a model, we provide strong evidence for a rate-limiting folding step of the aptamer domain being modulated through ligand binding, a feature that is important for regulation of the controlled gene. In the absence of ligand, the rate-determining step is dominated by the formation of long-range key tertiary contacts between peripheral stem-loop elements. In contrast, when the adenine ligand interacts with partially folded messenger RNAs, the aptamer requires specifically bound Mg²⁺ ions, as those observed in the crystal structure, to progress further towards the native form. Moreover, despite that the ligand-free and ligand-bound states are indistinguishable by FRET, their different stability against urea-induced denaturation allowed us to discriminate them, even when they coexist within a single FRET trajectory; a feature not accessible by existing methods.

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

Dynamic energy landscapes of riboswitches help interpret conformational rearrangements and function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way that sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we mimic transcription by computing energy landscapes which describe accessible secondary structures for a range of sequence lengths. Consistent with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: Low-barrier landscape with an ensemble of conformations in equilibrium before encountering a substrate; barrier-free landscape with a dominant “downhill” pathway to the minimum free energy structure; and barrier-dominated landscape with two isolated conformational states with different functions. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved between riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Role of ligand binding in structural organization of add A-riboswitch aptamer: a molecular dynamics simulation

The specific binding of ligands is the first step of gene expression or translation regulation by riboswitches. However, understanding the mechanism of the specific binding is still difficult because the tertiary structures of the riboswitch aptamers are available almost only for ligand-bound state at present. In this paper we hope to give some insights into this problem through the studies of the role of ligand-aptamer interaction in the structural organization of add A-riboswitch aptamer, based on the crystal structure of the ligand-bound aptamer. We use all-atom molecular dynamics to simulate the behaviors of the aptamer in ligand-bound, free and mutated states by Amber force field. The results show that the correct paring of the ligand adenine with the nucleotide U74 in the binding pocket is crucial to stabilizing the conformations of the ligand-bound aptamer, especially the helix P1 connecting the expression platform. Our results also suggest that both the nucleotide U74 and U51 may be the key sites of the ligand recognition but the former has much higher probability as the initial docking site. This is in agreement with previous experimental results.

Atomic-scale characterization of conformational changes in the preQ₁ riboswitch aptamer upon ligand binding

Riboswitches are mRNA structural elements that act as intracellular sensors of small-molecule metabolites. By undergoing conformational changes capable of modulating translation or terminating transcription, riboswitches are able to play a role in regulating the concentration of essential metabolites in the cell. Computer-guided fluorescence experiments were carried out to interrogate molecular dynamics and conformational changes in the minimal riboswitch aptamer that binds 7-aminomethyl-7-deazaguanine (preQ₁). Our combined experimental results and computational analysis suggest that the preQ₁ riboswitch apo form is structured but shows no evidence of a ligand-binding pocket. Simulations of the apo and bound forms indicate a large conformational change is triggered by the breaking of the Watson-Crick base pairing of nucleotides G11 and C31 upon preQ₁ removal, followed by collapse of the pocket due to interfering π-stacking. Computational predictions of local aptamer dynamics were validated by fluorescence experiments employing 2-aminopurine substitutions. In-line probing reactions confirmed that fluorophore-labeled riboswitches retain similar higher-order structural features as the unlabeled aptamer upon ligand binding, although their affinity for the ligand is reduced by the introduction of the fluorescent reporter.