1
|
Menon AP, Villanueva H, Meraviglia-Crivelli D, van Santen HM, Hellmeier J, Zheleva A, Nonateli F, Peters T, Wachsmann TL, Hernandez-Rueda M, Huppa JB, Schütz GJ, Sevcsik E, Moreno B, Pastor F. CD3 aptamers promote expansion and persistence of tumor-reactive T cells for adoptive T cell therapy in cancer. MOLECULAR THERAPY. NUCLEIC ACIDS 2024; 35:102198. [PMID: 38745854 PMCID: PMC11091522 DOI: 10.1016/j.omtn.2024.102198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 04/21/2024] [Indexed: 05/16/2024]
Abstract
The CD3/T cell receptor (TCR) complex is responsible for antigen-specific pathogen recognition by T cells, and initiates the signaling cascade necessary for activation of effector functions. CD3 agonistic antibodies are commonly used to expand T lymphocytes in a wide range of clinical applications, including in adoptive T cell therapy for cancer patients. A major drawback of expanding T cell populations ex vivo using CD3 agonistic antibodies is that they expand and activate T cells independent of their TCR antigen specificity. Therapeutic agents that facilitate expansion of T cells in an antigen-specific manner and reduce their threshold of T cell activation are therefore of great interest for adoptive T cell therapy protocols. To identify CD3-specific T cell agonists, several RNA aptamers were selected against CD3 using Systematic Evolution of Ligands by EXponential enrichment combined with high-throughput sequencing. The extent and specificity of aptamer binding to target CD3 were assessed through surface plasma resonance, P32 double-filter assays, and flow cytometry. Aptamer-mediated modulation of the threshold of T cell activation was observed in vitro and in preclinical transgenic TCR mouse models. The aptamers improved efficacy and persistence of adoptive T cell therapy by low-affinity TCR-reactive T lymphocytes in melanoma-bearing mice. Thus, CD3-specific aptamers can be applied as therapeutic agents which facilitate the expansion of tumor-reactive T lymphocytes while conserving their tumor specificity. Furthermore, selected CD3 aptamers also exhibit cross-reactivity to human CD3, expanding their potential for clinical translation and application in the future.
Collapse
Affiliation(s)
- Ashwathi Puravankara Menon
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Helena Villanueva
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Daniel Meraviglia-Crivelli
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Hisse M. van Santen
- Immune System Development and Function Unit, Centro Biología Molecular Severo Ochoa, Consejo Superior de Investigaciones Científicas/Universidad Autónoma de Madrid (CSIC/UAM), 28049 Madrid, Spain
| | - Joschka Hellmeier
- Institute of Applied Physics, TU Wien, Lehargasse 6, 1060 Vienna, Austria
| | - Angelina Zheleva
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Francesca Nonateli
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Timo Peters
- Institute for Hygiene and Applied Immunology, Center for Pathophysiology, Infectiology, Immunology, Medical University of Vienna, 1090 Vienna, Austria
| | | | - Mercedes Hernandez-Rueda
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Johannes B. Huppa
- Institute for Hygiene and Applied Immunology, Center for Pathophysiology, Infectiology, Immunology, Medical University of Vienna, 1090 Vienna, Austria
| | - Gerhard J. Schütz
- Institute of Applied Physics, TU Wien, Lehargasse 6, 1060 Vienna, Austria
| | - Eva Sevcsik
- Institute of Applied Physics, TU Wien, Lehargasse 6, 1060 Vienna, Austria
| | - Beatriz Moreno
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
| | - Fernando Pastor
- Molecular Therapeutics Program, Center for Applied Medical Research, CIMA, University of Navarra, 31008 Pamplona, Spain
- Navarra Institute for Health Research (IdiSNA), 31008 Pamplona, Spain
- Spanish Center for Biomedical Research Network in Oncology (CIBERONC), 28029 Madrid, Spain
| |
Collapse
|
2
|
Ota M, Takahashi H, Nogi Y, Kagotani Y, Saito-Tarashima N, Kondo J, Minakawa N. Synthesis and properties of fully-modified 4'-selenoRNA, an endonuclease-resistant RNA analog. Bioorg Med Chem 2022; 76:117093. [PMID: 36434923 DOI: 10.1016/j.bmc.2022.117093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Revised: 11/02/2022] [Accepted: 11/05/2022] [Indexed: 11/21/2022]
Abstract
A large number of chemically modified oligonucleotides (ONs) have been developed for RNA-based technologies. In most modified RNAs, the characteristic 2'-hydroxyl (2'-OH) groups are removed to enhance both nuclease resistance and hybridization ability. However, the importance of the 2'-OH group in RNA structure and function is well known. Here, we report the synthesis and properties of 4'-selenoRNA in which all four nucleoside units retain the 2'-OH groups but contain a selenium atom instead of an oxygen atom at the 4'-position of the furanose ring. 4'-SelenoRNA has enhanced ability to form duplexes with RNA, and high endonuclease resistance despite the presence of the 2'-OH groups. X-ray crystallography analysis showed that the 4'-selenoRNA duplex adopts an A-conformation, similar to natural RNA, although one 4'-selenocytidine residue has unusual South-type sugar puckering. Furthermore, preliminary studies using 4'-seleno-modified siRNAs suggest that 4'-selenoRNA may be applicable to RNA interference technology. Collectively, our results raise the possibility of a new class of modified RNA in which 2'-OH groups do not need to be masked.
Collapse
Affiliation(s)
- Masashi Ota
- Graduate School of Pharmaceutical Science, Tokushima University, Shomachi 1-78-1, Tokushima 770-8505, Japan
| | - Hiromi Takahashi
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, 7-1 Kioi-cho, Chiyoda-ku, 102-8554 Tokyo, Japan
| | - Yuhei Nogi
- Graduate School of Pharmaceutical Science, Tokushima University, Shomachi 1-78-1, Tokushima 770-8505, Japan
| | - Yuma Kagotani
- Graduate School of Pharmaceutical Science, Tokushima University, Shomachi 1-78-1, Tokushima 770-8505, Japan
| | - Noriko Saito-Tarashima
- Graduate School of Pharmaceutical Science, Tokushima University, Shomachi 1-78-1, Tokushima 770-8505, Japan
| | - Jiro Kondo
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, 7-1 Kioi-cho, Chiyoda-ku, 102-8554 Tokyo, Japan
| | - Noriaki Minakawa
- Graduate School of Pharmaceutical Science, Tokushima University, Shomachi 1-78-1, Tokushima 770-8505, Japan.
| |
Collapse
|
3
|
Zhang Y, Yan M, Huang T, Wang X. Understanding the Structural Elasticity of RNA and DNA: All‐Atom Molecular Dynamics. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Yingtong Zhang
- Department of Physics Wenzhou University Wenzhou 325035 China
| | - Miao Yan
- Department of Physics Wenzhou University Wenzhou 325035 China
| | - Tingting Huang
- Department of Mechanical Engineering Shanghai Techanical Institute of Electronics and Information Shanghai 201411 China
| | - Xianghong Wang
- Department of Physics Wenzhou University Wenzhou 325035 China
- Department of Mechanical Engineering Shanghai Techanical Institute of Electronics and Information Shanghai 201411 China
| |
Collapse
|
4
|
Harp JM, Lybrand TP, Pallan PS, Coates L, Sullivan B, Egli M. Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure. Nucleic Acids Res 2022; 50:7721-7738. [PMID: 35819202 PMCID: PMC9303348 DOI: 10.1093/nar/gkac577] [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: 03/21/2022] [Revised: 06/15/2022] [Accepted: 06/21/2022] [Indexed: 11/14/2022] Open
Abstract
The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.
Collapse
Affiliation(s)
- Joel M Harp
- Department of Biochemistry, Vanderbilt University, Nashville, TN 37232, USA.,Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Terry P Lybrand
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA.,Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
| | - Pradeep S Pallan
- Department of Biochemistry, Vanderbilt University, Nashville, TN 37232, USA.,Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Leighton Coates
- Neutron Scattering Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, TN 37831, USA
| | - Brendan Sullivan
- Neutron Scattering Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, TN 37831, USA
| | - Martin Egli
- Department of Biochemistry, Vanderbilt University, Nashville, TN 37232, USA.,Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| |
Collapse
|
5
|
Intermolecular interactions in microhydrated ribonucleoside and deoxyribonucleoside: A computational study. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113422] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
6
|
Liczner C, Duke K, Juneau G, Egli M, Wilds CJ. Beyond ribose and phosphate: Selected nucleic acid modifications for structure-function investigations and therapeutic applications. Beilstein J Org Chem 2021; 17:908-931. [PMID: 33981365 PMCID: PMC8093555 DOI: 10.3762/bjoc.17.76] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Accepted: 04/14/2021] [Indexed: 12/16/2022] Open
Abstract
Over the past 25 years, the acceleration of achievements in the development of oligonucleotide-based therapeutics has resulted in numerous new drugs making it to the market for the treatment of various diseases. Oligonucleotides with alterations to their scaffold, prepared with modified nucleosides and solid-phase synthesis, have yielded molecules with interesting biophysical properties that bind to their targets and are tolerated by the cellular machinery to elicit a therapeutic outcome. Structural techniques, such as crystallography, have provided insights to rationalize numerous properties including binding affinity, nuclease stability, and trends observed in the gene silencing. In this review, we discuss the chemistry, biophysical, and structural properties of a number of chemically modified oligonucleotides that have been explored for gene silencing.
Collapse
Affiliation(s)
- Christopher Liczner
- Department of Chemistry and Biochemistry, Concordia University, Montréal, Québec H4B 1R6, Canada
| | - Kieran Duke
- Department of Chemistry and Biochemistry, Concordia University, Montréal, Québec H4B 1R6, Canada
| | - Gabrielle Juneau
- Department of Chemistry and Biochemistry, Concordia University, Montréal, Québec H4B 1R6, Canada
| | - Martin Egli
- Department of Biochemistry, Vanderbilt Institute of Chemical Biology, and Center for Structural Biology, School of Medicine, Vanderbilt University, Nashville, Tennessee 37232, United States
| | - Christopher J Wilds
- Department of Chemistry and Biochemistry, Concordia University, Montréal, Québec H4B 1R6, Canada
| |
Collapse
|
7
|
Pallan PS, Lybrand TP, Schlegel MK, Harp JM, Jahns H, Manoharan M, Egli M. Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost. Biochemistry 2020; 59:4627-4637. [PMID: 33275419 DOI: 10.1021/acs.biochem.0c00685] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature Tm and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.
Collapse
Affiliation(s)
| | | | - Mark K Schlegel
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | | | - Hartmut Jahns
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | - Muthiah Manoharan
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | | |
Collapse
|
8
|
Wang Z, Liu J, Zhang Y, Qi J, Han X, Zhao X, Bai D, Zhao H, Chen Q. Intrinsic Contributions of 2'-Hydroxyl to the Hydration of Nucleosides at the Monomeric Level. Chemistry 2020; 26:17046-17055. [PMID: 32786015 DOI: 10.1002/chem.202002835] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Indexed: 02/05/2023]
Abstract
Although many reports have revealed structural differences between DNA and RNA at the polymeric level, there are no comparative studies with 2'-deoxyribonucleoside and ribonucleoside to explore the role of the 2'-OH group at the monomeric level under the same conditions. Inspired by this, herein, the intrinsic contributions of the 2'-OH group in the nucleoside have been systematically investigated by directly solving the single-crystal structures of 2'-deoxy-2-aminoadenosine (1), 2-aminoadenosine (2), and 2-aminoarabinofuranosyladenine (3) in water. The 2'-OH group not only influenced the conformation and base-pair pattern of the single-nucleoside molecule, but also played a fundamental role in the entire supramolecular structure. Interestingly, compound 1, which did not contain the 2'-OH group, displayed strong hydration, whereas 2 and 3 (with the 2'-OH group in the opposite direction) exhibited no hydration, which was completely different from that observed in nucleic acids. Meanwhile, compound 1 trapped water molecules to form unique trihydrol moieties, which further served as the backbone to construct the simplest double-chain DNA-like structures. To this end, to investigate the effect of the biological environment on these unique structures, the solvent was changed from water to phosphate-buffered saline (PBS). Surprisingly, such a subtle adjustment led to entirely different superstructures, consisting of 2D lamellar structures in water and 3D porous structures in PBS. These large morphological differences could be attributed to delicate ion hydration, which was also confirmed through variable-temperature X-ray analysis, SEM, and intermolecular interaction energy calculations. In summary, this study comprehensively investigated the intrinsic contributions of 2'-hydroxyl to the hydration of nucleosides at the monomeric level; this is helpful to further understand the differences in DNA/RNA and the impact of their surrounding environment.
Collapse
Affiliation(s)
- Zheng Wang
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Jiang Liu
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Yanan Zhang
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Jiajia Qi
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Xianglong Han
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Xuefeng Zhao
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Ding Bai
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Hang Zhao
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| | - Qianming Chen
- State Key Laboratory of Oral Diseases, National Clinical, Research Center for Oral Diseases, Chinese Academy of Medical Sciences, Research Unit of Oral Carcinogenesis and Management, West China, Hospital of Stomatology, Sichuan University, 610041, Chengdu, P.R. China
| |
Collapse
|
9
|
Salsbury AM, Lemkul JA. Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics. Curr Opin Struct Biol 2020; 67:9-17. [PMID: 32950937 DOI: 10.1016/j.sbi.2020.08.003] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 08/04/2020] [Accepted: 08/13/2020] [Indexed: 01/24/2023]
Abstract
Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States; Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, United States.
| |
Collapse
|
10
|
Palomino‐Hernandez O, Margreiter MA, Rossetti G. Challenges in RNA Regulation in Huntington's Disease: Insights from Computational Studies. Isr J Chem 2020. [DOI: 10.1002/ijch.202000021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Oscar Palomino‐Hernandez
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Faculty 1RWTH Aachen 52425 Aachen Germany
- Computation-based Science and Technology Research CenterThe Cyprus Institute Nicosia 2121 Cyprus
- Institute of Life ScienceThe Hebrew University of Jerusalem Jerusalem 91904 Israel
| | - Michael A. Margreiter
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Faculty 1RWTH Aachen 52425 Aachen Germany
| | - Giulia Rossetti
- Computational Biomedicine, Institute of Neuroscience and Medicine (INM-9)/Instute for advanced simulations (IAS-5)Forschungszentrum Juelich 52425 Jülich Germany
- Jülich Supercomputing Centre (JSC)Forschungszentrum Jülich 52425 Jülich Germany
- Department of Hematology, Oncology, Hemostaseology and Stem Cell Transplantation University Hospital AachenRWTH Aachen University Pauwelsstraße 30 52074 Aachen Germany
| |
Collapse
|
11
|
Light-controlled twister ribozyme with single-molecule detection resolves RNA function in time and space. Proc Natl Acad Sci U S A 2020; 117:12080-12086. [PMID: 32430319 DOI: 10.1073/pnas.2003425117] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Small ribozymes such as Oryza sativa twister spontaneously cleave their own RNA when the ribozyme folds into its active conformation. The coupling between twister folding and self-cleavage has been difficult to study, however, because the active ribozyme rapidly converts to product. Here, we describe the synthesis of a photocaged nucleotide that releases guanosine within microseconds upon photosolvolysis with blue light. Application of this tool to O. sativa twister achieved the spatial (75 µm) and temporal (≤30 ms) control required to resolve folding and self-cleavage events when combined with single-molecule fluorescence detection of the ribozyme folding pathway. Real-time observation of single ribozymes after photo-deprotection showed that the precleaved folded state is unstable and quickly unfolds if the RNA does not react. Kinetic analysis showed that Mg2+ and Mn2+ ions increase ribozyme efficiency by making transitions to the high energy active conformation more probable, rather than by stabilizing the folded ground state or the cleaved product. This tool for light-controlled single RNA folding should offer precise and rapid control of other nucleic acid systems.
Collapse
|
12
|
Lemkul JA. Same fold, different properties: polarizable molecular dynamics simulations of telomeric and TERRA G-quadruplexes. Nucleic Acids Res 2020; 48:561-575. [PMID: 31807754 PMCID: PMC6954416 DOI: 10.1093/nar/gkz1154] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 11/21/2019] [Accepted: 11/26/2019] [Indexed: 12/12/2022] Open
Abstract
DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, USA
| |
Collapse
|
13
|
Zhao J, Kennedy SD, Berger KD, Turner DH. Nuclear Magnetic Resonance of Single-Stranded RNAs and DNAs of CAAU and UCAAUC as Benchmarks for Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:1968-1984. [PMID: 31904966 DOI: 10.1021/acs.jctc.9b00912] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
RNA and DNA are rapidly emerging as targets for therapeutics and as potential frameworks for nanotechnology. Accurate methods for predicting and designing structures and dynamics of nucleic acids would accelerate progress in these and other applications. Suitable approximations for modeling nucleic acids are being developed but require validation against disparate experimental observations. Here, nuclear magnetic resonance spectra for RNA and DNA single strands, CAAU and UCAAUC, are used as benchmarks to test molecular dynamics simulations with AMBER force fields OL3 and ROC-RNA for RNA and BSC1 for DNA. A detailed scheme for making comparisons is also presented. The results reflect recent progress in approximations and reveal remaining challenges.
Collapse
Affiliation(s)
- Jianbo Zhao
- Department of Chemistry, University of Rochester, Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| | - Scott D Kennedy
- Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States.,Department of Biochemistry and Biophysics, School of Medicine & Dentistry, University of Rochester, Rochester, New York 14642, United States
| | - Kyle D Berger
- Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States.,Department of Biochemistry and Biophysics, School of Medicine & Dentistry, University of Rochester, Rochester, New York 14642, United States
| | - Douglas H Turner
- Department of Chemistry, University of Rochester, Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
14
|
Li Y, Bedi RK, Wiedmer L, Huang D, Śledź P, Caflisch A. Flexible Binding of m 6A Reader Protein YTHDC1 to Its Preferred RNA Motif. J Chem Theory Comput 2019; 15:7004-7014. [PMID: 31670957 DOI: 10.1021/acs.jctc.9b00987] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
N6-Methyladenosine (m6A) is the most prevalent chemical modification in human mRNAs. Its recognition by reader proteins enables many cellular functions, including splicing and translation of mRNAs. However, the binding mechanisms of m6A-containing RNAs to their readers are still elusive due to the unclear roles of m6A-flanking ribonucleotides. Here, we use a model system, YTHDC1 with its RNA motif 5'-G-2G-1(m6A)C+1U+2-3', to investigate the binding mechanisms by atomistic simulations, X-ray crystallography, and isothermal titration calorimetry. The experimental data and simulation results show that m6A is captured by an aromatic cage of YTHDC1 and the 3' terminus nucleotides are stabilized by cation-π-π interactions, while the 5' terminus remains flexible. Notably, simulations of unbound RNA motifs reveal that the methyl group of m6A and the 5' terminus shift the conformational preferences of the oligoribonucleotide to the bound-like conformation, thereby facilitating the association process. The binding mechanisms may help in the discovery of chemical probes against m6A reader proteins.
Collapse
Affiliation(s)
- Yaozong Li
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | - Rajiv Kumar Bedi
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | - Lars Wiedmer
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | - Danzhi Huang
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | - Paweł Śledź
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | - Amedeo Caflisch
- Department of Biochemistry , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| |
Collapse
|
15
|
Lemkul JA, MacKerell AD. Polarizable force field for RNA based on the classical drude oscillator. J Comput Chem 2018; 39:2624-2646. [PMID: 30515902 PMCID: PMC6284239 DOI: 10.1002/jcc.25709] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Revised: 08/01/2018] [Accepted: 09/23/2018] [Indexed: 12/15/2022]
Abstract
RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
16
|
Bao L, Zhang X, Shi YZ, Wu YY, Tan ZJ. Understanding the Relative Flexibility of RNA and DNA Duplexes: Stretching and Twist-Stretch Coupling. Biophys J 2017; 112:1094-1104. [PMID: 28355538 DOI: 10.1016/j.bpj.2017.02.022] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 02/02/2017] [Accepted: 02/21/2017] [Indexed: 01/16/2023] Open
Abstract
The flexibility of double-stranded (ds) RNA and dsDNA is crucial for their biological functions. Recent experiments have shown that the flexibility of dsRNA and dsDNA can be distinctively different in the aspects of stretching and twist-stretch coupling. Although various studies have been performed to understand the flexibility of dsRNA and dsDNA, there is still a lack of deep understanding of the distinctive differences in the flexibility of dsRNA and dsDNA helices as pertains to their stretching and twist-stretch coupling. In this work, we have explored the relative flexibility in stretching and twist-stretch coupling between dsRNA and dsDNA by all-atom molecular dynamics simulations. The calculated stretch modulus and twist-stretch coupling are in good accordance with the existing experiments. Our analyses show that the differences in stretching and twist-stretch coupling between dsRNA and dsDNA helices are mainly attributed to their different (A- and B-form) helical structures. Stronger basepair inclination and slide in dsRNA is responsible for the apparently weaker stretching rigidity versus that of dsDNA, and the opposite twist-stretch coupling for dsRNA and dsDNA is also attributed to the stronger basepair inclination in dsRNA than in dsDNA. Our calculated macroscopic elastic parameters and microscopic analyses are tested and validated by different force fields for both dsRNA and dsDNA.
Collapse
Affiliation(s)
- Lei Bao
- Center for Theoretical Physics and Key Laboratory of Artificial Micro- & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Xi Zhang
- Center for Theoretical 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
- Center for Theoretical Physics and Key Laboratory of Artificial Micro- & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China; Research Center of Nonlinear Science, School of Mathematics and Computer Science, Wuhan Textile University, Wuhan, China
| | - Yuan-Yan Wu
- Center for Theoretical Physics and Key Laboratory of Artificial Micro- & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China; College of Physical Science and Technology, Yangzhou University, Yangzhou, China
| | - Zhi-Jie Tan
- Center for Theoretical Physics and Key Laboratory of Artificial Micro- & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
17
|
Darré L, Ivani I, Dans PD, Gómez H, Hospital A, Orozco M. Small Details Matter: The 2′-Hydroxyl as a Conformational Switch in RNA. J Am Chem Soc 2016; 138:16355-16363. [DOI: 10.1021/jacs.6b09471] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Leonardo Darré
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Ivan Ivani
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Pablo D. Dans
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Hansel Gómez
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute
for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Joint
BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine, 08028 Barcelona, Spain
- Department
of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
18
|
Xu Y, Vanommeslaeghe K, Aleksandrov A, MacKerell AD, Nilsson L. Additive CHARMM force field for naturally occurring modified ribonucleotides. J Comput Chem 2016; 37:896-912. [PMID: 26841080 PMCID: PMC4801715 DOI: 10.1002/jcc.24307] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Revised: 11/16/2015] [Accepted: 01/06/2016] [Indexed: 01/13/2023]
Abstract
More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- You Xu
- Department of Biosciences and Nutrition, Karolinska Institutet, HUDDINGE, SE-141 83, Sweden
| | - Kenno Vanommeslaeghe
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland, 21201.,Department of Analytical Chemistry and Pharmaceutical Technology (FABI), Center for Pharmaceutical Research (CePhaR), Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, Brussels, B-1090, Belgium
| | - Alexey Aleksandrov
- Department of Biology, Ecole Polytechnique, Laboratoire De Biochimie (CNRS UMR7654), Palaiseau, F-91128, France
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland, 21201
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, HUDDINGE, SE-141 83, Sweden
| |
Collapse
|
19
|
Juneja A, Villa A, Nilsson L. Elucidating the Relation between Internal Motions and Dihedral Angles in an RNA Hairpin Using Molecular Dynamics. J Chem Theory Comput 2015; 10:3532-40. [PMID: 26588317 DOI: 10.1021/ct500203m] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Molecular dynamics simulations were performed to characterize the internal motions of the ribonucleic acid apical stem loop of human hepatitis B virus. The NMR relaxation rates calculated directly from the trajectory are in good agreement with the experiment. Calculated order parameters follow the experimental pattern. Order parameters lower than 0.8 are observed for nucleotides that are weakly hydrogen bonded to their base pair partner, unpaired, or part of the loop. These residues show slow decay of the internal correlation functions of their base and sugar C-H vectors. Concerted motions around backbone dihedral angles influence the amplitude of motion of the sugar and base C-H vectors. The order parameters for base C-H vectors are also affected by the fluctuation of the glycosidic dihedral angle.
Collapse
Affiliation(s)
- Alok Juneja
- Department of Biosciences and Nutrition, Center of Biosciences, Karolinska Institutet , SE-141 83 Huddinge, Sweden
| | - Alessandra Villa
- Department of Biosciences and Nutrition, Center of Biosciences, Karolinska Institutet , SE-141 83 Huddinge, Sweden
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Center of Biosciences, Karolinska Institutet , SE-141 83 Huddinge, Sweden
| |
Collapse
|
20
|
Vanommeslaeghe K, MacKerell AD. CHARMM additive and polarizable force fields for biophysics and computer-aided drug design. BIOCHIMICA ET BIOPHYSICA ACTA 2015; 1850:861-871. [PMID: 25149274 PMCID: PMC4334745 DOI: 10.1016/j.bbagen.2014.08.004] [Citation(s) in RCA: 174] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/07/2014] [Revised: 08/10/2014] [Accepted: 08/12/2014] [Indexed: 11/29/2022]
Abstract
BACKGROUND Molecular Mechanics (MM) is the method of choice for computational studies of biomolecular systems owing to its modest computational cost, which makes it possible to routinely perform molecular dynamics (MD) simulations on chemical systems of biophysical and biomedical relevance. SCOPE OF REVIEW As one of the main factors limiting the accuracy of MD results is the empirical force field used, the present paper offers a review of recent developments in the CHARMM additive force field, one of the most popular biomolecular force fields. Additionally, we present a detailed discussion of the CHARMM Drude polarizable force field, anticipating a growth in the importance and utilization of polarizable force fields in the near future. Throughout the discussion emphasis is placed on the force fields' parametrization philosophy and methodology. MAJOR CONCLUSIONS Recent improvements in the CHARMM additive force field are mostly related to newly found weaknesses in the previous generation of additive force fields. Beyond the additive approximation is the newly available CHARMM Drude polarizable force field, which allows for MD simulations of up to 1μs on proteins, DNA, lipids and carbohydrates. GENERAL SIGNIFICANCE Addressing the limitations ensures the reliability of the new CHARMM36 additive force field for the types of calculations that are presently coming into routine computational reach while the availability of the Drude polarizable force fields offers an inherently more accurate model of the underlying physical forces driving macromolecular structures and dynamics. This article is part of a Special Issue entitled "Recent developments of molecular dynamics".
Collapse
Affiliation(s)
- K Vanommeslaeghe
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| | - A D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA.
| |
Collapse
|
21
|
Unzippers, resolvers and sensors: a structural and functional biochemistry tale of RNA helicases. Int J Mol Sci 2015; 16:2269-93. [PMID: 25622248 PMCID: PMC4346836 DOI: 10.3390/ijms16022269] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2014] [Revised: 01/09/2015] [Accepted: 01/12/2015] [Indexed: 12/28/2022] Open
Abstract
The centrality of RNA within the biological world is an irrefutable fact that currently attracts increasing attention from the scientific community. The panoply of functional RNAs requires the existence of specific biological caretakers, RNA helicases, devoted to maintain the proper folding of those molecules, resolving unstable structures. However, evolution has taken advantage of the specific position and characteristics of RNA helicases to develop new functions for these proteins, which are at the interface of the basic processes for transference of information from DNA to proteins. RNA helicases are involved in many biologically relevant processes, not only as RNA chaperones, but also as signal transducers, scaffolds of molecular complexes, and regulatory elements. Structural biology studies during the last decade, founded in X-ray crystallography, have characterized in detail several RNA-helicases. This comprehensive review summarizes the structural knowledge accumulated in the last two decades within this family of proteins, with special emphasis on the structure-function relationships of the most widely-studied families of RNA helicases: the DEAD-box, RIG-I-like and viral NS3 classes.
Collapse
|
22
|
Sterpone F, Melchionna S, Tuffery P, Pasquali S, Mousseau N, Cragnolini T, Chebaro Y, St-Pierre JF, Kalimeri M, Barducci A, Laurin Y, Tek A, Baaden M, Nguyen PH, Derreumaux P. The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems. Chem Soc Rev 2014; 43:4871-93. [PMID: 24759934 PMCID: PMC4426487 DOI: 10.1039/c4cs00048j] [Citation(s) in RCA: 120] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.
Collapse
Affiliation(s)
- Fabio Sterpone
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, 13 rue Pierre et Marie Curie, 75005, Paris, France.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
23
|
Huang M, Giese TJ, Lee TS, York DM. Improvement of DNA and RNA Sugar Pucker Profiles from Semiempirical Quantum Methods. J Chem Theory Comput 2014; 10:1538-1545. [PMID: 24803866 PMCID: PMC3985690 DOI: 10.1021/ct401013s] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2013] [Indexed: 01/26/2023]
Abstract
Neglect of diatomic differential overlap (NDDO) and self-consistent density-functional tight-binding (SCC-DFTB) semiempirical models commonly employed in combined quantum mechanical/molecular mechanical simulations fail to adequately describe the deoxyribose and ribose sugar ring puckers. This failure limits the application of these methods to RNA and DNA systems. In this work, we provide benchmark ab initio gas-phase two-dimensional potential energy scans of the RNA and DNA sugar puckering. The benchmark calculations are compared with semiempirical models. Pucker corrections are introduced into the semiempirical models via B-spline interpolation of the potential energy difference surface relative to the benchmark data. The corrected semiempirical models are shown to well reproduce the ab initio puckering profiles. Furthermore, we demonstrate that the uncorrected semiempirical models do not usually produce a transition state between the A-form and B-form sugar puckers, but the ab initio transition state is reproduced when the B-spline correction is used.
Collapse
Affiliation(s)
- Ming Huang
- Center
for Integrative Proteomics Research, BioMaPS Institute and Department
of Chemistry and Chemical Biology, Rutgers
University, 174 Frelinghuysen
Road, Piscataway, New Jersey 08854−8076, United States
- Scientific
Computation, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455−0431, United States
| | - Timothy J. Giese
- Center
for Integrative Proteomics Research, BioMaPS Institute and Department
of Chemistry and Chemical Biology, Rutgers
University, 174 Frelinghuysen
Road, Piscataway, New Jersey 08854−8076, United States
| | - Tai-Sung Lee
- Center
for Integrative Proteomics Research, BioMaPS Institute and Department
of Chemistry and Chemical Biology, Rutgers
University, 174 Frelinghuysen
Road, Piscataway, New Jersey 08854−8076, United States
| | - Darrin M. York
- Center
for Integrative Proteomics Research, BioMaPS Institute and Department
of Chemistry and Chemical Biology, Rutgers
University, 174 Frelinghuysen
Road, Piscataway, New Jersey 08854−8076, United States
| |
Collapse
|
24
|
Mládek A, Banáš P, Jurečka P, Otyepka M, Zgarbová M, Šponer J. Energies and 2'-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations. J Chem Theory Comput 2013; 10:463-80. [PMID: 26579924 DOI: 10.1021/ct400837p] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.
Collapse
Affiliation(s)
- Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC, Central European Institute of Technology , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
25
|
Šponer J, Mládek A, Šponer JE, Svozil D, Zgarbová M, Banáš P, Jurečka P, Otyepka M. The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies. Phys Chem Chem Phys 2012; 14:15257-77. [PMID: 23072945 DOI: 10.1039/c2cp41987d] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic.
| | | | | | | | | | | | | | | |
Collapse
|