1
|
Walbrun A, Wang T, Matthies M, Šulc P, Simmel FC, Rief M. Single-molecule force spectroscopy of toehold-mediated strand displacement. Nat Commun 2024; 15:7564. [PMID: 39217165 PMCID: PMC11365964 DOI: 10.1038/s41467-024-51813-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Accepted: 08/16/2024] [Indexed: 09/04/2024] Open
Abstract
Toehold-mediated strand displacement (TMSD) is extensively utilized in dynamic DNA nanotechnology and for a wide range of DNA or RNA-based reaction circuits. Investigation of TMSD kinetics typically relies on bulk fluorescence measurements providing effective, bulk-averaged reaction rates. Information on individual molecules or even base pairs is scarce. In this work, we explore the dynamics of strand displacement processes at the single-molecule level using single-molecule force spectroscopy with a microfluidics-enhanced optical trap supported by state-of-the-art coarse-grained simulations. By applying force, we can trigger and observe TMSD in real-time with microsecond and nanometer resolution. We find TMSD proceeds very rapidly under load with single step times of 1 µs. Tuning invasion efficiency by introducing mismatches allows studying thousands of forward/backward invasion events on a single molecule and analyze the kinetics of the invasion process. Extrapolation to zero force reveals single step times for DNA invading DNA four times faster than for RNA invading RNA. We also study the kinetics of DNA invading RNA, a process that in the absence of force would rarely occur. Our results reveal the importance of sequence effects for the TMSD process and have relevance for a wide range of applications in nucleic acid nanotechnology and synthetic biology.
Collapse
Affiliation(s)
- Andreas Walbrun
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Center for Functional Protein Assemblies (CPA), Garching, Germany
| | - Tianhe Wang
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Garching, Germany
| | - Michael Matthies
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Garching, Germany
| | - Petr Šulc
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Garching, Germany
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, Arizona, USA
| | - Friedrich C Simmel
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Garching, Germany.
| | - Matthias Rief
- Technical University of Munich, TUM School of Natural Sciences, Department of Bioscience, Center for Functional Protein Assemblies (CPA), Garching, Germany.
| |
Collapse
|
2
|
Elonen A, Wimbes L, Mohammed A, Orponen P. DNAforge: a design tool for nucleic acid wireframe nanostructures. Nucleic Acids Res 2024; 52:W13-W18. [PMID: 38747339 PMCID: PMC11223811 DOI: 10.1093/nar/gkae367] [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] [Received: 02/17/2024] [Revised: 04/09/2024] [Accepted: 04/29/2024] [Indexed: 07/06/2024] Open
Abstract
DNAforge is an online tool that provides a unified, user-friendly interface to several recent design methods for DNA and RNA wireframe nanostructures, with the possibility of integrating additional methods into the same framework. Currently, DNAforge supports three design methods for DNA nanostructures and two for RNA nanostructures. The tool enables the design, visualisation and sequence generation for highly complex wireframe nanostructures with a simple fully automated process. DNAforge is freely accessible at https://dnaforge.org/.
Collapse
Affiliation(s)
- Antti Elonen
- Department of Computer Science, Aalto University, Finland
| | - Leon Wimbes
- Department of Computer Science, Aalto University, Finland
- Department of Mathematics and Computer Science, Philipps University of Marburg, Germany
| | | | - Pekka Orponen
- Department of Computer Science, Aalto University, Finland
| |
Collapse
|
3
|
Todisco M, Radakovic A, Szostak JW. RNA Complexes with Nicks and Gaps: Thermodynamic and Kinetic Effects of Coaxial Stacking and Dangling Ends. J Am Chem Soc 2024; 146:18083-18094. [PMID: 38904115 PMCID: PMC11229006 DOI: 10.1021/jacs.4c05115] [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: 04/14/2024] [Revised: 06/13/2024] [Accepted: 06/14/2024] [Indexed: 06/22/2024]
Abstract
Multiple RNA strands can interact in solution and assume a large variety of configurations dictated by their potential for base pairing. Although duplex formation from two complementary oligonucleotides has been studied in detail, we still lack a systematic characterization of the behavior of higher order complexes. Here, we focus on the thermodynamic and kinetic effects of an upstream oligonucleotide on the binding of a downstream oligonucleotide to a common template, as we vary the sequence and structure of the contact interface. We show that coaxial stacking in RNA is well correlated with but much more stabilizing than helix propagation over an analogous intact double helix step (median ΔΔG°37 °C ≈ 1.7 kcal/mol). Consequently, approximating coaxial stacking in RNA with the helix propagation term leads to large discrepancies between predictions and our experimentally determined melting temperatures, with an offset of ≈10 °C. Our kinetic study reveals that the hybridization of the downstream probe oligonucleotide is impaired (lower kon) by the presence of the upstream oligonucleotide, with the thermodynamic stabilization coming entirely from an extended lifetime (lower koff) of the bound downstream oligonucleotide, which can increase from seconds to months. Surprisingly, we show that the effect of nicks is dependent on the length of the stacking oligonucleotides, and we discuss the binding of ultrashort (1-4 nt) oligonucleotides that are relevant in the context of the origin of life. The thermodynamic and kinetic data obtained in this work allow for the prediction of the formation and stability of higher-order multistranded complexes.
Collapse
Affiliation(s)
- Marco Todisco
- Howard Hughes Medical Institute, Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Aleksandar Radakovic
- Howard Hughes Medical Institute, Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Harvard Medical School,Department of Genetics, Boston, Massachusetts 02115, United States
| | - Jack W Szostak
- Howard Hughes Medical Institute, Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
4
|
Bernard C, Postic G, Ghannay S, Tahi F. State-of-the-RNArt: benchmarking current methods for RNA 3D structure prediction. NAR Genom Bioinform 2024; 6:lqae048. [PMID: 38745991 PMCID: PMC11091930 DOI: 10.1093/nargab/lqae048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Revised: 04/05/2024] [Accepted: 05/08/2024] [Indexed: 05/16/2024] Open
Abstract
RNAs are essential molecules involved in numerous biological functions. Understanding RNA functions requires the knowledge of their 3D structures. Computational methods have been developed for over two decades to predict the 3D conformations from RNA sequences. These computational methods have been widely used and are usually categorised as either ab initio or template-based. The performances remain to be improved. Recently, the rise of deep learning has changed the sight of novel approaches. Deep learning methods are promising, but their adaptation to RNA 3D structure prediction remains difficult. In this paper, we give a brief review of the ab initio, template-based and novel deep learning approaches. We highlight the different available tools and provide a benchmark on nine methods using the RNA-Puzzles dataset. We provide an online dashboard that shows the predictions made by benchmarked methods, freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr/evryrna/state_of_the_rnart/.
Collapse
Affiliation(s)
- Clément Bernard
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Guillaume Postic
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Sahar Ghannay
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Fariza Tahi
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| |
Collapse
|
5
|
Tosti Guerra F, Poppleton E, Šulc P, Rovigatti L. ANNaMo: Coarse-grained modeling for folding and assembly of RNA and DNA systems. J Chem Phys 2024; 160:205102. [PMID: 38814009 DOI: 10.1063/5.0202829] [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: 02/07/2024] [Accepted: 05/04/2024] [Indexed: 05/31/2024] Open
Abstract
The folding of RNA and DNA strands plays crucial roles in biological systems and bionanotechnology. However, studying these processes with high-resolution numerical models is beyond current computational capabilities due to the timescales and system sizes involved. In this article, we present a new coarse-grained model for investigating the folding dynamics of nucleic acids. Our model represents three nucleotides with a patchy particle and is parameterized using well-established nearest-neighbor models. Thanks to the reduction of degrees of freedom and to a bond-swapping mechanism, our model allows for simulations at timescales and length scales that are currently inaccessible to more detailed models. To validate the performance of our model, we conducted extensive simulations of various systems: We examined the thermodynamics of DNA hairpins, capturing their stability and structural transitions, the folding of an MMTV pseudoknot, which is a complex RNA structure involved in viral replication, and also explored the folding of an RNA tile containing a k-type pseudoknot. Finally, we evaluated the performance of the new model in reproducing the melting temperatures of oligomers and the dependence on the toehold length of the displacement rate in toehold-mediated displacement processes, a key reaction used in molecular computing. All in all, the successful reproduction of experimental data and favorable comparisons with existing coarse-grained models validate the effectiveness of the new model.
Collapse
Affiliation(s)
- F Tosti Guerra
- Department of Physics, Sapienza University of Rome, Roma, Italy
| | - E Poppleton
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, Arizona 85281, USA
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - P Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, Arizona 85281, USA
- Department of Bioscience, School of Natural Sciences, Technical University Munich, Munich, Germany
| | - L Rovigatti
- Department of Physics, Sapienza University of Rome, Roma, Italy
| |
Collapse
|
6
|
Koksaldi I, Park D, Atilla A, Kang H, Kim J, Seker UOS. RNA-Based Sensor Systems for Affordable Diagnostics in the Age of Pandemics. ACS Synth Biol 2024; 13:1026-1037. [PMID: 38588603 PMCID: PMC11036506 DOI: 10.1021/acssynbio.3c00698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Revised: 03/25/2024] [Accepted: 03/25/2024] [Indexed: 04/10/2024]
Abstract
In the era of the COVID-19 pandemic, the significance of point-of-care (POC) diagnostic tools has become increasingly vital, driven by the need for quick and precise virus identification. RNA-based sensors, particularly toehold sensors, have emerged as promising candidates for POC detection systems due to their selectivity and sensitivity. Toehold sensors operate by employing an RNA switch that changes the conformation when it binds to a target RNA molecule, resulting in a detectable signal. This review focuses on the development and deployment of RNA-based sensors for POC viral RNA detection with a particular emphasis on toehold sensors. The benefits and limits of toehold sensors are explored, and obstacles and future directions for improving their performance within POC detection systems are presented. The use of RNA-based sensors as a technology for rapid and sensitive detection of viral RNA holds great potential for effectively managing (dealing/coping) with present and future pandemics in resource-constrained settings.
Collapse
Affiliation(s)
- Ilkay
Cisil Koksaldi
- UNAM
− Institute of Materials Science and Nanotechnology, National
Nanotechnology Research Center (UNAM), Bilkent
University, Ankara 06800, Turkey
| | - Dongwon Park
- Department
of Life Sciences, Pohang University of Science
and Technology, Pohang 37673, South Korea
| | - Abdurahman Atilla
- UNAM
− Institute of Materials Science and Nanotechnology, National
Nanotechnology Research Center (UNAM), Bilkent
University, Ankara 06800, Turkey
| | - Hansol Kang
- Department
of Life Sciences, Pohang University of Science
and Technology, Pohang 37673, South Korea
| | - Jongmin Kim
- Department
of Life Sciences, Pohang University of Science
and Technology, Pohang 37673, South Korea
| | - Urartu Ozgur Safak Seker
- UNAM
− Institute of Materials Science and Nanotechnology, National
Nanotechnology Research Center (UNAM), Bilkent
University, Ankara 06800, Turkey
| |
Collapse
|
7
|
Sinzger-D'Angelo M, Hanst M, Reinhardt F, Koeppl H. Effects of mRNA conformational switching on translational noise in gene circuits. J Chem Phys 2024; 160:134108. [PMID: 38573847 DOI: 10.1063/5.0186927] [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/09/2023] [Accepted: 03/08/2024] [Indexed: 04/06/2024] Open
Abstract
Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.
Collapse
Affiliation(s)
| | - Maleen Hanst
- Centre for Synthetic Biology, Technische Universität Darmstadt, Darmstadt, Germany
| | - Felix Reinhardt
- Centre for Synthetic Biology, Technische Universität Darmstadt, Darmstadt, Germany
| | - Heinz Koeppl
- Centre for Synthetic Biology, Technische Universität Darmstadt, Darmstadt, Germany
| |
Collapse
|
8
|
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
|
9
|
Rinaldi S, Moroni E, Rozza R, Magistrato A. Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation. J Chem Theory Comput 2024; 20:993-1018. [PMID: 38287883 DOI: 10.1021/acs.jctc.3c01239] [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: 01/31/2024]
Abstract
Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.
Collapse
Affiliation(s)
- Silvia Rinaldi
- National Research Council of Italy (CNR) - Institute of Chemistry of OrganoMetallic Compounds (ICCOM), c/o Area di Ricerca CNR di Firenze Via Madonna del Piano 10, 50019 Sesto Fiorentino, Florence, Italy
| | - Elisabetta Moroni
- National Research Council of Italy (CNR) - Institute of Chemical Sciences and Technologies (SCITEC), via Mario Bianco 9, 20131 Milano, Italy
| | - Riccardo Rozza
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| |
Collapse
|
10
|
Zhang S, Li J, Chen SJ. Machine learning in RNA structure prediction: Advances and challenges. Biophys J 2024:S0006-3495(24)00067-5. [PMID: 38297836 DOI: 10.1016/j.bpj.2024.01.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/08/2024] [Accepted: 01/24/2024] [Indexed: 02/02/2024] Open
Abstract
RNA molecules play a crucial role in various biological processes, with their functionality closely tied to their structures. The remarkable advancements in machine learning techniques for protein structure prediction have shown promise in the field of RNA structure prediction. In this perspective, we discuss the advances and challenges encountered in constructing machine learning-based models for RNA structure prediction. We explore topics including model building strategies, specific challenges involved in predicting RNA secondary (2D) and tertiary (3D) structures, and approaches to these challenges. In addition, we highlight the advantages and challenges of constructing RNA language models. Given the rapid advances of machine learning techniques, we anticipate that machine learning-based models will serve as important tools for predicting RNA structures, thereby enriching our understanding of RNA structures and their corresponding functions.
Collapse
Affiliation(s)
- Sicheng Zhang
- Department of Physics and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri
| | - Jun Li
- Department of Physics and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri
| | - Shi-Jie Chen
- Department of Physics and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri; Department of Biochemistry, University of Missouri, Columbia, Missouri.
| |
Collapse
|
11
|
Bernard C, Postic G, Ghannay S, Tahi F. RNAdvisor: a comprehensive benchmarking tool for the measure and prediction of RNA structural model quality. Brief Bioinform 2024; 25:bbae064. [PMID: 38436560 PMCID: PMC10939302 DOI: 10.1093/bib/bbae064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/30/2024] [Accepted: 02/02/2024] [Indexed: 03/05/2024] Open
Abstract
RNA is a complex macromolecule that plays central roles in the cell. While it is well known that its structure is directly related to its functions, understanding and predicting RNA structures is challenging. Assessing the real or predictive quality of a structure is also at stake with the complex 3D possible conformations of RNAs. Metrics have been developed to measure model quality while scoring functions aim at assigning quality to guide the discrimination of structures without a known and solved reference. Throughout the years, many metrics and scoring functions have been developed, and no unique assessment is used nowadays. Each developed assessment method has its specificity and might be complementary to understanding structure quality. Therefore, to evaluate RNA 3D structure predictions, it would be important to calculate different metrics and/or scoring functions. For this purpose, we developed RNAdvisor, a comprehensive automated software that integrates and enhances the accessibility of existing metrics and scoring functions. In this paper, we present our RNAdvisor tool, as well as state-of-the-art existing metrics, scoring functions and a set of benchmarks we conducted for evaluating them. Source code is freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr.
Collapse
Affiliation(s)
- Clement Bernard
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Guillaume Postic
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Sahar Ghannay
- LISN - CNRS/Université Paris-Saclay, France, 91400 Orsay, France
| | - Fariza Tahi
- Université Paris Saclay, Univ Evry, IBISC, 91020 Evry-Courcouronnes, France
| |
Collapse
|
12
|
Chou HY, Aksimentiev A. RNA regulates cohesiveness and porosity of a biological condensate. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.09.574811. [PMID: 38260307 PMCID: PMC10802450 DOI: 10.1101/2024.01.09.574811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Biological condensates have emerged as key elements of a biological cell function, concentrating disparate biomolecules to accomplish specific biological tasks. RNA was identified as a key ingredient of such condensates, however, its effect on the physical properties of the condensate was found to depend on the condensate's composition while its effect on the microstructure has remained elusive. Here, we characterize the physical properties and the microstructure of a protein-RNA condensate by means of large-scale coarse-grained (CG) molecular dynamics simulations. By developing a custom CG model of RNA compatible with a popular CG model of proteins, we systematically investigate the structural, thermodynamic, and kinetic properties of condensate droplets containing thousands of individual protein and RNA molecules over a range of temperatures. While we find RNA to increase the condensate's cohesiveness, its effect on the condensate's fluidity is more nuanced with longer molecules compacting the condensate and making it less fluid. We show that a biological condensate has a sponge-like morphology of interconnected channels of size that increases with temperature and decreases in the presence of RNA. Our results suggest that longer RNA form a dynamic scaffold within a condensate, regulating not only its fluidity but also permeability to intruder molecules.
Collapse
|
13
|
Wang F, Xia R, Su Y, Cai P, Xu X. Quantifying RNA structures and interactions with a unified reduced chain representation model. Int J Biol Macromol 2023; 253:127181. [PMID: 37793523 DOI: 10.1016/j.ijbiomac.2023.127181] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Revised: 08/30/2023] [Accepted: 09/25/2023] [Indexed: 10/06/2023]
Abstract
RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.
Collapse
Affiliation(s)
- Fengfei Wang
- Institute of Bioinformatics and Medical Engineering, School of Mathematics and Physics, Jiangsu University of Technology, Changzhou 213001, China
| | - Renjie Xia
- Institute of Bioinformatics and Medical Engineering, School of Mathematics and Physics, Jiangsu University of Technology, Changzhou 213001, China
| | - Yangyang Su
- Institute of Bioinformatics and Medical Engineering, School of Mathematics and Physics, Jiangsu University of Technology, Changzhou 213001, China
| | - Pinggen Cai
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, School of Mathematics and Physics, Jiangsu University of Technology, Changzhou 213001, China.
| |
Collapse
|
14
|
Dabin A, Stirnemann G. Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation. J Phys Chem B 2023. [PMID: 37389985 DOI: 10.1021/acs.jpcb.3c00908] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/02/2023]
Abstract
RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.
Collapse
Affiliation(s)
- Aimeric Dabin
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005, Paris, France
| |
Collapse
|
15
|
Mu ZC, Tan YL, Liu J, Zhang BG, Shi YZ. Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding. Molecules 2023; 28:4833. [PMID: 37375388 DOI: 10.3390/molecules28124833] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 06/11/2023] [Accepted: 06/14/2023] [Indexed: 06/29/2023] Open
Abstract
DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.
Collapse
Affiliation(s)
- Zi-Chun Mu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| |
Collapse
|
16
|
Li J, Chen SJ. RNAJP: enhanced RNA 3D structure predictions with non-canonical interactions and global topology sampling. Nucleic Acids Res 2023; 51:3341-3356. [PMID: 36864729 PMCID: PMC10123122 DOI: 10.1093/nar/gkad122] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 01/14/2023] [Accepted: 02/25/2023] [Indexed: 03/04/2023] Open
Abstract
RNA 3D structures are critical for understanding their functions. However, only a limited number of RNA structures have been experimentally solved, so computational prediction methods are highly desirable. Nevertheless, accurate prediction of RNA 3D structures, especially those containing multiway junctions, remains a significant challenge, mainly due to the complicated non-canonical base pairing and stacking interactions in the junction loops and the possible long-range interactions between loop structures. Here we present RNAJP ('RNA Junction Prediction'), a nucleotide- and helix-level coarse-grained model for the prediction of RNA 3D structures, particularly junction structures, from a given 2D structure. Through global sampling of the 3D arrangements of the helices in junctions using molecular dynamics simulations and in explicit consideration of non-canonical base pairing and base stacking interactions as well as long-range loop-loop interactions, the model can provide significantly improved predictions for multibranched junction structures than existing methods. Moreover, integrated with additional restraints from experiments, such as junction topology and long-range interactions, the model may serve as a useful structure generator for various applications.
Collapse
Affiliation(s)
- Jun Li
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, Columbia, MO 65211, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, Columbia, MO 65211, USA
| |
Collapse
|
17
|
Wang X, Tan YL, Yu S, Shi YZ, Tan ZJ. Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions. Biophys J 2023; 122:1503-1516. [PMID: 36924021 PMCID: PMC10147842 DOI: 10.1016/j.bpj.2023.03.017] [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: 10/24/2022] [Revised: 03/05/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.
Collapse
Affiliation(s)
- Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
18
|
Yang C, Song X, Feng Y, Zhao G, Liu Y. Stability of DNA and RNA hairpins: a comparative study based on ox-DNA. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:265101. [PMID: 36972608 DOI: 10.1088/1361-648x/acc7eb] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2022] [Accepted: 03/27/2023] [Indexed: 06/18/2023]
Abstract
Advances in single-molecule experiments on macromolecular crowding urgently need an efficient simulation method to resolve their discrepancies quantitatively. Ox-DNA model has been since reworked to treat the thermodynamics and mechanical properties of DNA/RNA hairpin at a stretching force. In hopping experiments, the critical forces of RNA hairpins at different temperatures are greater than those of DNA hairpins, in addition, the Gibbs free energy at a fixed temperature required to convert an RNA hairpin into a single-stranded molecule at zero force is obviously greater than that of DNA hairpin and gradually decreases by increasing the temperature. As far as force-ramping experiments are concerned, the first-rupture forces of RNA/DNA hairpins corresponding to the maximum probability density linearly pertain to the force-loading rate, with those of RNA hairpins being greater. The extended ox-DNA model could potentially identify the interaction between biologically inert polymer and RNA/DNA hairpins in crowded environments.
Collapse
Affiliation(s)
- Chao Yang
- College of Physics, Guizhou University, Guiyang, Guizhou 550025, People's Republic of China
| | - Xiaoya Song
- College of Physics, Guizhou University, Guiyang, Guizhou 550025, People's Republic of China
| | - Yuyu Feng
- College of Physics, Guizhou University, Guiyang, Guizhou 550025, People's Republic of China
| | - Guangju Zhao
- College of Physics, Guizhou University, Guiyang, Guizhou 550025, People's Republic of China
| | - Yanhui Liu
- College of Physics, Guizhou University, Guiyang, Guizhou 550025, People's Republic of China
- Breeding and Reproduction in the Plateau Mountainous Region, Ministry of Education (Guizhou University), Guiyang, Guizhou 550025, People's Republic of China
- State Key Laboratory of Public Big Data, Guizhou University, Guiyang 550025, People's Republic of China
- Kechuang Industrial Development Company Limited,Gui'an New Area, Guiyang 550025, People's Republic of China
| |
Collapse
|
19
|
Doye JPK, Fowler H, Prešern D, Bohlin J, Rovigatti L, Romano F, Šulc P, Wong CK, Louis AA, Schreck JS, Engel MC, Matthies M, Benson E, Poppleton E, Snodin BEK. The oxDNA Coarse-Grained Model as a Tool to Simulate DNA Origami. Methods Mol Biol 2023; 2639:93-112. [PMID: 37166713 DOI: 10.1007/978-1-0716-3028-0_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
This chapter introduces how to run molecular dynamics simulations for DNA origami using the oxDNA coarse-grained model.
Collapse
Affiliation(s)
- Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford, UK.
| | - Hannah Fowler
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford, UK
| | - Domen Prešern
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford, UK
| | - Joakim Bohlin
- Department of Physics, Clarendon Laboratory, University of Oxford, Oxford, UK
| | | | - Flavio Romano
- Dipartimento di Fisica, Sapienza Universitá di Roma, Rome, Italy
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, USA
| | - Chak Kui Wong
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford, UK
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Parks Road, Oxford, UK
| | - John S Schreck
- Computational and Information Systems Laboratory, National Center for Atmospheric Research (NCAR), Boulder, USA
| | - Megan C Engel
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| | - Michael Matthies
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, USA
| | - Erik Benson
- Department of Physics, Clarendon Laboratory, University of Oxford, Oxford, UK
| | - Erik Poppleton
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, USA
| | - Benedict E K Snodin
- Department of Philosophy, Future of Humanity Institute, University of Oxford, Oxford, UK
| |
Collapse
|
20
|
Poppleton E, Urbanek N, Chakraborty T, Griffo A, Monari L, Göpfrich K. RNA origami: design, simulation and application. RNA Biol 2023; 20:510-524. [PMID: 37498217 PMCID: PMC10376919 DOI: 10.1080/15476286.2023.2237719] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Revised: 06/20/2023] [Accepted: 07/12/2023] [Indexed: 07/28/2023] Open
Abstract
Design strategies for DNA and RNA nanostructures have developed along parallel lines for the past 30 years, from small structural motifs derived from biology to large 'origami' structures with thousands to tens of thousands of bases. With the recent publication of numerous RNA origami structures and improved design methods-even permitting co-transcriptional folding of kilobase-sized structures - the RNA nanotechnolgy field is at an inflection point. Here, we review the key achievements which inspired and enabled RNA origami design and draw comparisons with the development and applications of DNA origami structures. We further present the available computational tools for the design and the simulation, which will be key to the growth of the RNA origami community. Finally, we portray the transition from RNA origami structure to function. Several functional RNA origami structures exist already, their expression in cells has been demonstrated and first applications in cell biology have already been realized. Overall, we foresee that the fast-paced RNA origami field will provide new molecular hardware for biophysics, synthetic biology and biomedicine, complementing the DNA origami toolbox.
Collapse
Affiliation(s)
- Erik Poppleton
- Biophysical Engineering Group, Center for Molecular Biology of Heidelberg University (ZMBH), Heidelberg University, Heidelberg, Germany
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
- Molecular Biomechanics, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany
| | - Niklas Urbanek
- Biophysical Engineering Group, Center for Molecular Biology of Heidelberg University (ZMBH), Heidelberg University, Heidelberg, Germany
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Taniya Chakraborty
- Biophysical Engineering Group, Center for Molecular Biology of Heidelberg University (ZMBH), Heidelberg University, Heidelberg, Germany
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Alessandra Griffo
- Biophysical Engineering Group, Center for Molecular Biology of Heidelberg University (ZMBH), Heidelberg University, Heidelberg, Germany
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Luca Monari
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
- Institut de Science Et D’ingénierie Supramoléculaires (ISIS), Université de Strasbourg, Strasbourg, France
| | - Kerstin Göpfrich
- Biophysical Engineering Group, Center for Molecular Biology of Heidelberg University (ZMBH), Heidelberg University, Heidelberg, Germany
- Biophysical Engineering Group, Max Planck Institute for Medical Research, Heidelberg, Germany
| |
Collapse
|
21
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
22
|
Elonen A, Natarajan AK, Kawamata I, Oesinghaus L, Mohammed A, Seitsonen J, Suzuki Y, Simmel FC, Kuzyk A, Orponen P. Algorithmic Design of 3D Wireframe RNA Polyhedra. ACS NANO 2022; 16:16608-16616. [PMID: 36178116 PMCID: PMC9620399 DOI: 10.1021/acsnano.2c06035] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2022] [Accepted: 09/26/2022] [Indexed: 06/01/2023]
Abstract
We address the problem of de novo design and synthesis of nucleic acid nanostructures, a challenge that has been considered in the area of DNA nanotechnology since the 1980s and more recently in the area of RNA nanotechnology. Toward this goal, we introduce a general algorithmic design process and software pipeline for rendering 3D wireframe polyhedral nanostructures in single-stranded RNA. To initiate the pipeline, the user creates a model of the desired polyhedron using standard 3D graphic design software. As its output, the pipeline produces an RNA nucleotide sequence whose corresponding RNA primary structure can be transcribed from a DNA template and folded in the laboratory. As case examples, we design and characterize experimentally three 3D RNA nanostructures: a tetrahedron, a triangular bipyramid, and a triangular prism. The design software is openly available and also provides an export of the targeted 3D structure into the oxDNA molecular dynamics simulator for easy simulation and visualization.
Collapse
Affiliation(s)
- Antti Elonen
- Department
of Computer Science, Aalto University, 00076 Aalto, Finland
| | | | - Ibuki Kawamata
- Department
of Robotics, Graduate School of Engineering, Tohoku University, Sendai 980-8577, Japan
- Natural
Science Division, Faculty of Core Research, Ochanomizu University, Tokyo 112-8610, Japan
| | - Lukas Oesinghaus
- Physics
Department E14, Technical University Munich, 85748 Garching, Germany
| | - Abdulmelik Mohammed
- Department
of Computer Science, Aalto University, 00076 Aalto, Finland
- Department
of Biomedical Engineering, San José
State University, San José, California 95192, United States
| | - Jani Seitsonen
- Department
of Applied Physics and Nanomicroscopy Center, Aalto University, 00076 Aalto, Finland
| | - Yuki Suzuki
- Department
of Robotics, Graduate School of Engineering, Tohoku University, Sendai 980-8577, Japan
- Frontier
Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai 980-8577, Japan
- Division
of Chemistry for Materials, Graduate School of Engineering, Mie University, Tsu 514-8507, Japan
| | - Friedrich C. Simmel
- Physics
Department E14, Technical University Munich, 85748 Garching, Germany
| | - Anton Kuzyk
- Department
of Neuroscience and Biomedical Engineering, Aalto University, 00076 Aalto, Finland
| | - Pekka Orponen
- Department
of Computer Science, Aalto University, 00076 Aalto, Finland
| |
Collapse
|
23
|
Zhou L, Wang X, Yu S, Tan YL, Tan ZJ. FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures. Biophys J 2022; 121:3381-3392. [PMID: 35978551 PMCID: PMC9515226 DOI: 10.1016/j.bpj.2022.08.017] [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: 04/30/2022] [Revised: 07/19/2022] [Accepted: 08/15/2022] [Indexed: 11/23/2022] Open
Abstract
Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.
Collapse
Affiliation(s)
- Li Zhou
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan 430073, China.
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China.
| |
Collapse
|
24
|
Sinha S, Tam B, Wang SM. Applications of Molecular Dynamics Simulation in Protein Study. MEMBRANES 2022; 12:membranes12090844. [PMID: 36135863 PMCID: PMC9505860 DOI: 10.3390/membranes12090844] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2022] [Revised: 08/24/2022] [Accepted: 08/25/2022] [Indexed: 05/29/2023]
Abstract
Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.
Collapse
|
25
|
Russo J, Romano F, Kroc L, Sciortino F, Rovigatti L, Šulc P. SAT-assembly: a new approach for designing self-assembling systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:354002. [PMID: 35148521 DOI: 10.1088/1361-648x/ac5479] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Accepted: 02/11/2022] [Indexed: 06/14/2023]
Abstract
We propose a general framework for solving inverse self-assembly problems, i.e. designing interactions between elementary units such that they assemble spontaneously into a predetermined structure. Our approach uses patchy particles as building blocks, where the different units bind at specific interaction sites (the patches), and we exploit the possibility of having mixtures with several components. The interaction rules between the patches is determined by transforming the combinatorial problem into a Boolean satisfiability problem (SAT) which searches for solutions where all bonds are formed in the target structure. Additional conditions, such as the non-satisfiability of competing structures (e.g. metastable states) can be imposed, allowing to effectively design the assembly path in order to avoid kinetic traps. We demonstrate this approach by designing and numerically simulating a cubic diamond structure from four particle species that assembles without competition from other polymorphs, including the hexagonal structure.
Collapse
Affiliation(s)
- John Russo
- Dipartimento di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy
| | - Flavio Romano
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari di Venezia Campus Scientifico, Edificio Alfa, via Torino 155, 30170 Venezia Mestre, Italy
- European Centre for Living Technology (ECLT) Ca' Bottacin, 3911 Dorsoduro Calle Crosera, 30123 Venice, Italy
| | - Lukáš Kroc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, United States of America
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy
| | - Lorenzo Rovigatti
- Dipartimento di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy
- Institute for Complex Systems, Uos Sapienza, CNR, Rome 00185, Italy
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, United States of America
| |
Collapse
|
26
|
Design and simulation of DNA, RNA and hybrid protein-nucleic acid nanostructures with oxView. Nat Protoc 2022; 17:1762-1788. [PMID: 35668321 DOI: 10.1038/s41596-022-00688-5] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Accepted: 01/18/2022] [Indexed: 11/08/2022]
Abstract
Molecular simulation has become an integral part of the DNA/RNA nanotechnology research pipeline. In particular, understanding the dynamics of structures and single-molecule events has improved the precision of nanoscaffolds and diagnostic tools. Here we present oxView, a design tool for visualization, design, editing and analysis of simulations of DNA, RNA and nucleic acid-protein nanostructures. oxView provides an accessible software platform for designing novel structures, tweaking existing designs, preparing them for simulation in the oxDNA/RNA molecular simulation engine and creating visualizations of simulation results. In several examples, we present procedures for using the tool, including its advanced features that couple the design capabilities with a coarse-grained simulation engine and scripting interface that can programmatically edit structures and facilitate design of complex structures from multiple substructures. These procedures provide a practical basis from which researchers, including experimentalists with limited computational experience, can integrate simulation and 3D visualization into their existing research programs.
Collapse
|
27
|
Afanasyev AY, Onufriev AV. Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach. J Chem Theory Comput 2022; 18:3911-3920. [PMID: 35544776 DOI: 10.1021/acs.jctc.1c01221] [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
We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.
Collapse
Affiliation(s)
- Alexander Y Afanasyev
- Department of Biomedical Engineering and Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|
28
|
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
|
29
|
Guo ZH, Yuan L, Tan YL, Zhang BG, Shi YZ. RNAStat: An Integrated Tool for Statistical Analysis of RNA 3D Structures. FRONTIERS IN BIOINFORMATICS 2022; 1:809082. [PMID: 36303785 PMCID: PMC9580920 DOI: 10.3389/fbinf.2021.809082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 12/17/2021] [Indexed: 11/13/2022] Open
Abstract
The 3D architectures of RNAs are essential for understanding their cellular functions. While an accurate scoring function based on the statistics of known RNA structures is a key component for successful RNA structure prediction or evaluation, there are few tools or web servers that can be directly used to make comprehensive statistical analysis for RNA 3D structures. In this work, we developed RNAStat, an integrated tool for making statistics on RNA 3D structures. For given RNA structures, RNAStat automatically calculates RNA structural properties such as size and shape, and shows their distributions. Based on the RNA structure annotation from DSSR, RNAStat provides statistical information of RNA secondary structure motifs including canonical/non-canonical base pairs, stems, and various loops. In particular, the geometry of base-pairing/stacking can be calculated in RNAStat by constructing a local coordinate system for each base. In addition, RNAStat also supplies the distribution of distance between any atoms to the users to help build distance-based RNA statistical potentials. To test the usability of the tool, we established a non-redundant RNA 3D structure dataset, and based on the dataset, we made a comprehensive statistical analysis on RNA structures, which could have the guiding significance for RNA structure modeling. The python code of RNAStat, the dataset used in this work, and corresponding statistical data files are freely available at GitHub (https://github.com/RNA-folding-lab/RNAStat).
Collapse
Affiliation(s)
- Zhi-Hao Guo
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Li Yuan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- *Correspondence: Ya-Zhou Shi,
| |
Collapse
|
30
|
Poppleton E, Mallya A, Dey S, Joseph J, Šulc P. Nanobase.org: a repository for DNA and RNA nanostructures. Nucleic Acids Res 2022; 50:D246-D252. [PMID: 34747480 PMCID: PMC8728195 DOI: 10.1093/nar/gkab1000] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Revised: 10/07/2021] [Accepted: 10/11/2021] [Indexed: 12/17/2022] Open
Abstract
We introduce a new online database of nucleic acid nanostructures for the field of DNA and RNA nanotechnology. The database implements an upload interface, searching and database browsing. Each deposited nanostructures includes an image of the nanostructure, design file, an optional 3D view, and additional metadata such as experimental data, protocol or literature reference. The database accepts nanostructures in any preferred format used by the uploader for the nanostructure design. We further provide a set of conversion tools that encourage design file conversion into common formats (oxDNA and PDB) that can be used for setting up simulations, interactive editing or 3D visualization. The aim of the repository is to provide to the DNA/RNA nanotechnology community a resource for sharing their designs for further reuse in other systems and also to function as an archive of the designs that have been achieved in the field so far. Nanobase.org is available at https://nanobase.org/.
Collapse
Affiliation(s)
- Erik Poppleton
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Aatmik Mallya
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Swarup Dey
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
- Wyss Institute, Harvard University, Boston, MA 02115, USA
| | - Joel Joseph
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| |
Collapse
|
31
|
Poblete S, Božič A, Kanduč M, Podgornik R, Guzman HV. RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates. ACS OMEGA 2021; 6:32823-32831. [PMID: 34901632 PMCID: PMC8655909 DOI: 10.1021/acsomega.1c04774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 11/04/2021] [Indexed: 06/14/2023]
Abstract
RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto
de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Computational
Biology Lab, Fundación Ciencia &
Vida, Santiago 7780272, Chile
| | - Anže Božič
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Matej Kanduč
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Rudolf Podgornik
- School
of Physical Sciences and Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Institute
of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Wenzhou
Institute of the University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Horacio V. Guzman
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
32
|
Liu H, Hong F, Smith F, Goertz J, Ouldridge T, Stevens MM, Yan H, Šulc P. Kinetics of RNA and RNA:DNA Hybrid Strand Displacement. ACS Synth Biol 2021; 10:3066-3073. [PMID: 34752075 DOI: 10.1021/acssynbio.1c00336] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In nucleic acid nanotechnology, strand displacement is a widely used mechanism where one strand from a hybridized duplex is exchanged with an invading strand that binds to a toehold, a single-stranded region on the duplex. It is used to perform logic operations on a molecular level, initiate cascaded reactions, or even for in vivo diagnostics and treatments. While systematic experimental studies have been carried out to probe the kinetics of strand displacement in DNA with different toehold lengths, sequences, and mismatch positions, there has not been a comparable investigation of RNA or RNA-DNA hybrid systems. Here, we experimentally study how toehold length, toehold location (5' or 3' end of the strand), and mismatches influence the strand displacement kinetics. We observe reaction acceleration with increasing toehold length and placement of the toehold at the 5' end of the substrate. We find that mismatches closer to the interface of toehold and duplex slow down the reaction more than remote mismatches. A comparison of RNA and DNA displacement with hybrid displacement (RNA invading DNA or DNA invading RNA) is partly explainable by the thermodynamic stabilities of the respective toehold regions, but also suggests that the rearrangement from B-form to A-form helix in the case of RNA invading DNA might play a role in the kinetics.
Collapse
Affiliation(s)
- Hao Liu
- Center for Molecular Design and Biomimetics at the Biodesign Institute and School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| | - Fan Hong
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts 02115, United States
| | - Francesca Smith
- Department of Materials, Department of Bioengineering and Institute of Biomedical Engineering, Imperial College London, London, SW7 2AZ, U.K
| | - John Goertz
- Department of Materials, Department of Bioengineering and Institute of Biomedical Engineering, Imperial College London, London, SW7 2AZ, U.K
| | - Thomas Ouldridge
- Department of Bioengineering, Imperial College London, London, SW7 2AZ, U.K
| | - Molly M. Stevens
- Department of Materials, Department of Bioengineering and Institute of Biomedical Engineering, Imperial College London, London, SW7 2AZ, U.K
| | - Hao Yan
- Center for Molecular Design and Biomimetics at the Biodesign Institute and School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| | - Petr Šulc
- Center for Molecular Design and Biomimetics at the Biodesign Institute and School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| |
Collapse
|
33
|
KOH HEEYUEN, LEE JAEGYUNG, LEE JAEYOUNG, KIM RYAN, TABATA OSAMU, JIN-WOO KIM, KIM DONYUN. Design Approaches and Computational Tools for DNA Nanostructures. IEEE OPEN JOURNAL OF NANOTECHNOLOGY 2021; 2:86-100. [PMID: 35756857 PMCID: PMC9232119 DOI: 10.1109/ojnano.2021.3119913] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Designing a structure in nanoscale with desired shape and properties has been enabled by structural DNA nanotechnology. Design strategies in this research field have evolved to interpret various aspects of increasingly more complex nanoscale assembly and to realize molecular-level functionality by exploring static to dynamic characteristics of the target structure. Computational tools have naturally been of significant interest as they are essential to achieve a fine control over both shape and physicochemical properties of the structure. Here, we review the basic design principles of structural DNA nanotechnology together with its computational analysis and design tools.
Collapse
Affiliation(s)
- HEEYUEN KOH
- Institute of Advanced Machines and Design, Seoul National University, Seoul 08826, Republic of Korea
| | - JAE GYUNG LEE
- Department of Mechanical Engineering, Seoul National University, Seoul 08826, Republic of Korea
| | - JAE YOUNG LEE
- Institute of Advanced Machines and Design, Seoul National University, Seoul 08826, Republic of Korea
| | - RYAN KIM
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138 USA
- Bio/Nano Technology Group, Institute for Nanoscience and Engineering, University of Arkansas, Fayetteville, AR 72701 USA
| | - OSAMU TABATA
- Faculty of Engineering, Kyoto University of Advanced Science, Kyoto 621-8555, Japan
| | - KIM JIN-WOO
- Bio/Nano Technology Group, Institute for Nanoscience and Engineering, University of Arkansas, Fayetteville, AR 72701 USA
- Department of Biological and Agricultural Engineering, University of Arkansas, Fayetteville, AR 72701 USA
- Materials Science and Engineering Program, University of Arkansas, Fayetteville, AR 72701 USA
| | - DO-NYUN KIM
- Institute of Advanced Machines and Design, Seoul National University, Seoul 08826, Republic of Korea
- Department of Mechanical Engineering, Seoul National University, Seoul 08826, Republic of Korea
- Institute of Engineering Research, Seoul National University, Seoul 08826, Republic of Korea
| |
Collapse
|
34
|
Mazzanti L, Alferkh L, Frezza E, Pasquali S. Biasing RNA Coarse-Grained Folding Simulations with Small-Angle X-ray Scattering Data. J Chem Theory Comput 2021; 17:6509-6521. [PMID: 34506136 DOI: 10.1021/acs.jctc.1c00441] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.
Collapse
Affiliation(s)
- Liuba Mazzanti
- Laboratoire CiTCoM, CNRS UMR 8038, Université de Paris, 4 Avenue de l'observatoire, 75006 Paris, France
| | - Lina Alferkh
- Laboratoire CiTCoM, CNRS UMR 8038, Université de Paris, 4 Avenue de l'observatoire, 75006 Paris, France
| | - Elisa Frezza
- Laboratoire CiTCoM, CNRS UMR 8038, Université de Paris, 4 Avenue de l'observatoire, 75006 Paris, France
| | - Samuela Pasquali
- Laboratoire CiTCoM, CNRS UMR 8038, Université de Paris, 4 Avenue de l'observatoire, 75006 Paris, France
| |
Collapse
|
35
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|
36
|
Poblete S, Guzman HV. Structural 3D Domain Reconstruction of the RNA Genome from Viruses with Secondary Structure Models. Viruses 2021; 13:1555. [PMID: 34452420 PMCID: PMC8402887 DOI: 10.3390/v13081555] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 07/28/2021] [Accepted: 07/30/2021] [Indexed: 02/07/2023] Open
Abstract
Three-dimensional RNA domain reconstruction is important for the assembly, disassembly and delivery functionalities of a packed proteinaceus capsid. However, to date, the self-association of RNA molecules is still an open problem. Recent chemical probing reports provide, with high reliability, the secondary structure of diverse RNA ensembles, such as those of viral genomes. Here, we present a method for reconstructing the complete 3D structure of RNA genomes, which combines a coarse-grained model with a subdomain composition scheme to obtain the entire genome inside proteinaceus capsids based on secondary structures from experimental techniques. Despite the amount of sampling involved in the folded and also unfolded RNA molecules, advanced microscope techniques can provide points of anchoring, which enhance our model to include interactions between capsid pentamers and RNA subdomains. To test our method, we tackle the satellite tobacco mosaic virus (STMV) genome, which has been widely studied by both experimental and computational communities. We provide not only a methodology to structurally analyze the tertiary conformations of the RNA genome inside capsids, but a flexible platform that allows the easy implementation of features/descriptors coming from both theoretical and experimental approaches.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Chile and Computational Biology Lab, Fundación Ciencia & Vida, Santiago 7780272, Chile
| | - Horacio V. Guzman
- Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
37
|
Poppleton E, Romero R, Mallya A, Rovigatti L, Šulc P. OxDNA.org: a public webserver for coarse-grained simulations of DNA and RNA nanostructures. Nucleic Acids Res 2021; 49:W491-W498. [PMID: 34009383 PMCID: PMC8265093 DOI: 10.1093/nar/gkab324] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2021] [Revised: 04/09/2021] [Accepted: 04/19/2021] [Indexed: 12/15/2022] Open
Abstract
OxDNA and oxRNA are popular coarse-grained models used by the DNA/RNA nanotechnology community to prototype, analyze and rationalize designed DNA and RNA nanostructures. Here, we present oxDNA.org, a graphical web interface for running, visualizing and analyzing oxDNA and oxRNA molecular dynamics simulations on a GPU-enabled high performance computing server. OxDNA.org automatically generates simulation files, including a multi-step relaxation protocol for structures exported in non-physical states from DNA/RNA design tools. Once the simulation is complete, oxDNA.org provides an interactive visualization and analysis interface using the browser-based visualizer oxView to facilitate the understanding of simulation results for a user's specific structure. This online tool significantly lowers the entry barrier of integrating simulations in the nanostructure design pipeline for users who are not experts in the technical aspects of molecular simulation. The webserver is freely available at oxdna.org.
Collapse
Affiliation(s)
- Erik Poppleton
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Roger Romero
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Aatmik Mallya
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Lorenzo Rovigatti
- Department of Physics, Sapienza Università di Roma, P.le A. Moro 2 00185, Rome, Italy
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| |
Collapse
|
38
|
Li J, Chen SJ. RNA 3D Structure Prediction Using Coarse-Grained Models. Front Mol Biosci 2021; 8:720937. [PMID: 34277713 PMCID: PMC8283274 DOI: 10.3389/fmolb.2021.720937] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2021] [Accepted: 06/17/2021] [Indexed: 12/12/2022] Open
Abstract
The three-dimensional (3D) structures of Ribonucleic acid (RNA) molecules are essential to understanding their various and important biological functions. However, experimental determination of the atomic structures is laborious and technically difficult. The large gap between the number of sequences and the experimentally determined structures enables the thriving development of computational approaches to modeling RNAs. However, computational methods based on all-atom simulations are intractable for large RNA systems, which demand long time simulations. Facing such a challenge, many coarse-grained (CG) models have been developed. Here, we provide a review of CG models for modeling RNA 3D structures, compare the performance of the different models, and offer insights into potential future developments.
Collapse
Affiliation(s)
| | - Shi-Jie Chen
- Departments of Physics and Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, MO, United States
| |
Collapse
|
39
|
Sengar A, Ouldridge TE, Henrich O, Rovigatti L, Šulc P. A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results. Front Mol Biosci 2021; 8:693710. [PMID: 34235181 PMCID: PMC8256390 DOI: 10.3389/fmolb.2021.693710] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Accepted: 05/27/2021] [Indexed: 11/13/2022] Open
Abstract
The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.
Collapse
Affiliation(s)
- A. Sengar
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - T. E. Ouldridge
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - O. Henrich
- Department of Physics, SUPA, University of Strathclyde, Glasgow, United Kingdom
| | - L. Rovigatti
- Department of Physics, Sapienza University of Rome, Rome, Italy
- CNR Institute of Complex Systems, Sapienza University of Rome, Rome, Italy
| | - P. Šulc
- Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, United States
- School of Molecular Sciences, Arizona State University, Tempe, AZ, United States
| |
Collapse
|
40
|
RNA origami design tools enable cotranscriptional folding of kilobase-sized nanoscaffolds. Nat Chem 2021; 13:549-558. [PMID: 33972754 PMCID: PMC7610888 DOI: 10.1038/s41557-021-00679-1] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2019] [Accepted: 03/08/2021] [Indexed: 12/18/2022]
Abstract
RNA origami is a framework for the modular design of nanoscaffolds that can be folded from a single strand of RNA, and used to organize molecular components with nanoscale precision. Design of genetically expressible RNA origami, which must cotranscriptionally fold, requires modeling and design tools that simultaneously consider thermodynamics, folding pathway, sequence constraints, and pseudoknot optimization. Here, we describe RNA Origami Automated Design software (ROAD), which builds origami models from a library of structural modules, identifies potential folding barriers, and designs optimized sequences. Using ROAD, we extend the scale and functional diversity of RNA scaffolds, creating 32 designs of up to 2360 nucleotides, five that scaffold two proteins, and seven that scaffold two small molecules at precise distances. Micrographic and chromatographic comparison of optimized and nonoptimized structures validates that our principles for strand routing and sequence design substantially improve yield. By providing efficient design of RNA origami, ROAD may simplify construction of custom RNA scaffolds for nanomedicine and synthetic biology.
Collapse
|
41
|
Feng C, Tan YL, Cheng YX, Shi YZ, Tan ZJ. Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement. Front Mol Biosci 2021; 8:666369. [PMID: 33928126 PMCID: PMC8078894 DOI: 10.3389/fmolb.2021.666369] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Accepted: 03/17/2021] [Indexed: 12/27/2022] Open
Abstract
Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.
Collapse
Affiliation(s)
- Chenjie Feng
- Key Laboratory of Artificial Micro and Nano-structures of Ministry of Education, Center for Theoretical Physics, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Key Laboratory of Artificial Micro and Nano-structures of Ministry of Education, Center for Theoretical Physics, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Yu-Xuan Cheng
- Key Laboratory of Artificial Micro and Nano-structures of Ministry of Education, Center for Theoretical Physics, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematics and Computer Science, Wuhan Textile University, Wuhan, China
| | - Zhi-Jie Tan
- Key Laboratory of Artificial Micro and Nano-structures of Ministry of Education, Center for Theoretical Physics, School of Physics and Technology, Wuhan University, Wuhan, China
| |
Collapse
|
42
|
Procyk J, Poppleton E, Šulc P. Coarse-grained nucleic acid-protein model for hybrid nanotechnology. SOFT MATTER 2021; 17:3586-3593. [PMID: 33398312 DOI: 10.1039/d0sm01639j] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The emerging field of hybrid DNA-protein nanotechnology brings with it the potential for many novel materials which combine the addressability of DNA nanotechnology with the versatility of protein interactions. However, the design and computational study of these hybrid structures is difficult due to the system sizes involved. To aid in the design and in silico analysis process, we introduce here a coarse-grained DNA/RNA-protein model that extends the oxDNA/oxRNA models of DNA/RNA with a coarse-grained model of proteins based on an anisotropic network model representation. Fully equipped with analysis scripts and visualization, our model aims to facilitate hybrid nanomaterial design towards eventual experimental realization, as well as enabling study of biological complexes. We further demonstrate its usage by simulating DNA-protein nanocage, DNA wrapped around histones, and a nascent RNA in polymerase.
Collapse
Affiliation(s)
- Jonah Procyk
- 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.
| | | | | |
Collapse
|
43
|
Zhang D, Li J, Chen SJ. IsRNA1: De Novo Prediction and Blind Screening of RNA 3D Structures. J Chem Theory Comput 2021; 17:1842-1857. [PMID: 33560836 DOI: 10.1021/acs.jctc.0c01148] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Modeling structures and functions of large ribonucleic acid (RNAs) especially with complicated topologies is highly challenging due to the inefficiency of large conformational sampling and the presence of complicated tertiary interactions. To address this problem, one highly promising approach is coarse-grained modeling. Here, following an iterative simulated reference state approach to decipher the correlations between different structural parameters, we developed a potent coarse-grained RNA model named as IsRNA1 for RNA studies. Molecular dynamics simulations in the IsRNA1 can predict the native structures of small RNAs from a sequence and fold medium-sized RNAs into near-native tertiary structures with the assistance of secondary structure constraints. A large-scale benchmark test on RNA 3D structure prediction shows that IsRNA1 exhibits improved performance for relatively large RNAs of complicated topologies, such as large stem-loop structures and structures containing long-range tertiary interactions. The advantages of IsRNA1 include the consideration of the correlations between the different structural variables, the appropriate characterization of canonical base-pairing and base-stacking interactions, and the better sampling for the backbone conformations. Moreover, a blind screening protocol was developed based on IsRNA1 to identify good structural models from a pool of candidates without prior knowledge of the native structures.
Collapse
Affiliation(s)
- Dong Zhang
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| | - Jun Li
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| |
Collapse
|
44
|
Ferreira I, Amarante TD, Weber G. Salt dependent mesoscopic model for RNA at multiple strand concentrations. Biophys Chem 2021; 271:106551. [PMID: 33662903 DOI: 10.1016/j.bpc.2021.106551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Revised: 01/19/2021] [Accepted: 01/19/2021] [Indexed: 12/12/2022]
Abstract
Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.
Collapse
Affiliation(s)
- Izabela Ferreira
- Departamento de Física, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil; Programa Interunidades de Pós-Graduação em Bioinformática, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
| | - Tauanne D Amarante
- MRC Cancer Unit, University of Cambridge, Hutchison/MRC Research Centre, Cambridge Biomedical Campus, Cambridge, UK
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil.
| |
Collapse
|
45
|
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
|
46
|
Regy RM, Dignon GL, Zheng W, Kim YC, Mittal J. Sequence dependent phase separation of protein-polynucleotide mixtures elucidated using molecular simulations. Nucleic Acids Res 2020; 48:12593-12603. [PMID: 33264400 PMCID: PMC7736803 DOI: 10.1093/nar/gkaa1099] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 10/21/2020] [Accepted: 10/27/2020] [Indexed: 12/22/2022] Open
Abstract
Ribonucleoprotein (RNP) granules are membraneless organelles (MLOs), which majorly consist of RNA and RNA-binding proteins and are formed via liquid-liquid phase separation (LLPS). Experimental studies investigating the drivers of LLPS have shown that intrinsically disordered proteins (IDPs) and nucleic acids like RNA and other polynucleotides play a key role in modulating protein phase separation. There is currently a dearth of modelling techniques which allow one to delve deeper into how polynucleotides play the role of a modulator/promoter of LLPS in cells using computational methods. Here, we present a coarse-grained polynucleotide model developed to fill this gap, which together with our recently developed HPS model for protein LLPS, allows us to capture the factors driving protein-polynucleotide phase separation. We explore the capabilities of the modelling framework with the LAF-1 RGG system which has been well studied in experiments and also with the HPS model previously. Further taking advantage of the fact that the HPS model maintains sequence specificity we explore the role of charge patterning on controlling polynucleotide incorporation into condensates. With increased charge patterning we observe formation of structured or patterned condensates which suggests the possible roles of polynucleotides in not only shifting the phase boundaries but also introducing microscopic organization in MLOs.
Collapse
Affiliation(s)
- Roshan Mammen Regy
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, PA 18015, USA
| | - Gregory L Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, PA 18015, USA
| | - Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, AZ 85212, USA
| | - Young C Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, DC 20375, USA
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, PA 18015, USA
| |
Collapse
|
47
|
Halder A, Kumar S, Valsson O, Reddy G. Mg 2+ Sensing by an RNA Fragment: Role of Mg 2+-Coordinated Water Molecules. J Chem Theory Comput 2020; 16:6702-6715. [PMID: 32941038 DOI: 10.1021/acs.jctc.0c00589] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg2+, Ca2+, and K+) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg2+ solution. We show that the rigid network of Mg2+-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K+ and Ca2+ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.
Collapse
Affiliation(s)
- Antarip Halder
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Omar Valsson
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| |
Collapse
|
48
|
Badu S, Melnik R, Singh S. Mathematical and computational models of RNA nanoclusters and their applications in data-driven environments. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1804564] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Shyam Badu
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
| | - Roderick Melnik
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
- BCAM-Basque Center for Applied Mathematics, Bilbao, Spain
| | - Sundeep Singh
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
| |
Collapse
|
49
|
Torelli E, Kozyra J, Shirt-Ediss B, Piantanida L, Voïtchovsky K, Krasnogor N. Cotranscriptional Folding of a Bio-orthogonal Fluorescent Scaffolded RNA Origami. ACS Synth Biol 2020; 9:1682-1692. [PMID: 32470289 DOI: 10.1021/acssynbio.0c00009] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The scaffolded origami technique is an attractive tool for engineering nucleic acid nanostructures. This paper demonstrates scaffolded RNA origami folding in vitro in which, for the first time, all components are transcribed simultaneously in a single-pot reaction. Double-stranded DNA sequences are transcribed by T7 RNA polymerase into scaffold and staple strands able to correctly fold in a high synthesis yield into the nanoribbon. Synthesis is successfully confirmed by atomic force microscopy, and the unpurified transcription reaction mixture is analyzed by an in gel-imaging assay where the transcribed RNA nanoribbons are able to capture the specific dye through the reconstituted split Broccoli aptamer showing a clear green fluorescent band. Finally, we simulate the RNA origami in silico using the nucleotide-level coarse-grained model oxRNA to investigate the thermodynamic stability of the assembled nanostructure in isothermal conditions over a period of time. Our work suggests that the scaffolded origami technique is a viable, and potentially more powerful, assembly alternative to the single-stranded origami technique for future in vivo applications.
Collapse
Affiliation(s)
- Emanuela Torelli
- Interdisciplinary Computing and Complex BioSystems (ICOS), Centre for Synthetic Biology and Bioeconomy (CSBB), Devonshire Building, Newcastle University, Newcastle upon Tyne NE1 7RX, United Kingdom
| | - Jerzy Kozyra
- Interdisciplinary Computing and Complex BioSystems (ICOS), Centre for Synthetic Biology and Bioeconomy (CSBB), Devonshire Building, Newcastle University, Newcastle upon Tyne NE1 7RX, United Kingdom
| | - Ben Shirt-Ediss
- Interdisciplinary Computing and Complex BioSystems (ICOS), Centre for Synthetic Biology and Bioeconomy (CSBB), Devonshire Building, Newcastle University, Newcastle upon Tyne NE1 7RX, United Kingdom
| | - Luca Piantanida
- Department of Physics, Durham University, Durham DH1 3LE, United Kingdom
| | - Kislon Voïtchovsky
- Department of Physics, Durham University, Durham DH1 3LE, United Kingdom
| | - Natalio Krasnogor
- Interdisciplinary Computing and Complex BioSystems (ICOS), Centre for Synthetic Biology and Bioeconomy (CSBB), Devonshire Building, Newcastle University, Newcastle upon Tyne NE1 7RX, United Kingdom
| |
Collapse
|
50
|
Poppleton E, Bohlin J, Matthies M, Sharma S, Zhang F, Šulc P. Design, optimization and analysis of large DNA and RNA nanostructures through interactive visualization, editing and molecular simulation. Nucleic Acids Res 2020; 48:e72. [PMID: 32449920 PMCID: PMC7337935 DOI: 10.1093/nar/gkaa417] [Citation(s) in RCA: 73] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 04/22/2020] [Accepted: 05/07/2020] [Indexed: 12/27/2022] Open
Abstract
This work seeks to remedy two deficiencies in the current nucleic acid nanotechnology software environment: the lack of both a fast and user-friendly visualization tool and a standard for structural analyses of simulated systems. We introduce here oxView, a web browser-based visualizer that can load structures with over 1 million nucleotides, create videos from simulation trajectories, and allow users to perform basic edits to DNA and RNA designs. We additionally introduce open-source software tools for extracting common structural parameters to characterize large DNA/RNA nanostructures simulated using the coarse-grained modeling tool, oxDNA, which has grown in popularity in recent years and is frequently used to prototype new nucleic acid nanostructural designs, model biophysics of DNA/RNA processes, and rationalize experimental results. The newly introduced software tools facilitate the computational characterization of DNA/RNA designs by providing multiple analysis scripts, including mean structures and structure flexibility characterization, hydrogen bond fraying, and interduplex angles. The output of these tools can be loaded into oxView, allowing users to interact with the simulated structure in a 3D graphical environment and modify the structures to achieve the required properties. We demonstrate these newly developed tools by applying them to design and analysis of a range of DNA/RNA nanostructures.
Collapse
Affiliation(s)
- Erik Poppleton
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Joakim Bohlin
- Department of Physics, Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, UK
| | - Michael Matthies
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Shuchi Sharma
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| | - Fei Zhang
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
- Department of Chemistry, Rutgers University-Newark, 73 Warren St, Newark, NJ 07102, USA
| | - Petr Šulc
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, AZ 85281, USA
| |
Collapse
|