1
|
Coskuner-Weber O. Structures prediction and replica exchange molecular dynamics simulations of α-synuclein: A case study for intrinsically disordered proteins. Int J Biol Macromol 2024; 276:133813. [PMID: 38996889 DOI: 10.1016/j.ijbiomac.2024.133813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2024] [Revised: 07/08/2024] [Accepted: 07/09/2024] [Indexed: 07/14/2024]
Abstract
In recent years, a variety of three-dimensional structure prediction tools, including AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold, have been employed in the investigation of intrinsically disordered proteins. However, a comprehensive validation of these tools specifically for intrinsically disordered proteins has yet to be conducted. In this study, we utilize AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold to predict the structure of a model intrinsically disordered α-synuclein protein. Additionally, extensive replica exchange molecular dynamics simulations of the intrinsically disordered protein are conducted. The resulting structures from both structure prediction tools and replica exchange molecular dynamics simulations are analyzed for radius of gyration, secondary and tertiary structure properties, as well as Cα and Hα chemical shift values. A comparison of the obtained results with experimental data reveals that replica exchange molecular dynamics simulations provide results in excellent agreement with experimental observations. However, none of the structure prediction tools utilized in this study can fully capture the structural characteristics of the model intrinsically disordered protein. This study shows that a cluster of ensembles are required for intrinsically disordered proteins. Artificial-intelligence based structure prediction tools such as AlphaFold3 and C-I-TASSER could benefit from stochastic sampling or Monte Carlo simulations for generating an ensemble of structures for intrinsically disordered proteins.
Collapse
Affiliation(s)
- Orkid Coskuner-Weber
- Turkish-German University, Molecular Biotechnology, Sahinkaya Caddesi, No. 106, Beykoz, Istanbul 34820, Turkey.
| |
Collapse
|
2
|
Muscat S, Martino G, Manigrasso J, Marcia M, De Vivo M. On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules. J Chem Theory Comput 2024. [PMID: 39150960 DOI: 10.1021/acs.jctc.4c00773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2024]
Abstract
RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.
Collapse
Affiliation(s)
- Stefano Muscat
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Gianfranco Martino
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Jacopo Manigrasso
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, 431 50 Mölndal, Sweden
| | - Marco Marcia
- European Molecular Biology Laboratory Grenoble, 71 Avenue des Martyrs, 38042 Grenoble, France
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, 751 23 Uppsala, Sweden
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
3
|
Knappeová B, Mlýnský V, Pykal M, Šponer J, Banáš P, Otyepka M, Krepl M. Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes. J Chem Theory Comput 2024; 20:6917-6929. [PMID: 39012172 PMCID: PMC11325551 DOI: 10.1021/acs.jctc.4c00601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.
Collapse
Affiliation(s)
- Barbora Knappeová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Martin Pykal
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Pavel Banáš
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Michal Otyepka
- Czech Advanced Technology and Research Institute, CATRIN, Palacký University, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB-Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| |
Collapse
|
4
|
Svensson O, Bakker MJ, Skepö M. Deeper Insight of the Conformational Ensemble of Intrinsically Disordered Proteins. J Chem Inf Model 2024; 64:6105-6114. [PMID: 39056166 PMCID: PMC11323008 DOI: 10.1021/acs.jcim.4c00941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2024] [Revised: 07/17/2024] [Accepted: 07/17/2024] [Indexed: 07/28/2024]
Abstract
It is generally known that, unlike structured proteins, intrinsically disordered proteins, IDPs, exhibit various structures and conformers, the so-called conformational ensemble, CoE. This study aims to better understand the conformers that make up the IDP ensemble by decomposing the CoE into groups separated by their radius of gyration, Rg. A common approach to studying CoE for IDPs is to use low-resolution techniques, such as small-angle scattering, and combine those with computer simulations on different length scales. Herein, the well-studied antimicrobial saliva protein histatin 5 was utilized as a model peptide for an IDP; the average intensity curves were obtained from small-angle X-ray scattering; and compared with fully atomistic, explicit water, molecular dynamics simulations; then, the intensity curve was decomposed with respect to the different Rg values; and their secondary structure propensities were investigated. We foresee that this approach can provide important information on the CoE and the individual conformers within; in that case, it will serve as an additional tool for understanding the IDP structure-function relationship on a more detailed level.
Collapse
Affiliation(s)
- Oskar Svensson
- Division
of Computational Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
- NanoLund, Lund University, Box 118, 22100 Lund, Sweden
| | - Michael J. Bakker
- Faculty
of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic
| | - Marie Skepö
- Division
of Computational Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
- NanoLund, Lund University, Box 118, 22100 Lund, Sweden
| |
Collapse
|
5
|
McIvor JAP, Larsen DS, Mercadante D. Charge Relaying within a Phospho-Motif Rescue Binding Competency of a Disordered Transcription Factor. J Chem Inf Model 2024; 64:6041-6052. [PMID: 39074869 DOI: 10.1021/acs.jcim.4c00286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
Structural disorder in proteins is central to cellular signaling, where conformational plasticity equips molecules to promiscuously interact with different partners. By engaging with multiple binding partners via the rearrangement of its three helices, the nuclear coactivator binding domain (NCBD) of the CBP/p300 transcription factor is a paradigmatic example of promiscuity. Recently, molecular simulations and experiments revealed that, through the establishment of long-range electrostatic interactions, intended as salt-bridges formed between the post-translationally inserted phosphate and positively charged residues in helix H3 of NCBD, phosphorylation triggers NCBD compaction, lowering its affinity for binding partners. By means of extensive molecular simulations, we here investigated the effect of short-range electrostatics on the conformational ensemble of NCBD, by monitoring the interactions between a phosphorylated serine and conserved positively charged residues within the NCBD phospho-motif. We found that empowering proximal electrostatic interactions, as opposed to long-range electrostatics, can reshape the NCBD ensemble rescuing the binding competency of phosphorylated NCBD. Given the conservation of positive charges in phospho-motifs, proximal electrostatic interactions might dampen the effects of phosphorylation and act as a relay to regulate phosphorylated intrinsically disordered proteins, ultimately tuning the binding affinity for different cellular partners.
Collapse
Affiliation(s)
- Jordan A P McIvor
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Danaé S Larsen
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| |
Collapse
|
6
|
Wang P, Demaray J, Moroz S, Stuchebrukhov AA. Searching for proton transfer channels in respiratory complex I. Biophys J 2024:S0006-3495(24)00518-6. [PMID: 39095988 DOI: 10.1016/j.bpj.2024.07.041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 06/19/2024] [Accepted: 07/30/2024] [Indexed: 08/04/2024] Open
Abstract
We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.
Collapse
Affiliation(s)
- Panyue Wang
- Department of Chemistry, University of California at Davis, Davis, California
| | - Jackson Demaray
- Department of Chemistry, University of California at Davis, Davis, California
| | - Stanislav Moroz
- Department of Chemistry, University of California at Davis, Davis, California
| | | |
Collapse
|
7
|
Gorantla K, Krishnan A, Waheed SO, Varghese A, DiCastri I, LaRouche C, Paik M, Fields GB, Karabencheva-Christova TG. Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1. Biochemistry 2024; 63:1925-1940. [PMID: 38963231 PMCID: PMC11309001 DOI: 10.1021/acs.biochem.4c00076] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Revised: 06/14/2024] [Accepted: 06/24/2024] [Indexed: 07/05/2024]
Abstract
Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.
Collapse
Affiliation(s)
- Koteswara
Rao Gorantla
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Anandhu Krishnan
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sodiq O. Waheed
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Ann Varghese
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Isabella DiCastri
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Ciara LaRouche
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Meredith Paik
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Gregg B. Fields
- Department
of Chemistry and Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida 33458, United States
| | | |
Collapse
|
8
|
Sedano LF, Blazquez S, Vega C. Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant? J Chem Phys 2024; 161:044505. [PMID: 39046346 DOI: 10.1063/5.0211871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2024] [Accepted: 06/30/2024] [Indexed: 07/25/2024] Open
Abstract
The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.
Collapse
Affiliation(s)
- L F Sedano
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
9
|
Otteson L, Nagy G, Kunkel J, Kodis G, Zheng W, Bignon C, Longhi S, Grubmüller H, Vaiana AC, Vaiana SM. Transient Non-local Interactions Dominate the Dynamics of Measles Virus N TAIL. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.22.604679. [PMID: 39091801 PMCID: PMC11291014 DOI: 10.1101/2024.07.22.604679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
The RNA genome of measles virus is encapsidated by the nucleoprotein within a helical nucleocapsid that serves as template for both transcription and replication. The intrinsically disordered domain of the nucleoprotein (NTAIL), partly protruding outward from the nucleocapsid, is essential for binding the polymerase complex responsible for viral transcription and replication. As for many IDPs, binding of NTAIL occurs through a short molecular recognition element (MoRE) that folds upon binding, with the majority of NTAIL remaining disordered. Though NTAIL regions far from the MoRE influence the binding affinity, interactions between them and the MoRE have not been investigated in depth. Using an integrated approach, relying on photo-induced electron transfer (PET) experiments between tryptophan and cysteine pairs placed at different positions in the protein under varying salt and pH conditions, combined with simulations and analytical models, we identified transient interactions between two disordered regions distant in sequence, which dominate NTAIL dynamics, and regulate the conformational preferences of both the MoRE and the entire NTAIL domain. Co-evolutionary analysis corroborates our findings, and suggests an important functional role for the same intramolecular interactions. We propose mechanisms by which these non-local interactions may regulate binding to the phosphoprotein, polymerase recruitment, and ultimately viral transcription and replication. Our findings may be extended to other IDPs, where non-local intra-protein interactions affect the conformational preferences of intermolecular binding sites.
Collapse
Affiliation(s)
- Lillian Otteson
- Center for Biological Physics, Arizona State University, Tempe, AZ, USA
- Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Gabor Nagy
- Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
| | - John Kunkel
- Center for Biological Physics, Arizona State University, Tempe, AZ, USA
- Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Gerdenis Kodis
- Center for Biological Physics, Arizona State University, Tempe, AZ, USA
- Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, AZ 85212, USA
| | | | - Sonia Longhi
- Aix Marseille Univ, CNRS, AFMB, UMR 7257, Marseille, France
| | - Helmut Grubmüller
- Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
| | - Andrea C Vaiana
- Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
- Present address: Nature's Toolbox, Inc. (NTx), Rio Rancho, NM 87144, USA
| | - Sara M Vaiana
- Center for Biological Physics, Arizona State University, Tempe, AZ, USA
- Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| |
Collapse
|
10
|
Nüesch MF, Pietrek L, Holmstrom ED, Nettels D, von Roten V, Kronenberg-Tenga R, Medalia O, Hummer G, Schuler B. Nanosecond chain dynamics of single-stranded nucleic acids. Nat Commun 2024; 15:6010. [PMID: 39019880 PMCID: PMC11255343 DOI: 10.1038/s41467-024-50092-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Accepted: 07/01/2024] [Indexed: 07/19/2024] Open
Abstract
The conformational dynamics of single-stranded nucleic acids are fundamental for nucleic acid folding and function. However, their elementary chain dynamics have been difficult to resolve experimentally. Here we employ a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, and nanophotonic enhancement to determine the conformational ensembles and rapid chain dynamics of short single-stranded nucleic acids in solution. To interpret the experimental results in terms of end-to-end distance dynamics, we utilize the hierarchical chain growth approach, simple polymer models, and refinement with Bayesian inference to generate structural ensembles that closely align with the experimental data. The resulting chain reconfiguration times are exceedingly rapid, in the 10-ns range. Solvent viscosity-dependent measurements indicate that these dynamics of single-stranded nucleic acids exhibit negligible internal friction and are thus dominated by solvent friction. Our results provide a detailed view of the conformational distributions and rapid dynamics of single-stranded nucleic acids.
Collapse
Affiliation(s)
- Mark F Nüesch
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Lisa Pietrek
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438, Frankfurt am Main, Germany
| | - Erik D Holmstrom
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland.
- Department of Chemistry, University of Kansas, Lawrence, KS, USA.
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA.
| | - Daniel Nettels
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Valentin von Roten
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Rafael Kronenberg-Tenga
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Ohad Medalia
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438, Frankfurt am Main, Germany.
- Institute for Biophysics, Goethe University Frankfurt, 60438, Frankfurt am Main, Germany.
| | - Benjamin Schuler
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland.
- Department of Physics, University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland.
| |
Collapse
|
11
|
Xu Q, Yang M, Ji J, Weng J, Wang W, Xu X. Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis. J Chem Inf Model 2024; 64:5219-5231. [PMID: 38916177 DOI: 10.1021/acs.jcim.3c01833] [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/26/2024]
Abstract
Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.
Collapse
Affiliation(s)
- Qianjun Xu
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Maohua Yang
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Jie Ji
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Jingwei Weng
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Wenning Wang
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Xin Xu
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| |
Collapse
|
12
|
Boushehri S, Holey H, Brosz M, Gumbsch P, Pastewka L, Aponte-Santamaría C, Gräter F. O-glycans Expand Lubricin and Attenuate Its Viscosity and Shear Thinning. Biomacromolecules 2024; 25:3893-3908. [PMID: 38815979 PMCID: PMC11238335 DOI: 10.1021/acs.biomac.3c01348] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 04/26/2024] [Accepted: 04/26/2024] [Indexed: 06/01/2024]
Abstract
Lubricin, an intrinsically disordered glycoprotein, plays a pivotal role in facilitating smooth movement and ensuring the enduring functionality of synovial joints. The central domain of this protein serves as a source of this excellent lubrication and is characterized by its highly glycosylated, negatively charged, and disordered structure. However, the influence of O-glycans on the viscosity of lubricin remains unclear. In this study, we employ molecular dynamics simulations in the absence and presence of shear, along with continuum simulations, to elucidate the intricate interplay between O-glycans and lubricin and the impact of O-glycans on lubricin's conformational properties and viscosity. We found the presence of O-glycans to induce a more extended conformation in fragments of the disordered region of lubricin. These O-glycans contribute to a reduction in solution viscosity but at the same time weaken shear thinning at high shear rates, compared to nonglycosylated systems with the same density. This effect is attributed to the steric and electrostatic repulsion between the fragments, which prevents their conglomeration and structuring. Our computational study yields a mechanistic mechanism underlying previous experimental observations of lubricin and paves the way to a more rational understanding of its function in the synovial fluid.
Collapse
Affiliation(s)
- Saber Boushehri
- Heidelberg Institute for Theoretical Studies, Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany
- University of Heidelberg, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
- Karlsruhe Institute of Technology (KIT), Karlsruhe 76131, Germany
| | - Hannes Holey
- Karlsruhe Institute of Technology (KIT), Karlsruhe 76131, Germany
- Department of Microsystems Engineering, University of Freiburg, Georges-Köhler-Allee 103, Freiburg 79110, Germany
| | - Matthias Brosz
- Heidelberg Institute for Theoretical Studies, Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany
- University of Heidelberg, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| | - Peter Gumbsch
- Karlsruhe Institute of Technology (KIT), Karlsruhe 76131, Germany
- Fraunhofer IWM, Wöhlerstraße 11, Freiburg 79108, Germany
| | - Lars Pastewka
- Department of Microsystems Engineering, University of Freiburg, Georges-Köhler-Allee 103, Freiburg 79110, Germany
| | - Camilo Aponte-Santamaría
- Heidelberg Institute for Theoretical Studies, Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies, Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany
- University of Heidelberg, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| |
Collapse
|
13
|
Ramachandran V, Potoyan DA. Atomistic insights into the reentrant phase-transitions in polyuracil and polylysine mixtures. J Chem Phys 2024; 161:015101. [PMID: 38949285 DOI: 10.1063/5.0206190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Accepted: 06/10/2024] [Indexed: 07/02/2024] Open
Abstract
The phase separation of protein and RNA mixtures underpins the assembly and regulation of numerous membraneless organelles in cells. The ubiquity of protein-RNA condensates in cellular regulatory processes is in part due to their sensitivity to RNA concentration, which affects their physical properties and stability. Recent experiments with poly-cationic peptide-RNA mixtures have revealed closed-loop phase diagrams featuring lower and upper critical solution temperatures. These diagrams indicate reentrant phase transitions shaped by biomolecular interactions and entropic forces such as solvent and ion reorganization. We employed atomistic simulations to study mixtures with various RNA-polylysine stoichiometries and temperatures to elucidate the microscopic driving forces behind reentrant phase transitions in protein-RNA mixtures. Our findings reveal an intricate interplay between hydration, ion condensation, and specific RNA-polylysine hydrogen bonding, resulting in distinct stoichiometry-dependent phase equilibria governing stabilities and structures of the condensate phase. Our simulations show that reentrant transitions are accompanied by desolvation around the phosphate groups of RNA, with increased contacts between phosphate and lysine side chains. In RNA-rich systems at lower temperatures, RNA molecules can form an extensive pi-stacking and hydrogen bond network, leading to percolation. In protein-rich systems, no such percolation-induced transitions are observed. Furthermore, we assessed the performance of three prominent water force fields-Optimal Point Charge (OPC), TIP4P-2005, and TIP4P-D-in capturing reentrant phase transitions. OPC provided a superior balance of interactions, enabling effective capture of reentrant transitions and accurate characterization of changes in solvent reorganization. This study offers atomistic insights into the nature of reentrant phase transitions using simple model peptide and nucleotide mixtures. We believe that our results are broadly applicable to larger classes of peptide-RNA mixtures exhibiting reentrant phase transitions.
Collapse
Affiliation(s)
| | - Davit A Potoyan
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
14
|
Andrews B, Schweitzer-Stenner R, Urbanc B. Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data. J Phys Chem B 2024; 128:6217-6231. [PMID: 38877893 PMCID: PMC11215781 DOI: 10.1021/acs.jpcb.4c02278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/29/2024] [Accepted: 05/30/2024] [Indexed: 06/28/2024]
Abstract
Molecular dynamics (MD) is a great tool for elucidating conformational dynamics of proteins and peptides in water at the atomistic level that often surpasses the level of detail available experimentally. Structure predictions, however, are limited by the accuracy of the underlying MD force field. This limitation is particularly stark in the case of intrinsically disordered peptides and proteins, which are characterized by solvent-accessible and disordered peptide regions and domains. Recent studies show that most additive MD force fields, including CHARMM36m, do not reproduce the intrinsic conformational distributions of guest amino acid residues x in cationic GxG peptides in water in line with experimental data. Positing that a lack of polarizability in additive MD force fields may be the culprit for the reported discrepancies, we here examine the conformational dynamics of guest glycine and alanine residues in cationic GxG peptides in water using two polarizable MD force fields, CHARMM Drude and AMOEBA. Our results indicate that while AMOEBA captures the experimental data better than CHARMM Drude, neither of the two polarizable force fields offers an improvement of the Ramachandran distributions of glycine and alanine residues in cationic GGG and GAG peptides, respectively, over CHARMM36m.
Collapse
Affiliation(s)
- Brian Andrews
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | | | - Brigita Urbanc
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
15
|
Bickel D, Vranken W. Effects of Phosphorylation on Protein Backbone Dynamics and Conformational Preferences. J Chem Theory Comput 2024; 20:4998-5011. [PMID: 38830621 PMCID: PMC11210476 DOI: 10.1021/acs.jctc.4c00206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 05/22/2024] [Accepted: 05/22/2024] [Indexed: 06/05/2024]
Abstract
Phosphorylations are the most common and extensively studied post-translational modification (PTM) of proteins in eukaryotes. They constitute a major regulatory mechanism, modulating protein function, protein-protein interactions, as well as subcellular localization. Phosphorylation sites are preferably located in intrinsically disordered regions and have been shown to trigger structural rearrangements and order-to-disorder transitions. They can therefore have a significant effect on protein backbone dynamics or conformation, but only sparse experimental data are available. To obtain a more general description of how and when phosphorylations have a significant effect on protein behavior, molecular dynamics (MD) currently provides the only suitable framework to study these effects at a large scale in atomistic detail. This study develops a systematic MD simulation framework to explore the influence of phosphorylations on the local backbone dynamics and conformational propensities of proteins. Through a series of glycine-backbone peptides, we studied the effects of amino acid residues including the three most common phosphorylations (Ser, Thr, and Tyr), on local backbone dynamics and conformational propensities. We further extended our study to investigate the interactions of all such residues between position i to positions i + 1, i + 2, i + 3, and i + 4 in such peptides. The final data set comprises structural ensembles for 3393 sequences with more than 1 μs of sampling for each ensemble. To validate the relevance of the results, the structural and conformational properties extracted from the MD simulations are compared to NMR data from the Biological Magnetic Resonance Data Bank. The systematic nature of this study enables the projection of the gained knowledge onto any phosphorylation site in the proteome and provides a general framework for the study of further PTMs. The full data set is publicly available, as a training and reference set.
Collapse
Affiliation(s)
- David Bickel
- Interuniversity
Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
- Structural
Biology Brussels, Vrije Universiteit Brussels, 1050 Brussels, Belgium
| | - Wim Vranken
- Interuniversity
Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
- Structural
Biology Brussels, Vrije Universiteit Brussels, 1050 Brussels, Belgium
| |
Collapse
|
16
|
Duran C, Casadevall G, Osuna S. Harnessing conformational dynamics in enzyme catalysis to achieve nature-like catalytic efficiencies: the shortest path map tool for computational enzyme redesign. Faraday Discuss 2024. [PMID: 38910409 DOI: 10.1039/d3fd00156c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
Enzymes exhibit diverse conformations, as represented in the free energy landscape (FEL). Such conformational diversity provides enzymes with the ability to evolve towards novel functions. The challenge lies in identifying mutations that enhance specific conformational changes, especially if located in distal sites from the active site cavity. The shortest path map (SPM) method, which we developed to address this challenge, constructs a graph based on the distances and correlated motions of residues observed in nanosecond timescale molecular dynamics (MD) simulations. We recently introduced a template based AlphaFold2 (tAF2) approach coupled with 10 nanosecond MD simulations to quickly estimate the conformational landscape of enzymes and assess how the FEL is shifted after mutation. In this study, we evaluate the potential of SPM when coupled with tAF2-MD in estimating conformational heterogeneity and identifying key conformationally-relevant positions. The selected model system is the beta subunit of tryptophan synthase (TrpB). We compare how the SPM pathways differ when integrating tAF2 with different MD simulation lengths from as short as 10 ns until 50 ns and considering two distinct Amber forcefield and water models (ff14SB/TIP3P versus ff19SB/OPC). The new methodology can more effectively capture the distal mutations found in laboratory evolution, thus showcasing the efficacy of tAF2-MD-SPM in rapidly estimating enzyme dynamics and identifying the key conformationally relevant hotspots for computational enzyme engineering.
Collapse
Affiliation(s)
- Cristina Duran
- Departament de Química, Institut de Química Computacional i Catàlisi, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003, Girona, Spain.
| | - Guillem Casadevall
- Departament de Química, Institut de Química Computacional i Catàlisi, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003, Girona, Spain.
| | - Sílvia Osuna
- Departament de Química, Institut de Química Computacional i Catàlisi, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003, Girona, Spain.
- ICREA, Pg. Lluís Companys 23, 08010, Barcelona, Spain
| |
Collapse
|
17
|
Smardz P, Anila MM, Rogowski P, Li MS, Różycki B, Krupa P. A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta. Int J Mol Sci 2024; 25:6698. [PMID: 38928405 PMCID: PMC11204378 DOI: 10.3390/ijms25126698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2024] [Revised: 06/12/2024] [Accepted: 06/13/2024] [Indexed: 06/28/2024] Open
Abstract
Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-β (Aβ42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the β-content of Aβ42; however, it increases the tendency of the monomeric Aβ42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aβ42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aβ42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.
Collapse
Affiliation(s)
| | | | | | | | | | - Pawel Krupa
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland; (P.S.); (M.M.A.); (P.R.); (M.S.L.); (B.R.)
| |
Collapse
|
18
|
Wan K, He J, Shi X. Construction of High Accuracy Machine Learning Interatomic Potential for Surface/Interface of Nanomaterials-A Review. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2305758. [PMID: 37640376 DOI: 10.1002/adma.202305758] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Revised: 08/24/2023] [Indexed: 08/31/2023]
Abstract
The inherent discontinuity and unique dimensional attributes of nanomaterial surfaces and interfaces bestow them with various exceptional properties. These properties, however, also introduce difficulties for both experimental and computational studies. The advent of machine learning interatomic potential (MLIP) addresses some of the limitations associated with empirical force fields, presenting a valuable avenue for accurate simulations of these surfaces/interfaces of nanomaterials. Central to this approach is the idea of capturing the relationship between system configuration and potential energy, leveraging the proficiency of machine learning (ML) to precisely approximate high-dimensional functions. This review offers an in-depth examination of MLIP principles and their execution and elaborates on their applications in the realm of nanomaterial surface and interface systems. The prevailing challenges faced by this potent methodology are also discussed.
Collapse
Affiliation(s)
- Kaiwei Wan
- Laboratory of Theoretical and Computational Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing, 100190, China
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, 100049, China
| | - Jianxin He
- Laboratory of Theoretical and Computational Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing, 100190, China
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, 100049, China
| | - Xinghua Shi
- Laboratory of Theoretical and Computational Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing, 100190, China
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, 100049, China
| |
Collapse
|
19
|
Borsatto A, Gianquinto E, Rizzi V, Gervasio FL. SWISH-X, an Expanded Approach to Detect Cryptic Pockets in Proteins and at Protein-Protein Interfaces. J Chem Theory Comput 2024; 20:3335-3348. [PMID: 38563746 PMCID: PMC11044271 DOI: 10.1021/acs.jctc.3c01318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Revised: 03/22/2024] [Accepted: 03/23/2024] [Indexed: 04/04/2024]
Abstract
Protein-protein interactions mediate most molecular processes in the cell, offering a significant opportunity to expand the set of known druggable targets. Unfortunately, targeting these interactions can be challenging due to their typically flat and featureless interaction surfaces, which often change as the complex forms. Such surface changes may reveal hidden (cryptic) druggable pockets. Here, we analyze a set of well-characterized protein-protein interactions harboring cryptic pockets and investigate the predictive power of current computational methods. Based on our observations, we developed a new computational strategy, SWISH-X (SWISH Expanded), which combines the established cryptic pocket identification capabilities of SWISH with the rapid temperature range exploration of OPES MultiThermal. SWISH-X is able to reliably identify cryptic pockets at protein-protein interfaces while retaining its predictive power for revealing cryptic pockets in isolated proteins, such as TEM-1 β-lactamase.
Collapse
Affiliation(s)
- Alberto Borsatto
- School
of Pharmaceutical Sciences, University of
Geneva, 1205 Geneva, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland, University of Geneva, 1205 Geneva, Switzerland
- Swiss
Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Eleonora Gianquinto
- Department
of Drug Science and Technology, University
of Turin, 10125 Turin, Italy
| | - Valerio Rizzi
- School
of Pharmaceutical Sciences, University of
Geneva, 1205 Geneva, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland, University of Geneva, 1205 Geneva, Switzerland
- Swiss
Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Francesco Luigi Gervasio
- School
of Pharmaceutical Sciences, University of
Geneva, 1205 Geneva, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland, University of Geneva, 1205 Geneva, Switzerland
- Swiss
Institute of Bioinformatics, 1015 Lausanne, Switzerland
- Department
of Chemistry, University College London, WC1 H0AJ London, United Kingdom
- Institute
of Structural and Molecular Biology, University
College London, WC1E7JE London, United Kingdom
| |
Collapse
|
20
|
Alavizargar A, Gass M, Krahn MP, Heuer A. Elucidating the Membrane Binding Process of a Disordered Protein: Dynamic Interplay of Anionic Lipids and the Polybasic Region. ACS PHYSICAL CHEMISTRY AU 2024; 4:167-179. [PMID: 38560754 PMCID: PMC10979486 DOI: 10.1021/acsphyschemau.3c00051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/06/2023] [Accepted: 12/06/2023] [Indexed: 04/04/2024]
Abstract
Intrinsically disordered regions of proteins are responsible for many biological processes such as in the case of liver kinase B1 (LKB1)-a serine/threonine kinase relevant for cell proliferation and cell polarity. LKB1 becomes fully activated upon recruitment to the plasma membrane by binding of its disordered C-terminal polybasic motif consisting of eight lysines/arginines to phospholipids. Here, we present extensive molecular dynamics (MD) simulations of the polybasic motif interacting with a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleyl phosphatidic acid (PA) and cell culture experiments. Protein-membrane binding effects are due to the electrostatic interactions between the polybasic amino acids and PAs. For significant binding, the first three lysines turn out to be dispensable, which was also recapitulated in cell culture using transfected GFP-LKB1 variants. LKB1-membrane binding results in nonmonotonous changes in the structure of the protein as well as the membrane, in particular, accumulation of PAs and reduced thickness at the protein-membrane contact area. The protein-lipid binding turns out to be highly dynamic due to an interplay of PA-PA repulsion and protein-PA attraction. The thermodynamics of this interplay is captured by a statistical fluctuation model, which allows the estimation of both energies. Quantification of the significance of each polar amino acid in the polybasic provides detailed insights into the molecular mechanism of protein-membrane binding of LKB1. These results can likely be transferred to other proteins, which interact by intrinsically disordered polybasic regions with anionic membranes.
Collapse
Affiliation(s)
- Azadeh Alavizargar
- Institute
of Physical Chemistry, University of Münster, Corrensstr. 28/30, 48149 Münster, Germany
| | - Maximilian Gass
- Medical
Cell Biology, Medical Clinic D, University
Hospital of Münster, Albert-Schweitzer Campus 1-A14, 48149 Münster, Germany
| | - Michael P. Krahn
- Medical
Cell Biology, Medical Clinic D, University
Hospital of Münster, Albert-Schweitzer Campus 1-A14, 48149 Münster, Germany
| | - Andreas Heuer
- Institute
of Physical Chemistry, University of Münster, Corrensstr. 28/30, 48149 Münster, Germany
| |
Collapse
|
21
|
Hsueh SCC, Nijland M, Aina A, Plotkin SS. Cyclization Scaffolding for Improved Vaccine Immunogen Stability: Application to Tau Protein in Alzheimer's Disease. J Chem Inf Model 2024; 64:2035-2044. [PMID: 38427576 DOI: 10.1021/acs.jcim.3c01556] [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: 03/03/2024]
Abstract
Effective scaffolding of immunogens is crucial for generating conformationally selective antibodies through active immunization, particularly in the treatment of protein misfolding diseases such as Alzheimer's and Parkinson's disease. Previous computational work has revealed that a disorder-prone region of the tau protein, when in a stacked form, is predicted to structurally resemble a small, soluble protofibril, having conformational properties similar to those of experimental in vitro tau oligomers. Such an oligomeric structural mimic has the potential to serve as a vaccine immunogen design for Alzheimer's disease. In this study, we developed a cyclization scaffolding method in Rosetta, in which multiple cyclic peptides are stacked into a protofibril. Cyclization results in significant stabilization of protofibril-like structures by constraining the conformational space. Applying this method to the disorder-prone region of the tau fibril, we evaluated the metastability of the cyclized tau immunogen using molecular dynamics simulations, and we identified sequences of two cyclic constructs having high metastability in the protofibril. We then assessed their thermodynamic stability by computing the free energy required to separate a distal chain from the rest of the stacked structure. Our computational results, based on molecular dynamics simulations and free energy calculations, demonstrate that two cyclized constructs, cyclo-(VKSEKLDFKDRVQSKIFyN) and cyclo-(VKSEKLDFKDRVQSKIYvG) (lowercase letters indicate d-form amino acids), possess significantly increased thermodynamic stability in the protofibril over an uncyclized linear construct VKSEKLDFKDRVQSKI. The cyclization scaffolding approach proposed here holds promise as a means to effectively design immunogens for protein misfolding diseases, particularly those involving liposome-conjugated peptide constructs.
Collapse
Affiliation(s)
- Shawn C C Hsueh
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - Mark Nijland
- Laboratory of Physical Chemistry, Wageningen University, Wageningen 6708 WG, The Netherlands
| | - Adekunle Aina
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - Steven S Plotkin
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
- Genome Science and Technology Program, The University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
22
|
Montero de Hijes P, Dellago C, Jinnouchi R, Schmiedmayer B, Kresse G. Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks. J Chem Phys 2024; 160:114107. [PMID: 38506284 DOI: 10.1063/5.0197105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.
Collapse
Affiliation(s)
- Pablo Montero de Hijes
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- University of Vienna, Faculty of Earth Sciences, Geography and Astronomy, Josef-Holaubuek-Platz 2, 1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
| | - Ryosuke Jinnouchi
- Toyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | | | - Georg Kresse
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- VASP Software GmbH, Berggasse 21, A-1090 Vienna, Austria
| |
Collapse
|
23
|
Pietrek LM, Stelzl LS, Hummer G. Hierarchical Assembly of Single-Stranded RNA. J Chem Theory Comput 2024; 20:2246-2260. [PMID: 38361440 PMCID: PMC10938505 DOI: 10.1021/acs.jctc.3c01049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 12/09/2023] [Accepted: 01/25/2024] [Indexed: 02/17/2024]
Abstract
Single-stranded RNA (ssRNA) plays a major role in the flow of genetic information-most notably, in the form of messenger RNA (mRNA)-and in the regulation of biological processes. The highly dynamic nature of chains of unpaired nucleobases challenges structural characterizations of ssRNA by experiments or molecular dynamics (MD) simulations alike. Here, we use hierarchical chain growth (HCG) to construct ensembles of ssRNA chains. HCG assembles the structures of protein and nucleic acid chains from fragment libraries created by MD simulations. Applied to homo- and heteropolymeric ssRNAs of different lengths, we find that HCG produces structural ensembles that overall are in good agreement with diverse experiments, including nuclear magnetic resonance (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET). The agreement can be further improved by ensemble refinement using Bayesian inference of ensembles (BioEn). HCG can also be used to assemble RNA structures that combine base-paired and base-unpaired regions, as illustrated for the 5' untranslated region (UTR) of SARS-CoV-2 RNA.
Collapse
Affiliation(s)
- Lisa M. Pietrek
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Lukas S. Stelzl
- Faculty
of Biology, Johannes Gutenberg University
Mainz, Gresemundweg 2, 55128 Mainz, Germany
- KOMET
1, Institute of Physics, Johannes Gutenberg
University Mainz, 55099 Mainz, Germany
- Institute
of Molecular Biology (IMB), 55128 Mainz, Germany
| | - Gerhard Hummer
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
- Institute
for Biophysics, Goethe University, Max-von-Laue-Straße 9, 60438 Frankfurt am Main, Germany
| |
Collapse
|
24
|
Li S, Wu B, Luo YL, Han W. Simulations of Functional Motions of Super Large Biomolecules with a Mixed-Resolution Model. J Chem Theory Comput 2024; 20:2228-2245. [PMID: 38374639 PMCID: PMC10938502 DOI: 10.1021/acs.jctc.3c01046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 01/18/2024] [Accepted: 01/29/2024] [Indexed: 02/21/2024]
Abstract
Many large protein machines function through an interplay between large-scale movements and intricate conformational changes. Understanding functional motions of these proteins through simulations becomes challenging for both all-atom and coarse-grained (CG) modeling techniques because neither approach alone can readily capture the full details of these motions. In this study, we develop a multiscale model by employing the popular MARTINI CG model to represent a heterogeneous environment and structurally stable proteins and using the united-atom (UA) model PACE to describe proteins undergoing subtle conformational changes. PACE was previously developed to be compatible with the MARTINI solvent and membrane. Here, we couple the protein descriptions of the two models by directly mixing UA and CG interaction parameters to greatly simplify parameter determination. Through extensive validations with diverse protein systems in solution or membrane, we demonstrate that only additional parameter rescaling is needed to enable the resulting model to recover the stability of native structures of proteins under mixed representation. Moreover, we identify the optimal scaling factors that can be applied to various protein systems, rendering the model potentially transferable. To further demonstrate its applicability for realistic systems, we apply the model to a mechanosensitive ion channel Piezo1 that has peripheral arms for sensing membrane tension and a central pore for ion conductance. The model can reproduce the coupling between Piezo1's large-scale arm movement and subtle pore opening in response to membrane stress while consuming much less computational costs than all-atom models. Therefore, our model shows promise for studying functional motions of large protein machines.
Collapse
Affiliation(s)
- Shu Li
- Centre
for Artificial Intelligence Driven Drug Discovery, Faculty of Applied
Sciences, Macao Polytechnic University, Macao 999078, China
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Bohua Wu
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun Lyna Luo
- Department
of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Wei Han
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Department
of Chemistry, Faculty of Science, Hong Kong
Baptist University, Hong Kong SAR 999077, China
- Shenzhen
Bay Laboratory, Institute of Chemical Biology, Shenzhen 518132, China
| |
Collapse
|
25
|
Dohnalová H, Matoušková E, Lankaš F. Temperature-dependent elasticity of DNA, RNA, and hybrid double helices. Biophys J 2024; 123:572-583. [PMID: 38340722 PMCID: PMC10938081 DOI: 10.1016/j.bpj.2024.01.032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Revised: 12/19/2023] [Accepted: 01/26/2024] [Indexed: 02/12/2024] Open
Abstract
Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.
Collapse
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Praha 6, Czech Republic
| | - Eva Matoušková
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Praha 6, Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Praha 6, Czech Republic.
| |
Collapse
|
26
|
Holehouse AS, Kragelund BB. The molecular basis for cellular function of intrinsically disordered protein regions. Nat Rev Mol Cell Biol 2024; 25:187-211. [PMID: 37957331 DOI: 10.1038/s41580-023-00673-0] [Citation(s) in RCA: 48] [Impact Index Per Article: 48.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/26/2023] [Indexed: 11/15/2023]
Abstract
Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.
Collapse
Affiliation(s)
- Alex S Holehouse
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St Louis, MO, USA.
- Center for Biomolecular Condensates, Washington University in St Louis, St Louis, MO, USA.
| | - Birthe B Kragelund
- REPIN, Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
27
|
Latorraca NR, Sabaat S, Habrian C, Bleier J, Stanley C, Marqusee S, Isacoff EY. Domain coupling in activation of a family C GPCR. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.28.582567. [PMID: 38464305 PMCID: PMC10925283 DOI: 10.1101/2024.02.28.582567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
The G protein-coupled metabotropic glutamate receptors form homodimers and heterodimers with highly diverse responses to glutamate and varying physiological function. The molecular basis for this diversity remains poorly delineated. We employ molecular dynamics, single-molecule spectroscopy, and hydrogen-deuterium exchange to dissect the pathway of activation triggered by glutamate. We find that activation entails multiple loosely coupled steps and identify a novel pre-active intermediate whose transition to the active state forms dimer interactions that set signaling efficacy. Such subunit interactions generate functional diversity that differs across homodimers and heterodimers. The agonist-bound receptor is remarkably dynamic, with low occupancy of G protein-coupling conformations, providing considerable headroom for modulation of the landscape by allosteric ligands. Sites of sequence diversity within the dimerization interface and diverse coupling between activation rearrangements may contribute to precise decoding of glutamate signals and transients over broad spatial and temporal scales.
Collapse
Affiliation(s)
- Naomi R. Latorraca
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
| | - Sam Sabaat
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
| | - Chris Habrian
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
| | - Julia Bleier
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
| | - Cherise Stanley
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
| | - Susan Marqusee
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
- Department of Chemistry, University of California, Berkeley, California, 94720 USA
- Chan Zuckerberg Biohub, San Francisco, CA 94158
- California Institute for Quantitative Biosciences, University of California, Berkeley, California, 94720 USA
| | - Ehud Y. Isacoff
- Department of Molecular and Cell Biology, University of California, Berkeley, California, 94720, USA
- California Institute for Quantitative Biosciences, University of California, Berkeley, California, 94720 USA
- Molecular Biology & Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, USA
- Weill Neurohub, University of California, Berkeley, California, 94720 USA
- Molecular Biology & Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, USA
| |
Collapse
|
28
|
Sisk TR, Robustelli P. Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model. Proc Natl Acad Sci U S A 2024; 121:e2313360121. [PMID: 38294935 PMCID: PMC10861926 DOI: 10.1073/pnas.2313360121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 11/22/2023] [Indexed: 02/02/2024] Open
Abstract
A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.
Collapse
Affiliation(s)
- Thomas R. Sisk
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| | - Paul Robustelli
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| |
Collapse
|
29
|
Galano‐Frutos JJ, Sancho J. Energy, water, and protein folding: A molecular dynamics-based quantitative inventory of molecular interactions and forces that make proteins stable. Protein Sci 2024; 33:e4905. [PMID: 38284492 PMCID: PMC10804899 DOI: 10.1002/pro.4905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 12/12/2023] [Accepted: 01/05/2024] [Indexed: 01/30/2024]
Abstract
Protein folding energetics can be determined experimentally on a case-by-case basis but it is not understood in sufficient detail to provide deep control in protein design. The fundamentals of protein stability have been outlined by calorimetry, protein engineering, and biophysical modeling, but these approaches still face great difficulty in elucidating the specific contributions of the intervening molecules and physical interactions. Recently, we have shown that the enthalpy and heat capacity changes associated to the protein folding reaction can be calculated within experimental error using molecular dynamics simulations of native protein structures and their corresponding unfolded ensembles. Analyzing in depth molecular dynamics simulations of four model proteins (CI2, barnase, SNase, and apoflavodoxin), we dissect here the energy contributions to ΔH (a key component of protein stability) made by the molecular players (polypeptide and solvent molecules) and physical interactions (electrostatic, van der Waals, and bonded) involved. Although the proteins analyzed differ in length, isoelectric point and fold class, their folding energetics is governed by the same quantitative pattern. Relative to the unfolded ensemble, the native conformations are enthalpically stabilized by comparable contributions from protein-protein and solvent-solvent interactions, and almost equally destabilized by interactions between protein and solvent molecules. The native protein surface seems to interact better with water than the unfolded one, but this is outweighed by the unfolded surface being larger. From the perspective of physical interactions, the native conformations are stabilized by van de Waals and Coulomb interactions and destabilized by conformational strain arising from bonded interactions. Also common to the four proteins, the sign of the heat capacity change is set by interactions between protein and solvent molecules or, from the alternative perspective, by Coulomb interactions.
Collapse
Affiliation(s)
- Juan José Galano‐Frutos
- Biocomputation and Complex Systems Physics Institute (BIFI)‐Joint Unit GBsC‐CSICUniversity of ZaragozaZaragozaSpain
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de CienciasUniversity of ZaragozaZaragozaSpain
| | - Javier Sancho
- Biocomputation and Complex Systems Physics Institute (BIFI)‐Joint Unit GBsC‐CSICUniversity of ZaragozaZaragozaSpain
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de CienciasUniversity of ZaragozaZaragozaSpain
- Aragon Health Research Institute (IIS Aragón)ZaragozaSpain
| |
Collapse
|
30
|
Sreenivasan S, Heffren P, Suh K, Rodnin MV, Kosa E, Fenton AW, Ladokhin AS, Smith PE, Fontes JD, Swint‐Kruse L. The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm? Protein Sci 2024; 33:e4863. [PMID: 38073129 PMCID: PMC10806935 DOI: 10.1002/pro.4863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 12/04/2023] [Accepted: 12/07/2023] [Indexed: 01/27/2024]
Abstract
During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.
Collapse
Affiliation(s)
- Shwetha Sreenivasan
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Paul Heffren
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
- Present address:
Department of BiosciencesKansas City UniversityKansas CityMissouriUSA
| | - Kyung‐Shin Suh
- Department of ChemistryKansas State UniversityManhattanKansasUSA
| | - Mykola V. Rodnin
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Edina Kosa
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Aron W. Fenton
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Alexey S. Ladokhin
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Paul E. Smith
- Department of ChemistryKansas State UniversityManhattanKansasUSA
| | - Joseph D. Fontes
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| | - Liskin Swint‐Kruse
- Department of Biochemistry and Molecular BiologyUniversity of Kansas Medical CenterKansas CityKansasUSA
| |
Collapse
|
31
|
Dohnalová H, Seifert M, Matoušková E, Klein M, Papini FS, Lipfert J, Dulin D, Lankaš F. Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations. J Phys Chem B 2024; 128:664-675. [PMID: 38197365 PMCID: PMC10823466 DOI: 10.1021/acs.jpcb.3c06280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Revised: 12/19/2023] [Accepted: 12/20/2023] [Indexed: 01/11/2024]
Abstract
RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.
Collapse
Affiliation(s)
- Hana Dohnalová
- Department
of Informatics and Chemistry, University
of Chemistry and Technology Prague, Technická 5, 166 28 Praha
6, Czech Republic
| | - Mona Seifert
- Junior
Research Group 2, Interdisciplinary Center for Clinical Research, Friedrich-Alexander-University Erlangen-Nürnberg, Cauerstr. 3, Erlangen 91058, Germany
| | - Eva Matoušková
- Department
of Informatics and Chemistry, University
of Chemistry and Technology Prague, Technická 5, 166 28 Praha
6, Czech Republic
| | - Misha Klein
- Department
of Physics and Astronomy and LaserLaB Amsterdam, Vrije Universiteit Amsterdam, De Boelelaan 1081, Amsterdam 1081 HV, The Netherlands
| | - Flávia S. Papini
- Junior
Research Group 2, Interdisciplinary Center for Clinical Research, Friedrich-Alexander-University Erlangen-Nürnberg, Cauerstr. 3, Erlangen 91058, Germany
| | - Jan Lipfert
- Soft
Condensed Matter and Biophysics, Department of Physics and Debye Institute, Utrecht University, Utrecht 3584 CC, The Netherlands
| | - David Dulin
- Junior
Research Group 2, Interdisciplinary Center for Clinical Research, Friedrich-Alexander-University Erlangen-Nürnberg, Cauerstr. 3, Erlangen 91058, Germany
- Department
of Physics and Astronomy and LaserLaB Amsterdam, Vrije Universiteit Amsterdam, De Boelelaan 1081, Amsterdam 1081 HV, The Netherlands
| | - Filip Lankaš
- Department
of Informatics and Chemistry, University
of Chemistry and Technology Prague, Technická 5, 166 28 Praha
6, Czech Republic
| |
Collapse
|
32
|
Singh Y, Hocky GM. Improved Prediction of Molecular Response to Pulling by Combining Force Tempering with Replica Exchange Methods. J Phys Chem B 2024; 128:706-715. [PMID: 38230998 PMCID: PMC10823473 DOI: 10.1021/acs.jpcb.3c07081] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 01/02/2024] [Accepted: 01/03/2024] [Indexed: 01/18/2024]
Abstract
Small mechanical forces play important functional roles in many crucial cellular processes, including in the dynamic behavior of the cytoskeleton and in the regulation of osmotic pressure through membrane-bound proteins. Molecular simulations offer the promise of being able to design the behavior of proteins that sense and respond to these forces. However, it is difficult to predict and identify the effect of the relevant piconewton (pN) scale forces due to their small magnitude. Previously, we introduced the Infinite Switch Simulated Tempering in Force (FISST) method, which allows one to estimate the effect of a range of applied forces from a single molecular dynamics simulation, and also demonstrated that FISST additionally accelerates sampling of a molecule's conformational landscape. For some problems, we find that this acceleration is not sufficient to capture all relevant conformational fluctuations, and hence, here we demonstrate that FISST can be combined with either temperature replica exchange or solute tempering approaches to produce a hybrid method that enables more robust prediction of the effect of small forces on molecular systems.
Collapse
Affiliation(s)
- Yuvraj Singh
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Glen M. Hocky
- Department
of Chemistry, New York University, New York, New York 10003, United States
- Simons
Center for Computational Physical Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
33
|
Kota D, Prasad R, Zhou HX. Adenosine Triphosphate Mediates Phase Separation of Disordered Basic Proteins by Bridging Intermolecular Interaction Networks. J Am Chem Soc 2024; 146:1326-1336. [PMID: 38174879 PMCID: PMC10843746 DOI: 10.1021/jacs.3c09134] [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: 01/05/2024]
Abstract
Adenosine triphosphate (ATP) is an abundant molecule with crucial cellular roles as the energy currency and a building block of nucleic acids and for protein phosphorylation. Here we show that ATP mediates the phase separation of basic intrinsically disordered proteins (bIDPs). In the resulting condensates, ATP is highly concentrated (apparent partition coefficients up to 7700) and serves as bridges between bIDP chains. These liquid-like droplets have some of the lowest interfacial tension (∼25 pN/μm) but high zero-shear viscosities (1-15 Pa s) due to the bridged protein networks, and yet their fusion has some of the highest speeds (∼1 μm/ms). The rapid fusion manifests extreme shear thinning, where the apparent viscosity is lower than zero-shear viscosity by over 100-fold, made possible by fast reformation of the ATP bridges. At still higher concentrations, ATP does not dissolve bIDP droplets but results in aggregates and fibrils.
Collapse
Affiliation(s)
- Divya Kota
- Department of Chemistry, University of Illinois Chicago, Chicago IL 60607, USA
| | - Ramesh Prasad
- Department of Chemistry, University of Illinois Chicago, Chicago IL 60607, USA
| | - Huan-Xiang Zhou
- Department of Chemistry, University of Illinois Chicago, Chicago IL 60607, USA
- Department of Physics, University of Illinois Chicago, Chicago IL 60607, USA
| |
Collapse
|
34
|
Bass L, Elder LH, Folescu DE, Forouzesh N, Tolokh IS, Karpatne A, Onufriev AV. Improving the Accuracy of Physics-Based Hydration-Free Energy Predictions by Machine Learning the Remaining Error Relative to the Experiment. J Chem Theory Comput 2024; 20:396-410. [PMID: 38149593 PMCID: PMC10950260 DOI: 10.1021/acs.jctc.3c00981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
The accuracy of computational models of water is key to atomistic simulations of biomolecules. We propose a computationally efficient way to improve the accuracy of the prediction of hydration-free energies (HFEs) of small molecules: the remaining errors of the physics-based models relative to the experiment are predicted and mitigated by machine learning (ML) as a postprocessing step. Specifically, the trained graph convolutional neural network attempts to identify the "blind spots" in the physics-based model predictions, where the complex physics of aqueous solvation is poorly accounted for, and partially corrects for them. The strategy is explored for five classical solvent models representing various accuracy/speed trade-offs, from the fast analytical generalized Born (GB) to the popular TIP3P explicit solvent model; experimental HFEs of small neutral molecules from the FreeSolv set are used for the training and testing. For all of the models, the ML correction reduces the resulting root-mean-square error relative to the experiment for HFEs of small molecules, without significant overfitting and with negligible computational overhead. For example, on the test set, the relative accuracy improvement is 47% for the fast analytical GB, making it, after the ML correction, almost as accurate as uncorrected TIP3P. For the TIP3P model, the accuracy improvement is about 39%, bringing the ML-corrected model's accuracy below the 1 kcal/mol threshold. In general, the relative benefit of the ML corrections is smaller for more accurate physics-based models, reaching the lower limit of about 20% relative accuracy gain compared with that of the physics-based treatment alone. The proposed strategy of using ML to learn the remaining error of physics-based models offers a distinct advantage over training ML alone directly on reference HFEs: it preserves the correct overall trend, even well outside of the training set.
Collapse
Affiliation(s)
- Lewis Bass
- Department of Computer Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Luke H Elder
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Dan E Folescu
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Mathematics, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Negin Forouzesh
- Department of Computer Science, California State University, Los Angeles, California 90032, United States
| | - Igor S Tolokh
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Anuj Karpatne
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
35
|
Singh K, Reddy G. Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms. J Chem Theory Comput 2024; 20:421-435. [PMID: 38134376 DOI: 10.1021/acs.jctc.3c00999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2023]
Abstract
Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.
Collapse
Affiliation(s)
- Kushal Singh
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| |
Collapse
|
36
|
Lebedenko OO, Salikov VA, Izmailov SA, Podkorytov IS, Skrynnikov NR. Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4. Biophys J 2024; 123:80-100. [PMID: 37990496 PMCID: PMC10808029 DOI: 10.1016/j.bpj.2023.11.020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Revised: 10/28/2023] [Accepted: 11/17/2023] [Indexed: 11/23/2023] Open
Abstract
MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.
Collapse
Affiliation(s)
- Olga O Lebedenko
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Vladislav A Salikov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Sergei A Izmailov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Ivan S Podkorytov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Nikolai R Skrynnikov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia; Department of Chemistry, Purdue University, West Lafayette, Indiana.
| |
Collapse
|
37
|
Barbhuiya TK, Jayarathna DK, Gilmour R, Smet-Nocca C, Gandhi NS. Characterization of Posttranslationally Modified PHF-1 Tau Peptides Using Gaussian Accelerated Molecular Dynamics Simulation. Methods Mol Biol 2024; 2754:3-31. [PMID: 38512658 DOI: 10.1007/978-1-0716-3629-9_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
The microtubule-associated protein, Tau, is an intrinsically disordered protein that plays a crucial role in neurodegenerative diseases like Alzheimer's disease. The posttranslational modifications across the Tau protein domains are involved in regulating Tau protein's function and disease onset. Of the various posttranslational modifications at Ser, Thr, and Tyr sites, O-GlcNAcylation and phosphorylation are the most critical ones, playing a vital role in Tau aggregation and tauopathies. To understand the function, it is essential to characterize the structural changes associated with Tau modification. Previous experimental studies have focused on high-resolution nuclear magnetic resonance techniques to structurally characterize the effect of phosphorylation, O-GlcNAcylation, and combination of both PTMs on Tau conformation in small peptides centered on the PHF-1 epitope from amino acid 392 to 411. The structural characterization using atomistic molecular dynamics simulation of such disordered peptides requires long simulation time, proper sampling method, and utilization of appropriate force fields for accurate determination of conformational ensembles, resembling the experimental data. This chapter details the protocol for the structural characterization of modified Tau peptides using the CHARMM36m force field and enhanced sampling methods like Gaussian accelerated molecular dynamics (GaMD) simulation. We have focused on a detailed explanation of the GaMD method and analyses of molecular dynamics trajectories to explain the relationship between two modifications, phospho- and glyco-, at C-terminus of Tau protein and its stable conformation over the longer simulation timeframes. The analyses involve energetics reweighting, clustering of simulation trajectories, and characterization of secondary structure using circular dichroism data from the simulation. The reader can utilize this protocol to investigate the structures of complex proteins, especially the disordered ones.
Collapse
Affiliation(s)
- Tabassum Khair Barbhuiya
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
- Cancer and Ageing Research Program, Translational Research Institute, Woolloongabba, QLD, Australia
| | - Dulari K Jayarathna
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
| | - Raechelle Gilmour
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia
| | - Caroline Smet-Nocca
- University of Lille, Inserm, CHU Lille, Institut Pasteur de Lille, U1167 - RID-AGE - Risk Factors and Molecular Determinants of Aging-Related Diseases, Lille, France
- CNRS EMR9002 Integrative Structural Biology, Lille, France
| | - Neha S Gandhi
- School of Chemistry and Physics, Faculty of Science, Centre for Genomics and Personalised Health, Queensland University of Technology, Brisbane, QLD, Australia.
- Cancer and Ageing Research Program, Translational Research Institute, Woolloongabba, QLD, Australia.
| |
Collapse
|
38
|
Sharma B, Mattaparthi VSK. Prediction of interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins from dihedral angles. J Biomol Struct Dyn 2023:1-11. [PMID: 38116756 DOI: 10.1080/07391102.2023.2294837] [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: 04/21/2023] [Accepted: 12/06/2023] [Indexed: 12/21/2023]
Abstract
Intrinsically disordered proteins (IDPs) are proteins that do not form uniquely defined three-dimensional (3-D) structures. Experimental research on IDPs is difficult since they go against the traditional protein structure-function paradigm. Although there are several predictors of disorder based on amino acid sequences, but very limited based on the 3-D structures of proteins. Dihedral angles have a significant role in predicting protein structure because they establish a protein's backbone, which, coupled with its side chain, establishes its overall shape. Here, we have carried out atomistic Molecular Dynamics (MD) simulations on four different proteins: one ordered protein (Monellin), two partially disordered proteins (p53-TAD and Amyloid beta (Aβ1-42) peptide), and one completely disordered protein (Histatin 5). The MD simulation trajectories for the corresponding four proteins were used to conduct dihedral angle (ϕ and ѱ) analysis. Then, the average dihedral angles for each of the residues were calculated and plotted against the residue index. We noticed steep rises or falls in the average ϕ value at certain locations in the plot. These sudden shifts in the average ϕ value reflect the interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins. Using this method, the probable conformer of a protein with a higher degree of disorder can be found among the ensembles of structures sampled during the MD simulations. The results of our study offer new understandings on precisely identifying regions of various degrees of disorder in intrinsically disordered proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Babli Sharma
- Molecular Modelling and Simulation Laboratory, Department of Molecular Biology and Biotechnology, Tezpur University, Assam, India
| | - Venkata Satish Kumar Mattaparthi
- Molecular Modelling and Simulation Laboratory, Department of Molecular Biology and Biotechnology, Tezpur University, Assam, India
| |
Collapse
|
39
|
Yamada T, Miyazaki Y, Harada S, Kumar A, Vanni S, Shinoda W. Improved Protein Model in SPICA Force Field. J Chem Theory Comput 2023; 19:8967-8977. [PMID: 37989551 DOI: 10.1021/acs.jctc.3c01016] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.
Collapse
Affiliation(s)
- Teppei Yamada
- Graduate School of Natural Science and Technology, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Shogo Harada
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Ashutosh Kumar
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
40
|
Linse JB, Hub JS. Scrutinizing the protein hydration shell from molecular dynamics simulations against consensus small-angle scattering data. Commun Chem 2023; 6:272. [PMID: 38086909 PMCID: PMC10716392 DOI: 10.1038/s42004-023-01067-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 11/20/2023] [Indexed: 06/09/2024] Open
Abstract
Biological macromolecules in solution are surrounded by a hydration shell, whose structure differs from the structure of bulk solvent. While the importance of the hydration shell for numerous biological functions is widely acknowledged, it remains unknown how the hydration shell is regulated by macromolecular shape and surface composition, mainly because a quantitative probe of the hydration shell structure has been missing. We show that small-angle scattering in solution using X-rays (SAXS) or neutrons (SANS) provide a protein-specific probe of the protein hydration shell that enables quantitative comparison with molecular simulations. Using explicit-solvent SAXS/SANS predictions, we derived the effect of the hydration shell on the radii of gyration Rg of five proteins using 18 combinations of protein force field and water model. By comparing computed Rg values from SAXS relative to SANS in D2O with consensus SAXS/SANS data from a recent worldwide community effort, we found that several but not all force fields yield a hydration shell contrast in remarkable agreement with experiments. The hydration shell contrast captured by Rg values depends strongly on protein charge and geometric shape, thus providing a protein-specific footprint of protein-water interactions and a novel observable for scrutinizing atomistic hydration shell models against experimental data.
Collapse
Affiliation(s)
- Johanna-Barbara Linse
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany.
| |
Collapse
|
41
|
Bronkhorst AW, Lee CY, Möckel MM, Ruegenberg S, de Jesus Domingues AM, Sadouki S, Piccinno R, Sumiyoshi T, Siomi MC, Stelzl L, Luck K, Ketting RF. An extended Tudor domain within Vreteno interconnects Gtsf1L and Ago3 for piRNA biogenesis in Bombyx mori. EMBO J 2023; 42:e114072. [PMID: 37984437 PMCID: PMC10711660 DOI: 10.15252/embj.2023114072] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 10/25/2023] [Accepted: 10/30/2023] [Indexed: 11/22/2023] Open
Abstract
Piwi-interacting RNAs (piRNAs) direct PIWI proteins to transposons to silence them, thereby preserving genome integrity and fertility. The piRNA population can be expanded in the ping-pong amplification loop. Within this process, piRNA-associated PIWI proteins (piRISC) enter a membraneless organelle called nuage to cleave their target RNA, which is stimulated by Gtsf proteins. The resulting cleavage product gets loaded into an empty PIWI protein to form a new piRISC complex. However, for piRNA amplification to occur, the new RNA substrates, Gtsf-piRISC, and empty PIWI proteins have to be in physical proximity. In this study, we show that in silkworm cells, the Gtsf1 homolog BmGtsf1L binds to piRNA-loaded BmAgo3 and localizes to granules positive for BmAgo3 and BmVreteno. Biochemical assays further revealed that conserved residues within the unstructured tail of BmGtsf1L directly interact with BmVreteno. Using a combination of AlphaFold modeling, atomistic molecular dynamics simulations, and in vitro assays, we identified a novel binding interface on the BmVreteno-eTudor domain, which is required for BmGtsf1L binding. Our study reveals that a single eTudor domain within BmVreteno provides two binding interfaces and thereby interconnects piRNA-loaded BmAgo3 and BmGtsf1L.
Collapse
Affiliation(s)
| | - Chop Y Lee
- International PhD Programme on Gene Regulation, Epigenetics & Genome StabilityMainzGermany
- Integrative Systems Biology GroupInstitute of Molecular BiologyMainzGermany
| | - Martin M Möckel
- Protein Production Core FacilityInstitute of Molecular BiologyMainzGermany
| | - Sabine Ruegenberg
- Protein Production Core FacilityInstitute of Molecular BiologyMainzGermany
| | - Antonio M de Jesus Domingues
- Biology of Non‐coding RNA GroupInstitute of Molecular BiologyMainzGermany
- Present address:
Dewpoint Therapeutics GmbHDresdenGermany
| | - Shéraz Sadouki
- Biology of Non‐coding RNA GroupInstitute of Molecular BiologyMainzGermany
| | - Rossana Piccinno
- Microscopy Core FacilityInstitute of Molecular BiologyMainzGermany
| | - Tetsutaro Sumiyoshi
- Department of Biological Sciences, Graduate School of ScienceThe University of TokyoTokyoJapan
- Present address:
Department of Medical Innovations, Osaka Research Center for Drug DiscoveryOtsuka Pharmaceutical Co., Ltd.OsakaJapan
| | - Mikiko C Siomi
- Department of Biological Sciences, Graduate School of ScienceThe University of TokyoTokyoJapan
| | - Lukas Stelzl
- Faculty of BiologyJohannes Gutenberg University MainzMainzGermany
- KOMET 1, Institute of PhysicsJohannes Gutenberg University MainzMainzGermany
| | - Katja Luck
- Integrative Systems Biology GroupInstitute of Molecular BiologyMainzGermany
| | - René F Ketting
- Biology of Non‐coding RNA GroupInstitute of Molecular BiologyMainzGermany
- Institute of Developmental Biology and NeurobiologyJohannes Gutenberg UniversityMainzGermany
| |
Collapse
|
42
|
Ballabio F, Paissoni C, Bollati M, de Rosa M, Capelli R, Camilloni C. Accurate and Efficient SAXS/SANS Implementation Including Solvation Layer Effects Suitable for Molecular Simulations. J Chem Theory Comput 2023; 19:8401-8413. [PMID: 37923304 PMCID: PMC10687869 DOI: 10.1021/acs.jctc.3c00864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Revised: 10/11/2023] [Accepted: 10/24/2023] [Indexed: 11/07/2023]
Abstract
Small-angle X-ray and neutron scattering (SAXS/SANS) provide valuable insights into the structure and dynamics of biomolecules in solution, complementing a wide range of structural techniques, including molecular dynamics simulations. As contrast-based methods, they are sensitive not only to structural properties but also to solvent-solute interactions. Their use in molecular dynamics simulations requires a forward model that should be as fast and accurate as possible. In this work, we demonstrate the feasibility of calculating SAXS and SANS intensities using a coarse-grained representation consisting of one bead per amino acid and three beads per nucleic acid, with form factors that can be corrected on the fly to account for solvation effects at no additional computational cost. By coupling this forward model with molecular dynamics simulations restrained with SAS data, it is possible to determine conformational ensembles or refine the structure and dynamics of proteins and nucleic acids in agreement with the experimental results. To assess the robustness of this approach, we applied it to gelsolin, for which we acquired SAXS data on its closed state, and to a UP1-microRNA complex, for which we used previously collected measurements. Our hybrid-resolution small-angle scattering (hySAS) implementation, being distributed in PLUMED, can be used with atomistic and coarse-grained simulations using diverse restraining strategies.
Collapse
Affiliation(s)
- Federico Ballabio
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
| | - Cristina Paissoni
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
| | - Michela Bollati
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
- Istituto
di Biofisica, Consiglio Nazionale delle
Ricerche (IBF-CNR), via
Alfonso Corti 12, 20133 Milano, Italy
| | - Matteo de Rosa
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
- Istituto
di Biofisica, Consiglio Nazionale delle
Ricerche (IBF-CNR), via
Alfonso Corti 12, 20133 Milano, Italy
| | - Riccardo Capelli
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
| | - Carlo Camilloni
- Dipartimento
di Bioscienze, Università degli Studi
di Milano, via Celoria 26, 20133 Milano, Italy
| |
Collapse
|
43
|
Sun W, Lebedenko OO, Salguero NG, Shannon MD, Zandian M, Poirier MG, Skrynnikov NR, Jaroniec CP. Conformational and Interaction Landscape of Histone H4 Tails in Nucleosomes Probed by Paramagnetic NMR Spectroscopy. J Am Chem Soc 2023; 145:25478-25485. [PMID: 37943892 PMCID: PMC10719895 DOI: 10.1021/jacs.3c10340] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2023]
Abstract
The fundamental repeat unit of chromatin, the nucleosome, consists of approximately 147 base pairs of double-stranded DNA and a histone protein octamer containing two copies each of histones H2A, H2B, H3, and H4. Each histone possesses a dynamically disordered N-terminal tail domain, and it is well-established that the tails of histones H3 and H4 play key roles in chromatin compaction and regulation. Here we investigate the conformational ensemble and interactions of the H4 tail in nucleosomes by means of solution NMR measurements of paramagnetic relaxation enhancements (PREs) in recombinant samples reconstituted with 15N-enriched H4 and nitroxide spin-label tagged H3. The experimental PREs, which report on the proximities of individual H4 tail residues to the different H3 spin-label sites, are interpreted by using microsecond time-scale molecular dynamics simulations of the nucleosome core particle. Collectively, these data enable improved localization of histone H4 tails in nucleosomes and support the notion that H4 tails engage in a fuzzy complex interaction with nucleosomal DNA.
Collapse
Affiliation(s)
- Wenjun Sun
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Olga O. Lebedenko
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg 199034, Russia
| | - Nicole Gonzalez Salguero
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Matthew D. Shannon
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Mohamad Zandian
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Michael G. Poirier
- Department of Physics, The Ohio State University, Columbus, Ohio 43210, United States
| | - Nikolai R. Skrynnikov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg 199034, Russia
- Department of Chemistry, Purdue University, West Lafayette 47907, United States
| | - Christopher P. Jaroniec
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
44
|
Jana AK, Güven Ö, Yaşar F. The stability and dynamics of the Aβ40/Aβ42 interlaced mixed fibrils. J Biomol Struct Dyn 2023:1-14. [PMID: 37964619 DOI: 10.1080/07391102.2023.2280765] [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: 06/21/2023] [Accepted: 11/01/2023] [Indexed: 11/16/2023]
Abstract
The accumulation of fibrillar amyloid-β (Aβ) aggregates in the brain, predominantly comprising 40- and 42-residue amyloid-β (Aβ40 and Aβ42), is a major pathological hallmark of Alzheimer's disease (AD). Aβ40 and Aβ42 naturally coexist in the brain under normal physiological conditions, and their interplay is generally considered to be a critical factor in the progression of AD. In addition to forming homogeneous oligomers and fibrils, Aβ40 and Aβ42 are also reported to co-assemble into hetero-oligomers and interlaced mixed fibrils, as evidenced by solid-state nuclear magnetic resonance spectroscopy (NMR), high molecular weight mass spectrometry and cross-seeding experiments. However, the exact molecular mechanisms underlying these processes remain unclear. In this study, we have used a recently resolved structurally uniform 1:1 mixture of Aβ40/Aβ42 interlaced mixed fibril as a prototype to gain insights into the molecular-level interactions between Aβ40 and Aβ42. We employed fully atomistic molecular dynamics simulation and compared the results with a homogeneous U-shaped Aβ40 fibrillar model. Our simulations using two different force fields provide conclusive evidence that the Aβ40/Aβ42 interlaced mixed fibril is energetically more favorable than the homogeneous Aβ40 fibrillar model. Furthermore, we also show that the increase in stability observed in the mixed model stems primarily from the packing interfaces and the stacking interfaces between C-termini. Our simulation results provide valuable mechanistic insights that are not readily accessible in experiment and could have significant implications for both the pathogenesis of AD and the development of current therapeutic strategies.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Asis K Jana
- Department of Microbiology and Biotechnology, Sister Nivedita University, Kolkata, West Bengal, India
| | - Özgür Güven
- Department of Physics Engineering, Hacettepe University, Ankara, Türkiye
| | - Fatih Yaşar
- Department of Physics Engineering, Hacettepe University, Ankara, Türkiye
| |
Collapse
|
45
|
Wilson C, Karttunen M, de Groot BL, Gapsys V. Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy. J Chem Theory Comput 2023; 19:7833-7845. [PMID: 37820376 PMCID: PMC10653114 DOI: 10.1021/acs.jctc.3c00721] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Indexed: 10/13/2023]
Abstract
The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.
Collapse
Affiliation(s)
- Carter
J. Wilson
- Department
of Mathematics, The University of Western
Ontario, N6A 5B7 London, Canada
- Centre
for Advanced Materials and Biomaterials Research (CAMBR), The University of Western Ontario, N6A 5B7 London, Canada
| | - Mikko Karttunen
- Centre
for Advanced Materials and Biomaterials Research (CAMBR), The University of Western Ontario, N6A 5B7 London, Canada
- Department
of Physics & Astronomy, The University
of Western Ontario, N6A
5B7 London, Canada
- Department
of Chemistry, The University of Western
Ontario, N6A 5B7 London, Canada
| | - Bert L. de Groot
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Multidisciplinary
Sciences, 37077 Göttingen, Germany
| | - Vytautas Gapsys
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Multidisciplinary
Sciences, 37077 Göttingen, Germany
- Computational
Chemistry, Janssen Research & Development, Janssen Pharmaceutica N. V., Turnhoutseweg 30, B-2340 Beerse, Belgium
| |
Collapse
|
46
|
Herrera-Rodríguez AM, Dasanna AK, Daday C, Cruz-Chú ER, Aponte-Santamaría C, Schwarz US, Gräter F. The role of flow in the self-assembly of dragline spider silk proteins. Biophys J 2023; 122:4241-4253. [PMID: 37803828 PMCID: PMC10645567 DOI: 10.1016/j.bpj.2023.09.020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 07/14/2023] [Accepted: 09/29/2023] [Indexed: 10/08/2023] Open
Abstract
Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into β sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.
Collapse
Affiliation(s)
| | - Anil Kumar Dasanna
- BioQuant, Heidelberg University, Heidelberg, Germany; Institute for Theoretical Physics, Heidelberg University, Heidelberg, Germany
| | - Csaba Daday
- Heidelberg Institute for Theoretical Studies, Heidelberg, Germany
| | - Eduardo R Cruz-Chú
- Laboratorios de Investigación y Desarrollo, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | | | - Ulrich S Schwarz
- BioQuant, Heidelberg University, Heidelberg, Germany; Institute for Theoretical Physics, Heidelberg University, Heidelberg, Germany.
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies, Heidelberg, Germany; Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany.
| |
Collapse
|
47
|
Kumar S, Reddy G. Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch. J Phys Chem B 2023; 127:9267-9281. [PMID: 37851949 DOI: 10.1021/acs.jpcb.3c03941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.
Collapse
Affiliation(s)
- Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| |
Collapse
|
48
|
Pérez-Trejo I, Dominguez L. GaMD simulations as an alternative in the TFE-water mixture description. J Mol Model 2023; 29:352. [PMID: 37906368 PMCID: PMC10618327 DOI: 10.1007/s00894-023-05749-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 10/04/2023] [Indexed: 11/02/2023]
Abstract
CONTEXT 2,2,2-Trifluoroethanol has been widely used to study the structure and dynamic properties of intrinsically disordered proteins. Experimentally, it is known that TFE-water mixtures stabilize secondary structures of IDPs, and therefore, it allows the studying of conformational ensembles of these proteins. In the last decades, molecular dynamic simulations have helped study the IDPs' conformational ensemble. Unfortunately, conventional MD requires very long simulation times to describe the properties of IDPs. Therefore, a variety of accelerated sampling techniques have been developed and employed. The TFE-water mixture arrangement description through MD has faced substantial difficulties since emulating the TFE nanocrowding at certain TFE:H[Formula: see text]O ratios (around 15-40% of TFE). In this work, we determine the most suitable conditions that reproduce experimentally reported properties of TFE-water mixtures. We compared the employment of conventional MD and GaMD simulations and various water parameters. Our results show that the combination of parameters that better reproduce the experimental information is the combination of the TIP4PD water model and GaMD simulations. Therefore, these conditions help accurately describe the structural ensemble of IDPs in TFE-water mixtures. METHODS Conventional MD and GaMD simulations were performed under AMBER 18 software. The TFE and water molecules were described using GAFF2 and a variety of water models, such as TIP3P, TIP4P2005, TIP4PD, and TIP5P, respectively. The systems were simulated a 100 ns at 298 K.
Collapse
Affiliation(s)
- Itzel Pérez-Trejo
- Departamento de Fisicoquímica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City, 04510, Mexico
| | - Laura Dominguez
- Departamento de Fisicoquímica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City, 04510, Mexico.
| |
Collapse
|
49
|
Sridar J, Mafi A, Judge RA, Xu J, Kong KA, Wang JCK, Stoll VS, Koukos G, Simon RJ, Eaton D, Bratkowski M, Hao Q. Cryo-EM structure of human PAPP-A2 and mechanism of substrate recognition. Commun Chem 2023; 6:234. [PMID: 37898658 PMCID: PMC10613257 DOI: 10.1038/s42004-023-01032-y] [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: 05/05/2023] [Accepted: 10/18/2023] [Indexed: 10/30/2023] Open
Abstract
Pregnancy-Associated Plasma Protein A isoforms, PAPP-A and PAPP-A2, are metalloproteases that cleave insulin-like growth factor binding proteins (IGFBPs) to modulate insulin-like growth factor signaling. The structures of homodimeric PAPP-A in complex with IGFBP5 anchor peptide, and inhibitor proteins STC2 and proMBP have been recently reported. Here, we present the single-particle cryo-EM structure of the monomeric, N-terminal LG, MP, and the M1 domains (with the exception of LNR1/2) of human PAPP-A2 to 3.13 Å resolution. Our structure together with functional studies provides insight into a previously reported patient mutation that inactivates PAPP-A2 in a distal region of the protein. Using a combinational approach, we suggest that PAPP-A2 recognizes IGFBP5 in a similar manner as PAPP-A and show that PAPP-A2 cleaves IGFBP5 less efficiently due to differences in the M2 domain. Overall, our studies characterize the cleavage mechanism of IGFBP5 by PAPP-A2 and shed light onto key differences with its paralog PAPP-A.
Collapse
Affiliation(s)
- Janani Sridar
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | | | | | - Jun Xu
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | - Kailyn A Kong
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | - John C K Wang
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | | | - Georgios Koukos
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | - Reyna J Simon
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | - Dan Eaton
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA
| | | | - Qi Hao
- Calico Life Sciences LLC, South San Francisco, CA, 94080, USA.
| |
Collapse
|
50
|
Yildiz CB, Kundu T, Gehrmann J, Koesling J, Ravaei A, Wolff P, Kraft F, Maié T, Jakovcevski M, Pensold D, Zimmermann O, Rossetti G, Costa IG, Zimmer-Bensch G. EphrinA5 regulates cell motility by modulating Snhg15/DNA triplex-dependent targeting of DNMT1 to the Ncam1 promoter. Epigenetics Chromatin 2023; 16:42. [PMID: 37880732 PMCID: PMC10601256 DOI: 10.1186/s13072-023-00516-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Accepted: 10/13/2023] [Indexed: 10/27/2023] Open
Abstract
Cell-cell communication is mediated by membrane receptors and their ligands, such as the Eph/ephrin system, orchestrating cell migration during development and in diverse cancer types. Epigenetic mechanisms are key for integrating external "signals", e.g., from neighboring cells, into the transcriptome in health and disease. Previously, we reported ephrinA5 to trigger transcriptional changes of lncRNAs and protein-coding genes in cerebellar granule cells, a cell model for medulloblastoma. LncRNAs represent important adaptors for epigenetic writers through which they regulate gene expression. Here, we investigate a lncRNA-mediated targeting of DNMT1 to specific gene loci by the combined power of in silico modeling of RNA/DNA interactions and wet lab approaches, in the context of the clinically relevant use case of ephrinA5-dependent regulation of cellular motility of cerebellar granule cells. We provide evidence that Snhg15, a cancer-related lncRNA, recruits DNMT1 to the Ncam1 promoter through RNA/DNA triplex structure formation and the interaction with DNMT1. This mediates DNA methylation-dependent silencing of Ncam1, being abolished by ephrinA5 stimulation-triggered reduction of Snhg15 expression. Hence, we here propose a triple helix recognition mechanism, underlying cell motility regulation via lncRNA-targeted DNA methylation in a clinically relevant context.
Collapse
Affiliation(s)
- Can Bora Yildiz
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
- Research Training Group 2416 Multi Senses - Multi Scales, RWTH Aachen University, 52074, Aachen, Germany
| | - Tathagata Kundu
- Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| | - Julia Gehrmann
- Institute for Computational Genomics, RWTH Aachen University, Medical Faculty, 52074, Aachen, Germany
| | - Jannis Koesling
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Amin Ravaei
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
- Department of Neurosciences and Rehabilitation, Section of Medical Biochemistry, Molecular Biology and Genetics, University of Ferrara, Ferrara, Italy
| | - Philip Wolff
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Florian Kraft
- Institute for Human Genetics and Genomic Medicine, Medical Faculty, RWTH Aachen University, 52074, Aachen, Germany
| | - Tiago Maié
- Institute for Computational Genomics, RWTH Aachen University, Medical Faculty, 52074, Aachen, Germany
| | - Mira Jakovcevski
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Daniel Pensold
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Olav Zimmermann
- Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| | - Giulia Rossetti
- Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
- Department of Neurology, University Hospital Aachen, RWTH Aachen University, Aachen, Germany
- Institute of Neuroscience and Medicine (INM-9)/Institute of Advanced Simulations (IAS-5), Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| | - Ivan G Costa
- Institute for Computational Genomics, RWTH Aachen University, Medical Faculty, 52074, Aachen, Germany
| | - Geraldine Zimmer-Bensch
- Institute of Zoology (Biology 2), Division of Neuroepigenetics, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany.
- Research Training Group 2416 Multi Senses - Multi Scales, RWTH Aachen University, 52074, Aachen, Germany.
| |
Collapse
|