1
|
Muscat S, Martino G, Manigrasso J, Marcia M, De Vivo M. On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules. J Chem Theory Comput 2024. [PMID: 39150960 DOI: 10.1021/acs.jctc.4c00773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2024]
Abstract
RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.
Collapse
Affiliation(s)
- Stefano Muscat
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Gianfranco Martino
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Jacopo Manigrasso
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, 431 50 Mölndal, Sweden
| | - Marco Marcia
- European Molecular Biology Laboratory Grenoble, 71 Avenue des Martyrs, 38042 Grenoble, France
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, 751 23 Uppsala, Sweden
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
2
|
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; 20:6082-6097. [PMID: 38980289 DOI: 10.1021/acs.jctc.4c00497] [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: 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
|
3
|
Ratajczyk EJ, Šulc P, Turberfield AJ, Doye JPK, Louis AA. Coarse-grained modeling of DNA-RNA hybrids. J Chem Phys 2024; 160:115101. [PMID: 38497475 DOI: 10.1063/5.0199558] [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: 11/17/2023] [Accepted: 01/26/2024] [Indexed: 03/19/2024] Open
Abstract
We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.
Collapse
Affiliation(s)
- Eryk J Ratajczyk
- Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
- Kavli Institute for Nanoscience Discovery, University of Oxford, Dorothy Crowfoot Hodgkin Building, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, Arizona 85281, USA
- School of Natural Sciences, Department of Bioscience, Technical University Munich, 85748 Garching, Germany
| | - Andrew J Turberfield
- Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
- Kavli Institute for Nanoscience Discovery, University of Oxford, Dorothy Crowfoot Hodgkin Building, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| |
Collapse
|
4
|
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
|
5
|
Hamerla C, Mondal P, Hegger R, Burghardt I. Controlled destabilization of caged circularized DNA oligonucleotides predicted by replica exchange molecular dynamics simulations. Phys Chem Chem Phys 2023; 25:26132-26144. [PMID: 37740309 DOI: 10.1039/d3cp02961a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/24/2023]
Abstract
Spatiotemporal control is a critical issue in the design of strategies for the photoregulation of oligonucleotide activity. Efficient uncaging, i.e., activation by removal of photolabile protecting groups (PPGs), often necessitates multiple PPGs. An alternative approach is based on circularization strategies, exemplified by intrasequential circularization, also denoted photo-tethering, as introduced in [Seyfried et al., Angew. Chem., Int. Ed., 2017, 56, 359]. Here, we develop a computational protocol, relying on replica exchange molecular dynamics (REMD), in order to characterize the destabilization of a series of circularized, caged DNA oligonucleotides addressed in the aforementioned study. For these medium-sized (32 nt) oligonucleotides, melting temperatures are computed, whose trend is in good agreement with experiment, exhibiting a large destabilization and, hence, reduction of the melting temperature of the order of ΔTm ∼ 30 K as compared with the native species. The analysis of free energy landscapes confirms the destabilization pattern experienced by the circularized oligonucleotides. The present study underscores that computational protocols that capture controlled destabilization and uncaging of oligonucleotides are promising as predictive tools in the tailored photocontrol of nucleic acids.
Collapse
Affiliation(s)
- Carsten Hamerla
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt am Main, Germany.
| | - Padmabati Mondal
- Department of Chemistry and Center for Atomic, Molecular, and Optical Sciences and Technologies (CAMOST), Indian Institute of Science Education and Research (IISER) Tirupati, Panguru (G.P), Yerpedu Mandal, 517619 - Tirupati Dist., Andhra Pradesh, India
| | - Rainer Hegger
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt am Main, Germany.
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt am Main, Germany.
| |
Collapse
|
6
|
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
|
7
|
Schrenková V, Para Kkadan MS, Kessler J, Kapitán J, Bouř P. Molecular dynamics and Raman optical activity spectra reveal nucleotide conformation ratios in solution. Phys Chem Chem Phys 2023; 25:8198-8208. [PMID: 36880812 DOI: 10.1039/d2cp05756e] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
Nucleotide conformational flexibility affects their biological functions. Although the spectroscopy of Raman optical activity (ROA) is well suited to structural analyses in aqueous solutions, the link between the spectral shape and the nucleotide geometry is not fully understood. We recorded the Raman and ROA spectra of model nucleotides (rAMP, rGMP, rCMP, and dTMP) and interpreted them on the basis of molecular dynamics (MD) combined with density functional theory (DFT). The relation between the sugar puckering, base conformation and spectral intensities is discussed. Hydrogen bonds between the sugar's C3' hydroxyl and the phosphate groups were found to be important for the sugar puckering. The simulated spectra correlated well with the experimental data and provided an understanding of the dependence of the spectral shapes on conformational dynamics. Most of the strongest spectral bands could be assigned to vibrational molecular motions. Decomposition of the experimental spectra into calculated subspectra based on arbitrary maps of free energies provided experimental conformer populations, which could be used to verify and improve the MD predictions. The analyses indicate some flaws of common MD force fields, such as being unable to describe the fine conformer distribution. Also the accuracy of conformer populations obtained from the spectroscopic data depends on the simulations, improvement of which is desirable for gaining a more detailed insight in the future. Improvement of the spectroscopic and computational methodology for nucleotides also provides opportunities for its application to larger nucleic acids.
Collapse
Affiliation(s)
- Věra Schrenková
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic. .,Department of Analytical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague, Czech Republic
| | - Mohammed Siddhique Para Kkadan
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic. .,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16, Prague, Czech Republic
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| | - Josef Kapitán
- Department of Optics, Palacký University Olomouc, 17. listopadu 12, 77146, Olomouc, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic. .,Department of Analytical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague, Czech Republic
| |
Collapse
|
8
|
Janicki TD, Van Vleet MJ, Schmidt JR. Development and Implementation of Atomically Anisotropic First-Principles Force Fields: A Benzene Case Study. J Phys Chem A 2023; 127:1736-1749. [PMID: 36780209 DOI: 10.1021/acs.jpca.2c07244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
π-interactions are an important motif in chemical and biochemical systems. However, due to their anisotropic electron densities and complex balance of intermolecular interactions, aromatic molecules represent an ongoing challenge for accurate and transferable force field development. Historically, ab initio force fields for aromatics have not exhibited good accuracy with respect to bulk properties or have only been used to study gas-phase dimers. Using benzene as a proof of concept, herein we show how our own ab initio MASTIFF force field incorporates an atomically anisotropic description of intermolecular interactions to yield an accurate and robust model for aromatic interactions irrespective of phase. Compared to existing models, the MASTIFF benzene force field not only is accurate for liquid phase properties but also offers transferability to the gas and solid phases. Additionally, we introduce a computationally efficient OpenMM plugin which enables customizable anisotropic intermolecular functional forms and which can be generically used in any MD simulation where a model for nonspherical atomic features is required. Overall, our results demonstrate the importance of atomic-level anisotropy in enabling next-generation ab initio force field development.
Collapse
Affiliation(s)
- Tesia D Janicki
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Mary J Van Vleet
- Department of Chemistry and Biochemistry, Spelman College, 350 Spelman Ln SW, Atlanta, Georgia 30314, United States
| | - J R Schmidt
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
9
|
Oprzeska-Zingrebe EA, Smiatek J. Basket-type G-quadruplex with two tetrads in the presence of TMAO and urea: A molecular dynamics study. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2022.134375] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
10
|
Kensinger AH, Makowski JA, Pellegrene KA, Imperatore JA, Cunningham CL, Frye CJ, Lackey PE, Mihailescu MR, Evanseck JD. Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations. ACS PHYSICAL CHEMISTRY AU 2023; 3:30-43. [PMID: 36711027 PMCID: PMC9578647 DOI: 10.1021/acsphyschemau.2c00032] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 09/21/2022] [Accepted: 09/21/2022] [Indexed: 11/05/2022]
Abstract
The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.
Collapse
Affiliation(s)
- Adam H. Kensinger
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Joseph A. Makowski
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Kendy A. Pellegrene
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Joshua A. Imperatore
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Caylee L. Cunningham
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Caleb J. Frye
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Patrick E. Lackey
- Department
of Biochemistry and Chemistry, Westminster
College, New Wilmington, Pennsylvania16172, United States
| | - Mihaela Rita Mihailescu
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| | - Jeffrey D. Evanseck
- Department
of Chemistry and Biochemistry and Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania15282, United States
| |
Collapse
|
11
|
Olson WK, Li Y, Fenley MO. Insights into DNA solvation found in protein-DNA structures. Biophys J 2022; 121:4749-4758. [PMID: 36380591 PMCID: PMC9808563 DOI: 10.1016/j.bpj.2022.11.019] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Revised: 10/31/2022] [Accepted: 11/10/2022] [Indexed: 11/16/2022] Open
Abstract
The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.
Collapse
Affiliation(s)
- Wilma K Olson
- Department of Chemistry and Chemical Biology and Center for Quantitative Biology, Rutgers, the State University of New Jersey, Piscataway, New Jersey.
| | - Yun Li
- Department of Chemistry and Chemical Biology and Center for Quantitative Biology, Rutgers, the State University of New Jersey, Piscataway, New Jersey
| | - Marcia O Fenley
- Department of Chemistry and Chemical Biology and Center for Quantitative Biology, Rutgers, the State University of New Jersey, Piscataway, New Jersey; Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida
| |
Collapse
|
12
|
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
|
13
|
Base-specific RNA force field improving the dynamics conformation of nucleotide. Int J Biol Macromol 2022; 222:680-690. [PMID: 36167105 DOI: 10.1016/j.ijbiomac.2022.09.183] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 09/02/2022] [Accepted: 09/19/2022] [Indexed: 11/23/2022]
Abstract
RNA plays a key role in numerous biological processes. Traditional experimental methods have difficulties capturing the structure and dynamic conformation of RNA. Thus, Molecular dynamic simulations (MDs) has become an essential complementary for RNA experiment. However, state-of-the-art RNA force fields have two major limitations of overestimation base stacking propensity and generation of a high ratio of intercalated conformations. Therefore, a two-step strategy was used to optimize the parameters of ff99bsc0χOL3 (named BSFF1) to improve these limitations, which as well adjusted the unbonded parameters of nucleobase heavy atoms and added ζ/α grid-based energy correction map energy term with reweighting. MD simulations of tetranucleotides indicate that BSFF1 can significantly decrease the ratio of intercalated conformations. Tests of single-strand RNA and kink-turn show that BSFF1 force field can reproduce more accurate conformers than ff99bsc0χOL3 force field. BSFF1 can also stabilize the conformers of duplex and riboswitch. The successful ab initio folding of tetraloop further supports the performance of BSFF1. These findings confirm that the newly developed force field BSFF1 can improve the conformer sampling of RNA.
Collapse
|
14
|
Taghavi A, Riveros I, Wales DJ, Yildirim I. Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations. J Chem Theory Comput 2022; 18:3637-3653. [PMID: 35652685 DOI: 10.1021/acs.jctc.2c00178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.
Collapse
Affiliation(s)
- Amirhossein Taghavi
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States.,Department of Chemistry, Scripps Research Institute Florida, Jupiter, Florida 33458, United States
| | - Ivan Riveros
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| | - David J Wales
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Ilyas Yildirim
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| |
Collapse
|
15
|
A novel fluorescent aptasensor based on mesoporous silica nanoparticles for the selective detection of sulfadiazine in edible tissue. ARAB J CHEM 2022. [DOI: 10.1016/j.arabjc.2022.104067] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
|
16
|
Nicy, Chakraborty D, Wales DJ. Energy Landscapes for Base-Flipping in a Model DNA Duplex. J Phys Chem B 2022; 126:3012-3028. [PMID: 35427136 PMCID: PMC9098180 DOI: 10.1021/acs.jpcb.2c00340] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2022] [Revised: 03/24/2022] [Indexed: 12/31/2022]
Abstract
We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.
Collapse
Affiliation(s)
- Nicy
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K.
| | - Debayan Chakraborty
- Department
of Chemistry, The University of Texas at
Austin, Austin, Texas 78712, United States
| | - David J. Wales
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K.
| |
Collapse
|
17
|
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
|
18
|
Kumar S, Reddy G. TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg 2+ Concentration. J Phys Chem B 2022; 126:2369-2381. [PMID: 35298161 DOI: 10.1021/acs.jpcb.1c10794] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg2+. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg2+ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg2+ concentration ([Mg2+]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg2+], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg2+ aids in the formation of tertiary contacts. For [Mg2+] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.
Collapse
Affiliation(s)
- Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| |
Collapse
|
19
|
Esmaeeli R, Andal B, Perez A. Searching for Low Probability Opening Events in a DNA Sliding Clamp. Life (Basel) 2022; 12:life12020261. [PMID: 35207548 PMCID: PMC8876151 DOI: 10.3390/life12020261] [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: 12/30/2021] [Revised: 02/01/2022] [Accepted: 02/03/2022] [Indexed: 11/27/2022] Open
Abstract
The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.
Collapse
|
20
|
González-Olvera JC, Zamorano-Carrillo A, Arreola-Jardón G, Pless RC. Residue interactions affecting the deprotonation of internal guanine moieties in oligodeoxyribonucleotides, calculated by FMO methods. J Mol Model 2022; 28:43. [DOI: 10.1007/s00894-022-05033-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 01/18/2022] [Indexed: 11/24/2022]
|
21
|
Chen J, Liu H, Cui X, Li Z, Chen HF. RNA-Specific Force Field Optimization with CMAP and Reweighting. J Chem Inf Model 2022; 62:372-385. [PMID: 35021622 DOI: 10.1021/acs.jcim.1c01148] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
RNA plays a key role in a variety of cell activities. However, it is difficult to capture its structure dynamics by the traditional experimental methods because of the inherent limitations. Molecular dynamics simulation has become a valuable complement to the experimental methods. Previous studies have indicated that the current force fields cannot accurately reproduce the conformations and structural dynamics of RNA. Therefore, an RNA-specific force field was developed to improve the conformation sampling of RNA. The distribution of ζ/α dihedrals of tetranucleotides was optimized by a reweighting method, and the grid-based energy correction map (CMAP) term was first introduced into the Amber RNA force field of ff99bsc0χOL3, named ff99OL3_CMAP1. Extensive validations of tetranucleotides and tetraloops show that ff99OL3_CMAP1 can significantly decrease the population of an incorrect structure, increase the consistency between the simulation results and experimental values for tetranucleotides, and improve the stability of tetraloops. ff99OL3_CMAP1 can also precisely reproduce the conformation of a duplex and riboswitches. These findings confirm that the newly developed force field ff99OL3_CMAP1 can improve the conformer sampling of RNA.
Collapse
Affiliation(s)
- Jun Chen
- State Key Laboratory of Microbial Metabolism, 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, 20024 Shanghai, China
| | - Hao Liu
- Institute of Natural Sciences, Shanghai Jiao Tong University, 200240 Shanghai, China
| | - Xiaochen Cui
- State Key Laboratory of Microbial Metabolism, 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, 20024 Shanghai, China
| | - Zhengxin Li
- State Key Laboratory of Microbial Metabolism, 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, 20024 Shanghai, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, 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, 20024 Shanghai, China.,Shanghai Center for Bioinformation Technology, 200240 Shanghai, China
| |
Collapse
|
22
|
Muraru S, Muraru S, Nitu FR, Ionita M. Recent Efforts and Milestones for Simulating Nucleic Acid FRET Experiments through Computational Methods. J Chem Inf Model 2022; 62:232-239. [PMID: 35014791 DOI: 10.1021/acs.jcim.1c00957] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Computational methods can greatly aid nucleic acid fluorescence experiments by either offering fully detailed atomic insights into the conformations and interactions present in the studied system or by providing accurate simulations of the fundamental parameters. Fluorescence-based optical biosensors show great potential for clinical diagnosis of life-altering diseases with a very high specificity. Many of the designs for such rely on the concept of Förster resonance energy transfer (FRET). Currently, the methods used experimentally make use of theoretical assumptions which fundamentally affect the results. Having a detailed atomistic overview or significant simulated parameters could improve the understanding of the calculations and provide much more accurate outcomes. However, there are many challenges that need to be addressed before standardized computational protocols can be employed. This review is meant to highlight the progress made for computational methods used to simulate FRET experiments for nucleic acid probes. Recent advances have been made in computational tools, such as force field parametrizations and improved protocols. Complementary simulations to experimental data are also comprised in the this review.
Collapse
Affiliation(s)
- Sorin Muraru
- Faculty of Medical Engineering, University Politehnica of Bucharest, Gh. Polizu Street 1-7, 011061 Bucharest, Romania
| | - Sebastian Muraru
- Faculty of Medical Engineering, University Politehnica of Bucharest, Gh. Polizu Street 1-7, 011061 Bucharest, Romania
| | - Florentin Romeo Nitu
- Faculty of Medical Engineering, University Politehnica of Bucharest, Gh. Polizu Street 1-7, 011061 Bucharest, Romania
| | - Mariana Ionita
- Faculty of Medical Engineering, University Politehnica of Bucharest, Gh. Polizu Street 1-7, 011061 Bucharest, Romania.,Advanced Polymer Materials Group, University Polithenica of Bucharest, Gh. Polizu Street 1-7, 011061 Bucharest, Romania
| |
Collapse
|
23
|
Yang W, Miyazaki T, Chen P, Hong T, Naito M, Miyahara Y, Matsumoto A, Kataoka K, Miyata K, Cabral H. Block catiomer with flexible cationic segment enhances complexation with siRNA and the delivery performance in vitro. SCIENCE AND TECHNOLOGY OF ADVANCED MATERIALS 2021; 22:850-863. [PMID: 34658669 PMCID: PMC8519541 DOI: 10.1080/14686996.2021.1976055] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 08/29/2021] [Accepted: 08/30/2021] [Indexed: 05/14/2023]
Abstract
RNA interference (RNAi) by small interfering RNAs (siRNAs) is a promising therapeutic approach. Because siRNA has limited intracellular access and is rapidly cleared in vivo, the success of RNAi depends on efficient delivery technologies. Particularly, polyion complexation between block catiomers and siRNA is a versatile approach for constructing effective carriers, such as unit polyion complexes (uPIC), core-shell polyion complex (PIC) micelles and vesicular siRNAsomes, by engineering the structure of block catiomers. In this regard, the flexibility of block catiomers could be an important parameter in the formation of PIC nanostructures with siRNA, though its effect remains unknown. Here, we studied the influence of block catiomer flexibility on the assembly of PIC structures with siRNA using a complementary polymeric system, i.e. poly(ethylene glycol)-poly(L-lysine) (PEG-PLL) and PEG-poly(glycidylbutylamine) (PEG-PGBA), which has a relatively more flexible polycation segment than PEG-PLL. Mixing PEG-PGBA with siRNA at molar ratios of primary amines in polymer to phosphates in the siRNA (N/P ratios) higher than 1.5 promoted the multimolecular association of uPICs, whereas PEG-PLL formed uPIC at all N/P ratios higher than 1. Moreover, uPICs from PEG-PGBA were more stable against counter polyanion exchange than uPICs from PEG-PLL, probably due to a favorable complexation process, as suggested by computational studies of siRNA/block catiomer binding. In in vitro experiments, PEG-PGBA uPICs promoted effective intracellular delivery of siRNA and efficient gene knockdown. Our results indicate the significance of polycation flexibility on assembling PIC structures with siRNA, and its potential for developing innovative delivery systems.
Collapse
Affiliation(s)
- Wenqian Yang
- Department of Bioengineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
| | - Takuya Miyazaki
- Department of Bioengineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
- Kanagawa Institute of Industrial Science and Technology, Kanagawa, Japan
- Institute of Biomaterials and Bioengineering, Tokyo Medical and Dental University, Tokyo, Japan
| | - Pengwen Chen
- Department of Bioengineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
| | - Taehun Hong
- Department of Bioengineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
| | - Mitsuru Naito
- Center for Disease Biology and Integrative Medicine, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
| | - Yuji Miyahara
- Kanagawa Institute of Industrial Science and Technology, Kanagawa, Japan
| | - Akira Matsumoto
- Kanagawa Institute of Industrial Science and Technology, Kanagawa, Japan
- Institute of Biomaterials and Bioengineering, Tokyo Medical and Dental University, Tokyo, Japan
| | - Kazunori Kataoka
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
| | - Kanjiro Miyata
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
- Department of Materials Engineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
| | - Horacio Cabral
- Department of Bioengineering, Graduate School of Engineering, The University of Tokyo, Tokyo, Japan
- Innovation Center of NanoMedicine, Kawasaki Institute of Industrial Promotion, Kawasaki, Japan
| |
Collapse
|
24
|
Role and Perspective of Molecular Simulation-Based Investigation of RNA-Ligand Interaction: From Small Molecules and Peptides to Photoswitchable RNA Binding. Molecules 2021; 26:molecules26113384. [PMID: 34205049 PMCID: PMC8199858 DOI: 10.3390/molecules26113384] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 05/28/2021] [Accepted: 06/01/2021] [Indexed: 12/15/2022] Open
Abstract
Aberrant RNA–protein complexes are formed in a variety of diseases. Identifying the ligands that interfere with their formation is a valuable therapeutic strategy. Molecular simulation, validated against experimental data, has recently emerged as a powerful tool to predict both the pose and energetics of such ligands. Thus, the use of molecular simulation may provide insight into aberrant molecular interactions in diseases and, from a drug design perspective, may allow for the employment of less wet lab resources than traditional in vitro compound screening approaches. With regard to basic research questions, molecular simulation can support the understanding of the exact molecular interaction and binding mode. Here, we focus on examples targeting RNA–protein complexes in neurodegenerative diseases and viral infections. These examples illustrate that the strategy is rather general and could be applied to different pharmacologically relevant approaches. We close this study by outlining one of these approaches, namely the light-controllable association of small molecules with RNA, as an emerging approach in RNA-targeting therapy.
Collapse
|
25
|
Hurst T, Zhang D, Zhou Y, Chen SJ. A Bayes-inspired theory for optimally building an efficient coarse-grained folding force field. COMMUNICATIONS IN INFORMATION AND SYSTEMS 2021; 21:65-83. [PMID: 34354546 PMCID: PMC8336718 DOI: 10.4310/cis.2021.v21.n1.a4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Because of their potential utility in predicting conformational changes and assessing folding dynamics, coarse-grained (CG) RNA folding models are appealing for rapid characterization of RNA molecules. Previously, we reported the iterative simulated RNA reference state (IsRNA) method for parameterizing a CG force field for RNA folding, which consecutively updates the simulation force field to reflect marginal distributions of folding coordinates in the structure database and extract various energy terms. While the IsRNA model was validated by showing close agreement between the IsRNA-simulated and experimentally observed distributions, here, we expand our theoretical understanding of the model and, in doing so, improve the parameterization process to optimize the subset of included folding coordinates, which leads to accelerated simulations. Using statistical mechanical theory, we analyze the underlying, Bayesian concept that drives parameterization of the energy function, providing a general method for developing predictive, knowledge-based, polymer force fields on the basis of limited data. Furthermore, we propose an optimal parameterization procedure, based on the principal of maximum entropy.
Collapse
Affiliation(s)
- Travis Hurst
- Department of Physics, University of Missouri-Columbia, Columbia, MO 65211, USA
| | - Dong Zhang
- Department of Physics, University of Missouri-Columbia
| | - Yuanzhe Zhou
- Department of Physics, University of Missouri-Columbia, Columbia, MO 65211, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, MU Institute for Data Science and Informatics, University of Missouri-Columbia, Columbia, MO 65211, USA
| |
Collapse
|
26
|
Nicholson DA, Sengupta A, Nesbitt DJ. Chirality-Dependent Amino Acid Modulation of RNA Folding. J Phys Chem B 2020; 124:11561-11572. [PMID: 33296203 DOI: 10.1021/acs.jpcb.0c07420] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The preponderance of a specific d- or l-chirality in fats, sugars, amino acids, nucleic acids, and so on is ubiquitous in nature, yet the biological origin of such chiral dominance (i.e., with one enantiomer overwhelmingly present) remains an open question. One plausible proposal for the predominance of l-chirality in amino acids could be through evolutionary templating of chiral RNA-folding via chaperone activity. To help evaluate this possibility, single molecule fluorescence experiments have been performed that measure the chiral dependence of chaperone folding dynamics for the simple tetraloop-tetraloop receptor (TL-TLR) tertiary binding motif in the presence of a series of chiral amino acids. Specifically, d- vs l-arginine is found to accelerate the unfolding of this RNA motif in a chirally selective fashion, with temperature-dependent studies of the kinetics performed to extract free energy, enthalpy, and entropy landscapes for the underlying thermodynamics. Furthermore, all-atom molecular dynamics (MD) simulations are pursued to provide additional physical insight into this chiral sensitivity, which reveal enantiomer-specific sampling of nucleic acid surfaces by d- vs l-arginine and support a putative mechanism for chirally specific denaturation of RNA tertiary structure by arginine but not other amino acids.
Collapse
Affiliation(s)
- David A Nicholson
- JILA, National Institute of Standards and Technology and University of Colorado Boulder, Boulder, Colorado 80309 United States.,Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Abhigyan Sengupta
- Department of Physics, Technical University of Munich, Garching, Munich, Germany 85748
| | - David J Nesbitt
- JILA, National Institute of Standards and Technology and University of Colorado Boulder, Boulder, Colorado 80309 United States.,Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, United States.,Department of Physics, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
27
|
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
|
28
|
Kushwaha GS, Patra A, Bhavesh NS. Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins. Front Microbiol 2020; 11:575041. [PMID: 33224117 PMCID: PMC7674647 DOI: 10.3389/fmicb.2020.575041] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Accepted: 10/06/2020] [Indexed: 11/25/2022] Open
Abstract
(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.
Collapse
Affiliation(s)
- Gajraj Singh Kushwaha
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India.,KIIT Technology Business Incubator (KIIT-TBI), Kalinga Institute of Industrial Technology (KIIT) (Deemed to be University), Bhubaneswar, India
| | - Anupam Patra
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India
| | - Neel Sarovar Bhavesh
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India
| |
Collapse
|
29
|
The assessment of three dimensional modelling design for single strand DNA aptamers for computational chemistry application. Biophys Chem 2020; 267:106492. [PMID: 33035750 DOI: 10.1016/j.bpc.2020.106492] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Revised: 09/17/2020] [Accepted: 09/29/2020] [Indexed: 11/23/2022]
Abstract
Aptamers are oligonucleotides and peptides around 15-100 bases in length and are suitable as detection probes or as therapeutics molecules. There are growing interests in the aptamer screening approach through computational simulation methods. DNA and RNA modelling lacks of validation on their predicted 3D structures due to less number of validation tools, unlike protein structures. We suggest an approach to design the stem-loop/hairpin for the three dimensional structure of DNA aptamers through serial applications of computational prediction methods by comparing the simulated structures with the experimental data deposited in PDB Data bank, followed by MD simulations. The result shows minimal structural differences were observed between the designed and the original NMR aptamers, and the stem-loop conformational structures were also retained during the MD thus suggesting the proposed aptamers designing methods are able to synthesize a high quality molecular structure of hairpin aptamers, comparable to the NMR structures.
Collapse
|
30
|
Palomino‐Hernandez O, Margreiter MA, Rossetti G. Challenges in RNA Regulation in Huntington's Disease: Insights from Computational Studies. Isr J Chem 2020. [DOI: 10.1002/ijch.202000021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Oscar Palomino‐Hernandez
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Faculty 1RWTH Aachen 52425 Aachen Germany
- Computation-based Science and Technology Research CenterThe Cyprus Institute Nicosia 2121 Cyprus
- Institute of Life ScienceThe Hebrew University of Jerusalem Jerusalem 91904 Israel
| | - Michael A. Margreiter
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Faculty 1RWTH Aachen 52425 Aachen Germany
| | - Giulia Rossetti
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Jülich Supercomputing Centre (JSC)Forschungszentrum Jülich 52425 Jülich Germany
- Department of Hematology, Oncology, Hemostaseology and Stem Cell Transplantation University Hospital AachenRWTH Aachen University Pauwelsstraße 30 52074 Aachen Germany
| |
Collapse
|
31
|
Khoshbin Z, Housaindokht MR, Izadyar M, Bozorgmehr MR, Verdian A. Temperature and molecular crowding effects on the sensitivity of T30695 aptamer toward Pb2+ion: a joint molecular dynamics simulation and experimental study. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1751842] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Zahra Khoshbin
- Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran
| | | | - Mohammad Izadyar
- Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran
| | | | - Asma Verdian
- Department of Food Safety and Quality Control, Research Institute of Food Science and Technology (RIFST), Mashhad, Iran
| |
Collapse
|
32
|
Sereda YV, Ortoleva PJ. Temporally Coarse-Grained All-Atom Molecular Dynamics Achieved via Stochastic Padé Approximants. J Phys Chem B 2020; 124:1392-1410. [PMID: 31958947 DOI: 10.1021/acs.jpcb.9b10735] [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
A Padé approximant scheme for realizing the discrete-time evolution of the state of a many-atom system is introduced. This temporal coarse-graining scheme accounts for the underlying Newtonian physics and avoids the need for construction of spatially coarse-grained variables. Newtonian physics is incorporated through short molecular dynamics simulations at the beginning of each of the large coarse-grained timesteps. The balance between stochastic and coherent dynamics expressed by many-atom systems is captured via incorporation of the Ito formula into a Padé approximant for the time dependence of individual atom positions over large timesteps. Since the time for a many-atom system to express a characteristic ensemble of atomic velocity fluctuations is typically short relative to the characteristic time of large-scale atomic displacements, a computationally efficient and accurate temporal coarse-graining of the atom-resolved Newtonian dynamics is formulated, denoted all-atom Padé-Ito molecular dynamics (APIMD). Evolution of the system over a time step much longer than that required for standard molecular dynamics (MD) is achieved via incorporation of information from the short MD simulations into a Padé approximant extrapolation in time. The extrapolated atomic configuration is subjected to energy minimization and, when needed, thermal equilibration so as to avoid occasional unphysical close encounters deriving from the Padé approximant extrapolation and to represent configurations appropriate for the temperature of interest. APIMD is implemented and tested via comparison with traditional MD simulations of five phenomena: (1) pertussis toxin subunit deformation, (2) structural transition in a T = 1 capsid-like structure of HPV16 L1 protein, (3) coalescence of argon nanodroplets, and structural transitions in dialanine in (4) vacuum, and (5) water. Accuracy of APIMD is demonstrated using semimicroscopic descriptors (rmsd, radius of gyration, residue-residue contact maps, and densities) and the free energy. Significant computational acceleration relative to traditional molecular dynamics is illustrated.
Collapse
Affiliation(s)
- Yuriy V Sereda
- Department of Chemistry Indiana University Bloomington , Indiana 47405 , United States
| | - Peter J Ortoleva
- Department of Chemistry Indiana University Bloomington , Indiana 47405 , United States
| |
Collapse
|
33
|
Molecular Dynamics Simulations Suggest a Non-Doublet Decoding Model of -1 Frameshifting by tRNA Ser3. Biomolecules 2019; 9:biom9110745. [PMID: 31752208 PMCID: PMC6920855 DOI: 10.3390/biom9110745] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2019] [Revised: 11/14/2019] [Accepted: 11/14/2019] [Indexed: 12/28/2022] Open
Abstract
In-frame decoding in the ribosome occurs through canonical or wobble Watson-Crick pairing of three mRNA codon bases (a triplet) with a triplet of anticodon bases in tRNA. Departures from the triplet-triplet interaction can result in frameshifting, meaning downstream mRNA codons are then read in a different register. There are many mechanisms to induce frameshifting, and most are insufficiently understood. One previously proposed mechanism is doublet decoding, in which only codon bases 1 and 2 are read by anticodon bases 34 and 35, which would lead to -1 frameshifting. In E. coli, tRNASer3GCU can induce -1 frameshifting at alanine (GCA) codons. The logic of the doublet decoding model is that the Ala codon's GC could pair with the tRNASer3's GC, leaving the third anticodon residue U36 making no interactions with mRNA. Under that model, a U36C mutation would still induce -1 frameshifting, but experiments refute this. We perform all-atom simulations of wild-type tRNASer3, as well as a U36C mutant. Our simulations revealed a hydrogen bond between U36 of the anticodon and G1 of the codon. The U36C mutant cannot make this interaction, as it lacks the hydrogen-bond-donating H3. The simulation thus suggests a novel, non-doublet decoding mechanism for -1 frameshifting by tRNASer3 at Ala codons.
Collapse
|
34
|
Schmidt OP, Jurt S, Johannsen S, Karimi A, Sigel RKO, Luedtke NW. Concerted dynamics of metallo-base pairs in an A/B-form helical transition. Nat Commun 2019; 10:4818. [PMID: 31645548 PMCID: PMC6811676 DOI: 10.1038/s41467-019-12440-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Accepted: 09/05/2019] [Indexed: 01/31/2023] Open
Abstract
Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-HgII-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-HgII-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-HgII-(N4)C base pairs and greater A-form characteristics. Ten-fold different 1J coupling constants (15N,199Hg) are observed for (N3)C-HgII (114 Hz) versus (N4)C-HgII (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-HgII-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding HgII, demonstrating that C-HgII-T has unique switching potential in DNA-based materials and devices.
Collapse
Affiliation(s)
- Olivia P Schmidt
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Simon Jurt
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Silke Johannsen
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Ashkan Karimi
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Roland K O Sigel
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Nathan W Luedtke
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
35
|
Bag S, Maiti PK. Tuning molecular fluctuation to boost the conductance in DNA based molecular wires. Phys Chem Chem Phys 2019; 21:23514-23520. [PMID: 31617554 DOI: 10.1039/c9cp03589c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Inherent molecular fluctuations are known to have a significant influence on the charge transport properties of biomolecules like DNA, PNA and proteins. In this work, we show ways to control these fluctuations and further demonstrate their use to enhance the conductance of two widely studied molecular wires, namely dsDNA (DNA) and G4 Quadruplex (G4-Quad). We quantify the molecular fluctuation in terms of the root mean square deviation (RMSD) of the molecule. In the case of DNA, we use temperature to control the fluctuations, while in the case of G4-Quad the fluctuations are tuned by the ions inside the pore. The electronic coupling between the bases of dsDNA and G4-Quad, which measures the conductance of these molecular wires, shows a non-monotonic behaviour with the increase in fluctuation. We find values of fluctuation which give rise to maximum electronic coupling and hence high conductivity for both the cases. In the case of DNA, these optimal fluctuations (∼2.5 Å) are achieved at a temperature of 210 K, which gives rise to an electronic coupling of 0.135 eV between the DNA bases. The optimal fluctuations in G4-Quad are achieved (∼7 Å) in a 4 base pair long system with 2 Na+ ions inside the pore, giving rise to an electronic coupling of 0.09 eV.
Collapse
Affiliation(s)
- Saientan Bag
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, Karnataka, India.
| | | |
Collapse
|
36
|
Dodda LS, Cabeza de Vaca I, Tirado-Rives J, Jorgensen WL. LigParGen web server: an automatic OPLS-AA parameter generator for organic ligands. Nucleic Acids Res 2019; 45:W331-W336. [PMID: 28444340 PMCID: PMC5793816 DOI: 10.1093/nar/gkx312] [Citation(s) in RCA: 592] [Impact Index Per Article: 118.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Accepted: 04/20/2017] [Indexed: 11/24/2022] Open
Abstract
The accurate calculation of protein/nucleic acid–ligand interactions or condensed phase properties by force field-based methods require a precise description of the energetics of intermolecular interactions. Despite the progress made in force fields, small molecule parameterization remains an open problem due to the magnitude of the chemical space; the most critical issue is the estimation of a balanced set of atomic charges with the ability to reproduce experimental properties. The LigParGen web server provides an intuitive interface for generating OPLS-AA/1.14*CM1A(-LBCC) force field parameters for organic ligands, in the formats of commonly used molecular dynamics and Monte Carlo simulation packages. This server has high value for researchers interested in studying any phenomena based on intermolecular interactions with ligands via molecular mechanics simulations. It is free and open to all at jorgensenresearch.com/ligpargen, and has no login requirements.
Collapse
Affiliation(s)
- Leela S Dodda
- Department of Chemistry, Yale University, New Haven, CT 06520-8107, USA
| | | | | | | |
Collapse
|
37
|
Ekesan Ş, York DM. Framework for Conducting and Analyzing Crystal Simulations of Nucleic Acids to Aid in Modern Force Field Evaluation. J Phys Chem B 2019; 123:4611-4624. [PMID: 31002511 PMCID: PMC6614744 DOI: 10.1021/acs.jpcb.8b11923] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Crystal simulations provide useful tools, along with solution simulations, to test nucleic acid force fields, but should be interpreted with care owing to the difficulty of establishing the environmental conditions needed to reproduce experimental crystal packing. These challenges underscore the need to construct proper protocols for carrying out crystal simulations and analyzing results to identify the origin of deviations from crystallographic data. Toward this end, we introduce a novel framework for B-factor decomposition into additive intramolecular, rotational, and translational atomic fluctuation components and partitioning of each of these components into individual asymmetric unit and lattice contributions. We apply the framework to a benchmark set of A-DNA, Z-DNA, and B-DNA double helix systems of various chain lengths. Overall, the intramolecular deviations from the crystal were quite small (≤1.0 Å), suggesting high accuracy of the force field, whereas crystal packing was not well reproduced. The present work establishes a framework to conduct and analyze crystal simulations that ultimately take on issues of crystal packing and can provide insight into nucleic acid force fields.
Collapse
Affiliation(s)
- Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| |
Collapse
|
38
|
Kührová P, Mlýnský V, Zgarbová M, Krepl M, Bussi G, Best RB, Otyepka M, Šponer J, Banáš P. Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions. J Chem Theory Comput 2019; 15:3288-3305. [PMID: 30896943 PMCID: PMC7491206 DOI: 10.1021/acs.jctc.8b00955] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.
Collapse
Affiliation(s)
- 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
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Marie Zgarbová
- 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
| | - Miroslav Krepl
- 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
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Robert B. Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
| | - Michal Otyepka
- 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
- 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
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- 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
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
39
|
Kruse H, Banáš P, Šponer J. Investigations of Stacked DNA Base-Pair Steps: Highly Accurate Stacking Interaction Energies, Energy Decomposition, and Many-Body Stacking Effects. J Chem Theory Comput 2018; 15:95-115. [DOI: 10.1021/acs.jctc.8b00643] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, 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, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 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, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
40
|
Islam B, Stadlbauer P, Krepl M, Havrila M, Haider S, Sponer J. Structural Dynamics of Lateral and Diagonal Loops of Human Telomeric G-Quadruplexes in Extended MD Simulations. J Chem Theory Comput 2018; 14:5011-5026. [DOI: 10.1021/acs.jctc.8b00543] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Barira Islam
- 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
| | - Miroslav Krepl
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Marek Havrila
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Shozeb Haider
- UCL School of
Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, U.K
| | - Jiri Sponer
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
41
|
Krüger A, Zimbres FM, Kronenberger T, Wrenger C. Molecular Modeling Applied to Nucleic Acid-Based Molecule Development. Biomolecules 2018; 8:E83. [PMID: 30150587 PMCID: PMC6163985 DOI: 10.3390/biom8030083] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 08/12/2018] [Accepted: 08/16/2018] [Indexed: 12/15/2022] Open
Abstract
Molecular modeling by means of docking and molecular dynamics (MD) has become an integral part of early drug discovery projects, enabling the screening and enrichment of large libraries of small molecules. In the past decades, special emphasis was drawn to nucleic acid (NA)-based molecules in the fields of therapy, diagnosis, and drug delivery. Research has increased dramatically with the advent of the SELEX (systematic evolution of ligands by exponential enrichment) technique, which results in single-stranded DNA or RNA sequences that bind with high affinity and specificity to their targets. Herein, we discuss the role and contribution of docking and MD to the development and optimization of new nucleic acid-based molecules. This review focuses on the different approaches currently available for molecular modeling applied to NA interaction with proteins. We discuss topics ranging from structure prediction to docking and MD, highlighting their main advantages and limitations and the influence of flexibility on their calculations.
Collapse
Affiliation(s)
- Arne Krüger
- Unit for Drug Discovery, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, São Paulo, SP 05508-000, Brazil.
| | - Flávia M Zimbres
- Department of Biochemistry and Molecular Biology and Center for Tropical and Emerging Global Diseases, University of Georgia, Athens, GA 30602, USA.
| | - Thales Kronenberger
- Department of Internal Medicine VIII, University Hospital of Tübingen, 72076 Tübingen, Germany.
| | - Carsten Wrenger
- Unit for Drug Discovery, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, São Paulo, SP 05508-000, Brazil.
| |
Collapse
|
42
|
Abstract
The structure of RNA has been a natural subject for mathematical modeling, inviting many innovative computational frameworks. This single-stranded polynucleotide chain can fold upon itself in numerous ways to form hydrogen-bonded segments, imperfect with single-stranded loops. Illustrating these paired and non-paired interaction networks, known as RNA's secondary (2D) structure, using mathematical graph objects has been illuminating for RNA structure analysis. Building upon such seminal work from the 1970s and 1980s, graph models are now used to study not only RNA structure but also describe RNA's recurring modular units, sample the conformational space accessible to RNAs, predict RNA's three-dimensional folds, and apply the combined aspects to novel RNA design. In this article, we outline the development of the RNA-As-Graphs (or RAG) approach and highlight current applications to RNA structure prediction and design.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, 100 Washington Square East, Silver Building, New York University, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, NY 10012, USA; New York University ECNU - Center for Computational Chemistry at NYU Shanghai, 3663 North Zhongshan Road, Shanghai, 200062, China.
| |
Collapse
|
43
|
Kathuria P, Sharma P, Manderville RA, Wetmore SD. Molecular Dynamics Simulations of Mismatched DNA Duplexes Associated with the Major C8-Linked 2′-Deoxyguanosine Adduct of the Food Mutagen Ochratoxin A: Influence of Opposing Base, Adduct Ionization State, and Sequence on the Structure of Damaged DNA. Chem Res Toxicol 2018; 31:712-720. [DOI: 10.1021/acs.chemrestox.8b00064] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Preetleen Kathuria
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Richard A. Manderville
- Departments of Chemistry and Toxicology, University of Guelph, Guelph, Ontario N1G 2W1, Canada
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| |
Collapse
|
44
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
45
|
Maitarad A, Poomsuk N, Vilaivan C, Vilaivan T, Siriwong K. Insight into a conformation of the PNA-PNA duplex with (2′R,4′R)- and (2′R,4′S)-prolyl-(1S,2S)-2-aminocyclopentanecarboxylic acid backbones. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.03.024] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
46
|
New developments in force fields for biomolecular simulations. Curr Opin Struct Biol 2018; 49:129-138. [DOI: 10.1016/j.sbi.2018.02.002] [Citation(s) in RCA: 136] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Revised: 01/30/2018] [Accepted: 02/04/2018] [Indexed: 11/18/2022]
|
47
|
Poblete S, Bottaro S, Bussi G. A nucleobase-centered coarse-grained representation for structure prediction of RNA motifs. Nucleic Acids Res 2018; 46:1674-1683. [PMID: 29272539 PMCID: PMC5829650 DOI: 10.1093/nar/gkx1269] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Revised: 12/05/2017] [Accepted: 12/07/2017] [Indexed: 01/30/2023] Open
Abstract
We introduce the SPlit-and-conQueR (SPQR) model, a coarse-grained (CG) representation of RNA designed for structure prediction and refinement. In our approach, the representation of a nucleotide consists of a point particle for the phosphate group and an anisotropic particle for the nucleoside. The interactions are, in principle, knowledge-based potentials inspired by the $\mathcal {E}$SCORE function, a base-centered scoring function. However, a special treatment is given to base-pairing interactions and certain geometrical conformations which are lost in a raw knowledge-based model. This results in a representation able to describe planar canonical and non-canonical base pairs and base-phosphate interactions and to distinguish sugar puckers and glycosidic torsion conformations. The model is applied to the folding of several structures, including duplexes with internal loops of non-canonical base pairs, tetraloops, junctions and a pseudoknot. For the majority of these systems, experimental structures are correctly predicted at the level of individual contacts. We also propose a method for efficiently reintroducing atomistic detail from the CG representation.
Collapse
Affiliation(s)
- Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| | - Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
48
|
Smith LG, Zhao J, Mathews DH, Turner DH. Physics-based all-atom modeling of RNA energetics and structure. WILEY INTERDISCIPLINARY REVIEWS-RNA 2018; 8. [PMID: 28815951 DOI: 10.1002/wrna.1422] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Revised: 02/03/2017] [Accepted: 03/08/2017] [Indexed: 12/31/2022]
Abstract
The database of RNA sequences is exploding, but knowledge of energetics, structures, and dynamics lags behind. All-atom computational methods, such as molecular dynamics, hold promise for closing this gap. New algorithms and faster computers have accelerated progress in improving the reliability and accuracy of predictions. Currently, the methods can facilitate refinement of experimentally determined nuclear magnetic resonance and x-ray structures, but are 'unreliable' for predictions based only on sequence. Much remains to be discovered, however, about the many molecular interactions driving RNA folding and the best way to approximate them quantitatively. The large number of parameters required means that a wide variety of experimental results will be required to benchmark force fields and different approaches. As computational methods become more reliable and accessible, they will be used by an increasing number of biologists, much as x-ray crystallography has expanded. Thus, many fundamental physical principles underlying the computational methods are described. This review presents a summary of the current state of molecular dynamics as applied to RNA. It is designed to be helpful to students, postdoctoral fellows, and faculty who are considering or starting computational studies of RNA. WIREs RNA 2017, 8:e1422. doi: 10.1002/wrna.1422.
Collapse
Affiliation(s)
- Louis G Smith
- Department of Biochemistry and Biophysics and Center for RNA Biology, School of Medicine and Dentistry, University of Rochester, Rochester, NY, USA
| | - Jianbo Zhao
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, NY, USA
| | - David H Mathews
- Department of Biochemistry and Biophysics and Center for RNA Biology, School of Medicine and Dentistry, University of Rochester, Rochester, NY, USA
| | - Douglas H Turner
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, NY, USA
| |
Collapse
|
49
|
Šponer JE, Szabla R, Góra RW, Saitta AM, Pietrucci F, Saija F, Di Mauro E, Saladino R, Ferus M, Civiš S, Šponer J. Prebiotic synthesis of nucleic acids and their building blocks at the atomic level - merging models and mechanisms from advanced computations and experiments. Phys Chem Chem Phys 2018; 18:20047-66. [PMID: 27136968 DOI: 10.1039/c6cp00670a] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The origin of life on Earth is one of the most fascinating questions of contemporary science. Extensive research in the past decades furnished diverse experimental proposals for the emergence of first informational polymers that could form the basis of the early terrestrial life. Side by side with the experiments, the fast development of modern computational chemistry methods during the last 20 years facilitated the use of in silico modelling tools to complement the experiments. Modern computations can provide unique atomic-level insights into the structural and electronic aspects as well as the energetics of key prebiotic chemical reactions. Many of these insights are not directly obtainable from the experimental techniques and the computations are thus becoming indispensable for proper interpretation of many experiments and for qualified predictions. This review illustrates the synergy between experiment and theory in the origin of life research focusing on the prebiotic synthesis of various nucleic acid building blocks and on the self-assembly of nucleotides leading to the first functional oligonucleotides.
Collapse
Affiliation(s)
- Judit E Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| | - Rafał Szabla
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic.
| | - Robert W Góra
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - A Marco Saitta
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Fabio Pietrucci
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Franz Saija
- CNR-IPCF, Viale Ferdinando Stagno d'Alcontres 37, 98158 Messina, Italy
| | - Ernesto Di Mauro
- Dipartimento di Biologia e Biotecnologie "Charles Darwin", "Sapienza" Università di Roma, Piazzale Aldo Moro 5, Rome 00185, Italy
| | - Raffaele Saladino
- Dipartimento di Scienze Ecologiche e Biologiche Università della Tuscia, Via San Camillo De Lellis, 01100 Viterbo, Italy
| | - Martin Ferus
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Svatopluk Civiš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| |
Collapse
|
50
|
Abstract
The complex and often highly dynamic 3D structures of RNA molecules are central to their diverse cellular functions. Molecular dynamics (MD) simulations have played a major role in characterizing the structure and dynamics of proteins, but the physical models (“force fields”) used for simulating nucleic acids are substantially less accurate overall than those used in protein simulations, creating a major challenge for MD studies of RNA. Here, we report an RNA force field capable of describing the structural and thermodynamic properties of RNA molecules with accuracy comparable to state-of-the-art protein force fields. This force field should facilitate the use of MD simulation as a tool for the study of biologically significant RNA molecules and protein–RNA complexes. Molecular dynamics (MD) simulation has become a powerful tool for characterizing at an atomic level of detail the conformational changes undergone by proteins. The application of such simulations to RNA structures, however, has proven more challenging, due in large part to the fact that the physical models (“force fields”) available for MD simulations of RNA molecules are substantially less accurate in many respects than those currently available for proteins. Here, we introduce an extensive revision of a widely used RNA force field in which the parameters have been modified, based on quantum mechanical calculations and existing experimental information, to more accurately reflect the fundamental forces that stabilize RNA structures. We evaluate these revised parameters through long-timescale MD simulations of a set of RNA molecules that covers a wide range of structural complexity, including single-stranded RNAs, RNA duplexes, RNA hairpins, and riboswitches. The structural and thermodynamic properties measured in these simulations exhibited dramatically improved agreement with experimentally determined values. Based on the comparisons we performed, this RNA force field appears to achieve a level of accuracy comparable to that of state-of-the-art protein force fields, thus significantly advancing the utility of MD simulation as a tool for elucidating the structural dynamics and function of RNA molecules and RNA-containing biological assemblies.
Collapse
|