1
|
Chakraborty A, Samant D, Sarkar R, Sangeet S, Prusty S, Roy S. RNA's Dynamic Conformational Selection and Entropic Allosteric Mechanism in Controlling Cascade Protein Binding Events. J Phys Chem Lett 2024; 15:6115-6125. [PMID: 38830201 DOI: 10.1021/acs.jpclett.4c00740] [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: 06/05/2024]
Abstract
In the TAR RNA of immunodeficiency viruses, an allosteric communication exists between a distant loop and a bulge. The bulge interacts with the TAT protein vital for transactivating viral RNA, while the loop interacts with cyclin-T1, contingent on TAT binding. Through extensive atomistic and free energy simulations, we investigate TAR-TAT binding in nonpathogenic bovine immunodeficiency virus (BIV) and pathogenic human immunodeficiency virus (HIV). Thermodynamic analysis reveals enthalpically driven binding in BIV and entropically favored binding in HIV. The broader global basin in HIV is attributed to binding-induced loop fluctuation, corroborated by nuclear magnetic resonance (NMR), indicating classical entropic allostery onset. While this loop fluctuation affects the TAT binding affinity, it generates a binding-competent conformation that aids subsequent effector (cyclin-T1) binding. This study underscores how two structurally similar apo-RNA scaffolds adopt distinct conformational selection mechanisms to drive enthalpic and entropic allostery, influencing protein affinity in the signaling cascade.
Collapse
Affiliation(s)
- Amrita Chakraborty
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| | - Dibyamanjaree Samant
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| | - Raju Sarkar
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| | - Satyam Sangeet
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| | - Sangram Prusty
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Kolkata, West Bengal 741246, India
| |
Collapse
|
2
|
Levintov L, Vashisth H. Structural and computational studies of HIV-1 RNA. RNA Biol 2024; 21:1-32. [PMID: 38100535 PMCID: PMC10730233 DOI: 10.1080/15476286.2023.2289709] [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] [Accepted: 11/21/2023] [Indexed: 12/17/2023] Open
Abstract
Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering & Bioengineering, University of New Hampshire, Durham, USA
| | - Harish Vashisth
- Department of Chemical Engineering & Bioengineering, University of New Hampshire, Durham, USA
| |
Collapse
|
3
|
Maghsoud Y, Jayasinghe-Arachchige VM, Kumari P, Cisneros GA, Liu J. Leveraging QM/MM and Molecular Dynamics Simulations to Decipher the Reaction Mechanism of the Cas9 HNH Domain to Investigate Off-Target Effects. J Chem Inf Model 2023; 63:6834-6850. [PMID: 37877218 DOI: 10.1021/acs.jcim.3c01284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
The clustered regularly interspaced short palindromic repeats (CRISPR) technology is an RNA-guided targeted genome-editing tool using Cas family proteins. Two magnesium-dependent nuclease domains of the Cas9 enzyme, termed HNH and RuvC, are responsible for cleaving the target DNA (t-DNA) and nontarget DNA strands, respectively. The HNH domain is believed to determine the DNA cleavage activity of both endonuclease domains and is sensitive to complementary RNA-DNA base pairing. However, the underlying molecular mechanisms of CRISPR-Cas9, by which it rebukes or accepts mismatches, are poorly understood. Thus, investigation of the structure and dynamics of the catalytic state of Cas9 with either matched or mismatched t-DNA can provide insights into improving its specificity by reducing off-target cleavages. Here, we focus on a recently discovered catalytic-active form of the Streptococcus pyogenes Cas9 (SpCas9) and employ classical molecular dynamics and coupled quantum mechanics/molecular mechanics simulations to study two possible mechanisms of t-DNA cleavage reaction catalyzed by the HNH domain. Moreover, by designing a mismatched t-DNA structure called MM5 (C to G at the fifth position from the protospacer adjacent motif region), the impact of single-guide RNA (sgRNA) and t-DNA complementarity on the catalysis process was investigated. Based on these simulations, our calculated binding affinities, minimum energy paths, and analysis of catalytically important residues provide atomic-level details of the differences between matched and mismatched cleavage reactions. In addition, several residues exhibit significant differences in their catalytic roles for the two studied systems, including K253, K263, R820, K896, and K913.
Collapse
Affiliation(s)
- Yazdan Maghsoud
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Vindi M Jayasinghe-Arachchige
- Department of Pharmaceutical Sciences, University of North Texas System College of Pharmacy, University of North Texas Health Science Center, Fort Worth, Texas 76107, United States
| | - Pratibha Kumari
- Department of Pharmaceutical Sciences, University of North Texas System College of Pharmacy, University of North Texas Health Science Center, Fort Worth, Texas 76107, United States
| | - G Andrés Cisneros
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, United States
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Jin Liu
- Department of Pharmaceutical Sciences, University of North Texas System College of Pharmacy, University of North Texas Health Science Center, Fort Worth, Texas 76107, United States
| |
Collapse
|
4
|
Abstract
Being able to effectively target RNA with potent ligands will open up a large number of potential therapeutic options. The knowledge on how to achieve this is ever expanding but an important question that remains open is what chemical matter is suitable to achieve this goal. The high flexibility of an RNA as well as its more limited chemical diversity and featureless binding sites can be difficult to target selectively but can be addressed by well-designed cyclic peptides. In this review we will provide an overview of reported cyclic peptide ligands for therapeutically relevant RNA targets and discuss the methods used to discover them. We will also provide critical insights into the properties required for potent and selective interaction and suggestions on how to assess these parameters. The use of cyclic peptides to target RNA is still in its infancy but the lessons learned from past examples can be adopted for the development of novel potent and selective ligands.
Collapse
Affiliation(s)
- Sunit Pal
- Chemical Genomics Centre of the Max Planck Society, Max Planck Institute of Molecular Physiology, Dortmund, Germany
| | - Peter 't Hart
- Chemical Genomics Centre of the Max Planck Society, Max Planck Institute of Molecular Physiology, Dortmund, Germany
| |
Collapse
|
5
|
Naseem-Khan S, Berger MB, Leddin EM, Maghsoud Y, Cisneros GA. Impact of Remdesivir Incorporation along the Primer Strand on SARS-CoV-2 RNA-Dependent RNA Polymerase. J Chem Inf Model 2022; 62:2456-2465. [PMID: 35435671 DOI: 10.1021/acs.jcim.2c00201] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Remdesivir was the first antiviral drug that received emergency use authorization from the United States Food and Drug Administration and is now formally approved to treat COVID-19. Remdesivir is a nucleotide analogue that targets the RNA-dependent RNA polymerase (RdRp) of coronaviruses, including SARS-CoV-2. The solution of multiple RdRp structures has been one of the main axes of research in the race against the SARS-CoV-2 virus. Several hypotheses of the mechanism of inhibition of RdRp by remdesivir have been proposed, although open questions remain. This work uses molecular dynamics simulations to explore the impact of remdesivir and two analogues as incoming nucleotides and of up to four incorporations of remdesivir along the primer strand on RdRp. The simulation results suggest that the overall structure and the dynamical behavior of RdRp are destabilized by remdesivir and the two analogues in the incoming position. The incorporation of remdesivir along the primer strand impacts specific non-bonded interactions between the nascent RNA and the polymerase subunit, as well as the overall dynamical networks on RdRp. The strongest impact on the structure and dynamics are observed after three incorporations, when remdesivir is located at position -A3, in agreement with previously reported experimental and computational results. Our results provide atomic-level details of the role played by remdesivir on the disruption of RNA synthesis by RdRp and the main drivers of these disruptions.
Collapse
Affiliation(s)
- Sehr Naseem-Khan
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Madison B Berger
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Emmett M Leddin
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Yazdan Maghsoud
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| |
Collapse
|
6
|
Levintov L, Vashisth H. Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides. Biophys J 2021; 120:5060-5073. [PMID: 34710377 PMCID: PMC8633718 DOI: 10.1016/j.bpj.2021.10.007] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 08/17/2021] [Accepted: 10/06/2021] [Indexed: 12/14/2022] Open
Abstract
Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering, University of New Hampshire, Durham, New Hampshire
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, New Hampshire.
| |
Collapse
|
7
|
Kumar A, Vashisth H. Conformational dynamics and energetics of viral RNA recognition by lab-evolved proteins. Phys Chem Chem Phys 2021; 23:24773-24779. [PMID: 34714308 PMCID: PMC8579469 DOI: 10.1039/d1cp03822b] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2021] [Accepted: 10/14/2021] [Indexed: 12/01/2022]
Abstract
The conserved and structured elements in viral RNA genomes interact with proteins to regulate various events in the viral life cycle and have become key targets for developing novel therapeutic approaches. We probe physical interactions between lab-evolved proteins and a viral RNA element from the HIV-1 genome. Specifically, we study the role of an arginine-rich loop in recognition of designed proteins by the viral RNA element. We report free energy calculations to quantitatively estimate the protein/RNA binding energetics, focusing on the mutations of arginine residues involved in recognition of the major groove of RNA by proteins.
Collapse
Affiliation(s)
- Amit Kumar
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, USA.
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, USA.
| |
Collapse
|
8
|
Gong W, Liu Y, Zhao Y, Wang S, Han Z, Li C. Equally Weighted Multiscale Elastic Network Model and Its Comparison with Traditional and Parameter-Free Models. J Chem Inf Model 2021; 61:921-937. [PMID: 33496590 DOI: 10.1021/acs.jcim.0c01178] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dynamical properties of proteins play an essential role in their function exertion. The elastic network model (ENM) is an effective and efficient tool in characterizing the intrinsic dynamical properties encoded in biomacromolecule structures. The Gaussian network model (GNM) and anisotropic network model (ANM) are the two often-used ENM models. Here, we introduce an equally weighted multiscale ENM (equally weighted mENM) based on the original mENM (denoted as mENM), in which fitting weights of Kirchhoff/Hessian matrixes in mENM are removed since they neglect the details of pairwise interactions. Then, we perform its comparison with the mENM, traditional ENM, and parameter-free ENM (pfENM) in reproducing dynamical properties for the six representative proteins whose molecular dynamics (MD) trajectories are available in http://mmb.pcb.ub.es/MoDEL/. In the results, for B-factor prediction, mENM performs best, while the equally weighted mENM performs also well, better than the traditional ENM and pfENM models. As to the dynamical cross-correlation map calculation, mENM performs worst, while the results produced from the equally weighted mENM and pfENM models are close to those from MD trajectories with the latter a little better than the former. Furthermore, encouragingly, the equally weighted mANM displays the best performance in capturing the functional motional modes, followed by pfANM and traditional ANM models, while the mANM fails in all the cases. This work is helpful for strengthening the understanding of the elastic network model and provides a valuable guide for researchers to utilize the model to explore protein dynamics.
Collapse
Affiliation(s)
- Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yang Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yanpeng Zhao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Shihao Wang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
9
|
Collias D, Leenay RT, Slotkowski RA, Zuo Z, Collins SP, McGirr BA, Liu J, Beisel CL. A positive, growth-based PAM screen identifies noncanonical motifs recognized by the S. pyogenes Cas9. SCIENCE ADVANCES 2020; 6:eabb4054. [PMID: 32832642 PMCID: PMC7439565 DOI: 10.1126/sciadv.abb4054] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 05/28/2020] [Indexed: 05/16/2023]
Abstract
CRISPR technologies have overwhelmingly relied on the Streptococcus pyogenes Cas9 (SpyCas9), with its consensus NGG and less preferred NAG and NGA protospacer-adjacent motifs (PAMs). Here, we report that SpyCas9 also recognizes sequences within an N(A/C/T)GG motif. These sequences were identified on the basis of preferential enrichment in a growth-based screen in Escherichia coli. DNA binding, cleavage, and editing assays in bacteria and human cells validated recognition, with activities paralleling those for NAG(A/C/T) PAMs and dependent on the first two PAM positions. Molecular-dynamics simulations and plasmid-clearance assays with mismatch-intolerant variants supported induced-fit recognition of an extended PAM by SpyCas9 rather than recognition of NGG with a bulged R-loop. Last, the editing location for SpyCas9-derived base editors could be shifted by one nucleotide by selecting between (C/T)GG and adjacent N(C/T)GG PAMs. SpyCas9 and its enhanced variants thus recognize a larger repertoire of PAMs, with implications for precise editing, off-target predictions, and CRISPR-based immunity.
Collapse
Affiliation(s)
- D. Collias
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
| | - R. T. Leenay
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
| | - R. A. Slotkowski
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
| | - Z. Zuo
- College of Chemical Engineering, Shanghai University of Engineering Science, Shanghai 201620, China
- Department of Pharmaceutical Sciences, University of North Texas System College of Pharmacy, University of North Texas Health Science Center, Fort Worth, TX 76129, USA
| | - S. P. Collins
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
| | - B. A. McGirr
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
| | - J. Liu
- Department of Pharmaceutical Sciences, University of North Texas System College of Pharmacy, University of North Texas Health Science Center, Fort Worth, TX 76129, USA
| | - C. L. Beisel
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA
- Helmholtz Institute of RNA-based Infection Research (HIRI), Helmholtz Centre for Infection Research (HZI), 97080 Würzburg, Germany
- Faculty of Medicine, University of Würzburg, 97078 Würzburg, Germany
| |
Collapse
|
10
|
Moghimi S, Morsali A, Heravi MM, Beyramabadi SA. Quantum‐Chemical Modeling of Cyclic Peptide‐Selenium Nanoparticle as an Anticancer Drug Nanocarrier. B KOREAN CHEM SOC 2019. [DOI: 10.1002/bkcs.11912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Sara Moghimi
- Department of Chemistry, Mashhad BranchIslamic Azad University Mashhad Iran
| | - Ali Morsali
- Department of Chemistry, Mashhad BranchIslamic Azad University Mashhad Iran
- Research Center for Animal Development Applied Biology, Mashhad BranchIslamic Azad University Mashhad 917568 Iran
| | - Mohammad M. Heravi
- Department of Chemistry, Mashhad BranchIslamic Azad University Mashhad Iran
| | | |
Collapse
|
11
|
Zuo Z, Zolekar A, Babu K, Lin VJT, Hayatshahi HS, Rajan R, Wang YC, Liu J. Structural and functional insights into the bona fide catalytic state of Streptococcus pyogenes Cas9 HNH nuclease domain. eLife 2019; 8:e46500. [PMID: 31361218 PMCID: PMC6706240 DOI: 10.7554/elife.46500] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Accepted: 07/21/2019] [Indexed: 12/21/2022] Open
Abstract
The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (SpyCas9), along with a programmable single-guide RNA (sgRNA), has been exploited as a significant genome-editing tool. Despite the recent advances in determining the SpyCas9 structures and DNA cleavage mechanism, the cleavage-competent conformation of the catalytic HNH nuclease domain of SpyCas9 remains largely elusive and debatable. By integrating computational and experimental approaches, we unveiled and validated the activated Cas9-sgRNA-DNA ternary complex in which the HNH domain is neatly poised for cleaving the target DNA strand. In this catalysis model, the HNH employs the catalytic triad of D839-H840-N863 for cleavage catalysis, rather than previously implicated D839-H840-D861, D837-D839-H840, or D839-H840-D861-N863. Our study contributes critical information to defining the catalytic conformation of the HNH domain and advances the knowledge about the conformational activation underlying Cas9-mediated DNA cleavage.
Collapse
Affiliation(s)
- Zhicheng Zuo
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
- College of Chemistry and Chemical EngineeringShanghai University of Engineering ScienceShanghaiChina
| | - Ashwini Zolekar
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
| | - Kesavan Babu
- Department of Chemistry and Biochemistry, Price Family Foundation Institute of Structural Biology, Stephenson Life Sciences Research CenterUniversity of OklahomaNormanUnited States
| | - Victor JT Lin
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
| | - Hamed S Hayatshahi
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
| | - Rakhi Rajan
- Department of Chemistry and Biochemistry, Price Family Foundation Institute of Structural Biology, Stephenson Life Sciences Research CenterUniversity of OklahomaNormanUnited States
| | - Yu-Chieh Wang
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
| | - Jin Liu
- Department of Pharmaceutical SciencesUNT System College of Pharmacy, University of North Texas Health Science CenterFort WorthUnited States
| |
Collapse
|
12
|
Zuo Z, Liu J. Assessing the Performance of the Nonbonded Mg 2+ Models in a Two-Metal-Dependent Ribonuclease. J Chem Inf Model 2018; 59:399-408. [PMID: 30521334 DOI: 10.1021/acs.jcim.8b00627] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Magnesium ions (Mg2+), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg2+ has motivated continuous development and improvement of various Mg2+ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg2+ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg2+ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg2+-containing biomolecular systems.
Collapse
Affiliation(s)
- Zhicheng Zuo
- Department of Pharmaceutical Sciences , University of North Texas System College of Pharmacy, University of North Texas Health Science Center , Fort Worth , Texas 76107 , United States
| | - Jin Liu
- Department of Pharmaceutical Sciences , University of North Texas System College of Pharmacy, University of North Texas Health Science Center , Fort Worth , Texas 76107 , United States
| |
Collapse
|
13
|
Identification of a unique Ca 2+-binding site in rat acid-sensing ion channel 3. Nat Commun 2018; 9:2082. [PMID: 29802295 PMCID: PMC5970173 DOI: 10.1038/s41467-018-04424-0] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2017] [Accepted: 04/17/2018] [Indexed: 12/20/2022] Open
Abstract
Acid-sensing ion channels (ASICs) evolved to sense changes in extracellular acidity with the divalent cation calcium (Ca2+) as an allosteric modulator and channel blocker. The channel-blocking activity is most apparent in ASIC3, as removing Ca2+ results in channel opening, with the site's location remaining unresolved. Here we show that a ring of rat ASIC3 (rASIC3) glutamates (Glu435), located above the channel gate, modulates proton sensitivity and contributes to the formation of the elusive Ca2+ block site. Mutation of this residue to glycine, the equivalent residue in chicken ASIC1, diminished the rASIC3 Ca2+ block effect. Atomistic molecular dynamic simulations corroborate the involvement of this acidic residue in forming a high-affinity Ca2+ site atop the channel pore. Furthermore, the reported observations provide clarity for past controversies regarding ASIC channel gating. Our findings enhance understanding of ASIC gating mechanisms and provide structural and energetic insights into this unique calcium-binding site.
Collapse
|
14
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
15
|
Peddi SR, Sivan SK, Manga V. Molecular dynamics and MM/GBSA-integrated protocol probing the correlation between biological activities and binding free energies of HIV-1 TAR RNA inhibitors. J Biomol Struct Dyn 2017; 36:486-503. [PMID: 28081678 DOI: 10.1080/07391102.2017.1281762] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The interaction of HIV-1 transactivator protein Tat with its cognate transactivation response (TAR) RNA has emerged as a promising target for developing antiviral compounds and treating HIV infection, since it is a crucial step for efficient transcription and replication. In the present study, molecular dynamics (MD) simulations and MM/GBSA calculations have been performed on a series of neamine derivatives in order to estimate appropriate MD simulation time for acceptable correlation between ΔGbind and experimental pIC50 values. Initially, all inhibitors were docked into the active site of HIV-1 TAR RNA. Later to explore various conformations and examine the docking results, MD simulations were carried out. Finally, binding free energies were calculated using MM/GBSA method and were correlated with experimental pIC50 values at different time scales (0-1 to 0-10 ns). From this study, it is clear that in case of neamine derivatives as simulation time increased the correlation between binding free energy and experimental pIC50 values increased correspondingly. Therefore, the binding energies which can be interpreted at longer simulation times can be used to predict the bioactivity of new neamine derivatives. Moreover, in this work, we have identified some plausible critical nucleotide interactions with neamine derivatives that are responsible for potent inhibitory activity. Furthermore, we also provide some insights into a new class of oxadiazole-based back bone cyclic peptides designed by incorporating the structural features of neamine derivatives. On the whole, this approach can provide a valuable guidance for designing new potent inhibitors and modify the existing compounds targeting HIV-1 TAR RNA.
Collapse
Affiliation(s)
- Saikiran Reddy Peddi
- a Molecular Modeling and Medicinal Chemistry Group, Department of Chemistry , University College of Science, Osmania University , Hyderabad 500 007 , Telangana , India
| | - Sree Kanth Sivan
- b Department of Chemistry , Nizam College, Osmania University , Hyderabad 500 001 , Telangana , India
| | - Vijjulatha Manga
- a Molecular Modeling and Medicinal Chemistry Group, Department of Chemistry , University College of Science, Osmania University , Hyderabad 500 007 , Telangana , India
| |
Collapse
|
16
|
Xiao X, Zhao B, Agris PF, Hall CK. Simulation study of the ability of a computationally-designed peptide to recognize target tRNA Lys3 and other decoy tRNAs. Protein Sci 2016; 25:2243-2255. [PMID: 27680513 DOI: 10.1002/pro.3056] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2016] [Accepted: 09/25/2016] [Indexed: 11/06/2022]
Abstract
In this paper, we investigate the ability of our computationally-designed peptide, Pept10 (PNWNGNRWLNNCLRG), to recognize the anticodon stem and loop (ASL) domain of the hypermodified tRNALys3 (mcm5 s2 U34 ,ms2 t6 A37 ), a reverse transcription primer of HIV replication. Five other ASLs, the singly modified ASLLys3 (ms2 t6 A37 ), ASLLys3 (s2 U34 ), ASLLys3 (Ψ39 ), ASLLys1,2 (t6 A37 ), and ASLGlu (s2 U34 ), were used as decoys. Explicit-solvent atomistic molecular dynamics simulations were performed to examine the process of binding of Pept10 with the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) and the decoy ASLs. Simulation results demonstrated that Pept10 is capable of recognizing the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) as well as one of the decoys, ASLLys3 (Ψ39 ), but screens out the other four decoy ASLs. The interchain van der Waals (VDW) and charge-charge (ELE + EGB) energies for the two best complexes were evaluated to shed light on the molecular recognition mechanism between Pept10 and ASLs. The results indicated that Pept10 recognizes and binds to the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) through residues W3 and R7 which interact with the nucleotides mcm5 s2 U34 , U35 , and ms2 t6 A37 via the interchain VDW energy. Pept10 also recognizes the decoy ASLLys3 (Ψ39 ) through residue R14 which contacts the nucleotide U36 via the interchain VDW energy. Regardless of the type of ASL, the positively charged arginines on Pept10 are attracted to the negatively charged phosphate linkages on the ASL via the interchain ELE + EGB energy, thereby enhancing the binding affinity.
Collapse
Affiliation(s)
- Xingqing Xiao
- Chemical and Biomolecular Engineering Department, North Carolina State University, Raleigh, North Carolina, 95-2767905
| | - Binwu Zhao
- Chemical and Biomolecular Engineering Department, North Carolina State University, Raleigh, North Carolina, 95-2767905
| | - Paul F Agris
- The RNA Institute, University at Albany, State University of New York, Albany, New York, 12222
| | - Carol K Hall
- Chemical and Biomolecular Engineering Department, North Carolina State University, Raleigh, North Carolina, 95-2767905
| |
Collapse
|
17
|
Zuo Z, Weng J, Wang W. Insights into the Inhibitory Mechanism of D13-9001 to the Multidrug Transporter AcrB through Molecular Dynamics Simulations. J Phys Chem B 2016; 120:2145-54. [PMID: 26900716 DOI: 10.1021/acs.jpcb.5b11942] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The resistance-nodulation-cell division transporter AcrB is responsible for energy transduction and substrate recognition in the tripartite AcrAB-TolC efflux system in Escherichia coli. Despite a broad substrate specificity, only a few compounds have been cocrystallized with AcrB inside the distal binding pocket (DBP), including doxorubicin (DOX) and D13-9001. D13-9001 is a promising efflux pump inhibitor that potentiates the efficacy of a wide variety of antibiotics. To understand its inhibition effect under the framework of functional rotating mechanism, we performed targeted and steered molecular dynamics simulations to compare the binding and extrusion processes of this inhibitor and the substrate DOX in AcrB. The results demonstrate that, with respect to DOX, the interaction of D13-9001 with the hydrophobic trap results in delayed disassociation from the DBP. Notably, the detachment of D13-9001 is tightly correlated with the side-chain reorientation of Phe628 and large-scale displacement of Tyr327. Furthermore, the inhibitor induces much more significant conformational changes at the exit gate than DOX does, thereby causing higher energy cost for extrusion and contributing to the inhibitory effect in addition to the tight binding at DBP.
Collapse
Affiliation(s)
- Zhicheng Zuo
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, and Institutes of Biomedical Sciences, Fudan University , Shanghai 200433, People's Republic of China
| | - Jingwei Weng
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, and Institutes of Biomedical Sciences, Fudan University , Shanghai 200433, People's Republic of China
| | - Wenning Wang
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, and Institutes of Biomedical Sciences, Fudan University , Shanghai 200433, People's Republic of China
| |
Collapse
|
18
|
Post-Transcriptional Modifications of RNA: Impact on RNA Function and Human Health. MODIFIED NUCLEIC ACIDS IN BIOLOGY AND MEDICINE 2016. [DOI: 10.1007/978-3-319-34175-0_5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
|
19
|
Stefaniak F, Chudyk EI, Bodkin M, Dawson WK, Bujnicki JM. Modeling of ribonucleic acid-ligand interactions. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2015. [DOI: 10.1002/wcms.1226] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Filip Stefaniak
- Laboratory of Bioinformatics and Protein Engineering; International Institute of Molecular and Cell Biology; Warsaw Poland
| | - Ewa I. Chudyk
- Research Informatics; Evotec (UK) Ltd; Milton Park UK
| | | | - Wayne K. Dawson
- Laboratory of Bioinformatics and Protein Engineering; International Institute of Molecular and Cell Biology; Warsaw Poland
| | - Janusz M. Bujnicki
- Laboratory of Bioinformatics and Protein Engineering; International Institute of Molecular and Cell Biology; Warsaw Poland
- Laboratory of Bioinformatics, Institute of Molecular Biology and Biotechnology, Faculty of Biology; Adam Mickiewicz University; Poznan Poland
| |
Collapse
|
20
|
Xiao X, Agris PF, Hall CK. Molecular recognition mechanism of peptide chain bound to the tRNA(Lys3) anticodon loop in silico. J Biomol Struct Dyn 2014; 33:14-27. [PMID: 24417415 DOI: 10.1080/07391102.2013.869660] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
The mechanism by which proteins recognize and bind the post-transcriptional modifications of RNAs is unknown, yet these interactions play important functions in biology. Atomistic molecular dynamics simulations were performed to examine the folding of the model peptide chain -RVTHHAFLGAHRTVG- and the complex formed by the folded peptide with the native anticodon stem and loop of the human tRNA(Lys3) (hASL(Lys3)) in order to explore the binding mechanism. By analyzing and comparing two folded conformations of this peptide obtained from the folding simulation, we found that the van der Waals (VDW) energy is necessary for the thermal stability of the peptide, and the charge-charge (ELE + EGB) energy is crucial for determining the three-dimensional folded structure of the peptide backbone. Subsequently, two conformations of the peptide were employed to investigate their binding behaviors to hASL(Lys3). The metastable folded peptide was found to bind to hASL(Lys3) much easier than the stable folded peptide in the binding simulations. An energetic analysis reveals that the VDW energy favors the binding, whereas the ELE + EGB energies disfavor the binding. Arginines on the peptide preferentially attract the phosphate backbone via the inter-chain ELE + EGB interaction, significantly contributing to the binding affinity. The hydrophobic phenylalanine interacts with the anticodon loop of hASL(Lys3) via the inter-chain VDW interaction, significantly contributing to the binding specificity.
Collapse
Affiliation(s)
- Xingqing Xiao
- a Chemical and Biomolecular Engineering Department , North Carolina State University , Engineering Building I, 911 Partners Way, Raleigh , NC 27695-7905 , USA
| | | | | |
Collapse
|