1
|
Hien EDM, Chauvier A, St-Pierre P, Lafontaine DA. Structural Characterization of the Cotranscriptional Folding of the Thiamin Pyrophosphate Sensing thiC Riboswitch in Escherichia coli. Biochemistry 2024; 63:1608-1620. [PMID: 38864595 DOI: 10.1021/acs.biochem.3c00665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Riboswitches are RNA-regulating elements that mostly rely on structural changes to modulate gene expression at various levels. Recent studies have revealed that riboswitches may control several regulatory mechanisms cotranscriptionally, i.e., during the transcription elongation of the riboswitch or early in the coding region of the regulated gene. Here, we study the structure of the nascent thiamin pyrophosphate (TPP)-sensing thiC riboswitch in Escherichia coli by using biochemical and enzymatic conventional probing approaches. Our chemical (in-line and lead probing) and enzymatic (nucleases S1, A, T1, and RNase H) probing data provide a comprehensive model of how TPP binding modulates the structure of the thiC riboswitch. Furthermore, by using transcriptional roadblocks along the riboswitch sequence, we find that a certain portion of nascent RNA is needed to sense TPP that coincides with the formation of the P5 stem loop. Together, our data suggest that conventional techniques may readily be used to study cotranscriptional folding of nascent RNAs.
Collapse
Affiliation(s)
- Elsa D M Hien
- Department of Biology, Faculty of Science, RNA Group, Université de Sherbrooke, Sherbrooke, Quebec J1K 2R1, Canada
| | - Adrien Chauvier
- Department of Biology, Faculty of Science, RNA Group, Université de Sherbrooke, Sherbrooke, Quebec J1K 2R1, Canada
| | - Patrick St-Pierre
- Department of Biology, Faculty of Science, RNA Group, Université de Sherbrooke, Sherbrooke, Quebec J1K 2R1, Canada
| | - Daniel A Lafontaine
- Department of Biology, Faculty of Science, RNA Group, Université de Sherbrooke, Sherbrooke, Quebec J1K 2R1, Canada
| |
Collapse
|
2
|
Zheng L, Song Q, Xu X, Shen X, Li C, Li H, Chen H, Ren A. Structure-based insights into recognition and regulation of SAM-sensing riboswitches. SCIENCE CHINA. LIFE SCIENCES 2023; 66:31-50. [PMID: 36459353 DOI: 10.1007/s11427-022-2188-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2022] [Accepted: 07/17/2022] [Indexed: 12/03/2022]
Abstract
Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.
Collapse
Affiliation(s)
- Luqian Zheng
- Department of Gastroenterology, Department of Cardiology of the Second Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou, 310058, China.,Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China.,The Eighth Affiliated Hospital, Sun Yat-sen University, Shenzhen, 518033, China
| | - Qianqian Song
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Xiaochen Xu
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Xin Shen
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Chunyan Li
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Hongcheng Li
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Hao Chen
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Aiming Ren
- Department of Gastroenterology, Department of Cardiology of the Second Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou, 310058, China. .,Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China.
| |
Collapse
|
3
|
Yan S, Ilgu M, Nilsen-Hamilton M, Lamm MH. Computational Modeling of RNA Aptamers: Structure Prediction of the Apo State. J Phys Chem B 2022; 126:7114-7125. [PMID: 36097649 PMCID: PMC9512008 DOI: 10.1021/acs.jpcb.2c04649] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 08/23/2022] [Indexed: 11/28/2022]
Abstract
RNA aptamers are single-stranded oligonucleotides that bind to specific molecular targets with high affinity and specificity. To design aptamers for new applications, it is critical to understand the ligand binding mechanism in terms of the structure and dynamics of the ligand-bound and apo states. The problem is that most of the NMR or X-ray crystal structures available for RNA aptamers are for ligand-bound states. Available apo state structures, mostly characterized by crystallization under nonphysiological conditions or probed by low resolution techniques, might fail to represent the diverse structural variations of the apo state in solution. Here, we develop an approach to obtain a representative ensemble of apo structures that are based on in silico RNA 3D structure prediction and in vitro experiments that characterize base stacking. Using the neomycin-B aptamer as a case study, an ensemble of structures for the aptamer in the apo (unbound) state are validated and then used to investigate the ligand-binding mechanism for the aptamer in complex with neomycin-B.
Collapse
Affiliation(s)
- Shuting Yan
- Iowa
State University, Ames, Iowa 50011, United States
| | - Muslum Ilgu
- Iowa
State University, Ames, Iowa 50011, United States
- Ames
National Laboratory, Ames, Iowa 50011, United States
- Aptalogic
Inc., Ames, Iowa 50014, United States
| | - Marit Nilsen-Hamilton
- Iowa
State University, Ames, Iowa 50011, United States
- Ames
National Laboratory, Ames, Iowa 50011, United States
- Aptalogic
Inc., Ames, Iowa 50014, United States
| | | |
Collapse
|
4
|
Kumar S, Reddy G. TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg 2+ Concentration. J Phys Chem B 2022; 126:2369-2381. [PMID: 35298161 DOI: 10.1021/acs.jpcb.1c10794] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg2+. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg2+ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg2+ concentration ([Mg2+]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg2+], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg2+ aids in the formation of tertiary contacts. For [Mg2+] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.
Collapse
Affiliation(s)
- Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| |
Collapse
|
5
|
Negi I, Mahmi AS, Seelam Prabhakar P, Sharma P. Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands. J Chem Inf Model 2021; 61:5243-5255. [PMID: 34609872 DOI: 10.1021/acs.jcim.1c01022] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- Indu Negi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Amanpreet Singh Mahmi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Preethi Seelam Prabhakar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| |
Collapse
|
6
|
Ding J, Swain M, Yu P, Stagno JR, Wang YX. Conformational flexibility of adenine riboswitch aptamer in apo and bound states using NMR and an X-ray free electron laser. JOURNAL OF BIOMOLECULAR NMR 2019; 73:509-518. [PMID: 31606878 PMCID: PMC6817744 DOI: 10.1007/s10858-019-00278-w] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 09/12/2019] [Indexed: 06/10/2023]
Abstract
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 kex ~ 319 s-1, based on imino signals. In the presence of ligand, the P1 duplex adopts a highly ordered conformation, with kex~ 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.
Collapse
Affiliation(s)
- Jienv Ding
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, National Institute of Health, Frederick, MD, 21702, USA.
| | - Monalisa Swain
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, National Institute of Health, Frederick, MD, 21702, USA
| | - Ping Yu
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, National Institute of Health, Frederick, MD, 21702, USA
| | - Jason R Stagno
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, National Institute of Health, Frederick, MD, 21702, USA
| | - Yun-Xing Wang
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, National Institute of Health, Frederick, MD, 21702, USA.
| |
Collapse
|
7
|
Yonkunas MJ, Baird NJ. A highly ordered, nonprotective MALAT1 ENE structure is adopted prior to triplex formation. RNA (NEW YORK, N.Y.) 2019; 25:975-984. [PMID: 31113838 PMCID: PMC6633196 DOI: 10.1261/rna.069906.118] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2018] [Accepted: 05/13/2019] [Indexed: 05/06/2023]
Abstract
The 3' end of the ∼7 kb lncRNA MALAT1 contains an evolutionarily and structurally conserved element for nuclear expression (ENE) which confers protection from cellular degradation pathways. Formation of an ENE triple helix is required to support transcript accumulation, leading to persistent oncogenic activity of MALAT1 in multiple cancer types. Though the specific mechanism of triplex-mediated protection remains unknown, the MALAT1 ENE triplex has been identified as a promising target for therapeutic intervention. Interestingly, a maturation step of the nascent lncRNA 3' end is required prior to triplex formation. We hypothesize that disruption of the maturation or folding process may be a viable mechanism of inhibition. To assess putative cotranscriptional ENE conformations prior to triplex formation, we perform microsecond MD simulations of a partially folded ENE conformation and the ENE triplex. We identify a highly ordered ENE structure prior to triplex formation. Extensive formation of U•U base pairs within the large U-rich internal loops produces a global rod-like architecture. We present a three-dimensional structure of the isolated ENE motif, the global features of which are consistent with small angle X-ray scattering (SAXS) experiments. Our structural model represents a nonprotective conformation of the MALAT1 ENE, providing a molecular description useful for future mechanistic and inhibition studies. We anticipate that targeting stretches of U•U pairs within the ENE motif will prove advantageous for the design of therapeutics targeting this oncogenic lncRNA.
Collapse
Affiliation(s)
- Michael J Yonkunas
- Department of Chemistry & Biochemistry, University of the Sciences, Philadelphia, Pennsylvania 19104, USA
| | - Nathan J Baird
- Department of Chemistry & Biochemistry, University of the Sciences, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
8
|
Padhi S, Pradhan M, Bung N, Roy A, Bulusu G. TPP riboswitch aptamer: Role of Mg 2+ ions, ligand unbinding, and allostery. J Mol Graph Model 2019; 88:282-291. [PMID: 30818079 DOI: 10.1016/j.jmgm.2019.01.015] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Revised: 01/18/2019] [Accepted: 01/18/2019] [Indexed: 01/23/2023]
Abstract
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. Mg2+ 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.
Collapse
Affiliation(s)
- Siladitya Padhi
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited, Hyderabad, 500081, India
| | - Meenakshi Pradhan
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited, Hyderabad, 500081, India
| | - Navneet Bung
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited, Hyderabad, 500081, India
| | - Arijit Roy
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited, Hyderabad, 500081, India
| | - Gopalakrishnan Bulusu
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited, Hyderabad, 500081, India.
| |
Collapse
|
9
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|