1
|
Hu G, Zhang Y, Yu Z, Cui T, Cui W. Dynamical characterization and multiple unbinding paths of two PreQ 1 ligands in one pocket. Phys Chem Chem Phys 2023; 25:24004-24015. [PMID: 37646322 DOI: 10.1039/d3cp03142j] [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: 09/01/2023]
Abstract
Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China.
- Laoling People's Hospital, Dezhou 253600, China
| | | | - Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Tiejun Cui
- Laoling People's Hospital, Dezhou 253600, China
| | - Wanling Cui
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China.
- Laoling People's Hospital, Dezhou 253600, China
| |
Collapse
|
2
|
Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch. Commun Biol 2022; 5:1120. [PMID: 36273041 PMCID: PMC9588036 DOI: 10.1038/s42003-022-04096-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 10/11/2022] [Indexed: 11/09/2022] Open
Abstract
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 Mg2+ 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 (CJ12-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 CJ12-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.
Collapse
|
3
|
Zhou Y, Jiang Y, Chen SJ. RNA-ligand molecular docking: advances and challenges. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 12:e1571. [PMID: 37293430 PMCID: PMC10250017 DOI: 10.1002/wcms.1571] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/20/2021] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yuanzhe Zhou
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Yangwei Jiang
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Shi-Jie Chen
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| |
Collapse
|
4
|
Zhao J, Kennedy SD, Turner DH. Nuclear Magnetic Resonance Spectra and AMBER OL3 and ROC-RNA Simulations of UCUCGU Reveal Force Field Strengths and Weaknesses for Single-Stranded RNA. J Chem Theory Comput 2022; 18:1241-1254. [PMID: 34990548 DOI: 10.1021/acs.jctc.1c00643] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Single-stranded regions of RNA are important for folding of sequences into 3D structures and for design of therapeutics targeting RNA. Prediction of ensembles of 3D structures for single-stranded regions often involves classical mechanical approximations of interactions defined by quantum mechanical calculations on small model systems. Nuclear magnetic resonance (NMR) spectra and molecular dynamics (MD) simulations of short single strands provide tests for how well the approximations model many of the interactions. Here, the NMR spectra for UCUCGU at 2, 15, and 30 °C are compared to simulations with the AMBER force fields, OL3 and ROC-RNA. This is the first such comparison to an oligoribonucleotide containing an internal guanosine nucleotide (G). G is particularly interesting because of its many H-bonding groups, large dipole moment, and proclivity for both syn and anti conformations. Results reveal formation of a G amino to phosphate non-bridging oxygen H-bond. The results also demonstrate dramatic differences in details of the predicted structures. The variations emphasize the dependence of predictions on individual parameters and their balance with the rest of the force field. The NMR data can serve as a benchmark for future force fields.
Collapse
|
5
|
Hu G, Zhou HX. Binding free energy decomposition and multiple unbinding paths of buried ligands in a PreQ1 riboswitch. PLoS Comput Biol 2021; 17:e1009603. [PMID: 34767553 PMCID: PMC8612554 DOI: 10.1371/journal.pcbi.1009603] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 11/24/2021] [Accepted: 10/31/2021] [Indexed: 11/18/2022] Open
Abstract
Riboswitches are naturally occurring RNA elements that control bacterial gene expression by binding to specific small molecules. They serve as important models for RNA-small molecule recognition and have also become a novel class of targets for developing antibiotics. Here, we carried out conventional and enhanced-sampling molecular dynamics (MD) simulations, totaling 153.5 μs, to characterize the determinants of binding free energies and unbinding paths for the cognate and synthetic ligands of a PreQ1 riboswitch. Binding free energy analysis showed that two triplets of nucleotides, U6-C15-A29 and G5-G11-C16, contribute the most to the binding of the cognate ligands, by hydrogen bonding and by base stacking, respectively. Mg2+ ions are essential in stabilizing the binding pocket. For the synthetic ligands, the hydrogen-bonding contributions of the U6-C15-A29 triplet are significantly compromised, and the bound state resembles the apo state in several respects, including the disengagement of the C15-A14-A13 and A32-G33 base stacks. The bulkier synthetic ligands lead to significantly loosening of the binding pocket, including extrusion of the C15 nucleobase and a widening of the C15-C30 groove. Enhanced-sampling simulations further revealed that the cognate and synthetic ligands unbind in almost opposite directions. Our work offers new insight for designing riboswitch ligands. Riboswitches are bacterial RNA elements that change structures upon binding a cognate ligand. They are of great interest not only for understanding gene regulation but also as targets for designing small-molecule antibiotics and chemical tools. Understanding the molecular determinants for ligand affinity and selectivity is thus crucial for designing synthetic ligands. Here we carried out extensive molecular dynamics simulations of a PreQ1 riboswitch bound to either cognate or synthetic ligands. By comparing and contrasting these two groups of ligands, we learn how the chemical (e.g., number of hydrogen bond donors and acceptors) and physical (e.g., molecular size) features of ligands affect binding affinity and ligand exit paths. While the number of hydrogen bond donors and acceptors is a key determinant for RNA binding affinity, the ligand size affects the rigidity of the binding pocket and thereby regulates the unbinding of the ligand. These lessons provide guidance for designing riboswitch ligands.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou, China
- Department of Chemistry, University of Illinois at Chicago, Chicago, Illinois, United States of America
| | - Huan-Xiang Zhou
- Department of Chemistry, University of Illinois at Chicago, Chicago, Illinois, United States of America
- Department of Physics, University of Illinois at Chicago, Chicago, Illinois, United States of America
- * E-mail:
| |
Collapse
|
6
|
Chen J, Li N, Wang X, Chen J, Zhang JZH, Zhu T. Molecular mechanism related to the binding of fluorophores to Mango-II revealed by multiple-replica molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:10636-10649. [PMID: 33904542 DOI: 10.1039/d0cp06438f] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Recently, RNA aptamers activating small-molecule fluorophores have been successfully applied to tag and track RNAs in vivo. It is of significance to investigate the molecular mechanism of the fluorophore-RNA aptamer bindings at the atomic level to seek a possible pathway to enhance the fluorescence efficiency of fluorophores. In this work, multiple replica molecular dynamics (MRMD) simulations, essential dynamics (ED) analysis, and hierarchical clustering analysis were coupled to probe the effect of A22U mutation on the binding of two fluorophores, TO1-Biotin (TO1) and TO3-Biotin (TO3), to the Mango-II RNA aptamer (Mango-II). ED analysis reveals that A22U induces alterations in the binding pocket and sites of TO1 and TO3 to the Mango-II, which in turn tunes the fluorophore-RNA interface and changes the interactions of TO1 and TO3 with separate nucleotides of Mango-II. Dynamics analyses also uncover that A22U exerts the opposite impact on the molecular surface areas of the Mango-II and sugar puckers of nucleotides 22 and 23 in Mango-II complexed with TO1 and TO3. Moreover, the calculations of binding free energies suggest that A22U strengthens the binding ability of TO1 to the mutated Mango-II but weakens TO3 to the mutated Mango-II when compared with WT. These findings imply that point mutation in nucleotides possibly tune the fluorescence of fluorophores binding to RNA aptamers, providing a possible scheme to enhance the fluorescence of fluorophores.
Collapse
Affiliation(s)
- Junxiao Chen
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and School of Chemistry and Pharmaceutical Engineering, Qilu University of Technology, Jinan, 250353, People's Republic of China
| | - Na Li
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China.
| | - Xingyu Wang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, People's Republic of China.
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China and Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
7
|
Ma B, Bai G, Nussinov R, Ding J, Wang YX. Conformational Ensemble of TteAdoCbl Riboswitch Provides Stable Structural Elements for Conformation Selection and Population Shift in Cobalamin Recognition. J Phys Chem B 2021; 125:2589-2596. [PMID: 33683130 PMCID: PMC9272747 DOI: 10.1021/acs.jpcb.1c00038] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cobalamin riboswitch is a cis-regulatory element widely found in the 5'-UTRs of the vitamin B12-associated genes in bacteria, resulting in modulation and production of a particular protein. Thermoanaerobacter tengcongensis (Tte) AdoCbl riboswitches are the largest of the known riboswitches with 210 nucleotides, partially due to its long peripheral P6-extension, which enable high affinity of AdoCbl. Two structural elements, T-loop/T-looplike motif and kissing loop are key to the global folding of the RNA. While the structure of the TteAdoCbl riboswitch complex is known, we still do not understand the structure and conformation before AdoCbl ligand recognition. In order to delineate the conformational changes and the stabilities of long-range interactions, we have performed extensive all-atom replica-exchange molecular dynamics simulations of the TteAdoCbl riboswitch with a total simulation time of 2296 ns. We found that both the T-loop/T-looplike motif and kissing loop are very stable with ligand binding. The gating conformation changes of P6-extension allow the ligand to bind to the preorganized kissing loop binding pocket. The T-loop/T-looplike motif has much more hydrogen bonds than observed in TteAdoCbl riboswitch complex crystal structure, indicating an allosteric response of the T-loop/T-looplike motif. Our study demonstrated that the conformational ensemble of TteAdoCbl riboswitch provides stable structural elements for conformation selection and population shift in cobalamin recognition.
Collapse
Affiliation(s)
- Buyong Ma
- Engineering Research Center of Cell and Therapeutic Antibody, MOE, School of Pharmacy, Shanghai Jiao Tong University, Shanghai 200240, China
- Basic Science Program, Leidos Biomedical Research, Inc. Laboratory of Cancer ImmunoMetabolism, National Cancer Institute, Frederick, Maryland 21702, United States
| | - Ganggang Bai
- Engineering Research Center of Cell and Therapeutic Antibody, MOE, School of Pharmacy, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Laboratory of Cancer ImmunoMetabolism, National Cancer Institute, Frederick, Maryland 21702, United States
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Jienyu Ding
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702, United States
| | - Yun-Xing Wang
- Protein-Nucleic Acid Interaction Section, Structural Biophysics Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702, United States
| |
Collapse
|
8
|
Hu G, Li H, Xu S, Wang J. Ligand Binding Mechanism and Its Relationship with Conformational Changes in Adenine Riboswitch. Int J Mol Sci 2020; 21:ijms21061926. [PMID: 32168940 PMCID: PMC7139962 DOI: 10.3390/ijms21061926] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2020] [Revised: 03/01/2020] [Accepted: 03/09/2020] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Guodong Hu
- Correspondence: (G.H.); (J.W.); Tel.: +86-534-8987536 (G.H.); +86-534-8985933 (J.W.)
| | | | | | - Jihua Wang
- Correspondence: (G.H.); (J.W.); Tel.: +86-534-8987536 (G.H.); +86-534-8985933 (J.W.)
| |
Collapse
|
9
|
Bao L, Wang J, Xiao Y. Molecular dynamics simulation of the binding process of ligands to the add adenine riboswitch aptamer. Phys Rev E 2020; 100:022412. [PMID: 31574664 DOI: 10.1103/physreve.100.022412] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2019] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Lei Bao
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Jun Wang
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
10
|
Chen J, Wang X, Pang L, Zhang JZH, Zhu T. Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations. Nucleic Acids Res 2020; 47:6618-6631. [PMID: 31173143 PMCID: PMC6649850 DOI: 10.1093/nar/gkz499] [Citation(s) in RCA: 118] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 05/22/2019] [Accepted: 05/29/2019] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Xingyu Wang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Tong Zhu
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|
11
|
Common Secondary and Tertiary Structural Features of Aptamer-Ligand Interaction Shared by RNA Aptamers with Different Primary Sequences. Molecules 2019; 24:molecules24244535. [PMID: 31835789 PMCID: PMC6943582 DOI: 10.3390/molecules24244535] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Revised: 11/30/2019] [Accepted: 12/05/2019] [Indexed: 11/16/2022] Open
Abstract
Aptamer selection can yield many oligonucleotides with different sequences and affinities for the target molecule. Here, we have combined computational and experimental approaches to understand if aptamers with different sequences but the same molecular target share structural and dynamical features. NEO1A, with a known NMR-solved structure, displays a flexible loop that interacts differently with individual aminoglycosides, its ligand affinities and specificities are responsive to ionic strength, and it possesses an adenosine in the loop that is critical for high-affinity ligand binding. NEO2A was obtained from the same selection and, although they are only 43% identical in overall sequence, NEO1A and NEO2A share similar loop sequences. Experimental analysis by 1D NMR and 2-aminopurine reporters combined with molecular dynamics modeling revealed similar structural and dynamical characteristics in both aptamers. These results are consistent with the hypothesis that the target ligand drives aptamer structure and also selects relevant dynamical characteristics for high-affinity aptamer-ligand interaction. Furthermore, they suggest that it might be possible to “migrate” structural and dynamical features between aptamer group members with different primary sequences but with the same target ligand.
Collapse
|
12
|
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
|
13
|
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
|
14
|
Diggins P, Liu C, Deserno M, Potestio R. Optimal Coarse-Grained Site Selection in Elastic Network Models of Biomolecules. J Chem Theory Comput 2018; 15:648-664. [PMID: 30514085 PMCID: PMC6391041 DOI: 10.1021/acs.jctc.8b00654] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Elastic network models, simple structure-based representations of biomolecules where atoms interact via short-range harmonic potentials, provide great insight into a molecule's internal dynamics and mechanical properties at extremely low computational cost. Their efficiency and effectiveness have made them a pivotal instrument in the computer-aided study of proteins and, since a few years, also of nucleic acids. In general, the coarse-grained sites, i.e. those effective force centers onto which the all-atom structure is mapped, are constructed based on intuitive rules: a typical choice for proteins is to retain only the C α atoms of each amino acid. However, a mapping strategy relying only on the atom type and not the local properties of its embedding can be suboptimal compared to a more careful selection. Here, we present a strategy in which the subset of atoms, each of which is mapped onto a unique coarse-grained site of the model, is selected in a stochastic search aimed at optimizing a cost function. The latter is taken to be a simple measure of the consistency between the harmonic approximation of an elastic network model and the harmonic model obtained through exact integration of the discarded degrees of freedom. The method is applied to two representatives of structurally very different types of biomolecules: the protein adenylate kinase and the RNA molecule adenine riboswitch. Our analysis quantifies the substantial impact that an algorithm-driven selection of coarse-grained sites can have on a model's properties.
Collapse
Affiliation(s)
- Patrick Diggins
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Changjiang Liu
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States.,Department of Biophysics , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Markus Deserno
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Raffaello Potestio
- Physics Department , University of Trento , via Sommarive, 14 I-38123 Trento , Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications , I-38123 Trento , Italy
| |
Collapse
|
15
|
Thermodynamic investigation of kissing-loop interactions. Biochimie 2018; 157:177-183. [PMID: 30502370 DOI: 10.1016/j.biochi.2018.11.012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2018] [Accepted: 11/20/2018] [Indexed: 01/01/2023]
Abstract
Kissing loop interactions (KLIs) are a common motif that is critical in retroviral dimerization, viroid replication, mRNA, and riboswitches. In addition, KLIs are currently used in a variety of biotechnology applications, such as in aptamer sensors, RNA scaffolds and to stabilize vaccines for therapeutics. Here we describe the thermodynamics of a basic intramolecular DNA capable of engaging in a KLI, consisting of two hairpins connected by a flexible linker. Each hairpin loop has a five-nucleotide complementary sequence theoretically capable of engaging in a KLI. On either side of each loop is two thymines which will not engage in kissing but are present to provide more flexibility and optimal KLI positioning. Our results suggest that the KLI occurs even at physiological salt levels, and that the KLI does not alter the thermodynamics and stability of the two stem structures. The KLI does not involve all five nucleotides, or at least each base-pair stack is not making full contact. Adding a second strand complementary to the bottom of the kissing complex removes flexibility and causes destabilization of the stems. The KLI of this less flexible complex is maintained but the TM is reduced, indicating an entopic penalty to its formation.
Collapse
|
16
|
Kesherwani M, N H V K, Velmurugan D. Conformational Dynamics of thiM Riboswitch To Understand the Gene Regulation Mechanism Using Markov State Modeling and the Residual Fluctuation Network Approach. J Chem Inf Model 2018; 58:1638-1651. [PMID: 29939019 DOI: 10.1021/acs.jcim.8b00155] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Manish Kesherwani
- Centre for Advanced Study in Crystallography and Biophysics , University of Madras, Guindy Campus , Chennai - 600025 India
| | - Kutumbarao N H V
- Centre for Advanced Study in Crystallography and Biophysics , University of Madras, Guindy Campus , Chennai - 600025 India
| | - Devadasan Velmurugan
- Centre for Advanced Study in Crystallography and Biophysics , University of Madras, Guindy Campus , Chennai - 600025 India
| |
Collapse
|
17
|
Š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
|
18
|
Gong S, Wang Y, Wang Z, Zhang W. Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches. Molecules 2017; 22:molecules22071169. [PMID: 28703767 PMCID: PMC6152003 DOI: 10.3390/molecules22071169] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2017] [Revised: 07/07/2017] [Accepted: 07/09/2017] [Indexed: 11/16/2022] Open
Abstract
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, SMK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the riboswitches.
Collapse
Affiliation(s)
- Sha Gong
- Hubei Key Laboratory of Economic Forest Germplasm Improvement and Resources Comprehensive Utilization, Hubei Collaborative Innovation Center for the Characteristic Resources Exploitation of Dabie Mountains, Huanggang Normal University, Huanggang 438000, Hubei, China.
| | - Yanli Wang
- Department of Physics, Wuhan University, Wuhan 430072, Hubei, China.
| | - Zhen Wang
- Department of Physics, Wuhan University, Wuhan 430072, Hubei, China.
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan 430072, Hubei, China.
| |
Collapse
|
19
|
Bergonzo C, Cheatham TE. Mg 2+ Binding Promotes SLV as a Scaffold in Varkud Satellite Ribozyme SLI-SLV Kissing Loop Junction. Biophys J 2017; 113:313-320. [PMID: 28669407 PMCID: PMC5529310 DOI: 10.1016/j.bpj.2017.06.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Revised: 06/02/2017] [Accepted: 06/06/2017] [Indexed: 11/07/2022] Open
Abstract
Though the structure of the substrate stem loop I (SLI)-stem loop V (SLV) kissing loop junction of the Varkud Satellite ribozyme has been experimentally characterized, the dynamics of this Mg2+-dependent loop-loop interaction have been elusive. Specifically, each hairpin loop contains a U-turn motif, but only SLV shows a conformational shift triggered by Mg2+ ion association. Here, we use molecular dynamics simulations to analyze the binding and dynamics of this kissing loop junction. We show that SLV acts as a scaffold, providing stability to the junction. Mg2+ ions associate with SLV when it is part of the junction in a manner similar to when it is unbound, but there is no specificity in Mg2+ binding for the SLI loop. This suggests that the entropic penalty of ordering the larger SLI is too high, allowing SLV to act as a scaffold for multiple substrate loop sequences.
Collapse
|
20
|
Hanke CA, Gohlke H. Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch. PLoS One 2017; 12:e0179271. [PMID: 28640851 PMCID: PMC5480868 DOI: 10.1371/journal.pone.0179271] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2016] [Accepted: 05/27/2017] [Indexed: 12/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Christian A. Hanke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
21
|
Hayatshahi HS, Bergonzo C, Cheatham III TE. Investigating the ion dependence of the first unfolding step of GTPase-Associating Center ribosomal RNA. J Biomol Struct Dyn 2017; 36:243-253. [DOI: 10.1080/07391102.2016.1274272] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Hamed S. Hayatshahi
- Department of Medicinal Chemistry, College of Pharmacy, The University of Utah, 2000 East 30 South Skaggs 307, Salt Lake City, Utah 84112-5820, USA
| | - Christina Bergonzo
- Department of Medicinal Chemistry, College of Pharmacy, The University of Utah, 2000 East 30 South Skaggs 307, Salt Lake City, Utah 84112-5820, USA
| | - Thomas E. Cheatham III
- Department of Medicinal Chemistry, College of Pharmacy, The University of Utah, 2000 East 30 South Skaggs 307, Salt Lake City, Utah 84112-5820, USA
| |
Collapse
|
22
|
Hu G, Ma A, Wang J. Ligand Selectivity Mechanism and Conformational Changes in Guanine Riboswitch by Molecular Dynamics Simulations and Free Energy Calculations. J Chem Inf Model 2017; 57:918-928. [PMID: 28345904 DOI: 10.1021/acs.jcim.7b00139] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| | - Aijing Ma
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| |
Collapse
|
23
|
Li C, Lv D, Zhang L, Yang F, Wang C, Su J, Zhang Y. Approach to the unfolding and folding dynamics of add A-riboswitch upon adenine dissociation using a coarse-grained elastic network model. J Chem Phys 2017; 145:014104. [PMID: 27394096 DOI: 10.1063/1.4954992] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Chunhua Li
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Dashuai Lv
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Lei Zhang
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Feng Yang
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Cunxin Wang
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Jiguo Su
- College of Science, Yanshan University, Qinhuangdao 066004, China
| | - Yang Zhang
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 45108, USA
| |
Collapse
|
24
|
Chandra V, Hannan Z, Xu H, Mandal M. Single-molecule analysis reveals multi-state folding of a guanine riboswitch. Nat Chem Biol 2016; 13:194-201. [PMID: 27941758 DOI: 10.1038/nchembio.2252] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2016] [Accepted: 10/03/2016] [Indexed: 12/29/2022]
Abstract
Guanine-responsive riboswitches undergo ligand-dependent structural rearrangements to control gene expression by transcription termination. While the molecular basis for ligand recognition is well established, the associated structural rearrangements and the kinetics involved in the formation of the aptamer domain are less well understood. Using high-resolution optical tweezers, we followed the folding trajectories of a single molecule of the xpt-pbuX guanine aptamer from Bacillus subtilis. We report a rapid six-state conformational rearrangement, in which three of the states are guanine dependent, during the transition from the linear to the native receptor conformation. The folding completes in <1 s. The force-dependent equilibrium kinetics and the mutational data indicated that the flexible J2-J3 junction undergoes a ligand-dependent conformational switching, which triggers the formation of the long-range tertiary interactions and the P1 helix. In the absence of the right ligand, the junction failed to initiate the series of conformational rearrangements required for the riboswitch activities.
Collapse
Affiliation(s)
- Vishnu Chandra
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| | - Zain Hannan
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| | - Huizhong Xu
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| | - Maumita Mandal
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA.,Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA.,Department of Biomedical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| |
Collapse
|
25
|
Šponer J, Krepl M, Banáš P, Kührová P, Zgarbová M, Jurečka P, Havrila M, Otyepka M. How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes? WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27863061 DOI: 10.1002/wrna.1405] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 01/01/2023]
Abstract
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| |
Collapse
|
26
|
Wang J, Xiao Y. Types and concentrations of metal ions affect local structure and dynamics of RNA. Phys Rev E 2016; 94:040401. [PMID: 27841650 DOI: 10.1103/physreve.94.040401] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Indexed: 01/01/2023]
Abstract
The roles that metal ions play in the structure and dynamics of RNA molecules are long-standing problems that have been studied extensively but are still not well understood. Here we show that metal ions have distributions around RNA molecules that strongly depend on the types and concentrations of the metal ions and also the electrostatic surface of the molecule. In particular, the ion distributions may not balance all the local electronegativity of the molecule. These ion distributions do not only greatly affect local structures but also lead to different local dynamics of RNA. We studied the effects of different ion solutions on the structure and dynamics of RNA by taking the preQ_{1} riboswitch aptamer domain as an illustrative example and using molecular dynamics simulations. Since the local structures and dynamics of RNAs are important to their functions, our results also indicate that the selection of proper ion conditions is necessary to model them correctly, in contrast to the use of diverse ion solutions in current molecular dynamics simulations.
Collapse
Affiliation(s)
- Jun Wang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
27
|
Suresh G, Srinivasan H, Nanda S, Priyakumar UD. Ligand-Induced Stabilization of a Duplex-like Architecture Is Crucial for the Switching Mechanism of the SAM-III Riboswitch. Biochemistry 2016; 55:3349-60. [PMID: 27249101 DOI: 10.1021/acs.biochem.5b00973] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Gorle Suresh
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Hyderabad 500 032, India
| | - Harini Srinivasan
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Hyderabad 500 032, India
| | - Shivani Nanda
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Hyderabad 500 032, India
| | - U Deva Priyakumar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Hyderabad 500 032, India
| |
Collapse
|
28
|
Durand G, Dausse E, Goux E, Fiore E, Peyrin E, Ravelet C, Toulmé JJ. A combinatorial approach to the repertoire of RNA kissing motifs; towards multiplex detection by switching hairpin aptamers. Nucleic Acids Res 2016; 44:4450-9. [PMID: 27067541 PMCID: PMC4872101 DOI: 10.1093/nar/gkw206] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2015] [Accepted: 03/15/2016] [Indexed: 01/22/2023] Open
Abstract
Loop–loop (also known as kissing) interactions between RNA hairpins are involved in several mechanisms in both prokaryotes and eukaryotes such as the regulation of the plasmid copy number or the dimerization of retroviral genomes. The stability of kissing complexes relies on loop parameters (base composition, sequence and size) and base combination at the loop–loop helix - stem junctions. In order to identify kissing partners that could be used as regulatory elements or building blocks of RNA scaffolds, we analysed a pool of 5.2 × 106 RNA hairpins with randomized loops. We identified more than 50 pairs of kissing RNA hairpins. Two kissing motifs, 5′CCNY and 5′RYRY, generate highly stable complexes with KDs in the low nanomolar range. Such motifs were introduced in the apical loop of hairpin aptamers that switch between unfolded and folded state upon binding to their cognate target molecule, hence their name aptaswitch. The aptaswitch–ligand complex is specifically recognized by a second RNA hairpin named aptakiss through loop–loop interaction. Taking advantage of our kissing motif repertoire we engineered aptaswitch–aptakiss modules for purine derivatives, namely adenosine, GTP and theophylline and demonstrated that these molecules can be specifically and simultaneously detected by surface plasmon resonance or by fluorescence anisotropy.
Collapse
Affiliation(s)
- Guillaume Durand
- University of Bordeaux, ARNA Laboratory, 146 rue Léo Saignat, 33076 Bordeaux, France Inserm U1212, 146 rue Léo Saignat, 33076 Bordeaux, France CNRS UMR5320, 146 rue Léo Saignat, 33076 Bordeaux, France
| | - Eric Dausse
- University of Bordeaux, ARNA Laboratory, 146 rue Léo Saignat, 33076 Bordeaux, France Inserm U1212, 146 rue Léo Saignat, 33076 Bordeaux, France CNRS UMR5320, 146 rue Léo Saignat, 33076 Bordeaux, France
| | - Emma Goux
- University Grenoble Alpes, Département de Pharmacochimie Moléculaire, CNRS UMR5063, 38400 St Martin d'Hères, France
| | - Emmanuelle Fiore
- University Grenoble Alpes, Département de Pharmacochimie Moléculaire, CNRS UMR5063, 38400 St Martin d'Hères, France
| | - Eric Peyrin
- University Grenoble Alpes, Département de Pharmacochimie Moléculaire, CNRS UMR5063, 38400 St Martin d'Hères, France
| | - Corinne Ravelet
- University Grenoble Alpes, Département de Pharmacochimie Moléculaire, CNRS UMR5063, 38400 St Martin d'Hères, France
| | - Jean-Jacques Toulmé
- University of Bordeaux, ARNA Laboratory, 146 rue Léo Saignat, 33076 Bordeaux, France Inserm U1212, 146 rue Léo Saignat, 33076 Bordeaux, France CNRS UMR5320, 146 rue Léo Saignat, 33076 Bordeaux, France
| |
Collapse
|
29
|
Pinamonti G, Bottaro S, Micheletti C, Bussi G. Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments. Nucleic Acids Res 2015; 43:7260-9. [PMID: 26187990 PMCID: PMC4551938 DOI: 10.1093/nar/gkv708] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2015] [Accepted: 06/30/2015] [Indexed: 11/23/2022] Open
Abstract
Elastic network models (ENMs) are valuable and efficient tools for characterizing the collective internal dynamics of proteins based on the knowledge of their native structures. The increasing evidence that the biological functionality of RNAs is often linked to their innate internal motions poses the question of whether ENM approaches can be successfully extended to this class of biomolecules. This issue is tackled here by considering various families of elastic networks of increasing complexity applied to a representative set of RNAs. The fluctuations predicted by the alternative ENMs are stringently validated by comparison against extensive molecular dynamics simulations and SHAPE experiments. We find that simulations and experimental data are systematically best reproduced by either an all-atom or a three-beads-per-nucleotide representation (sugar-base-phosphate), with the latter arguably providing the best balance of accuracy and computational complexity.
Collapse
Affiliation(s)
- Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Cristian Micheletti
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
30
|
Shaw E, St-Pierre P, McCluskey K, Lafontaine DA, Penedo JC. Using sm-FRET and denaturants to reveal folding landscapes. Methods Enzymol 2015; 549:313-41. [PMID: 25432755 DOI: 10.1016/b978-0-12-801122-5.00014-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
RNA folding studies aim to clarify the relationship among sequence, tridimensional structure, and biological function. In the last decade, the application of single-molecule fluorescence resonance energy transfer (sm-FRET) techniques to investigate RNA structure and folding has revealed the details of conformational changes and timescale of the process leading to the formation of biologically active RNA structures with subnanometer resolution on millisecond timescales. In this review, we initially summarize the first wave of single-molecule FRET-based RNA techniques that focused on analyzing the influence of mono- and divalent metal ions on RNA function, and how these studies have provided very valuable information about folding pathways and the presence of intermediate and low-populated states. Next, we describe a second generation of single-molecule techniques that combine sm-FRET with the use of chemical denaturants as an emerging powerful approach to reveal information about the dynamics and energetics of RNA folding that remains hidden using conventional sm-FRET approaches. The main advantages of using the competing interplay between folding agents such as metal ions and denaturants to observe and manipulate the dynamics of RNA folding and RNA-ligand interactions is discussed in the context of the adenine riboswitch aptamer.
Collapse
Affiliation(s)
- Euan Shaw
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom
| | - Patrick St-Pierre
- RNA Group, Department of Biology, Faculty of Science, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Kaley McCluskey
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom
| | - Daniel A Lafontaine
- RNA Group, Department of Biology, Faculty of Science, Université de Sherbrooke, Sherbrooke, Québec, Canada.
| | - J Carlos Penedo
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom; Biomedical Sciences Research Complex, University of St. Andrews, St. Andrews, Fife, United Kingdom.
| |
Collapse
|
31
|
Condon D, Kennedy SD, Mort BC, Kierzek R, Yildirim I, Turner DH. Stacking in RNA: NMR of Four Tetramers Benchmark Molecular Dynamics. J Chem Theory Comput 2015; 11:2729-2742. [PMID: 26082675 PMCID: PMC4463549 DOI: 10.1021/ct501025q] [Citation(s) in RCA: 84] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2014] [Indexed: 12/31/2022]
Abstract
Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.
Collapse
Affiliation(s)
- David
E. Condon
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
| | - Scott D. Kennedy
- Department
of Biochemistry and Biophysics, University
of Rochester, Rochester, New York 14642, United States
| | - Brendan C. Mort
- University
of Rochester Center for Integrated Research Computing, Rochester, New York 14627, United States
| | - Ryszard Kierzek
- Institute
of Bioorganic Chemistry, Polish Academy
of Sciences, 60-704 Poznan, Poland
| | - Ilyas Yildirim
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
- Department
of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Douglas H. Turner
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
32
|
Di Palma F, Bottaro S, Bussi G. Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations. BMC Bioinformatics 2015; 16 Suppl 9:S6. [PMID: 26051557 PMCID: PMC4464220 DOI: 10.1186/1471-2105-16-s9-s6] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
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.
Collapse
|
33
|
Abstract
Riboswitches, RNA elements found in the untranslated region, regulate gene expression by binding to target metaboloites with exquisite specificity. Binding of metabolites to the conserved aptamer domain allosterically alters the conformation in the downstream expression platform. The fate of gene expression is determined by the changes in the downstream RNA sequence. As the metabolite-dependent cotranscriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate both the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolite. Single molecule force experiments that decipher the free energy landscape of riboswitches from their mechanical responses, theoretical and computational studies have recently shed light on the distinct mechanism of folding dynamics in different classes of riboswitches. Here, we first discuss the dynamics of water around riboswitch, highlighting that water dynamics can enhance the fluctuation of nucleic acid structure. To go beyond native state fluctuations, we used the Self-Organized Polymer model to predict the dynamics of add adenine riboswitch under mechanical forces. In addition to quantitatively predicting the folding landscape of add-riboswitch, our simulations also explain the difference in the dynamics between pbuE adenine- and add adenine-riboswitches. In order to probe the function in vivo, we use the folding landscape to propose a system level kinetic network model to quantitatively predict how gene expression is regulated for riboswitches that are under kinetic control.
Collapse
|
34
|
Gong S, Wang Y, Zhang W. Kinetic regulation mechanism of pbuE riboswitch. J Chem Phys 2015; 142:015103. [DOI: 10.1063/1.4905214] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Sha Gong
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, People’s Republic of China
| | - Yujie Wang
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, People’s Republic of China
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, People’s Republic of China
| |
Collapse
|
35
|
Perez-Gonzalez DC, Penedo JC. Single-Molecule Strategies for DNA and RNA Diagnostics. RNA TECHNOLOGIES 2015. [DOI: 10.1007/978-3-319-17305-4_15] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
36
|
Sund J, Lind C, Åqvist J. Binding site preorganization and ligand discrimination in the purine riboswitch. J Phys Chem B 2014; 119:773-82. [PMID: 25014157 DOI: 10.1021/jp5052358] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
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.
Collapse
Affiliation(s)
- Johan Sund
- Department of Cell and Molecular Biology, Uppsala University , Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden
| | | | | |
Collapse
|
37
|
DFT and MD study of the divalent-cation-mediated interaction of ochratoxin A with DNA nucleosides. J Mol Model 2014; 20:2274. [PMID: 24863531 DOI: 10.1007/s00894-014-2274-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2014] [Accepted: 04/24/2014] [Indexed: 01/01/2023]
Abstract
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.
Collapse
|
38
|
Disney MD, Yildirim I, Childs-Disney JL. Methods to enable the design of bioactive small molecules targeting RNA. Org Biomol Chem 2014; 12:1029-39. [PMID: 24357181 PMCID: PMC4020623 DOI: 10.1039/c3ob42023j] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
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.
Collapse
Affiliation(s)
- Matthew D Disney
- The Department of Chemistry, The Scripps Research Institute, 130 Scripps Way #3A1, Jupiter, FL 33458, USA.
| | | | | |
Collapse
|
39
|
Jain SS, Sonavane UB, Uppuladinne MVN, McLaughlin EC, Wang W, Black S, Joshi RR. Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study. J Biomol Struct Dyn 2014; 33:234-43. [PMID: 24404773 DOI: 10.1080/07391102.2013.870930] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.
Collapse
Affiliation(s)
- Swapan S Jain
- a Department of Chemistry , Bard College , 30 Campus Rd, Annandale-on-Hudson, NY 12504 , USA
| | | | | | | | | | | | | |
Collapse
|
40
|
Di Palma F, Colizzi F, Bussi G. Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch. RNA (NEW YORK, N.Y.) 2013; 19:1517-1524. [PMID: 24051105 PMCID: PMC3851719 DOI: 10.1261/rna.040493.113] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2013] [Accepted: 08/02/2013] [Indexed: 05/28/2023]
Abstract
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.
Collapse
|
41
|
Lin JC, Thirumalai D. Kinetics of allosteric transitions in S-adenosylmethionine riboswitch are accurately predicted from the folding landscape. J Am Chem Soc 2013; 135:16641-50. [PMID: 24087850 DOI: 10.1021/ja408595e] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Riboswitches are RNA elements that allosterically regulate gene expression by binding cellular metabolites. The SAM-III riboswitch, one of several classes that binds S-adenosylmethionine (SAM), represses translation upon binding SAM (OFF state) by encrypting the ribosome binding sequence. We have carried out simulations of the RNA by applying mechanical force (f) to the ends of SAM-III, with and without SAM, to get quantitative insights into the f-dependent structural changes. Force-extension (z) curves (FECs) for the apo (ON) state, obtained in simulations in which f is increased at a constant loading rate, show three intermediates, with the first one being the rupture of SAM binding region, which is greatly stabilized in the OFF state. Force-dependent free energy profiles, G(z,f), as a function of z, obtained in equilibrium constant force simulations, reveal the intermediates observed in FECs. The predicted stability difference between the ON and OFF states using G(z,f) is in excellent agreement with experiments. Remarkably, using G(z,f)s and estimate of an effective diffusion constant at a single value of f allows us to predict the f-dependent transition rates using theory of first passage times for both the apo and holo states. To resolve the kinetics of assembly of SAM-III riboswitch in structural terms, we use force stretch-quench pulse sequences in which the force on RNA is maintained at a low (fq) value starting from a high value for a time period tq. Variation of tq over a wide range results in resolution of elusive states involved in the SAM binding pocket and leads to accurate determination of folding times down to fq = 0. Quantitative measure of the folding kinetics, obtained from the folding landscape, allows us to propose that, in contrast to riboswitches regulating transcription, SAM-III functions under thermodynamic control provided the basal concentration of SAM exceeds a small critical value. All of the predictions are amenable to tests in single molecule pulling experiments.
Collapse
Affiliation(s)
- Jong-Chin Lin
- Department of Chemistry and Biochemistry, Biophysics Program, Institute for Physical Sciences and Technology, University of Maryland , College Park, Maryland 20742, United States
| | | |
Collapse
|