1
|
Li Z, Song G, Zhu J, Mu J, Sun Y, Hong X, Choi T, Cui X, Chen HF. Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field. J Chem Theory Comput 2024. [PMID: 38980289 DOI: 10.1021/acs.jctc.4c00497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.
Collapse
Affiliation(s)
- Zhengxin Li
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ge Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Junjie Zhu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Junxi Mu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yutong Sun
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Xiaokun Hong
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
- College of Biological Science and Engineering, Fuzhou University, Fuzhou, Fujian 350116, China
| | - Taeyoung Choi
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Xiaochen Cui
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic Developmental Sciences, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
2
|
Forget S, Juillé M, Duboué-Dijon E, Stirnemann G. Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme. J Chem Theory Comput 2024. [PMID: 38958594 DOI: 10.1021/acs.jctc.4c00294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.
Collapse
Affiliation(s)
- Sélène Forget
- PASTEUR, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| | - Marie Juillé
- PASTEUR, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
- Université Paris Cité, CNRS, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Elise Duboué-Dijon
- Université Paris Cité, CNRS, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Guillaume Stirnemann
- PASTEUR, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
3
|
Bernard C, Postic G, Ghannay S, Tahi F. State-of-the-RNArt: benchmarking current methods for RNA 3D structure prediction. NAR Genom Bioinform 2024; 6:lqae048. [PMID: 38745991 PMCID: PMC11091930 DOI: 10.1093/nargab/lqae048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Revised: 04/05/2024] [Accepted: 05/08/2024] [Indexed: 05/16/2024] Open
Abstract
RNAs are essential molecules involved in numerous biological functions. Understanding RNA functions requires the knowledge of their 3D structures. Computational methods have been developed for over two decades to predict the 3D conformations from RNA sequences. These computational methods have been widely used and are usually categorised as either ab initio or template-based. The performances remain to be improved. Recently, the rise of deep learning has changed the sight of novel approaches. Deep learning methods are promising, but their adaptation to RNA 3D structure prediction remains difficult. In this paper, we give a brief review of the ab initio, template-based and novel deep learning approaches. We highlight the different available tools and provide a benchmark on nine methods using the RNA-Puzzles dataset. We provide an online dashboard that shows the predictions made by benchmarked methods, freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr/evryrna/state_of_the_rnart/.
Collapse
Affiliation(s)
- Clément Bernard
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Guillaume Postic
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Sahar Ghannay
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Fariza Tahi
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| |
Collapse
|
4
|
Michel HM, Lemkul JA. Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction. Biophys J 2024; 123:1129-1138. [PMID: 38576161 PMCID: PMC11079942 DOI: 10.1016/j.bpj.2024.03.042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Revised: 03/10/2024] [Accepted: 03/29/2024] [Indexed: 04/06/2024] Open
Abstract
G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.
Collapse
Affiliation(s)
- Haley M Michel
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia; Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia.
| |
Collapse
|
5
|
Choi T, Li Z, Song G, Chen HF. Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields. J Chem Theory Comput 2024; 20:2676-2688. [PMID: 38447040 DOI: 10.1021/acs.jctc.4c00066] [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: 03/08/2024]
Abstract
Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.
Collapse
Affiliation(s)
- Taeyoung Choi
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Zhengxin Li
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ge Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
6
|
Pokorná P, Mlýnský V, Bussi G, Šponer J, Stadlbauer P. Molecular dynamics simulations reveal the parallel stranded d(GGGA) 3GGG DNA quadruplex folds via multiple paths from a coil-like ensemble. Int J Biol Macromol 2024; 261:129712. [PMID: 38286387 DOI: 10.1016/j.ijbiomac.2024.129712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Revised: 01/10/2024] [Accepted: 01/16/2024] [Indexed: 01/31/2024]
Abstract
G-quadruplexes (G4s) are non-canonical nucleic acid structures that fold through complex processes. Characterization of the G4 folding landscape may help to elucidate biological roles of G4s but is challenging both experimentally and computationally. Here, we achieved complete folding of a three-quartet parallel DNA G4 with (GGGA)3GGG sequence using all-atom explicit-solvent enhanced-sampling molecular dynamics (MD) simulations. The simulations suggested early formation of guanine stacks in the G-tracts, which behave as semi-rigid blocks in the folding process. The folding continues via the formation of a collapsed compact coil-like ensemble. Structuring of the G4 from the coil then proceeds via various cross-like, hairpin, slip-stranded and two-quartet ensembles and can bypass the G-triplex structure. Folding of the parallel G4 does not appear to involve any salient intermediates and is a multi-pathway process. We also carried out an extended set of simulations of parallel G-hairpins. While parallel G-hairpins are extremely unstable when isolated, they are more stable inside the coil structure. On the methodology side, we show that the AMBER DNA force field predicts the folded G4 to be less stable than the unfolded ensemble, uncovering substantial force-field issues. Overall, we provide unique atomistic insights into the folding landscape of parallel-stranded G4 but also reveal limitations of current state-of-the-art MD techniques.
Collapse
Affiliation(s)
- Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, Trieste 34136, Italy
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic.
| |
Collapse
|
7
|
Beghennou A, Rondot O, Corcé V, Botuha C. 1 H-1,2,3-triazolyl-1,6-naphthyridin-7(6 H)-ones as Potential Fluorescent Nucleoside Analogues: Synthesis and Optical Properties. Molecules 2024; 29:687. [PMID: 38338431 PMCID: PMC10856630 DOI: 10.3390/molecules29030687] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Revised: 01/24/2024] [Accepted: 01/30/2024] [Indexed: 02/12/2024] Open
Abstract
In this article, we present the synthesis and the optical properties of three original molecules as potential fluorescent ribonucleoside analogues incorporating a 1,6-naphthyridin-7(6H)-one scaffold as a fluorescent nucleobase and a 1,2,3-triazole as a linkage. The nucleosides were prepared via a Cu alkyne-azide cycloaddition (CuAAC) reaction between a ribofuranosyl azide and a 4-ethynylpyridine partner. Construction of substituted 1,6-naphthyridin-7(6H)-ones was achieved through two additional steps. Optical property studies were investigated on nucleoside analogues. Powerful fluorescence properties have been evidenced with a remarkable change of emissivity depending on the polarity of the solvent, making these molecules suitable as a new class of artificial fluorescent nucleosides for investigating enzyme binding sites as well as probing nucleic acids. In addition, we are convinced that such analogues could be of great interest in the search for new antiviral or antitumoral drugs based on nucleosides.
Collapse
Affiliation(s)
| | | | - Vincent Corcé
- Institut Parisien de Chimie Moléculaire, CNRS UMR 9232, Sorbonne Université, F-75252 Paris, France; (A.B.); (O.R.)
| | - Candice Botuha
- Institut Parisien de Chimie Moléculaire, CNRS UMR 9232, Sorbonne Université, F-75252 Paris, France; (A.B.); (O.R.)
| |
Collapse
|
8
|
Bernard C, Postic G, Ghannay S, Tahi F. RNAdvisor: a comprehensive benchmarking tool for the measure and prediction of RNA structural model quality. Brief Bioinform 2024; 25:bbae064. [PMID: 38436560 PMCID: PMC10939302 DOI: 10.1093/bib/bbae064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/30/2024] [Accepted: 02/02/2024] [Indexed: 03/05/2024] Open
Abstract
RNA is a complex macromolecule that plays central roles in the cell. While it is well known that its structure is directly related to its functions, understanding and predicting RNA structures is challenging. Assessing the real or predictive quality of a structure is also at stake with the complex 3D possible conformations of RNAs. Metrics have been developed to measure model quality while scoring functions aim at assigning quality to guide the discrimination of structures without a known and solved reference. Throughout the years, many metrics and scoring functions have been developed, and no unique assessment is used nowadays. Each developed assessment method has its specificity and might be complementary to understanding structure quality. Therefore, to evaluate RNA 3D structure predictions, it would be important to calculate different metrics and/or scoring functions. For this purpose, we developed RNAdvisor, a comprehensive automated software that integrates and enhances the accessibility of existing metrics and scoring functions. In this paper, we present our RNAdvisor tool, as well as state-of-the-art existing metrics, scoring functions and a set of benchmarks we conducted for evaluating them. Source code is freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr.
Collapse
Affiliation(s)
- Clement Bernard
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Guillaume Postic
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Sahar Ghannay
- LISN - CNRS/Université Paris-Saclay, France, 91400 Orsay, France
| | - Fariza Tahi
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| |
Collapse
|
9
|
Ashwood B, Jones MS, Lee Y, Sachleben JR, Ferguson AL, Tokmakoff A. Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics. Biophys J 2024; 123:118-133. [PMID: 38006207 PMCID: PMC10808028 DOI: 10.1016/j.bpj.2023.11.022] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Revised: 10/19/2023] [Accepted: 11/22/2023] [Indexed: 11/26/2023] Open
Abstract
Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Michael S Jones
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | - Yumin Lee
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Joseph R Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, Illinois
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois.
| | - Andrei Tokmakoff
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
10
|
Singh K, Reddy G. Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms. J Chem Theory Comput 2024; 20:421-435. [PMID: 38134376 DOI: 10.1021/acs.jctc.3c00999] [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: 12/24/2023]
Abstract
Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.
Collapse
Affiliation(s)
- Kushal Singh
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| |
Collapse
|
11
|
Thiel BC, Poblete S, Hofacker IL. The Multiscale Ernwin/SPQR RNA Structure Prediction Pipeline. Methods Mol Biol 2024; 2726:377-399. [PMID: 38780739 DOI: 10.1007/978-1-0716-3519-3_15] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
Aside from the well-known role in protein synthesis, RNA can perform catalytic, regulatory, and other essential biological functions which are determined by its three-dimensional structure. In this regard, a great effort has been made during the past decade to develop computational tools for the prediction of the structure of RNAs from the knowledge of their sequence, incorporating experimental data to refine or guide the modeling process. Nevertheless, this task can become exceptionally challenging when dealing with long noncoding RNAs, constituted by more than 200 nucleotides, due to their large size and the specific interactions involved. In this chapter, we describe a multiscale approach to predict such structures, incorporating SAXS experimental data into a hierarchical procedure which couples two coarse-grained representations: Ernwin, a helix-based approach, which deals with the global arrangement of secondary structure elements, and SPQR, a nucleotide-centered coarse-grained model, which corrects and refines the structures predicted at the coarser level.We describe the methodology through its application on the Braveheart long noncoding RNA, starting from the SAXS and secondary structure data to propose a refined, all-atom structure.
Collapse
Affiliation(s)
- Bernhard C Thiel
- Department of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Vienna, Austria
| | - Simón Poblete
- Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia, Chile
- Computational Biology Lab, Fundación Ciencia & Vida, Santiago, Chile
- Facultad de Ingeniería, Arquitectura y Diseño, Universidad SanSebastián, Santiago, Chile
| | - Ivo L Hofacker
- Department of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Vienna, Austria.
- Research Group Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Vienna, Austria.
| |
Collapse
|
12
|
Grille L, Gallego D, Darré L, da Rosa G, Battistini F, Orozco M, Dans PD. The pseudotorsional space of RNA. RNA (NEW YORK, N.Y.) 2023; 29:1896-1909. [PMID: 37793790 PMCID: PMC10653382 DOI: 10.1261/rna.079821.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 09/01/2023] [Indexed: 10/06/2023]
Abstract
The characterization of the conformational landscape of the RNA backbone is rather complex due to the ability of RNA to assume a large variety of conformations. These backbone conformations can be depicted by pseudotorsional angles linking RNA backbone atoms, from which Ramachandran-like plots can be built. We explore here different definitions of these pseudotorsional angles, finding that the most accurate ones are the traditional η (eta) and θ (theta) angles, which represent the relative position of RNA backbone atoms P and C4'. We explore the distribution of η - θ in known experimental structures, comparing the pseudotorsional space generated with structures determined exclusively by one experimental technique. We found that the complete picture only appears when combining data from different sources. The maps provide a quite comprehensive representation of the RNA accessible space, which can be used in RNA-structural predictions. Finally, our results highlight that protein interactions lead to significant changes in the population of the η - θ space, pointing toward the role of induced-fit mechanisms in protein-RNA recognition.
Collapse
Affiliation(s)
- Leandro Grille
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
| | - Diego Gallego
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Leonardo Darré
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Gabriela da Rosa
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
| | - Federica Battistini
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Modesto Orozco
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Pablo D Dans
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| |
Collapse
|
13
|
Mlýnský V, Kührová P, Stadlbauer P, Krepl M, Otyepka M, Banáš P, Šponer J. Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations. J Chem Theory Comput 2023; 19:8423-8433. [PMID: 37944118 PMCID: PMC10687871 DOI: 10.1021/acs.jctc.3c00990] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2023] [Revised: 10/18/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Petra Kührová
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Petr Stadlbauer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Michal Otyepka
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| |
Collapse
|
14
|
Kim H, Kim E, Pak Y. Computational Probing of the Folding Mechanism of Human Telomeric G-Quadruplex DNA. J Chem Inf Model 2023; 63:6366-6375. [PMID: 37782649 DOI: 10.1021/acs.jcim.3c01257] [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: 10/04/2023]
Abstract
The human telomeric (htel) sequences in the terminal regions of human telomeres form diverse G-quadruplex (GQ) structures. Despite much experimental efforts to elucidate the folding pathways of htel GQ, no comprehensive model of htel GQ folding has been presented. Here, we describe folding pathways of the htel GQ determined by state-of-the-art enhanced sampling molecular dynamics simulation at the all-atom level. Briefly, GQ folding is initiated by the formation of a single-hairpin and then followed by the formation of double-hairpins, which then branch via distinct folding pathways to produce different GQ topologies (antiparallel chair, antiparallel basket, hybrids 1 and 2, and parallel propeller). In addition to these double-hairpin states, three-triad and two-tetrad structures in antiparallel backbone alignment serve as key intermediates that connect the GQ folding and transition between two different GQs.
Collapse
Affiliation(s)
- Hyeonjun Kim
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 46241, South Korea
| | - Eunae Kim
- College of Pharmacy, Chosun University, Gwangju 61452, South Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 46241, South Korea
| |
Collapse
|
15
|
Malhotra S, Mulvaney T, Cragnolini T, Sidhu H, Joseph A, Beton J, Topf M. RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures. Nucleic Acids Res 2023; 51:9567-9575. [PMID: 37670532 PMCID: PMC10570027 DOI: 10.1093/nar/gkad721] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Revised: 08/10/2023] [Accepted: 08/21/2023] [Indexed: 09/07/2023] Open
Abstract
Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.
Collapse
Affiliation(s)
- Sony Malhotra
- Science and Technology Facilities Council, Scientific Computing, Research Complex at Harwell, Didcot OX11 0FA, UK
| | - Thomas Mulvaney
- Leibniz Institute of Virology, Hamburg 20251, Germany
- Centre for Structural Systems Biology, Hamburg D-22607, Germany
- Universitätsklinikum Hamburg Eppendorf (UKE), Hamburg 20246, Germany
| | - Tristan Cragnolini
- Leibniz Institute of Virology, Hamburg 20251, Germany
- Institute of Structural and Molecular Biology, Department of Biological Sciences, Birkbeck College, University of London, London WC1E 7HX, UK
| | - Haneesh Sidhu
- Institute of Structural and Molecular Biology, Department of Biological Sciences, Birkbeck College, University of London, London WC1E 7HX, UK
| | - Agnel P Joseph
- Science and Technology Facilities Council, Scientific Computing, Research Complex at Harwell, Didcot OX11 0FA, UK
| | - Joseph G Beton
- Leibniz Institute of Virology, Hamburg 20251, Germany
- Centre for Structural Systems Biology, Hamburg D-22607, Germany
| | - Maya Topf
- Leibniz Institute of Virology, Hamburg 20251, Germany
- Centre for Structural Systems Biology, Hamburg D-22607, Germany
- Universitätsklinikum Hamburg Eppendorf (UKE), Hamburg 20246, Germany
| |
Collapse
|
16
|
Viegas RG, Sanches MN, Chen AA, Paulovich FV, Garcia AE, Leite VBP. Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop. J Chem Inf Model 2023; 63:5641-5649. [PMID: 37606640 DOI: 10.1021/acs.jcim.3c00426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.
Collapse
Affiliation(s)
- Rafael G Viegas
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, São Paulo 15.808-305, Brazil
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Murilo N Sanches
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Alan A Chen
- Department of Chemistry and the RNA Institute, University at Albany, Albany, New York 12222, United States
| | - Fernando V Paulovich
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, Eindhoven 5600 MB, the Netherlands
| | - Angel E Garcia
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Physics Division, National Science Foundation, 2415 Eisenhower Ave, Alexandria, Virginia 22314, United States
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
17
|
Ashwood B, Jones MS, Radakovic A, Khanna S, Lee Y, Sachleben JR, Szostak JW, Ferguson AL, Tokmakoff A. Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs. Biophys J 2023; 122:3323-3339. [PMID: 37469144 PMCID: PMC10465710 DOI: 10.1016/j.bpj.2023.07.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Revised: 06/27/2023] [Accepted: 07/17/2023] [Indexed: 07/21/2023] Open
Abstract
Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Michael S Jones
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | | | - Smayan Khanna
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | - Yumin Lee
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Joseph R Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, Illinois
| | - Jack W Szostak
- Department of Chemistry, The University of Chicago, Chicago, Illinois
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois
| | - Andrei Tokmakoff
- Department of Chemistry, The University of Chicago, Chicago, Illinois; The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
18
|
Ashwood B, Jones MS, Lee Y, Sachleben JR, Ferguson AL, Tokmakoff A. Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.22.550182. [PMID: 37546925 PMCID: PMC10401965 DOI: 10.1101/2023.07.22.550182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| | - Michael S. Jones
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Yumin Lee
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| | - Joseph R. Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, IL 60637, United States
| | - Andrew L. Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Andrei Tokmakoff
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| |
Collapse
|
19
|
Oxenfarth A, Kümmerer F, Bottaro S, Schnieders R, Pinter G, Jonker HRA, Fürtig B, Richter C, Blackledge M, Lindorff-Larsen K, Schwalbe H. Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop. J Am Chem Soc 2023. [PMID: 37479220 PMCID: PMC10401711 DOI: 10.1021/jacs.3c03578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.
Collapse
Affiliation(s)
- Andreas Oxenfarth
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Felix Kümmerer
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- IRCCS Humanitas Research Hospital, Department of Biomedical Sciences, Humanitas University, Milan 20089, Italy
| | - Robbin Schnieders
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - György Pinter
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Hendrik R A Jonker
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Boris Fürtig
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Christian Richter
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| | - Martin Blackledge
- Institut de Biologie Structurale (IBS), CEA, CNRS, University Grenoble Alpes, Grenoble 38000, France
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Harald Schwalbe
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance (BMRZ), Goethe University Frankfurt am Main, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Hessen, Germany
| |
Collapse
|
20
|
Maksudov F, Kliuchnikov E, Pierson D, Ujwal M, Marx KA, Chanda A, Barsegov V. Therapeutic phosphorodiamidate morpholino oligonucleotides: Physical properties, solution structures, and folding thermodynamics. MOLECULAR THERAPY. NUCLEIC ACIDS 2023; 31:631-647. [PMID: 36910708 PMCID: PMC9996446 DOI: 10.1016/j.omtn.2023.02.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 02/10/2023] [Indexed: 02/16/2023]
Abstract
Elucidating the structure-function relationships for therapeutic RNA mimicking phosphorodiamidate morpholino oligonucleotides (PMOs) is challenging due to the lack of information about their structures. While PMOs have been approved by the US Food and Drug Administration for treatment of Duchenne muscular dystrophy, no structural information on these unique, charge-neutral, and stable molecules is available. We performed circular dichroism and solution viscosity measurements combined with molecular dynamics simulations and machine learning to resolve solution structures of 22-mer, 25-mer, and 30-mer length PMOs. The PMO conformational dynamics are defined by the competition between non-polar nucleobases and uncharged phosphorodiamidate groups for shielding from solvent exposure. PMO molecules form non-canonical, partially helical, stable folded structures with a small 1.4- to 1.7-nm radius of gyration, low count of three to six base pairs and six to nine base stacks, characterized by -34 to -51 kcal/mol free energy, -57 to -103 kcal/mol enthalpy, and -23 to -53 kcal/mol entropy for folding. The 4.5- to 6.2-cm3/g intrinsic viscosity and Huggins constant of 4.5-9.9 are indicative of extended and aggregating systems. The results obtained highlight the importance of the conformational ensemble view of PMO solution structures, thermodynamic stability of their non-canonical structures, and concentration-dependent viscosity properties. These principles form a paradigm to understand the structure-properties-function relationship for therapeutic PMOs to advance the design of new RNA-mimic-based drugs.
Collapse
Affiliation(s)
- Farkhad Maksudov
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
| | | | - Daniel Pierson
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
| | | | - Kenneth A. Marx
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
- Inciton, Inc., Andover, MA 01854, USA
| | - Arani Chanda
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
- Corresponding author: Arani Chanda, Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA.
| | - Valeri Barsegov
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
- Inciton, Inc., Andover, MA 01854, USA
- Corresponding author: Valeri Barsegov, Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA.
| |
Collapse
|
21
|
Serrano-Chacón I, Mir B, Cupellini L, Colizzi F, Orozco M, Escaja N, González C. pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs. J Am Chem Soc 2023; 145:3696-3705. [PMID: 36745195 PMCID: PMC9936585 DOI: 10.1021/jacs.2c13043] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 02/07/2023]
Abstract
We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.
Collapse
Affiliation(s)
- Israel Serrano-Chacón
- Instituto
de Química Física ”Rocasolano”, CSIC, Serrano 119, 28006Madrid, Spain
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028Barcelona, Spain
| | - Bartomeu Mir
- Instituto
de Química Física ”Rocasolano”, CSIC, Serrano 119, 28006Madrid, Spain
- Inorganic
and Organic Chemistry Department, Organic Chemistry Section, and IBUB, University of Barcelona, Martí i Franquès 1-11, 08028Barcelona, Spain
| | - Lorenzo Cupellini
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028Barcelona, Spain
| | - Francesco Colizzi
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028Barcelona, Spain
| | - Modesto Orozco
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028Barcelona, Spain
- Departament
de Bioquímica i Biomedicina. Facultat de Biologia, Universitat de Barcelona, 08028Barcelona, Spain
| | - Núria Escaja
- Inorganic
and Organic Chemistry Department, Organic Chemistry Section, and IBUB, University of Barcelona, Martí i Franquès 1-11, 08028Barcelona, Spain
- BIOESTRAN
Associated Unit UB-CSIC, 08028Barcelona, Spain
| | - Carlos González
- Instituto
de Química Física ”Rocasolano”, CSIC, Serrano 119, 28006Madrid, Spain
- BIOESTRAN
Associated Unit UB-CSIC, 08028Barcelona, Spain
| |
Collapse
|
22
|
NARall: a novel tool for reconstruction of the all-atom structure of nucleic acids from heavily coarse-grained model. CHEMICAL PAPERS 2022. [DOI: 10.1007/s11696-022-02634-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
AbstractNucleic acids are one of the most important cellular components. These molecules have been studied both experimentally and theoretically. As all-atom simulations are still limited to short time scales, coarse-grain modeling allows to study of those molecules on a longer time scale. Nucleic-Acid united RESidue (NARES-2P) is a low-resolution coarse-grained model with two centers of interaction per repeating unit. It has been successfully applied to study DNA self-assembly and telomeric properties. This force field enables the study of nucleic acids Behavior on a long time scale but lacks atomistic details. In this article, we present new software to reconstruct atomistic details from the NARES-2P model. It has been applied to RNA pseudoknot, nucleic acid four-way junction, G-quadruplex and DNA duplex converted to NARES-2P model and back. Moreover, it was applied to DNA structure folded and self-assembled with NARES-2P.
Collapse
|
23
|
Krepl M, Pokorná P, Mlýnský V, Stadlbauer P, Šponer J. Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field. Nucleic Acids Res 2022; 50:12480-12496. [PMID: 36454011 PMCID: PMC9757038 DOI: 10.1093/nar/gkac1106] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 10/25/2022] [Accepted: 11/04/2022] [Indexed: 12/03/2022] Open
Abstract
Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.
Collapse
Affiliation(s)
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
24
|
Kim H, Pak Y. Improving All-Atom Force Field to Accurately Describe DNA G-Quadruplex Loops. J Phys Chem B 2022; 126:6199-6209. [PMID: 35951994 DOI: 10.1021/acs.jpcb.2c04256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The DNA G-quadruplex (GQ) displays structural polymorphisms, and interactions between its loops and flanking sequences critically determine which of the diverse GQ conformers is adopted. All-atom molecular dynamics (MD) simulations of GQs are computationally challenging due to slow folding times and force field (ff) artifacts. In an earlier study, a direct folding simulation of the simplest DNA GQ (TBA15) was first reported using a modified version of the AMBER bsc1 ff (bsc1_vdW ff). Despite this successful folding simulation, it was later found that the bsc1_vdW ff is somewhat limited in terms of describing loop structures of GQs, which is problematic because GQ loop regions play key roles in ligand binding to modulate GQ activities. In this study, we further modified the bsc1_vdW ff to enhance the GQ loop prediction by fine-tuning a limited number of van der Waals (vdW) parameters of the standard AMBER bsc1 ff to improve the GQ loop distribution of a target GQ system (three-layered antiparallel GQ; mHtel21). Test simulations of this newly generated ff (bsc1_vdWL ff) on DNA GQs with diverse topologies (hybrid1, hybrid2, and parallel propeller) revealed that loop structures were predicted more accurately than by the bsc1_vdW ff. We consider that enhanced sampling MD simulation methods in combination with bsc1_vdWL provide useful simulation protocols for resolving outstanding issues of DNA GQ folding and GQ/ligand binding at the all-atom level.
Collapse
Affiliation(s)
- Hyeonjun Kim
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busandaehak-ro 63beon-gil, Geumjeong-gu, Busan 46241, S. Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busandaehak-ro 63beon-gil, Geumjeong-gu, Busan 46241, S. Korea
| |
Collapse
|
25
|
Length-dependent motions of SARS-CoV-2 frameshifting RNA pseudoknot and alternative conformations suggest avenues for frameshifting suppression. Nat Commun 2022; 13:4284. [PMID: 35879278 PMCID: PMC9310368 DOI: 10.1038/s41467-022-31353-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 06/10/2022] [Indexed: 12/16/2022] Open
Abstract
The SARS-CoV-2 frameshifting element (FSE), a highly conserved mRNA region required for correct translation of viral polyproteins, defines an excellent therapeutic target against Covid-19. As discovered by our prior graph-theory analysis with SHAPE experiments, the FSE adopts a heterogeneous, length-dependent conformational landscape consisting of an assumed 3-stem H-type pseudoknot (graph motif 3_6), and two alternative motifs (3_3 and 3_5). Here, for the first time, we build and simulate, by microsecond molecular dynamics, 30 models for all three motifs plus motif-stabilizing mutants at different lengths. Our 3_6 pseudoknot systems, which agree with experimental structures, reveal interconvertible L and linear conformations likely related to ribosomal pausing and frameshifting. The 3_6 mutant inhibits this transformation and could hamper frameshifting. Our 3_3 systems exhibit length-dependent stem interactions that point to a potential transition pathway connecting the three motifs during ribosomal elongation. Together, our observations provide new insights into frameshifting mechanisms and anti-viral strategies.
Collapse
|
26
|
Bergonzo C, Grishaev A, Bottaro S. Conformational heterogeneity of UCAAUC RNA oligonucleotide from molecular dynamics simulations, SAXS, and NMR experiments. RNA (NEW YORK, N.Y.) 2022; 28:937-946. [PMID: 35483823 PMCID: PMC9202585 DOI: 10.1261/rna.078888.121] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Accepted: 03/17/2022] [Indexed: 06/14/2023]
Abstract
We describe the conformational ensemble of the single-stranded r(UCAAUC) oligonucleotide obtained using extensive molecular dynamics (MD) simulations and Rosetta's FARFAR2 algorithm. The conformations observed in MD consist of A-form-like structures and variations thereof. These structures are not present in the pool generated using FARFAR2. By comparing with available nuclear magnetic resonance (NMR) measurements, we show that the presence of both A-form-like and other extended conformations is necessary to quantitatively explain experimental data. To further validate our results, we measure solution X-ray scattering (SAXS) data on the RNA hexamer and find that simulations result in more compact structures than observed from these experiments. The integration of simulations with NMR via a maximum entropy approach shows that small modifications to the MD ensemble lead to an improved description of the conformational ensemble. Nevertheless, we identify persisting discrepancies in matching experimental SAXS data.
Collapse
Affiliation(s)
- Christina Bergonzo
- National Institute of Standards and Technology and Institute for Bioscience and Biotechnology Research, Rockville, Maryland 20850, USA
| | - Alexander Grishaev
- National Institute of Standards and Technology and Institute for Bioscience and Biotechnology Research, Rockville, Maryland 20850, USA
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- Department of Biomedical Sciences, Humanitas University, 20090 Pieve Emanuele, Italy
| |
Collapse
|
27
|
Fröhlking T, Mlýnský V, Janeček M, Kührová P, Krepl M, Banáš P, Šponer J, Bussi G. Automatic Learning of Hydrogen-Bond Fixes in the AMBER RNA Force Field. J Chem Theory Comput 2022; 18:4490-4502. [PMID: 35699952 PMCID: PMC9281393 DOI: 10.1021/acs.jctc.2c00200] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
![]()
The
capability of
current force fields to reproduce RNA structural
dynamics is limited. Several methods have been developed to take advantage
of experimental data in order to enforce agreement with experiments.
Here, we extend an existing framework which allows arbitrarily chosen
force-field correction terms to be fitted by quantification of the
discrepancy between observables back-calculated from simulation and
corresponding experiments. We apply a robust regularization protocol
to avoid overfitting and additionally introduce and compare a number
of different regularization strategies, namely, L1, L2, Kish size,
relative Kish size, and relative entropy penalties. The training set
includes a GACC tetramer as well as more challenging systems, namely,
gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen
bonds in the AMBER RNA force field are corrected with automatically
determined parameters that we call gHBfixopt. A validation
involving a separate simulation of a system present in the training
set (gcUUCGgc) and new systems not seen during training (CAAU and
UUUU tetramers) displays improvements regarding the native population
of the tetraloop as well as good agreement with NMR experiments for
tetramers when using the new parameters. Then, we simulate folded
RNAs (a kink–turn and L1 stalk rRNA) including hydrogen bond
types not sufficiently present in the training set. This allows a
final modification of the parameter set which is named gHBfix21 and
is suggested to be applicable to a wider range of RNA systems.
Collapse
Affiliation(s)
- Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, Brno 612 65, Czech Republic
| | - Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, Olomouc 771 46, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - 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, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, Brno 612 65, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
| |
Collapse
|
28
|
Mlýnský V, Janeček M, Kührová P, Fröhlking T, Otyepka M, Bussi G, Banáš P, Šponer J. Toward Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Modifications. J Chem Theory Comput 2022; 18:2642-2656. [PMID: 35363478 DOI: 10.1021/acs.jctc.1c01222] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔGfold° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic.,IT4Innovations, VSB─Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava-Poruba, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| |
Collapse
|
29
|
The CAR-mRNA Interaction Surface Is a Zipper Extension of the Ribosome A Site. Int J Mol Sci 2022; 23:ijms23031417. [PMID: 35163343 PMCID: PMC8835751 DOI: 10.3390/ijms23031417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 01/21/2022] [Accepted: 01/22/2022] [Indexed: 12/10/2022] Open
Abstract
The ribosome CAR interaction surface behaves as an extension of the decoding center A site and has H-bond interactions with the +1 codon, which is next in line to enter the A site. Through molecular dynamic simulations, we investigated the codon sequence specificity of this CAR–mRNA interaction and discovered a strong preference for GCN codons, suggesting that there may be a sequence-dependent layer of translational regulation dependent on the CAR interaction surface. Dissection of the CAR–mRNA interaction through nucleotide substitution experiments showed that the first nucleotide of the +1 codon dominates over the second nucleotide position, consistent with an energetically favorable zipper-like activity that emanates from the A site through the CAR–mRNA interface. Moreover, the CAR/+1 codon interaction is affected by the identity of nucleotide 3 of +1 GCN codons, which influences the stacking of G and C. Clustering analysis suggests that the A-site decoding center adopts different neighborhood substates that depend on the identity of the +1 codon.
Collapse
|
30
|
Guo ZH, Yuan L, Tan YL, Zhang BG, Shi YZ. RNAStat: An Integrated Tool for Statistical Analysis of RNA 3D Structures. FRONTIERS IN BIOINFORMATICS 2022; 1:809082. [PMID: 36303785 PMCID: PMC9580920 DOI: 10.3389/fbinf.2021.809082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 12/17/2021] [Indexed: 11/13/2022] Open
Abstract
The 3D architectures of RNAs are essential for understanding their cellular functions. While an accurate scoring function based on the statistics of known RNA structures is a key component for successful RNA structure prediction or evaluation, there are few tools or web servers that can be directly used to make comprehensive statistical analysis for RNA 3D structures. In this work, we developed RNAStat, an integrated tool for making statistics on RNA 3D structures. For given RNA structures, RNAStat automatically calculates RNA structural properties such as size and shape, and shows their distributions. Based on the RNA structure annotation from DSSR, RNAStat provides statistical information of RNA secondary structure motifs including canonical/non-canonical base pairs, stems, and various loops. In particular, the geometry of base-pairing/stacking can be calculated in RNAStat by constructing a local coordinate system for each base. In addition, RNAStat also supplies the distribution of distance between any atoms to the users to help build distance-based RNA statistical potentials. To test the usability of the tool, we established a non-redundant RNA 3D structure dataset, and based on the dataset, we made a comprehensive statistical analysis on RNA structures, which could have the guiding significance for RNA structure modeling. The python code of RNAStat, the dataset used in this work, and corresponding statistical data files are freely available at GitHub (https://github.com/RNA-folding-lab/RNAStat).
Collapse
Affiliation(s)
- Zhi-Hao Guo
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Li Yuan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- *Correspondence: Ya-Zhou Shi,
| |
Collapse
|
31
|
Poblete S, Božič A, Kanduč M, Podgornik R, Guzman HV. RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates. ACS OMEGA 2021; 6:32823-32831. [PMID: 34901632 PMCID: PMC8655909 DOI: 10.1021/acsomega.1c04774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 11/04/2021] [Indexed: 06/14/2023]
Abstract
RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto
de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Computational
Biology Lab, Fundación Ciencia &
Vida, Santiago 7780272, Chile
| | - Anže Božič
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Matej Kanduč
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Rudolf Podgornik
- School
of Physical Sciences and Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Institute
of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Wenzhou
Institute of the University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Horacio V. Guzman
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
32
|
Paloncýová M, Čechová P, Šrejber M, Kührová P, Otyepka M. Role of Ionizable Lipids in SARS-CoV-2 Vaccines As Revealed by Molecular Dynamics Simulations: From Membrane Structure to Interaction with mRNA Fragments. J Phys Chem Lett 2021; 12:11199-11205. [PMID: 34761943 PMCID: PMC8609914 DOI: 10.1021/acs.jpclett.1c03109] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Accepted: 11/01/2021] [Indexed: 06/13/2023]
Abstract
Recent advances in RNA-based medicine have provided new opportunities for the global current challenge, i.e., the COVID-19 pandemic. Novel vaccines are based on a messenger RNA (mRNA) motif with a lipid nanoparticle (LNP) vector, consisting of high content of unique pH-sensitive ionizable lipids (ILs). Here we provide molecular insights into the role of the ILs and lipid mixtures used in current mRNA vaccines. We observed that the lipid mixtures adopted a nonlamellar organization, with ILs separating into a very disordered, pH-sensitive phase. We describe structural differences of the two ILs leading to their different congregation, with implications for the vaccine stability. Finally, as RNA interacts preferentially with IL-rich phases located at the regions with high curvature of lipid phase, local changes in RNA flexibility and base pairing are induced by lipids. A proper atomistic understanding of RNA-lipid interactions may enable rational tailoring of LNP composition for efficient RNA delivery.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Petra Čechová
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Martin Šrejber
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Petra Kührová
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
- IT4Innovations, VSB—Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava-Poruba, Czech Republic
| |
Collapse
|
33
|
Tulsiyan K, Jena S, González-Viegas M, Kar RK, Biswal HS. Structural Dynamics of RNA in the Presence of Choline Amino Acid Based Ionic Liquid: A Spectroscopic and Computational Outlook. ACS CENTRAL SCIENCE 2021; 7:1688-1697. [PMID: 34729412 PMCID: PMC8554839 DOI: 10.1021/acscentsci.1c00768] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Indexed: 05/03/2023]
Abstract
Ribonucleic acid (RNA) is exceedingly sensitive to degradation compared to DNA. The current protocol for storage of purified RNA requires freezing conditions below -20 °C. Recent advancements in biological chemistry have identified amino acid-based ionic liquids as suitable preservation media for RNA, even in the presence of degrading enzymes. However, the mechanistic insight into the interaction between ILs and RNA is unclear. To the best of our knowledge, no attempts are made so far to provide a molecular view. This work aims to establish a detailed understanding of how ILs enable structural stability to RNA sourced from Torula yeast. Herein, we manifest the hypothesis of multimodal binding of IL and its minimal perturbation to the macromolecular structure, with several spectroscopic techniques such as time-resolved fluorescence and fluorescence correlation spectroscopy (FCS) aided with molecular dynamics at microsecond time scales. Relevant structural and thermodynamic details from biophysical experiments confirm that even long-term RNA preservation with ILs is a possible alternative devoid of any structural deformation. These results establish a unifying mechanism of how ILs are maintaining conformational integrity and thermal stability. The atomistic insights are transferable for their potential applications in drug delivery and biomaterials by considering the advantages of having maximum structural retention and minimum toxicity.
Collapse
Affiliation(s)
- Kiran
Devi Tulsiyan
- School
of Chemical Sciences, National Institute
of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, District, Khurda, 752050, Bhubaneswar, India
- Homi
Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Subhrakant Jena
- School
of Chemical Sciences, National Institute
of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, District, Khurda, 752050, Bhubaneswar, India
- Homi
Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - María González-Viegas
- Institut
für Biologie, Humboldt Universität zu Berlin, Invalidenstraße 42, 10115 Berlin, Germany
| | - Rajiv K. Kar
- Faculty
II-Mathematics and Natural Sciences, Technische
Universität Berlin, Sekr. PC 14, Strasse des 17, Juni 135, D-10623 Berlin, Germany
| | - Himansu S. Biswal
- School
of Chemical Sciences, National Institute
of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, District, Khurda, 752050, Bhubaneswar, India
- Homi
Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
- . Phone: +91-674-2494 185/186
| |
Collapse
|
34
|
Bernetti M, Hall KB, Bussi G. Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles. Nucleic Acids Res 2021; 49:e84. [PMID: 34107023 PMCID: PMC8373061 DOI: 10.1093/nar/gkab459] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 05/07/2021] [Accepted: 05/16/2021] [Indexed: 01/03/2023] Open
Abstract
Small-angle X-ray scattering (SAXS) experiments are increasingly used to probe RNA structure. A number of forward models that relate measured SAXS intensities and structural features, and that are suitable to model either explicit-solvent effects or solute dynamics, have been proposed in the past years. Here, we introduce an approach that integrates atomistic molecular dynamics simulations and SAXS experiments to reconstruct RNA structural ensembles while simultaneously accounting for both RNA conformational dynamics and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward models and enhanced sampling methods to sample an heterogenous ensemble of structures, with no information towards the experiments provided on-the-fly. The generated structural ensemble is then reweighted through the maximum entropy principle so as to match reference SAXS experimental data at multiple ionic conditions. Importantly, accurate explicit-solvent forward models are used at this reweighting stage. We apply this framework to the GTPase-associated center, a relevant RNA molecule involved in protein translation, in order to elucidate its ion-dependent conformational ensembles. We show that (a) both solvent and dynamics are crucial to reproduce experimental SAXS data and (b) the resulting dynamical ensembles contain an ion-dependent fraction of extended structures.
Collapse
Affiliation(s)
- Mattia Bernetti
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| | - Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| |
Collapse
|
35
|
Poblete S, Guzman HV. Structural 3D Domain Reconstruction of the RNA Genome from Viruses with Secondary Structure Models. Viruses 2021; 13:1555. [PMID: 34452420 PMCID: PMC8402887 DOI: 10.3390/v13081555] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 07/28/2021] [Accepted: 07/30/2021] [Indexed: 02/07/2023] Open
Abstract
Three-dimensional RNA domain reconstruction is important for the assembly, disassembly and delivery functionalities of a packed proteinaceus capsid. However, to date, the self-association of RNA molecules is still an open problem. Recent chemical probing reports provide, with high reliability, the secondary structure of diverse RNA ensembles, such as those of viral genomes. Here, we present a method for reconstructing the complete 3D structure of RNA genomes, which combines a coarse-grained model with a subdomain composition scheme to obtain the entire genome inside proteinaceus capsids based on secondary structures from experimental techniques. Despite the amount of sampling involved in the folded and also unfolded RNA molecules, advanced microscope techniques can provide points of anchoring, which enhance our model to include interactions between capsid pentamers and RNA subdomains. To test our method, we tackle the satellite tobacco mosaic virus (STMV) genome, which has been widely studied by both experimental and computational communities. We provide not only a methodology to structurally analyze the tertiary conformations of the RNA genome inside capsids, but a flexible platform that allows the easy implementation of features/descriptors coming from both theoretical and experimental approaches.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Chile and Computational Biology Lab, Fundación Ciencia & Vida, Santiago 7780272, Chile
| | - Horacio V. Guzman
- Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
36
|
Bottaro S, Bussi G, Lindorff-Larsen K. Conformational Ensembles of Noncoding Elements in the SARS-CoV-2 Genome from Molecular Dynamics Simulations. J Am Chem Soc 2021; 143:8333-8343. [PMID: 34039006 PMCID: PMC8188756 DOI: 10.1021/jacs.1c01094] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2021] [Indexed: 12/17/2022]
Abstract
The 5' untranslated region (UTR) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome is a conserved, functional and structured genomic region consisting of several RNA stem-loop elements. While the secondary structure of such elements has been determined experimentally, their three-dimensional structures are not known yet. Here, we predict structure and dynamics of five RNA stem loops in the 5'-UTR of SARS-CoV-2 by extensive atomistic molecular dynamics simulations, more than 0.5 ms of aggregate simulation time, in combination with enhanced sampling techniques. We compare simulations with available experimental data, describe the resulting conformational ensembles, and identify the presence of specific structural rearrangements in apical and internal loops that may be functionally relevant. Our atomic-detailed structural predictions reveal a rich dynamics in these RNA molecules, could help the experimental characterization of these systems, and provide putative three-dimensional models for structure-based drug design studies.
Collapse
Affiliation(s)
- Sandro Bottaro
- Structural
Biology and NMR Laboratory & Linderstrøm-Lang Centre for
Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Giovanni Bussi
- Scuola
Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136, Trieste, Italy
| | - Kresten Lindorff-Larsen
- Structural
Biology and NMR Laboratory & Linderstrøm-Lang Centre for
Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
37
|
Szabla R, Zdrowowicz M, Spisz P, Green NJ, Stadlbauer P, Kruse H, Šponer J, Rak J. 2,6-diaminopurine promotes repair of DNA lesions under prebiotic conditions. Nat Commun 2021; 12:3018. [PMID: 34021158 PMCID: PMC8139960 DOI: 10.1038/s41467-021-23300-y] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 04/20/2021] [Indexed: 01/04/2023] Open
Abstract
High-yielding and selective prebiotic syntheses of RNA and DNA nucleotides involve UV irradiation to promote the key reaction steps and eradicate biologically irrelevant isomers. While these syntheses were likely enabled by UV-rich prebiotic environment, UV-induced formation of photodamages in polymeric nucleic acids, such as cyclobutane pyrimidine dimers (CPDs), remains the key unresolved issue for the origins of RNA and DNA on Earth. Here, we demonstrate that substitution of adenine with 2,6-diaminopurine enables repair of CPDs with yields reaching 92%. This substantial self-repairing activity originates from excellent electron donating properties of 2,6-diaminopurine in nucleic acid strands. We also show that the deoxyribonucleosides of 2,6-diaminopurine and adenine can be formed under the same prebiotic conditions. Considering that 2,6-diaminopurine was previously shown to increase the rate of nonenzymatic RNA replication, this nucleobase could have played critical roles in the formation of functional and photostable RNA/DNA oligomers in UV-rich prebiotic environments. UV-induced photodamage that likely occurred during the prebiotic synthesis of DNA and RNA is still an untackled issue for their origin on early Earth. Here, the authors show that substitution of 2,6-diaminopurine for adenine enables repair of cyclobutane pyrimidine dimers with high yields, and demonstrate that both 2,6-diaminopurine and adenine nucleosides can be formed under the same prebiotic conditions.
Collapse
Affiliation(s)
- Rafał Szabla
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, UK. .,Institute of Physics, Polish Academy of Sciences, Warsaw, Poland.
| | | | - Paulina Spisz
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | | | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Janusz Rak
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| |
Collapse
|
38
|
Mráziková K, Mlýnský V, Kührová P, Pokorná P, Kruse H, Krepl M, Otyepka M, Banáš P, Šponer J. UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations. J Chem Theory Comput 2020; 16:7601-7617. [PMID: 33215915 DOI: 10.1021/acs.jctc.0c00801] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
39
|
Niu X, Liu Q, Xu Z, Chen Z, Xu L, Xu L, Li J, Fang X. Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs). Nat Commun 2020; 11:5496. [PMID: 33127896 PMCID: PMC7603331 DOI: 10.1038/s41467-020-19260-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Accepted: 10/07/2020] [Indexed: 02/06/2023] Open
Abstract
Mechanical anisotropy is an essential property for many biomolecules to assume their structures, functions and applications, however, the mechanisms for their direction-dependent mechanical responses remain elusive. Herein, by using a single-molecule nanopore sensing technique, we explore the mechanisms of directional mechanical stability of the xrRNA1 RNA from ZIKA virus (ZIKV), which forms a complex ring-like architecture. We reveal extreme mechanical anisotropy in ZIKV xrRNA1 which highly depends on Mg2+ and the key tertiary interactions. The absence of Mg2+ and disruption of the key tertiary interactions strongly affect the structural integrity and attenuate mechanical anisotropy. The significance of ring structures in RNA mechanical anisotropy is further supported by steered molecular dynamics simulations in combination with force distribution analysis. We anticipate the ring structures can be used as key elements to build RNA-based nanostructures with controllable mechanical anisotropy for biomaterial and biomedical applications.
Collapse
Affiliation(s)
- Xiaolin Niu
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Qiuhan Liu
- Key Laboratory of Bioorganic Phosphorus Chemistry and Chemical Biology, Department of Chemistry, Tsinghua University, Beijing, 100084, China
| | - Zhonghe Xu
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Zhifeng Chen
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Linghui Xu
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Lilei Xu
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Jinghong Li
- Key Laboratory of Bioorganic Phosphorus Chemistry and Chemical Biology, Department of Chemistry, Tsinghua University, Beijing, 100084, China.
| | - Xianyang Fang
- Beijing Advanced Innovation Center for Structfural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China.
| |
Collapse
|
40
|
Halder A, Kumar S, Valsson O, Reddy G. Mg 2+ Sensing by an RNA Fragment: Role of Mg 2+-Coordinated Water Molecules. J Chem Theory Comput 2020; 16:6702-6715. [PMID: 32941038 DOI: 10.1021/acs.jctc.0c00589] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg2+, Ca2+, and K+) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg2+ solution. We show that the rigid network of Mg2+-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K+ and Ca2+ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.
Collapse
Affiliation(s)
- Antarip Halder
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Omar Valsson
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| |
Collapse
|
41
|
Zhang H, Zhang H, Chen C. Simulation Study of the Plasticity of k-Turn Motif in Different Environments. Biophys J 2020; 119:1416-1426. [PMID: 32918889 DOI: 10.1016/j.bpj.2020.08.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Revised: 07/15/2020] [Accepted: 08/12/2020] [Indexed: 10/23/2022] Open
Abstract
The k-turn is a widespread and important motif in RNA. According to the internal hydrogen bond network, it has two stable states, called N1 and N3. The relative stability between the states changes with the environment. It is able to accept different conformations in different environments. This is called the "plasticity" of a molecule. In this work, we study the plasticity of k-turn by the mixing REMD method in explicit solvent. The results are concluded as follows. First, N1 and N3 are almost equally stable when k-turn is in the solvent alone. The molecule is quite flexible as a hinge. However, after binding to different proteins, such as the proteins L7Ae and L24e, k-turn falls into one global minimum. The preferred state could be either N1 or N3. On the contrary, the other nonpreferred state becomes unstable with a weaker binding affinity to the protein. It reveals that RNA-binding protein is able to modulate the representative state of k-turn at equilibrium. This is in agreement with the findings in experiments. Moreover, free energy calculations show that the free energy barrier between the N1 and N3 states of k-turn increases in the complexes. The state-to-state transition is greatly impeded. We also give a deep discussion on the mechanism of the high plasticity of k-turn in different environments.
Collapse
Affiliation(s)
- Haomiao Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Haozhe Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China.
| |
Collapse
|
42
|
Cabaj MK, Dominiak PM. Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals. Nucleic Acids Res 2020; 48:8302-8319. [PMID: 32725210 PMCID: PMC7470937 DOI: 10.1093/nar/gkaa629] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 07/10/2020] [Accepted: 07/16/2020] [Indexed: 11/16/2022] Open
Abstract
We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.
Collapse
Affiliation(s)
- Małgorzata Katarzyna Cabaj
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
43
|
Duan Q, Tao P, Wang J, Xiao Y. Molecular dynamics study of ways of RNA base-pair formation. Phys Rev E 2020; 102:032403. [PMID: 33076020 DOI: 10.1103/physreve.102.032403] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Accepted: 08/21/2020] [Indexed: 06/11/2023]
Abstract
Base pairing is a fundamental phenomenon in RNA structure and function. However, although there have been considerable recent advances, some important aspects of base-pair formation are still unknown, including the ways of base-pair formation and detailed roles of metal ions. Here we show that base pairs can form through four different ways: stabilizing, bridging, rotating, and shifting. Among them the stabilizing and bridging ways involve direct binding of metal ions while the rotating and shifting ways do not in most cases. Furthermore, we find that the formations of base pairs in different positions of the hairpin stem may adopt different ways.
Collapse
Affiliation(s)
- Qiangqiang Duan
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Peng Tao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Jun Wang
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
44
|
Bian Y, Song F, Zhang J, Yu J, Wang J, Wang W. Insights into the Kinetic Partitioning Folding Dynamics of the Human Telomeric G-Quadruplex from Molecular Simulations and Machine Learning. J Chem Theory Comput 2020; 16:5936-5947. [PMID: 32794754 DOI: 10.1021/acs.jctc.0c00340] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The human telomeric DNA G-quadruplex follows a kinetic partitioning folding mechanism. The underlying folding landscape potentially has many minima separated by high free-energy barriers. However, using current theoretical models to characterize this complex folding landscape has remained a challenging problem. In this study, by developing a hybrid atomistic structure-based model that merges structural information on the hybrid-1, hybrid-2, and chair-type G-quadruplex topologies, we investigated a kinetic partitioning folding process of human telomeric DNA involving three native folds. The model was validated as it reproduced the experimental observation that the hybrid-1 conformation is the major fold and the hybrid-2 conformation is kinetically more accessible. A three-step mechanism was revealed for the formation of the hybrid-1 conformation, while a two-step mechanism was demonstrated for the formation of hybrid-2 and chair-type conformations. Likewise, a class of state in which structures adopted inappropriate combinations of syn/anti guanine nucleotides was found to greatly slow down the folding process. In addition, by employing the XGBoost machine learning algorithm, three interatom distances and six dihedral angles were identified as essential internal coordinates to represent the low-dimensional folding landscape. The strategy of coupling the multibasin model and the machine learning algorithm may be useful to investigate the conformational dynamics of other multistate biomolecules.
Collapse
Affiliation(s)
- Yunqiang Bian
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.,National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Feng Song
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| | - Jian Zhang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Jiafeng Yu
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| | - Wei Wang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| |
Collapse
|
45
|
Gallego D, Darré L, Dans PD, Orozco M. VeriNA3d: an R package for nucleic acids data mining. Bioinformatics 2020; 35:5334-5336. [PMID: 31286135 DOI: 10.1093/bioinformatics/btz553] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 06/20/2019] [Accepted: 07/06/2019] [Indexed: 11/12/2022] Open
Abstract
SUMMARY veriNA3d is an R package for the analysis of nucleic acids structural data, with an emphasis in complex RNA structures. In addition to single-structure analyses, veriNA3d also implements functions to handle whole datasets of mmCIF/PDB structures that could be retrieved from public/local repositories. Our package aims to fill a gap in the data mining of nucleic acids structures to produce flexible and high throughput analysis of structural databases. AVAILABILITY AND IMPLEMENTATION http://mmb.irbbarcelona.org/gitlab/dgallego/veriNA3d. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Diego Gallego
- Computational Biology Node, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology.,Department of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, Barcelona, Spain
| | - Leonardo Darré
- Computational Biology Node, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology.,Functional Genomics Laboratory and Biomolecular Simulations Laboratory, Institute Pasteur of Montevideo, Montevideo, Uruguay
| | - Pablo D Dans
- Computational Biology Node, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology
| | - Modesto Orozco
- Computational Biology Node, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology.,Department of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, Barcelona, Spain
| |
Collapse
|
46
|
Bottaro S, Nichols PJ, Vögeli B, Parrinello M, Lindorff-Larsen K. Integrating NMR and simulations reveals motions in the UUCG tetraloop. Nucleic Acids Res 2020; 48:5839-5848. [PMID: 32427326 PMCID: PMC7293013 DOI: 10.1093/nar/gkaa399] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2020] [Revised: 04/03/2020] [Accepted: 05/17/2020] [Indexed: 12/21/2022] Open
Abstract
We provide an atomic-level description of the structure and dynamics of the UUCG RNA stem-loop by combining molecular dynamics simulations with experimental data. The integration of simulations with exact nuclear Overhauser enhancements data allowed us to characterize two distinct states of this molecule. The most stable conformation corresponds to the consensus three-dimensional structure. The second state is characterized by the absence of the peculiar non-Watson-Crick interactions in the loop region. By using machine learning techniques we identify a set of experimental measurements that are most sensitive to the presence of non-native states. We find that although our MD ensemble, as well as the consensus UUCG tetraloop structures, are in good agreement with experiments, there are remaining discrepancies. Together, our results show that (i) the MD simulation overstabilize a non-native loop conformation, (ii) eNOE data support its presence with a population of ≈10% and (iii) the structural interpretation of experimental data for dynamic RNAs is highly complex, even for a simple model system such as the UUCG tetraloop.
Collapse
Affiliation(s)
- Sandro Bottaro
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Parker J Nichols
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Michele Parrinello
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
47
|
Mlýnský V, Kührová P, Kühr T, Otyepka M, Bussi G, Banáš P, Šponer J. Fine-Tuning of the AMBER RNA Force Field with a New Term Adjusting Interactions of Terminal Nucleotides. J Chem Theory Comput 2020; 16:3936-3946. [PMID: 32384244 DOI: 10.1021/acs.jctc.0c00228] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Determination of RNA structural-dynamic properties is challenging for experimental methods. Thus, atomistic molecular dynamics (MD) simulations represent a helpful technique complementary to experiments. However, contemporary MD methods still suffer from limitations of force fields (ffs), including imbalances in the nonbonded ff terms. We have recently demonstrated that some improvement of state-of-the-art AMBER RNA ff can be achieved by adding a new term for H-bonding called gHBfix, which increases tuning flexibility and reduces risk of side-effects. Still, the first gHBfix version did not fully correct simulations of short RNA tetranucleotides (TNs). TNs are key benchmark systems due to availability of unique NMR data, although giving too much weight on improving TN simulations can easily lead to overfitting to A-form RNA. Here we combine the gHBfix version with another term called tHBfix, which separately treats H-bond interactions formed by terminal nucleotides. This allows to refine simulations of RNA TNs without affecting simulations of other RNAs. The approach is in line with adopted strategy of current RNA ffs, where the terminal nucleotides possess different parameters for terminal atoms than the internal nucleotides. Combination of gHBfix with tHBfix significantly improves the behavior of RNA TNs during well-converged enhanced-sampling simulations using replica exchange with solute tempering. TNs mostly populate canonical A-form like states while spurious intercalated structures are largely suppressed. Still, simulations of r(AAAA) and r(UUUU) TNs show some residual discrepancies with primary NMR data which suggests that future tuning of some other ff terms might be useful. Nevertheless, the tHBfix has a clear potential to improve modeling of key biochemical processes, where interactions of RNA single stranded ends are involved.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Tomáš Kühr
- Department of Computer Science, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
48
|
Reißer S, Zucchelli S, Gustincich S, Bussi G. Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data. Nucleic Acids Res 2020; 48:1164-1174. [PMID: 31889193 PMCID: PMC7026608 DOI: 10.1093/nar/gkz1184] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Revised: 12/05/2019] [Accepted: 12/09/2019] [Indexed: 01/12/2023] Open
Abstract
Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.
Collapse
Affiliation(s)
- Sabine Reißer
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Silvia Zucchelli
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- Department of Health Sciences, Center for Autoimmune and Allergic Diseases (CAAD) and Interdisciplinary Research Center of Autoimmune Diseases (IRCAD), University of Piemonte Orientale, Novara, Italy
| | - Stefano Gustincich
- Central RNA Laboratory and Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia (IIT), 16163 Genova, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
49
|
Liao Q. Enhanced sampling and free energy calculations for protein simulations. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:177-213. [PMID: 32145945 DOI: 10.1016/bs.pmbts.2020.01.006] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Molecular dynamics simulation is a powerful computational technique to study biomolecular systems, which complements experiments by providing insights into the structural dynamics relevant to biological functions at atomic scale. It can also be used to calculate the free energy landscapes of the conformational transitions to better understand the functions of the biomolecules. However, the sampling of biomolecular configurations is limited by the free energy barriers that need to be overcome, leading to considerable gaps between the timescales reached by MD simulation and those governing biological processes. To address this issue, many enhanced sampling methodologies have been developed to increase the sampling efficiency of molecular dynamics simulations and free energy calculations. Usually, enhanced sampling algorithms can be classified into methods based on collective variables (CV-based) and approaches which do not require predefined CVs (CV-free). In this chapter, the theoretical basis of free energy estimation is briefly reviewed first, followed by the reviews of the most common CV-based and CV-free methods including the presentation of some examples and recent developments. Finally, the combination of different enhanced sampling methods is discussed.
Collapse
Affiliation(s)
- Qinghua Liao
- Science for Life Laboratory, Department of Chemistry-BMC, Uppsala University, Uppsala, Sweden.
| |
Collapse
|
50
|
Zhang H, Gong Q, Zhang H, Chen C. FSATOOL: A useful tool to do the conformational sampling and trajectory analysis work for biomolecules. J Comput Chem 2020; 41:156-164. [PMID: 31603251 DOI: 10.1002/jcc.26083] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 09/10/2019] [Accepted: 09/12/2019] [Indexed: 12/27/2022]
Abstract
Reliable conformational sampling and trajectory analysis are always important to the study of the folding or binding mechanisms of biomolecules. Generally, one has to prepare many complicated parameters and follow a lot of steps to obtain the final data. The whole process is too complicated to new users. In this article, we provide a convenient and user-friendly tool that is compatible to AMBER, called fast sampling and analysis tool (FSATOOL). FSATOOL has some useful features. First and the most important, the whole work is extremely simplified into two steps, one is the fast sampling procedure and the other is the trajectory analysis procedure. Second, it contains several powerful sampling methods for the simulation on graphics process unit, including our previous mixing replica exchange molecular dynamics method. The method combines the advantages of the biased and unbiased simulations. Finally, it extracts the dominant transition pathways automatically from the folding network by Markov state model. Users do not need to do the tedious intermediate steps by hand. To illustrate the usage of FSATOOL in practice, we perform one simulation for a RNA hairpin in explicit solvent. All the results are presented. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Haomiao Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Qiankun Gong
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Haozhe Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| |
Collapse
|