1
|
Li T, Motta S, He Y. Deciphering the Mystery in p300 Taz2-p53 TAD2 Recognition. J Chem Theory Comput 2024. [PMID: 39141804 DOI: 10.1021/acs.jctc.4c00541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/16/2024]
Abstract
Intrinsically disordered proteins (IDPs) engage in various fundamental biological activities, and their behavior is of particular importance for a better understanding of the verbose but well-organized signal transduction in cells. IDPs exhibit uniquely paradoxical features with low affinity but simultaneously high specificity in recognizing their binding targets. The transcription factor p53 plays a crucial role in cancer suppression, carrying out some of its biological functions using its disordered regions, such as N-terminal transactivation domain 2 (TAD2). Exploration of the binding and unbinding processes between proteins is challenging, and the inherently disordered properties of these regions further complicate the issue. Computer simulations are a powerful tool to complement the experiments to fill gaps to explore the binding/unbinding processes between proteins. Here, we investigated the binding mechanism between p300 Taz2 and p53 TAD2 through extensive molecular dynamics (MD) simulations using the physics-based UNited RESidue (UNRES) force field with additional Go̅-like potentials. Distance restraints extracted from the NMR-resolved structures were imposed on intermolecular residue pairs to accelerate binding simulations, in which Taz2 was immobilized in a native-like conformation and disordered TAD2 was fully free. Starting from six structures with TAD2 placed at different positions around Taz2, we observed a metastable intermediate state in which the middle helical segment of TAD2 is anchored in the binding pocket, highlighting the significance of the TAD2 helix in directing protein recognition. Physics-based binding simulations show that successful binding is achieved after a series of stages, including (1) protein collisions to initiate the formation of encounter complexes, (2) partial attachment of TAD2, and finally (3) full attachment of TAD2 to the correct binding pocket of Taz2. Furthermore, machine-learning-based PathDetect-SOM was used to identify two binding pathways, the encounter complexes, and the intermediate states.
Collapse
Affiliation(s)
- Tongtong Li
- Department of Chemistry & Chemical Biology, The University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Stefano Motta
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Milan 20126, Italy
| | - Yi He
- Department of Chemistry & Chemical Biology, The University of New Mexico, Albuquerque, New Mexico 87131, United States
- Translational Informatics Division, Department of Internal Medicine, The University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
2
|
Tripathi A, Del Galdo S, Chandramouli B, Kumar N. Distinct dynamical features of plasmodial and human HSP70-HSP110 highlight the divergence in their chaperone-assisted protein folding. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2023; 1871:140942. [PMID: 37516289 DOI: 10.1016/j.bbapap.2023.140942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2023] [Revised: 07/24/2023] [Accepted: 07/26/2023] [Indexed: 07/31/2023]
Abstract
HSP70 and its evolutionarily diverged co-chaperone HSP110, forms an important node in protein folding cascade. How these proteins maintain the aggregation-prone proteome of malaria parasite in functional state remains underexplored, in contrast to its human orthologs. In this study, we have probed into conformational dynamics of plasmodial HSP70 and HSP110 through multiple μs MD-simulations (ATP-state) and compared with their respective human counterparts. Simulations covered sampling of 3.4 and 2.8 μs for HSP70 and HSP110, respectively, for parasite and human orthologs. We provide a comprehensive description of the dynamic behaviors that characterize the systems and also introduce a parameter for quantifying protein rigidity. For HSP70, the interspecies comparison reveals enhanced flexibility in IA and IB subdomain within the conserved NBD, lesser solvent accessibility of the interdomain linker and distinct dynamics of the SBDβ of Pf HSP70 in comparison to Hs HSP70. In the case of HSP110, notable contrast in the dynamics of NBD, SBDβ and SBDα was observed between parasite and human ortholog. Although HSP70 and HSP110 are members of the same superfamily, we identified specific differences in the subdomain contacts in NBD, linker properties and interdomain movements in their human and parasite orthologs. Our study suggests that differences in conformational dynamics may translate into species-specific differences in the chaperoning activities of HSP70-HSP110 in the parasite and human, respectively. Dynamical features of Pf HSP70-HSP110 may contribute to the maintenance of proteostasis in the parasite during its intracellular survival in the host.
Collapse
Affiliation(s)
- Aradhya Tripathi
- Department of Molecular Microbiology and Immunology, CSIR-Central Drug Research Institute (CSIR-CDRI), Sector 10, Jankipuram Extension, Lucknow 226031, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Sara Del Galdo
- Science Department, University of Roma Tre, Via della Vasca Navale 84, Rome, Italy
| | | | - Niti Kumar
- Department of Molecular Microbiology and Immunology, CSIR-Central Drug Research Institute (CSIR-CDRI), Sector 10, Jankipuram Extension, Lucknow 226031, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India.
| |
Collapse
|
3
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
4
|
Danielsson A, Samsonov SA, Liwo A, Sieradzan AK. Extension of the SUGRES-1P Coarse-Grained Model of Polysaccharides to Heparin. J Chem Theory Comput 2023; 19:6023-6036. [PMID: 37587433 PMCID: PMC10500997 DOI: 10.1021/acs.jctc.3c00511] [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/17/2023] [Indexed: 08/18/2023]
Abstract
Heparin is an unbranched periodic polysaccharide composed of negatively charged monomers and involved in key biological processes, including anticoagulation, angiogenesis, and inflammation. Its structure and dynamics have been studied extensively using experimental as well as theoretical approaches. The conventional approach of computational chemistry applied to the analysis of biomolecules is all-atom molecular dynamics, which captures the interactions of individual atoms by solving Newton's equation of motion. An alternative is molecular dynamics simulations using coarse-grained models of biomacromolecules, which offer a reduction of the representation and consequently enable us to extend the time and size scale of simulations by orders of magnitude. In this work, we extend the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules developed in our laboratory to heparin. We carried out extensive tests to estimate the optimal weights of energy terms of the effective energy function as well as the optimal Debye-Hückel screening factor for electrostatic interactions. We applied the model to study unbound heparin molecules of polymerization degree ranging from 6 to 68 residues. We compare the obtained coarse-grained heparin conformations with models obtained from X-ray diffraction studies of heparin. The SUGRES-1P force field was able to accurately predict the general shape and global characteristics of heparin molecules.
Collapse
Affiliation(s)
- Annemarie Danielsson
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Sergey A. Samsonov
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam K. Sieradzan
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
5
|
Li T, Hendrix E, He Y. Simple and Effective Conformational Sampling Strategy for Intrinsically Disordered Proteins Using the UNRES Web Server. J Phys Chem B 2023; 127:2177-2186. [PMID: 36827446 DOI: 10.1021/acs.jpcb.2c08945] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/26/2023]
Abstract
Intrinsically disordered proteins (IDPs) contain more charged amino acids than folded proteins, resulting in a lack of hydrophobic core(s) and a tendency to adopt rapidly interconverting structures rather than well-defined structures. The structural heterogeneity of IDPs, encoded by the amino acid sequence, is closely related to their unique roles in biological pathways, which require them to interact with different binding partners. Recently, Robustelli and co-workers have demonstrated that a balanced all-atom force field can be used to sample heterogeneous structures of disordered proteins ( Proc. Natl. Acad. Sci. U.S.A. 2018, 115, E4758-E4766). However, such a solution requires extensive computational resources, such as Anton supercomputers. Here, we propose a simple and effective solution to sample the conformational space of IDPs using a publicly available web server, namely, the UNited-RESidue (UNRES) web server. Our proposed solution requires no investment in computational resources and no prior knowledge of UNRES. UNRES Replica Exchange Molecular Dynamics (REMD) simulations were carried out on a set of eight disordered proteins at temperatures spanning from 270 to 430 K. Utilizing the latest UNRES force field designed for structured proteins, with proper selections of temperatures, we were able to produce comparable results to all-atom force fields as reported in work done by Robustelli and co-workers. In addition, NMR observables and the radius of gyration calculated from UNRES ensembles were directly compared with the experimental data to further evaluate the accuracy of the UNRES model at all temperatures. Our results suggest that carrying out the UNRES simulations at optimal temperatures using the UNRES web server can be a good alternative to sample heterogeneous structures of IDPs.
Collapse
Affiliation(s)
- Tongtong Li
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Emily Hendrix
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yi He
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States.,Translational Informatics Division, Department of Internal Medicine, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
6
|
Lipska AG, Sieradzan AK, Czaplewski C, Lipińska AD, Ocetkiewicz KM, Proficz J, Czarnul P, Krawczyk H, Liwo A. Long-time scale simulations of virus-like particles from three human-norovirus strains. J Comput Chem 2023; 44:1470-1483. [PMID: 36799410 DOI: 10.1002/jcc.27087] [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: 11/07/2022] [Revised: 12/22/2022] [Accepted: 01/29/2023] [Indexed: 02/18/2023]
Abstract
The dynamics of the virus like particles (VLPs) corresponding to the GII.4 Houston, GII.2 SMV, and GI.1 Norwalk strains of human noroviruses (HuNoV) that cause gastroenteritis was investigated by means of long-time (about 30 μs in the laboratory timescale) molecular dynamics simulations with the coarse-grained UNRES force field. The main motion of VLP units turned out to be the bending at the junction between the P1 subdomain (that sits in the VLP shell) and the P2 subdomain (that protrudes outside) of the major VP1 protein, this resulting in a correlated wagging motion of the P2 subdomains with respect to the VLP surface. The fluctuations of the P2 subdomain were found to be more pronounced and the P2 domain made a greater angle with the normal to the VLP surface for the GII.2 strain, which could explain the inability of this strain to bind the histo-blood group antigens (HBGAs).
Collapse
Affiliation(s)
- Agnieszka G Lipska
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam K Sieradzan
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Cezary Czaplewski
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Andrea D Lipińska
- Laboratory of Virus Molecular Biology, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Krzysztof M Ocetkiewicz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Jerzy Proficz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Paweł Czarnul
- Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Henryk Krawczyk
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam Liwo
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| |
Collapse
|
7
|
Liu Y, Li M, Yu J, Ma A, Wang J, Yun DJ, Xu ZY. Plasma membrane-localized Hsp40/DNAJ chaperone protein facilitates OsSUVH7-OsBAG4-OsMYB106 transcriptional complex formation for OsHKT1;5 activation. JOURNAL OF INTEGRATIVE PLANT BIOLOGY 2023; 65:265-279. [PMID: 36349953 DOI: 10.1111/jipb.13403] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Accepted: 11/07/2022] [Indexed: 06/16/2023]
Abstract
The salinization of irrigated land affects agricultural productivity. HIGH-AFFINITY POTASSIUM (K+ ) TRANSPORTER 1;5 (OsHKT1;5)-dependent sodium (Na+ ) transport is a key salt tolerance mechanism during rice growth and development. Using a previously generated high-throughput activation tagging-based T-DNA insertion mutant pool, we isolated a mutant exhibiting salt stress-sensitive phenotype, caused by a reduction in OsHKT1;5 transcripts. The salt stress-sensitive phenotype of this mutant results from the loss of function of OsDNAJ15, which encodes plasma membrane-localized heat shock protein 40 (Hsp40). osdnaj15 loss-of-function mutants show decreased plant height, increased leaf angle, and reduced grain number caused by shorter panicle length and fewer branches. On the other h'and, OsDNAJ15-overexpression plants showed salt stress-tolerant phenotypes. Intriguingly, salt stress facilitates the nuclear relocation of OsDNAJ15 so that it can interact with OsBAG4, and OsDNAJ15 and OsBAG4 synergistically facilitate the DNA-binding activity of OsMYB106 to positively regulate the expression of OsHKT1;5. Overall, our results reveal a novel function of plasma membrane-localized Hsp40 protein in modulating, alongside chaperon regulator OsBAG4, transcriptional regulation under salinity stress tolerance.
Collapse
Affiliation(s)
- Yutong Liu
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| | - Mengting Li
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| | - Jinlei Yu
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| | - Ao Ma
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| | - Jie Wang
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| | - Dae-Jin Yun
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
- Department of Biomedical Science and Engineering, Konkuk University, Seoul, South Korea
| | - Zheng-Yi Xu
- Key Laboratory of Molecular Epigenetics of the Ministry of Education (MOE), Northeast Normal University, Changchun, 130024, China
| |
Collapse
|
8
|
Hendrix E, Motta S, Gahl RF, He Y. Insight into the Initial Stages of the Folding Process in Onconase Revealed by UNRES. J Phys Chem B 2022; 126:7934-7942. [PMID: 36179061 DOI: 10.1021/acs.jpcb.2c04770] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The unfolded state of proteins presents many challenges to elucidate the structural basis for biological function. This state is characterized by a large degree of structural heterogeneity which makes it difficult to generate structural models. However, recent experiments into the initial folding events of the 104-residue ribonuclease homologue onconase (ONC) were able to identify the regions in the protein that participate in the initial folding of this protein. Therefore, to gain additional structural insight into the unfolded state of proteins, this study utilized molecular dynamics simulations using the UNited-RESidue (UNRES) force field to evaluate whether there is a good agreement between the experimentally determined initial structures and the structures identified by computer simulations along a folding pathway. Indeed, these UNRES simulations accurately identified the two regions experimentally observed to form the initial native structure along the folding pathway of ONC. In addition, these regions are determined to be chain folding initiation sites (CFIS) according to methods developed previously. Subsequent self-organization maps (SOM) analysis has revealed key structural states involved in these early folding events.
Collapse
Affiliation(s)
- Emily Hendrix
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico87131, United States
| | - Stefano Motta
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Milan20126, Italy
| | - Robert F Gahl
- Division of Extramural Activities, National Cancer Institute, National Institutes of Health, Bethesda, Maryland20850, United States
| | - Yi He
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico87131, United States.,Translational Informatics Division, Department of Internal Medicine, University of New Mexico, Albuquerque, New Mexico87131, United States
| |
Collapse
|
9
|
Nordquist EB, Clerico EM, Chen J, Gierasch LM. Computationally-Aided Modeling of Hsp70-Client Interactions: Past, Present, and Future. J Phys Chem B 2022; 126:6780-6791. [PMID: 36040440 PMCID: PMC10309085 DOI: 10.1021/acs.jpcb.2c03806] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Hsp70 molecular chaperones play central roles in maintaining a healthy cellular proteome. Hsp70s function by binding to short peptide sequences in incompletely folded client proteins, thus preventing them from misfolding and/or aggregating, and in many cases holding them in a state that is competent for subsequent processes like translocation across membranes. There is considerable interest in predicting the sites where Hsp70s may bind their clients, as the ability to do so sheds light on the cellular functions of the chaperone. In addition, the capacity of the Hsp70 chaperone family to bind to a broad array of clients and to identify accessible sequences that enable discrimination of those that are folded from those that are not fully folded, which is essential to their cellular roles, is a fascinating puzzle in molecular recognition. In this article we discuss efforts to harness computational modeling with input from experimental data to develop a predictive understanding of the promiscuous yet selective binding of Hsp70 molecular chaperones to accessible sequences within their client proteins. We trace how an increasing understanding of the complexities of Hsp70-client interactions has led computational modeling to new underlying assumptions and design features. We describe the trend from purely data-driven analysis toward increased reliance on physics-based modeling that deeply integrates structural information and sequence-based functional data with physics-based binding energies. Notably, new experimental insights are adding to our understanding of the molecular origins of "selective promiscuity" in substrate binding by Hsp70 chaperones and challenging the underlying assumptions and design used in earlier predictive models. Taking the new experimental findings together with exciting progress in computational modeling of protein structures leads us to foresee a bright future for a predictive understanding of selective-yet-promiscuous binding exploited by Hsp70 molecular chaperones; the resulting new insights will also apply to substrate binding by other chaperones and by signaling proteins.
Collapse
Affiliation(s)
- Erik B. Nordquist
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts, 01003, United States
| | - Eugenia M. Clerico
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts, 01003, United States
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts, 01003, United States
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts, 01003, United States
| | - Lila M. Gierasch
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts, 01003, United States
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts, 01003, United States
| |
Collapse
|
10
|
Guzzo A, Delarue P, Rojas A, Nicolaï A, Maisuradze GG, Senet P. Wild-Type α-Synuclein and Variants Occur in Different Disordered Dimers and Pre-Fibrillar Conformations in Early Stage of Aggregation. Front Mol Biosci 2022; 9:910104. [PMID: 35836937 PMCID: PMC9273784 DOI: 10.3389/fmolb.2022.910104] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 05/16/2022] [Indexed: 12/15/2022] Open
Abstract
α-Synuclein is a 140 amino-acid intrinsically disordered protein mainly found in the brain. Toxic α-synuclein aggregates are the molecular hallmarks of Parkinson’s disease. In vitro studies showed that α-synuclein aggregates in oligomeric structures of several 10th of monomers and into cylindrical structures (fibrils), comprising hundred to thousands of proteins, with polymorphic cross-β-sheet conformations. Oligomeric species, formed at the early stage of aggregation remain, however, poorly understood and are hypothezised to be the most toxic aggregates. Here, we studied the formation of wild-type (WT) and mutant (A30P, A53T, and E46K) dimers of α-synuclein using coarse-grained molecular dynamics. We identified two principal segments of the sequence with a higher propensity to aggregate in the early stage of dimerization: residues 36–55 and residues 66–95. The transient α-helices (residues 53–65 and 73–82) of α-synuclein monomers are destabilized by A53T and E46K mutations, which favors the formation of fibril native contacts in the N-terminal region, whereas the helix 53–65 prevents the propagation of fibril native contacts along the sequence for the WT in the early stages of dimerization. The present results indicate that dimers do not adopt the Greek key motif of the monomer fold in fibrils but form a majority of disordered aggregates and a minority (9–15%) of pre-fibrillar dimers both with intra-molecular and intermolecular β-sheets. The percentage of residues in parallel β-sheets is by increasing order monomer < disordered dimers < pre-fibrillar dimers. Native fibril contacts between the two monomers are present in the NAC domain for WT, A30P, and A53T and in the N-domain for A53T and E46K. Structural properties of pre-fibrillar dimers agree with rupture-force atomic force microscopy and single-molecule Förster resonance energy transfer available data. This suggests that the pre-fibrillar dimers might correspond to the smallest type B toxic oligomers. The probability density of the dimer gyration radius is multi-peaks with an average radius that is 10 Å larger than the one of the monomers for all proteins. The present results indicate that even the elementary α-synuclein aggregation step, the dimerization, is a complicated phenomenon that does not only involve the NAC region.
Collapse
Affiliation(s)
- Adrien Guzzo
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Patrice Delarue
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Ana Rojas
- Schrödinger, Inc., New York, NY, United States
| | - Adrien Nicolaï
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Gia G. Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, United States
| | - Patrick Senet
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, United States
- *Correspondence: Patrick Senet,
| |
Collapse
|
11
|
Modeling the Structure, Dynamics, and Transformations of Proteins with the UNRES Force Field. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2022; 2376:399-416. [PMID: 34845623 DOI: 10.1007/978-1-0716-1716-8_23] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The physics-based united-residue (UNRES) model of proteins ( www.unres.pl ) has been designed to carry out large-scale simulations of protein folding. The force field has been derived and parameterized based on the principles of statistical-mechanics, which makes it independent of structural databases and applicable to treat nonstandard situations such as, proteins that contain D-amino-acid residues. Powered by Langevin dynamics and its replica-exchange extensions, UNRES has found a variety of applications, including ab initio and database-assisted protein-structure prediction, simulating protein-folding pathways, exploring protein free-energy landscapes, and solving biological problems. This chapter provides a summary of UNRES and a guide for potential users regarding the application of the UNRES package in a variety of research tasks.
Collapse
|
12
|
Li J, Xie J, Godec A, Weninger KR, Liu C, Smith JC, Hong L. Non-ergodicity of a globular protein extending beyond its functional timescale. Chem Sci 2022; 13:9668-9677. [PMID: 36091909 PMCID: PMC9400594 DOI: 10.1039/d2sc03069a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2022] [Accepted: 07/18/2022] [Indexed: 11/21/2022] Open
Abstract
Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble. Here, by performing single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular dynamics (MD) simulations of a multi-domain globular protein, cytoplasmic protein-tyrosine phosphatase (SHP2), we demonstrate that the functional inter-domain motion is observationally non-ergodic over the time spans 10−12 to 10−7 s and 10−1 to 102 s. The difference between observational non-ergodicity and simple non-convergence is discussed. In comparison, a single-strand DNA of similar size behaves ergodically with an energy landscape resembling a one-dimensional linear chain. The observed non-ergodicity results from the hierarchical connectivity of the high-dimensional energy landscape of the protein molecule. As the characteristic time for the protein to conduct its dephosphorylation function is ∼10 s, our findings suggest that, due to the non-ergodicity, individual, seemingly identical protein molecules can be dynamically and functionally different. Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble.![]()
Collapse
Affiliation(s)
- Jun Li
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
| | - JingFei Xie
- Interdisciplinary Research Center on Biology and Chemistry, Center for Excellence in Molecular Synthesis, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai 201203, China
- University of the Chinese Academy of Sciences, Beijing 100049, China
| | - Aljaž Godec
- Mathematical BioPhysics Group, Max Planck Institute for Biophysical Chemistry, Göttingen 37077, Germany
| | - Keith R. Weninger
- Department of Physics, North Carolina State University, Raleigh, NC 27695, USA
| | - Cong Liu
- Interdisciplinary Research Center on Biology and Chemistry, Center for Excellence in Molecular Synthesis, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai 201203, China
| | - Jeremy C. Smith
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Liang Hong
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
13
|
Guzzo A, Delarue P, Rojas A, Nicolaï A, Maisuradze GG, Senet P. Missense Mutations Modify the Conformational Ensemble of the α-Synuclein Monomer Which Exhibits a Two-Phase Characteristic. Front Mol Biosci 2021; 8:786123. [PMID: 34912851 PMCID: PMC8667727 DOI: 10.3389/fmolb.2021.786123] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 10/25/2021] [Indexed: 12/15/2022] Open
Abstract
α-Synuclein is an intrinsically disordered protein occurring in different conformations and prone to aggregate in β-sheet structures, which are the hallmark of the Parkinson disease. Missense mutations are associated with familial forms of this neuropathy. How these single amino-acid substitutions modify the conformations of wild-type α-synuclein is unclear. Here, using coarse-grained molecular dynamics simulations, we sampled the conformational space of the wild type and mutants (A30P, A53P, and E46K) of α-synuclein monomers for an effective time scale of 29.7 ms. To characterize the structures, we developed an algorithm, CUTABI (CUrvature and Torsion based of Alpha-helix and Beta-sheet Identification), to identify residues in the α-helix and β-sheet from Cα-coordinates. CUTABI was built from the results of the analysis of 14,652 selected protein structures using the Dictionary of Secondary Structure of Proteins (DSSP) algorithm. DSSP results are reproduced with 93% of success for 10 times lower computational cost. A two-dimensional probability density map of α-synuclein as a function of the number of residues in the α-helix and β-sheet is computed for wild-type and mutated proteins from molecular dynamics trajectories. The density of conformational states reveals a two-phase characteristic with a homogeneous phase (state B, β-sheets) and a heterogeneous phase (state HB, mixture of α-helices and β-sheets). The B state represents 40% of the conformations for the wild-type, A30P, and E46K and only 25% for A53T. The density of conformational states of the B state for A53T and A30P mutants differs from the wild-type one. In addition, the mutant A53T has a larger propensity to form helices than the others. These findings indicate that the equilibrium between the different conformations of the α-synuclein monomer is modified by the missense mutations in a subtle way. The α-helix and β-sheet contents are promising order parameters for intrinsically disordered proteins, whereas other structural properties such as average gyration radius, Rg, or probability distribution of Rg cannot discriminate significantly the conformational ensembles of the wild type and mutants. When separated in states B and HB, the distributions of Rg are more significantly different, indicating that global structural parameters alone are insufficient to characterize the conformational ensembles of the α-synuclein monomer.
Collapse
Affiliation(s)
- Adrien Guzzo
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Patrice Delarue
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Ana Rojas
- Schrödinger, Inc., New York, NY, United States
| | - Adrien Nicolaï
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, United States
| | - Patrick Senet
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, Dijon, France.,Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, United States
| |
Collapse
|
14
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|
15
|
Giulini M, Rigoli M, Mattiotti G, Menichetti R, Tarenzi T, Fiorentini R, Potestio R. From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules. Front Mol Biosci 2021; 8:676976. [PMID: 34164432 PMCID: PMC8215203 DOI: 10.3389/fmolb.2021.676976] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 05/06/2021] [Indexed: 12/18/2022] Open
Abstract
The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.
Collapse
Affiliation(s)
- Marco Giulini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Marta Rigoli
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Giovanni Mattiotti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Roberto Menichetti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Thomas Tarenzi
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaele Fiorentini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaello Potestio
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| |
Collapse
|
16
|
Grindle MP, Carter B, Alao JP, Connors K, Tehver R, Kravats AN. Structural Communication between the E. coli Chaperones DnaK and Hsp90. Int J Mol Sci 2021; 22:ijms22042200. [PMID: 33672263 PMCID: PMC7926864 DOI: 10.3390/ijms22042200] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Revised: 02/10/2021] [Accepted: 02/12/2021] [Indexed: 01/03/2023] Open
Abstract
The 70 kDa and 90 kDa heat shock proteins Hsp70 and Hsp90 are two abundant and highly conserved ATP-dependent molecular chaperones that participate in the maintenance of cellular homeostasis. In Escherichia coli, Hsp90 (Hsp90Ec) and Hsp70 (DnaK) directly interact and collaborate in protein remodeling. Previous work has produced a model of the direct interaction of both chaperones. The locations of the residues involved have been confirmed and the model has been validated. In this study, we investigate the allosteric communication between Hsp90Ec and DnaK and how the chaperones couple their conformational cycles. Using elastic network models (ENM), normal mode analysis (NMA), and a structural perturbation method (SPM) of asymmetric and symmetric DnaK-Hsp90Ec, we extract biologically relevant vibrations and identify residues involved in allosteric signaling. When one DnaK is bound, the dominant normal modes favor biological motions that orient a substrate protein bound to DnaK within the substrate/client binding site of Hsp90Ec and release the substrate from the DnaK substrate binding domain. The presence of one DnaK molecule stabilizes the entire Hsp90Ec protomer to which it is bound. Conversely, the symmetric model of DnaK binding results in steric clashes of DnaK molecules and suggests that the Hsp90Ec and DnaK chaperone cycles operate independently. Together, this data supports an asymmetric binding of DnaK to Hsp90Ec.
Collapse
Affiliation(s)
- Matthew P. Grindle
- Department of Chemistry & Biochemistry, Miami University, Oxford, OH 45056, USA; (M.P.G.); (J.P.A.); (K.C.)
| | - Ben Carter
- Department of Physics, Denison University, Granville, OH 43023, USA; (B.C.); (R.T.)
| | - John Paul Alao
- Department of Chemistry & Biochemistry, Miami University, Oxford, OH 45056, USA; (M.P.G.); (J.P.A.); (K.C.)
| | - Katherine Connors
- Department of Chemistry & Biochemistry, Miami University, Oxford, OH 45056, USA; (M.P.G.); (J.P.A.); (K.C.)
| | - Riina Tehver
- Department of Physics, Denison University, Granville, OH 43023, USA; (B.C.); (R.T.)
| | - Andrea N. Kravats
- Department of Chemistry & Biochemistry, Miami University, Oxford, OH 45056, USA; (M.P.G.); (J.P.A.); (K.C.)
- Correspondence:
| |
Collapse
|
17
|
Stevens AO, He Y. Residue-Level Contact Reveals Modular Domain Interactions of PICK1 Are Driven by Both Electrostatic and Hydrophobic Forces. Front Mol Biosci 2021; 7:616135. [PMID: 33585564 PMCID: PMC7873044 DOI: 10.3389/fmolb.2020.616135] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2020] [Accepted: 12/15/2020] [Indexed: 11/13/2022] Open
Abstract
PICK1 is a multi-domain scaffolding protein that is uniquely comprised of both a PDZ domain and a BAR domain. While previous experiments have shown that the PDZ domain and the linker positively regulate the BAR domain and the C-terminus negatively regulates the BAR domain, the details of internal regulation mechanisms are unknown. Molecular dynamics (MD) simulations have been proven to be a useful tool in revealing the intramolecular interactions at atomic-level resolution. PICK1 performs its biological functions in a dimeric form which is extremely computationally demanding to simulate with an all-atom force field. Here, we use coarse-grained MD simulations to expose the key residues and driving forces in the internal regulations of PICK1. While the PDZ and BAR domains do not form a stable complex, our simulations show the PDZ domain preferentially interacting with the concave surface of the BAR domain over other BAR domain regions. Furthermore, our simulations show that the short helix in the linker region can form interactions with the PDZ domain. Our results reveal that the surface of the βB-βC loop, βC strand, and αA-βD loop of the PDZ domain can form a group of hydrophobic interactions surrounding the linker helix. These interactions are driven by hydrophobic forces. In contrast, our simulations reveal a very dynamic C-terminus that most often resides on the convex surface of the BAR domain rather than the previously suspected concave surface. These interactions are driven by a combination of electrostatic and hydrophobic interactions.
Collapse
Affiliation(s)
- Amy O Stevens
- Department of Chemistry and Chemical Biology, The University of New Mexico, Albuquerque, NM, United States
| | - Yi He
- Department of Chemistry and Chemical Biology, The University of New Mexico, Albuquerque, NM, United States
| |
Collapse
|
18
|
Liwo A, Czaplewski C, Sieradzan AK, Lubecka EA, Lipska AG, Golon Ł, Karczyńska A, Krupa P, Mozolewska MA, Makowski M, Ganzynkowicz R, Giełdoń A, Maciejczyk M. Scale-consistent approach to the derivation of coarse-grained force fields for simulating structure, dynamics, and thermodynamics of biopolymers. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:73-122. [PMID: 32145953 DOI: 10.1016/bs.pmbts.2019.12.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this chapter the scale-consistent approach to the derivation of coarse-grained force fields developed in our laboratory is presented, in which the effective energy function originates from the potential of mean force of the system under consideration and embeds atomistically detailed interactions in the resulting energy terms through use of Kubo's cluster-cumulant expansion, appropriate selection of the major degrees of freedom to be averaged out in the derivation of analytical approximations to the energy terms, and appropriate expression of the interaction energies at the all-atom level in these degrees of freedom. Our approach enables the developers to find correct functional forms of the effective coarse-grained energy terms, without having to import them from all-atom force fields or deriving them on a heuristic basis. In particular, the energy terms derived in such a way exhibit correct dependence on coarse-grained geometry, in particular on site orientation. Moreover, analytical formulas for the multibody (correlation) terms, which appear to be crucial for coarse-grained modeling of many of the regular structures such as, e.g., protein α-helices and β-sheets, can be derived in a systematic way. Implementation of the developed theory to the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules, which consists of the UNRES (for proteins), NARES-2P (for nucleic acids), and SUGRES-1P (for polysaccharides) components, and is being developed in our laboratory is described. Successful applications of UNICORN to the prediction of protein structure, simulating the folding and stability of proteins and nucleic acids, and solving biological problems are discussed.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland; School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea.
| | | | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland; School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea
| | - Emilia A Lubecka
- Institute of Informatics, Faculty of Mathematics, Physics, and Informatics, University of Gdańsk, Gdańsk, Poland
| | | | - Łukasz Golon
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | | | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| | | | | | | | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | - Maciej Maciejczyk
- Department of Physics and Biophysics, Faculty of Food Science, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
| |
Collapse
|
19
|
Ziȩba K, Ślusarz M, Ślusarz R, Liwo A, Czaplewski C, Sieradzan AK. Extension of the UNRES Coarse-Grained Force Field to Membrane Proteins in the Lipid Bilayer. J Phys Chem B 2019; 123:7829-7839. [PMID: 31454484 DOI: 10.1021/acs.jpcb.9b06700] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The physics-based UNRES coarse-grained force field for the simulations of protein structure and dynamics has been extended to treat membrane proteins. The lipid bilayer has been modeled by introducing a continuous nonpolar phase with the water-interface region of appropriate thickness. The potentials for average electrostatic and correlation interactions of the peptide groups have been rescaled to account for the reduction of the dielectric permittivity compared to the water phase and new potentials for protein side-chain-side-chain interactions inside and across the lipid phase have been introduced. The model was implemented in the UNRES package for coarse-grained simulations of proteins, and the package with the new functionality was tested for total energy conservation and thermostat behavior in microcanonical and canonical molecular dynamics simulations runs, respectively. The method was validated by running unrestricted ab initio blind-prediction tests of 10 short α-helical membrane proteins, all runs started from the extended structures. The modified UNRES force field was able to predict correctly the overall folds of the membrane proteins studied.
Collapse
Affiliation(s)
- Karolina Ziȩba
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Magdalena Ślusarz
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Rafał Ślusarz
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Adam Liwo
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Cezary Czaplewski
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Adam K Sieradzan
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| |
Collapse
|
20
|
Xue YL, Zhang Q, Sun Y, Zhou X, Hurley IP, Jones GW, Song Y. Using steered molecular dynamics to study the interaction between ADP and the nucleotide-binding domain of yeast Hsp70 protein Ssa1. J Comput Aided Mol Des 2018; 32:1217-1227. [PMID: 30392073 DOI: 10.1007/s10822-018-0136-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Accepted: 07/19/2018] [Indexed: 11/25/2022]
Abstract
Genetics experiments have identified six mutations located in the subdomain IA (A17V, R23H, G32D, G32S, R34K, V372I) of Ssa1 that influence propagation of the yeast [PSI+] prion. However, the underlining molecular mechanisms of these mutations are still unclear. The six mutation sites are present in the IA subdomain of the nucleotide-binding domain (NBD). The ATPase subdomain IA is a critical mediator of inter-domain allostery in Hsp70 molecular chaperones, so the mutation and changes in this subdomain may influence the function of the substrate-binding domain. In addition, ADP release is a rate-limiting step of the ATPase cycle and dysregulation of the ATPase cycle influences the propagation of the yeast [PSI+] prion. In this work, steered molecular dynamics (SMD) simulations were performed to explore the interaction between ADP and NBD. Results suggest that during the SMD simulations, hydrophobic interactions are predominant and variations in the binding state of ADP within the mutants is a potential reason for in vivo effects on yeast [PSI+] prion propagation. Additionally, we identify the primary residues in the ATPase domain that directly constitute the main hydrophobic interaction network and directly influence the ADP interaction state with the NBD of Ssa1. Furthermore, this in silico analysis reaffirms the importance of previously experimentally-determined residues in the Hsp70 ATPase domain involved in ADP binding and also identifies new residues potentially involved in this process.
Collapse
Affiliation(s)
- You-Lin Xue
- School of Environmental Science, College of Environment, Liaoning University, No. 66 Chongshan Middle Road, Huanggu District, Shenyang, 110036, Liaoning, China.,College of Light Industry, Liaoning University, Shenyang, 110036, China
| | - Qiaoshi Zhang
- School of Environmental Science, College of Environment, Liaoning University, No. 66 Chongshan Middle Road, Huanggu District, Shenyang, 110036, Liaoning, China
| | - Yuna Sun
- School of Environmental Science, College of Environment, Liaoning University, No. 66 Chongshan Middle Road, Huanggu District, Shenyang, 110036, Liaoning, China
| | - Xiaohong Zhou
- School of Environmental Science, College of Environment, Liaoning University, No. 66 Chongshan Middle Road, Huanggu District, Shenyang, 110036, Liaoning, China
| | - Ian P Hurley
- Centre for Biomedical Science Research, School of Clinical and Applied Sciences, Faculty of Health and Social Sciences, Leeds Beckett University, Leeds, LS13HE, UK
| | - Gary W Jones
- Centre for Biomedical Science Research, School of Clinical and Applied Sciences, Faculty of Health and Social Sciences, Leeds Beckett University, Leeds, LS13HE, UK
| | - Youtao Song
- School of Environmental Science, College of Environment, Liaoning University, No. 66 Chongshan Middle Road, Huanggu District, Shenyang, 110036, Liaoning, China.
| |
Collapse
|
21
|
Ferraro M, D’Annessa I, Moroni E, Morra G, Paladino A, Rinaldi S, Compostella F, Colombo G. Allosteric Modulators of HSP90 and HSP70: Dynamics Meets Function through Structure-Based Drug Design. J Med Chem 2018; 62:60-87. [DOI: 10.1021/acs.jmedchem.8b00825] [Citation(s) in RCA: 63] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Mariarosaria Ferraro
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | - Ilda D’Annessa
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | | | - Giulia Morra
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | - Antonella Paladino
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | - Silvia Rinaldi
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | - Federica Compostella
- Dipartimento di Biotecnologie Mediche e Medicina Traslazionale, Università degli Studi di Milano, Via Saldini, 50, 20133 Milano, Italy
| | - Giorgio Colombo
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
- Dipartimento di Chimica, Università di Pavia, V.le Taramelli 12, 27100 Pavia, Italy
| |
Collapse
|
22
|
Lubecka EA, Liwo A. A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. II. Backbone-local potentials of coarse-grained O1→4-bonded polyglucose chains. J Chem Phys 2017; 147:115101. [DOI: 10.1063/1.4994130] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Emilia A. Lubecka
- Institute of Informatics, University of Gdańsk, Wita Stwosza 57, 80-308 Gdańsk, Poland
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
23
|
Krupa P, Hałabis A, Żmudzińska W, Ołdziej S, Scheraga HA, Liwo A. Maximum Likelihood Calibration of the UNRES Force Field for Simulation of Protein Structure and Dynamics. J Chem Inf Model 2017; 57:2364-2377. [DOI: 10.1021/acs.jcim.7b00254] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland
| | - Anna Hałabis
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Wioletta Żmudzińska
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Stanisław Ołdziej
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Harold A. Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| |
Collapse
|
24
|
Penkler D, Sensoy Ö, Atilgan C, Tastan Bishop Ö. Perturbation-Response Scanning Reveals Key Residues for Allosteric Control in Hsp70. J Chem Inf Model 2017; 57:1359-1374. [PMID: 28505454 DOI: 10.1021/acs.jcim.6b00775] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Hsp70 molecular chaperones play an important role in maintaining cellular homeostasis, and are implicated in a wide array of cellular processes, including protein recovery from aggregates, cross-membrane protein translocation, and protein biogenesis. Hsp70 consists of two domains, a nucleotide binding domain (NBD) and a substrate binding domain (SBD), each of which communicates via an allosteric mechanism such that the protein interconverts between two functional states, an ATP-bound open conformation and an ADP-bound closed conformation. The exact mechanism for interstate conversion is not as yet fully understood. However, the ligand-bound states of the NBD and SBD as well as interactions with cochaperones such as DnaJ and nucleotide exchange factor are thought to play crucial regulatory roles. In this study, we apply the perturbation-response scanning (PRS) method in combination with molecular dynamics simulations as a computational tool for the identification of allosteric hot residues in the large multidomain Hsp70 protein. We find evidence in support of the hypothesis that substrate binding triggers ATP hydrolysis and that the ADP-substrate complex favors interstate conversion to the closed state. Furthermore, our data are in agreement with the proposal that there is an allosterically active intermediate state between the open and closed states and vice versa, as we find evidence that ATP binding to the closed structure and peptide binding to the open structure allosterically "activate" the respective complexes. We conclude our analysis by showing how our PRS data fit the current opinion on the Hsp70 conformational cycle and present several allosteric hot residues that may provide a platform for further studies to gain additional insight into Hsp70 allostery.
Collapse
Affiliation(s)
- David Penkler
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown 6140, South Africa
| | - Özge Sensoy
- School of Engineering and Natural Sciences, Istanbul Medipol University , Beykoz 34810, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University , Tuzla 34956, Istanbul, Turkey
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown 6140, South Africa
| |
Collapse
|
25
|
Nanomechanics of the substrate binding domain of Hsp70 determine its allosteric ATP-induced conformational change. Proc Natl Acad Sci U S A 2017; 114:6040-6045. [PMID: 28533394 DOI: 10.1073/pnas.1619843114] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Owing to the cooperativity of protein structures, it is often almost impossible to identify independent subunits, flexible regions, or hinges simply by visual inspection of static snapshots. Here, we use single-molecule force experiments and simulations to apply tension across the substrate binding domain (SBD) of heat shock protein 70 (Hsp70) to pinpoint mechanical units and flexible hinges. The SBD consists of two nanomechanical units matching 3D structural parts, called the α- and β-subdomain. We identified a flexible region within the rigid β-subdomain that gives way under load, thus opening up the α/β interface. In exactly this region, structural changes occur in the ATP-induced opening of Hsp70 to allow substrate exchange. Our results show that the SBD's ability to undergo large conformational changes is already encoded by passive mechanics of the individual elements.
Collapse
|
26
|
Stetz G, Verkhivker GM. Computational Analysis of Residue Interaction Networks and Coevolutionary Relationships in the Hsp70 Chaperones: A Community-Hopping Model of Allosteric Regulation and Communication. PLoS Comput Biol 2017; 13:e1005299. [PMID: 28095400 PMCID: PMC5240922 DOI: 10.1371/journal.pcbi.1005299] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2016] [Accepted: 12/06/2016] [Indexed: 12/28/2022] Open
Abstract
Allosteric interactions in the Hsp70 proteins are linked with their regulatory mechanisms and cellular functions. Despite significant progress in structural and functional characterization of the Hsp70 proteins fundamental questions concerning modularity of the allosteric interaction networks and hierarchy of signaling pathways in the Hsp70 chaperones remained largely unexplored and poorly understood. In this work, we proposed an integrated computational strategy that combined atomistic and coarse-grained simulations with coevolutionary analysis and network modeling of the residue interactions. A novel aspect of this work is the incorporation of dynamic residue correlations and coevolutionary residue dependencies in the construction of allosteric interaction networks and signaling pathways. We found that functional sites involved in allosteric regulation of Hsp70 may be characterized by structural stability, proximity to global hinge centers and local structural environment that is enriched by highly coevolving flexible residues. These specific characteristics may be necessary for regulation of allosteric structural transitions and could distinguish regulatory sites from nonfunctional conserved residues. The observed confluence of dynamics correlations and coevolutionary residue couplings with global networking features may determine modular organization of allosteric interactions and dictate localization of key mediating sites. Community analysis of the residue interaction networks revealed that concerted rearrangements of local interacting modules at the inter-domain interface may be responsible for global structural changes and a population shift in the DnaK chaperone. The inter-domain communities in the Hsp70 structures harbor the majority of regulatory residues involved in allosteric signaling, suggesting that these sites could be integral to the network organization and coordination of structural changes. Using a network-based formalism of allostery, we introduced a community-hopping model of allosteric communication. Atomistic reconstruction of signaling pathways in the DnaK structures captured a direction-specific mechanism and molecular details of signal transmission that are fully consistent with the mutagenesis experiments. The results of our study reconciled structural and functional experiments from a network-centric perspective by showing that global properties of the residue interaction networks and coevolutionary signatures may be linked with specificity and diversity of allosteric regulation mechanisms. The diversity of allosteric mechanisms in the Hsp70 proteins could range from modulation of the inter-domain interactions and conformational dynamics to fine-tuning of the Hsp70 interactions with co-chaperones. The goal of this study is to present a systematic computational analysis of the dynamic and evolutionary factors underlying allosteric structural transformations of the Hsp70 proteins. We investigated the relationship between functional dynamics, residue coevolution, and network organization of residue interactions in the Hsp70 proteins. The results of this study revealed that conformational dynamics of the Hsp70 proteins may be linked with coevolutionary propensities and mutual information dependencies of the protein residues. Modularity and connectivity of allosteric interactions in the Hsp70 chaperones are coordinated by stable functional sites that feature unique coevolutionary signatures and high network centrality. The emergence of the inter-domain communities that are coordinated by functional centers and include highly coevolving residues could facilitate structural transitions through cooperative reorganization of the local interacting modules. We determined that the differences in the modularity of the residue interactions and organization of coevolutionary networks in DnaK may be associated with variations in their allosteric mechanisms. The network signatures of the DnaK structures are characteristic of a population-shift allostery that allows for coordinated structural rearrangements of local communities. A dislocation of mediating centers and insufficient coevolutionary coupling between functional regions may render a reduced cooperativity and promote a limited entropy-driven allostery in the Sse1 chaperone that occurs without structural changes. The results of this study showed that a network-centric framework and a community-hopping model of allosteric communication pathways may provide novel insights into molecular and evolutionary principles of allosteric regulation in the Hsp70 proteins.
Collapse
Affiliation(s)
- Gabrielle Stetz
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
| | - Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
- Chapman University School of Pharmacy, Irvine, California, United States of America
- * E-mail:
| |
Collapse
|
27
|
Mozolewska MA, Krupa P, Zaborowski B, Liwo A, Lee J, Joo K, Czaplewski C. Use of Restraints from Consensus Fragments of Multiple Server Models To Enhance Protein-Structure Prediction Capability of the UNRES Force Field. J Chem Inf Model 2016; 56:2263-2279. [DOI: 10.1021/acs.jcim.6b00189] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
| | - Paweł Krupa
- Faculty
of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | | | - Adam Liwo
- Faculty
of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Jooyoung Lee
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Keehyoung Joo
- Center
for Advanced Computation, Korea Institute for Advanced Study, 85
Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Cezary Czaplewski
- Faculty
of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
28
|
Mozolewska MA, Sieradzan AK, Niadzvedstki A, Czaplewski C, Liwo A, Krupa P. Role of the sulfur to α-carbon thioether bridges in thurincin H. J Biomol Struct Dyn 2016; 35:2868-2879. [PMID: 27615507 DOI: 10.1080/07391102.2016.1234414] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Thurincin H is a small protein produced by Bacillus thuringiensis SF361 with gram-positive antimicrobial properties. The toxins produced by B. thuringiensis are widely used in the agriculture as, e.g. natural preservatives in dairy products. The structure of thurincin H possesses four covalent sulfur to [Formula: see text]-carbon bonds that involve the cysteine side-chains; these bonds are probably responsible for the shape and stability of the protein and, thereby, for its antimicrobial properties. To examine the influence of the formation of the sulfur-carbon bonds on the folding pathways and stability of the protein, a series of canonical and multiplexed replica-exchange simulations with the coarse-grained UNRES force field was carried out without and with distance restraints imposed on selected S-C[Formula: see text] atom pairs. It was found that the order of the formation and breaking of the S-C[Formula: see text] thioether bonds significantly impacts on the foldability and stability of the thurincin H. It was also observed that thioether bridges play a major role in stabilizing the global fold of the protein, although it significantly diminishes the entropy of the system. The maximum foldability of thurincin H was observed in the presence of the optimal set of three out of four thioether bridges. Thus, the results suggest that the presence of ThnB enzyme and other agents that catalyze the formation of thioether bridges can be essential for correct folding of thurincin H and that the formation of the fourth bridge does not seem to facilitate folding; instead, it seems to rigidify the loop and prevent proteolysis.
Collapse
Affiliation(s)
| | - Adam K Sieradzan
- a Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk , Poland
| | - Andrei Niadzvedstki
- a Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk , Poland
| | - Cezary Czaplewski
- a Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk , Poland
| | - Adam Liwo
- a Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk , Poland
| | - Paweł Krupa
- a Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk , Poland
| |
Collapse
|
29
|
Alderson TR, Kim JH, Markley JL. Dynamical Structures of Hsp70 and Hsp70-Hsp40 Complexes. Structure 2016; 24:1014-30. [PMID: 27345933 DOI: 10.1016/j.str.2016.05.011] [Citation(s) in RCA: 77] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Revised: 05/05/2016] [Accepted: 05/10/2016] [Indexed: 12/25/2022]
Abstract
Protein misfolding and aggregation are pathological events that place a significant amount of stress on the maintenance of protein homeostasis (proteostasis). For prevention and repair of protein misfolding and aggregation, cells are equipped with robust mechanisms that mainly rely on molecular chaperones. Two classes of molecular chaperones, heat shock protein 70 kDa (Hsp70) and Hsp40, recognize and bind to misfolded proteins, preventing their toxic biomolecular aggregation and enabling refolding or targeted degradation. Here, we review the current state of structural biology of Hsp70 and Hsp40-Hsp70 complexes and examine the link between their structures, dynamics, and functions. We highlight the power of nuclear magnetic resonance spectroscopy to untangle complex relationships behind molecular chaperones and their mechanism(s) of action.
Collapse
Affiliation(s)
- Thomas Reid Alderson
- Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3TA, UK; Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jin Hae Kim
- National Magnetic Resonance Facility at Madison, Biochemistry Department, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - John Lute Markley
- National Magnetic Resonance Facility at Madison, Biochemistry Department, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
30
|
Kmiecik S, Gront D, Kolinski M, Wieteska L, Dawid AE, Kolinski A. Coarse-Grained Protein Models and Their Applications. Chem Rev 2016; 116:7898-936. [DOI: 10.1021/acs.chemrev.6b00163] [Citation(s) in RCA: 555] [Impact Index Per Article: 69.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Sebastian Kmiecik
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dominik Gront
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michal Kolinski
- Bioinformatics
Laboratory, Mossakowski Medical Research Center of the Polish Academy of Sciences, Pawinskiego 5, 02-106 Warsaw, Poland
| | - Lukasz Wieteska
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
- Department
of Medical Biochemistry, Medical University of Lodz, Mazowiecka 6/8, 92-215 Lodz, Poland
| | | | - Andrzej Kolinski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
31
|
Lipska AG, Seidman SR, Sieradzan AK, Giełdoń A, Liwo A, Scheraga HA. Molecular dynamics of protein A and a WW domain with a united-residue model including hydrodynamic interaction. J Chem Phys 2016; 144:184110. [PMID: 27179474 PMCID: PMC4866947 DOI: 10.1063/1.4948710] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2016] [Accepted: 04/25/2016] [Indexed: 01/01/2023] Open
Abstract
The folding of the N-terminal part of the B-domain of staphylococcal protein A (PDB ID: 1BDD, a 46-residue three-α-helix bundle) and the formin-binding protein 28 WW domain (PDB ID: 1E0L, a 37-residue three-stranded anti-parallel β protein) was studied by means of Langevin dynamics with the coarse-grained UNRES force field to assess the influence of hydrodynamic interactions on protein-folding pathways and kinetics. The unfolded, intermediate, and native-like structures were identified by cluster analysis, and multi-exponential functions were fitted to the time dependence of the fractions of native and intermediate structures, respectively, to determine bulk kinetics. It was found that introducing hydrodynamic interactions slows down both the formation of an intermediate state and the transition from the collapsed structures to the final native-like structures by creating multiple kinetic traps. Therefore, introducing hydrodynamic interactions considerably slows the folding, as opposed to the results obtained from earlier studies with the use of Gō-like models.
Collapse
Affiliation(s)
- Agnieszka G Lipska
- Laboratory of Molecular Modeling, Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Steven R Seidman
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | - Adam K Sieradzan
- Laboratory of Molecular Modeling, Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Artur Giełdoń
- Laboratory of Molecular Modeling, Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Laboratory of Molecular Modeling, Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| |
Collapse
|
32
|
Elengoe A, Naser MA, Hamdan S. Molecular dynamics simulation and docking studies on novel mutants (T11V, T12P and D364S) of the nucleotide-binding domain of human heat shock 70 kDa protein. Biologia (Bratisl) 2015. [DOI: 10.1515/biolog-2015-0194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
33
|
Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins. PLoS One 2015; 10:e0143752. [PMID: 26619280 PMCID: PMC4664246 DOI: 10.1371/journal.pone.0143752] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2015] [Accepted: 11/09/2015] [Indexed: 01/04/2023] Open
Abstract
Hsp70 and Hsp110 chaperones play an important role in regulating cellular processes that involve protein folding and stabilization, which are essential for the integrity of signaling networks. Although many aspects of allosteric regulatory mechanisms in Hsp70 and Hsp110 chaperones have been extensively studied and significantly advanced in recent experimental studies, the atomistic picture of signal propagation and energetics of dynamics-based communication still remain unresolved. In this work, we have combined molecular dynamics simulations and protein stability analysis of the chaperone structures with the network modeling of residue interaction networks to characterize molecular determinants of allosteric mechanisms. We have shown that allosteric mechanisms of Hsp70 and Hsp110 chaperones may be primarily determined by nucleotide-induced redistribution of local conformational ensembles in the inter-domain regions and the substrate binding domain. Conformational dynamics and energetics of the peptide substrate binding with the Hsp70 structures has been analyzed using free energy calculations, revealing allosteric hotspots that control negative cooperativity between regulatory sites. The results have indicated that cooperative interactions may promote a population-shift mechanism in Hsp70, in which functional residues are organized in a broad and robust allosteric network that can link the nucleotide-binding site and the substrate-binding regions. A smaller allosteric network in Hsp110 structures may elicit an entropy-driven allostery that occurs in the absence of global structural changes. We have found that global mediating residues with high network centrality may be organized in stable local communities that are indispensable for structural stability and efficient allosteric communications. The network-centric analysis of allosteric interactions has also established that centrality of functional residues could correlate with their sensitivity to mutations across diverse chaperone functions. This study reconciles a wide spectrum of structural and functional experiments by demonstrating how integration of molecular simulations and network-centric modeling may explain thermodynamic and mechanistic aspects of allosteric regulation in chaperones.
Collapse
|
34
|
Zaborowski B, Jagieła D, Czaplewski C, Hałabis A, Lewandowska A, Żmudzińska W, Ołdziej S, Karczyńska A, Omieczynski C, Wirecki T, Liwo A. A Maximum-Likelihood Approach to Force-Field Calibration. J Chem Inf Model 2015; 55:2050-70. [DOI: 10.1021/acs.jcim.5b00395] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Bartłomiej Zaborowski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Dawid Jagieła
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Anna Hałabis
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Agnieszka Lewandowska
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Wioletta Żmudzińska
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Stanisław Ołdziej
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Agnieszka Karczyńska
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Christian Omieczynski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Tomasz Wirecki
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 87 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| |
Collapse
|
35
|
Mozolewska MA, Krupa P, Scheraga HA, Liwo A. Molecular modeling of the binding modes of the iron-sulfur protein to the Jac1 co-chaperone from Saccharomyces cerevisiae by all-atom and coarse-grained approaches. Proteins 2015; 83:1414-26. [PMID: 25973573 DOI: 10.1002/prot.24824] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2014] [Revised: 04/18/2015] [Accepted: 04/29/2015] [Indexed: 11/06/2022]
Abstract
The iron-sulfur protein 1 (Isu1) and the J-type co-chaperone Jac1 from yeast are part of a huge ATP-dependent system, and both interact with Hsp70 chaperones. Interaction of Isu1 and Jac1 is a part of the iron-sulfur cluster biogenesis system in mitochondria. In this study, the structure and dynamics of the yeast Isu1-Jac1 complex has been modeled. First, the complete structure of Isu1 was obtained by homology modeling using the I-TASSER server and YASARA software and thereafter tested for stability in the all-atom force field AMBER. Then, the known experimental structure of Jac1 was adopted to obtain initial models of the Isu1-Jac1 complex by using the ZDOCK server for global and local docking and the AutoDock software for local docking. Three most probable models were subsequently subjected to the coarse-grained molecular dynamics simulations with the UNRES force field to obtain the final structures of the complex. In the most probable model, Isu1 binds to the left face of the Γ-shaped Jac1 molecule by the β-sheet section of Isu1. Residues L105 , L109 , and Y163 of Jac1 have been assessed by mutation studies to be essential for binding (Ciesielski et al., J Mol Biol 2012; 417:1-12). These residues were also found, by UNRES/molecular dynamics simulations, to be involved in strong interactions between Isu1 and Jac1 in the complex. Moreover, N(95), T(98), P(102), H(112), V(159), L(167), and A(170) of Jac1, not yet tested experimentally, were also found to be important in binding.
Collapse
Affiliation(s)
- Magdalena A Mozolewska
- Faculty of Chemistry, University of Gdansk, Gdansk, 80-308, Poland.,Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301
| | - Paweł Krupa
- Faculty of Chemistry, University of Gdansk, Gdansk, 80-308, Poland.,Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Gdansk, 80-308, Poland
| |
Collapse
|
36
|
Abstract
This is a tour of a physical chemist through 65 years of protein chemistry from the time when emphasis was placed on the determination of the size and shape of the protein molecule as a colloidal particle, with an early breakthrough by James Sumner, followed by Linus Pauling and Fred Sanger, that a protein was a real molecule, albeit a macromolecule. It deals with the recognition of the nature and importance of hydrogen bonds and hydrophobic interactions in determining the structure, properties, and biological function of proteins until the present acquisition of an understanding of the structure, thermodynamics, and folding pathways from a linear array of amino acids to a biological entity. Along the way, with a combination of experiment and theoretical interpretation, a mechanism was elucidated for the thrombin-induced conversion of fibrinogen to a fibrin blood clot and for the oxidative-folding pathways of ribonuclease A. Before the atomic structure of a protein molecule was determined by x-ray diffraction or nuclear magnetic resonance spectroscopy, experimental studies of the fundamental interactions underlying protein structure led to several distance constraints which motivated the theoretical approach to determine protein structure, and culminated in the Empirical Conformational Energy Program for Peptides (ECEPP), an all-atom force field, with which the structures of fibrous collagen-like proteins and the 46-residue globular staphylococcal protein A were determined. To undertake the study of larger globular proteins, a physics-based coarse-grained UNited-RESidue (UNRES) force field was developed, and applied to the protein-folding problem in terms of structure, thermodynamics, dynamics, and folding pathways. Initially, single-chain and, ultimately, multiple-chain proteins were examined, and the methodology was extended to protein-protein interactions and to nucleic acids and to protein-nucleic acid interactions. The ultimate results led to an understanding of a variety of biological processes underlying natural and disease phenomena.
Collapse
|
37
|
Lipska AG, Sieradzan AK, Krupa P, Mozolewska MA, D’Auria S, Liwo A. Studies of conformational changes of an arginine-binding protein from Thermotoga maritima in the presence and absence of ligand via molecular dynamics simulations with the coarse-grained UNRES force field. J Mol Model 2015; 21:64. [DOI: 10.1007/s00894-015-2609-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Accepted: 02/08/2015] [Indexed: 11/30/2022]
|
38
|
Sieradzan AK, Krupa P, Scheraga HA, Liwo A, Czaplewski C. Physics-based potentials for the coupling between backbone- and side-chain-local conformational states in the UNited RESidue (UNRES) force field for protein simulations. J Chem Theory Comput 2015; 11:817-31. [PMID: 25691834 PMCID: PMC4327884 DOI: 10.1021/ct500736a] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The UNited RESidue (UNRES) model of polypeptide chains is a coarse-grained model in which each amino-acid residue is reduced to two interaction sites, namely, a united peptide group (p) located halfway between the two neighboring α-carbon atoms (Cαs), which serve only as geometrical points, and a united side chain (SC) attached to the respective Cα. Owing to this simplification, millisecond molecular dynamics simulations of large systems can be performed. While UNRES predicts overall folds well, it reproduces the details of local chain conformation with lower accuracy. Recently, we implemented new knowledge-based torsional potentials (Krupa et al. J. Chem. Theory Comput. 2013, 9, 4620–4632) that depend on the virtual-bond dihedral angles involving side chains: Cα···Cα···Cα···SC (τ(1)), SC···Cα···Cα···Cα (τ(2)), and SC···Cα···Cα···SC (τ(3)) in the UNRES force field. These potentials resulted in significant improvement of the simulated structures, especially in the loop regions. In this work, we introduce the physics-based counterparts of these potentials, which we derived from the all-atom energy surfaces of terminally blocked amino-acid residues by Boltzmann integration over the angles λ(1) and λ(2) for rotation about the Cα···Cα virtual-bond angles and over the side-chain angles χ. The energy surfaces were, in turn, calculated by using the semiempirical AM1 method of molecular quantum mechanics. Entropy contribution was evaluated with use of the harmonic approximation from Hessian matrices. One-dimensional Fourier series in the respective virtual-bond-dihedral angles were fitted to the calculated potentials, and these expressions have been implemented in the UNRES force field. Basic calibration of the UNRES force field with the new potentials was carried out with eight training proteins, by selecting the optimal weight of the new energy terms and reducing the weight of the regular torsional terms. The force field was subsequently benchmarked with a set of 22 proteins not used in the calibration. The new potentials result in a decrease of the root-mean-square deviation of the average conformation from the respective experimental structure by 0.86 Å on average; however, improvement of up to 5 Å was observed for some proteins.
Collapse
Affiliation(s)
- Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Harold A. Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
| |
Collapse
|
39
|
Gołaś EI, Czaplewski C, Scheraga HA, Liwo A. Common functionally important motions of the nucleotide-binding domain of Hsp70. Proteins 2014; 83:282-99. [PMID: 25412765 DOI: 10.1002/prot.24731] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2014] [Revised: 11/05/2014] [Accepted: 11/13/2014] [Indexed: 01/14/2023]
Abstract
The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate-binding domain (SBD) that binds client substrates, and the nucleotide-binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure-function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone from Bos taurus (PDB 3C7N:B) by all-atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP- and ATP-unique classes, which reflect conformational trends that are unique to either the ADP- or ATP-bound states, respectively. "Mutual" class motions generally describe "in-plane" and/or "out-of-plane" (scissor-like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The "unique" class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the "unique" type, regions of enhanced mobility can be identified; these are termed "hot spots," and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide-binding pocket was also found to influence the dynamics of the NBD significantly.
Collapse
Affiliation(s)
- Ewa I Gołaś
- Faculty of Chemistry, University of Gdańsk, 80-308, Gdańsk, Poland
| | | | | | | |
Collapse
|
40
|
Liwo A, Baranowski M, Czaplewski C, Gołaś E, He Y, Jagieła D, Krupa P, Maciejczyk M, Makowski M, Mozolewska MA, Niadzvedtski A, Ołdziej S, Scheraga HA, Sieradzan AK, Slusarz R, Wirecki T, Yin Y, Zaborowski B. A unified coarse-grained model of biological macromolecules based on mean-field multipole-multipole interactions. J Mol Model 2014; 20:2306. [PMID: 25024008 PMCID: PMC4139597 DOI: 10.1007/s00894-014-2306-5] [Citation(s) in RCA: 111] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2014] [Accepted: 05/12/2014] [Indexed: 12/14/2022]
Abstract
A unified coarse-grained model of three major classes of biological molecules—proteins, nucleic acids, and polysaccharides—has been developed. It is based on the observations that the repeated units of biopolymers (peptide groups, nucleic acid bases, sugar rings) are highly polar and their charge distributions can be represented crudely as point multipoles. The model is an extension of the united residue (UNRES) coarse-grained model of proteins developed previously in our laboratory. The respective force fields are defined as the potentials of mean force of biomacromolecules immersed in water, where all degrees of freedom not considered in the model have been averaged out. Reducing the representation to one center per polar interaction site leads to the representation of average site–site interactions as mean-field dipole–dipole interactions. Further expansion of the potentials of mean force of biopolymer chains into Kubo’s cluster-cumulant series leads to the appearance of mean-field dipole–dipole interactions, averaged in the context of local interactions within a biopolymer unit. These mean-field interactions account for the formation of regular structures encountered in biomacromolecules, e.g., α-helices and β-sheets in proteins, double helices in nucleic acids, and helicoidally packed structures in polysaccharides, which enables us to use a greatly reduced number of interacting sites without sacrificing the ability to reproduce the correct architecture. This reduction results in an extension of the simulation timescale by more than four orders of magnitude compared to the all-atom representation. Examples of the performance of the model are presented. Components of the Unified Coarse Grained Model (UCGM) of biological macromolecules ![]()
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308, Gdańsk, Poland,
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
41
|
Ha-Duong T. Coarse-grained models of the proteins backbone conformational dynamics. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 805:157-69. [PMID: 24446361 DOI: 10.1007/978-3-319-02970-2_7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Coarse-grained models are more and more frequently used in the studies of the proteins structural and dynamic properties, since the reduced number of degrees of freedom allows to enhance the conformational space exploration. This chapter attempts to provide an overview of the various coarse-grained models that were applied to study the functional conformational changes of the polypeptides main chain around their native state. It will more specifically discuss the methods used to represent the protein backbone flexibility and to account for the physico-chemical interactions that stabilize the secondary structure elements.
Collapse
Affiliation(s)
- Tap Ha-Duong
- BIOCIS - UMR CNRS 8076, Faculté de Pharmacie - Université Paris Sud, 5 rue Jean-Baptiste Clément, 92296, Châtenay-Malabry, France,
| |
Collapse
|
42
|
General IJ, Liu Y, Blackburn ME, Mao W, Gierasch LM, Bahar I. ATPase subdomain IA is a mediator of interdomain allostery in Hsp70 molecular chaperones. PLoS Comput Biol 2014; 10:e1003624. [PMID: 24831085 PMCID: PMC4022485 DOI: 10.1371/journal.pcbi.1003624] [Citation(s) in RCA: 79] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2013] [Accepted: 03/31/2014] [Indexed: 11/18/2022] Open
Abstract
The versatile functions of the heat shock protein 70 (Hsp70) family of molecular chaperones rely on allosteric interactions between their nucleotide-binding and substrate-binding domains, NBD and SBD. Understanding the mechanism of interdomain allostery is essential to rational design of Hsp70 modulators. Yet, despite significant progress in recent years, how the two Hsp70 domains regulate each other's activity remains elusive. Covariance data from experiments and computations emerged in recent years as valuable sources of information towards gaining insights into the molecular events that mediate allostery. In the present study, conservation and covariance properties derived from both sequence and structural dynamics data are integrated with results from Perturbation Response Scanning and in vivo functional assays, so as to establish the dynamical basis of interdomain signal transduction in Hsp70s. Our study highlights the critical roles of SBD residues D481 and T417 in mediating the coupled motions of the two domains, as well as that of G506 in enabling the movements of the α-helical lid with respect to the β-sandwich. It also draws attention to the distinctive role of the NBD subdomains: Subdomain IA acts as a key mediator of signal transduction between the ATP- and substrate-binding sites, this function being achieved by a cascade of interactions predominantly involving conserved residues such as V139, D148, R167 and K155. Subdomain IIA, on the other hand, is distinguished by strong coevolutionary signals (with the SBD) exhibited by a series of residues (D211, E217, L219, T383) implicated in DnaJ recognition. The occurrence of coevolving residues at the DnaJ recognition region parallels the behavior recently observed at the nucleotide-exchange-factor recognition region of subdomain IIB. These findings suggest that Hsp70 tends to adapt to co-chaperone recognition and activity via coevolving residues, whereas interdomain allostery, critical to chaperoning, is robustly enabled by conserved interactions.
Collapse
Affiliation(s)
- Ignacio J. General
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Ying Liu
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Mandy E. Blackburn
- Department of Biochemistry & Molecular Biology, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
| | - Wenzhi Mao
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Department of Pharmacology, Tsinghua University, Beijing, China
| | - Lila M. Gierasch
- Department of Biochemistry & Molecular Biology, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- * E-mail:
| |
Collapse
|
43
|
Krokhotin A, Liwo A, Maisuradze GG, Niemi AJ, Scheraga HA. Kinks, loops, and protein folding, with protein A as an example. J Chem Phys 2014; 140:025101. [PMID: 24437917 PMCID: PMC3899063 DOI: 10.1063/1.4855735] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2013] [Accepted: 12/10/2013] [Indexed: 02/02/2023] Open
Abstract
The dynamics and energetics of formation of loops in the 46-residue N-terminal fragment of the B-domain of staphylococcal protein A has been studied. Numerical simulations have been performed using coarse-grained molecular dynamics with the united-residue (UNRES) force field. The results have been analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger (DNLS) equation. In the case of proteins, the DNLS equation arises from a C(α)-trace-based energy function. Three individual kink profiles were identified in the experimental three-α-helix structure of protein A, in the range of the Glu16-Asn29, Leu20-Asn29, and Gln33-Asn44 residues, respectively; these correspond to two loops in the native structure. UNRES simulations were started from the full right-handed α-helix to obtain a clear picture of kink formation, which would otherwise be blurred by helix formation. All three kinks emerged during coarse-grained simulations. It was found that the formation of each is accompanied by a local free energy increase; this is expressed as the change of UNRES energy which has the physical sense of the potential of mean force of a polypeptide chain. The increase is about 7 kcal/mol. This value can thus be considered as the free energy barrier to kink formation in full α-helical segments of polypeptide chains. During the simulations, the kinks emerge, disappear, propagate, and annihilate each other many times. It was found that the formation of a kink is initiated by an abrupt change in the orientation of a pair of consecutive side chains in the loop region. This resembles the formation of a Bloch wall along a spin chain, where the C(α) backbone corresponds to the chain, and the amino acid side chains are interpreted as the spin variables. This observation suggests that nearest-neighbor side chain-side chain interactions are responsible for initiation of loop formation. It was also found that the individual kinks are reflected as clear peaks in the principal modes of the analyzed trajectory of protein A, the shapes of which resemble the directional derivatives of the kinks along the chain. These observations suggest that the kinks of the DNLS equation determine the functionally important motions of proteins.
Collapse
Affiliation(s)
- Andrey Krokhotin
- Department of Physics and Astronomy and Science for Life Laboratory, Uppsala University, P.O. Box 803, S-75108 Uppsala, Sweden
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, ul. Sobieskiego 18, 80-952 Gdansk, Poland
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | - Antti J Niemi
- Department of Physics and Astronomy and Science for Life Laboratory, Uppsala University, P.O. Box 803, S-75108 Uppsala, Sweden
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| |
Collapse
|
44
|
Low-Frequency, Functional, Modes of Proteins: All-Atom and Coarse-Grained Normal Mode Analysis. COMPUTATIONAL METHODS TO STUDY THE STRUCTURE AND DYNAMICS OF BIOMOLECULES AND BIOMOLECULAR PROCESSES 2014. [DOI: 10.1007/978-3-642-28554-7_15] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
45
|
Hills RD. Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics. Methods Mol Biol 2014; 1084:123-140. [PMID: 24061919 DOI: 10.1007/978-1-62703-658-0_7] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Pharmaceutical Sciences, University of New England, Portland, ME, USA
| |
Collapse
|
46
|
Spiga E, Degiacomi MT, Dal Peraro M. New Strategies for Integrative Dynamic Modeling of Macromolecular Assembly. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 96:77-111. [DOI: 10.1016/bs.apcsb.2014.06.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
47
|
Decipher the mechanisms of protein conformational changes induced by nucleotide binding through free-energy landscape analysis: ATP binding to Hsp70. PLoS Comput Biol 2013; 9:e1003379. [PMID: 24348227 PMCID: PMC3861046 DOI: 10.1371/journal.pcbi.1003379] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Accepted: 10/15/2013] [Indexed: 11/30/2022] Open
Abstract
ATP regulates the function of many proteins in the cell by transducing its binding and hydrolysis energies into protein conformational changes by mechanisms which are challenging to identify at the atomic scale. Based on molecular dynamics (MD) simulations, a method is proposed to analyze the structural changes induced by ATP binding to a protein by computing the effective free-energy landscape (FEL) of a subset of its coordinates along its amino-acid sequence. The method is applied to characterize the mechanism by which the binding of ATP to the nucleotide-binding domain (NBD) of Hsp70 propagates a signal to its substrate-binding domain (SBD). Unbiased MD simulations were performed for Hsp70-DnaK chaperone in nucleotide-free, ADP-bound and ATP-bound states. The simulations revealed that the SBD does not interact with the NBD for DnaK in its nucleotide-free and ADP-bound states whereas the docking of the SBD was found in the ATP-bound state. The docked state induced by ATP binding found in MD is an intermediate state between the initial nucleotide-free and final ATP-bound states of Hsp70. The analysis of the FEL projected along the amino-acid sequence permitted to identify a subset of 27 protein internal coordinates corresponding to a network of 91 key residues involved in the conformational change induced by ATP binding. Among the 91 residues, 26 are identified for the first time, whereas the others were shown relevant for the allosteric communication of Hsp70 s in several experiments and bioinformatics analysis. The FEL analysis revealed also the origin of the ATP-induced structural modifications of the SBD recently measured by Electron Paramagnetic Resonance. The pathway between the nucleotide-free and the intermediate state of DnaK was extracted by applying principal component analysis to the subset of internal coordinates describing the transition. The methodology proposed is general and could be applied to analyze allosteric communication in other proteins. The precise biophysical characterization of the mechanisms of the protein conformational changes controlled by a nucleotide remains a challenge in biology. Molecular dynamics simulations of proteins in different nucleotide-binding states contain information on the nucleotide-dependent conformational dynamics. However, it is difficult to extract relevant information about the conformation-induced mechanism from the raw molecular dynamics data. Herein, we addressed this issue for the major ATP-dependent molecular chaperones Hsp70 s, which contribute to crucial cellular processes and are involved in several neurodegenerative diseases and in cancer. To function, Hsp70 undergoes several conformational changes controlled by the state of its nucleotide-binding domain. We demonstrated that the analysis of the effective free-energy landscape of the protein projected along the amino-acid sequence and computed from the molecular dynamics simulations of Hsp70 in different nucleotide-binding states, holds the key to identify the key residues of the conformational induced pathway. Identification of the key residues involved in the propagation of the structural changes induced by ATP binding offer alternative druggable specific sites other than the ligand binding clefts. The methodology developed for Hsp70 is general and can be adapted to any ligand induced conformational change in proteins.
Collapse
|
48
|
Identification of key hinge residues important for nucleotide-dependent allostery in E. coli Hsp70/DnaK. PLoS Comput Biol 2013; 9:e1003279. [PMID: 24277995 PMCID: PMC3836694 DOI: 10.1371/journal.pcbi.1003279] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2013] [Accepted: 08/29/2013] [Indexed: 01/20/2023] Open
Abstract
DnaK is a molecular chaperone that has important roles in protein folding. The hydrolysis of ATP is essential to this activity, and the effects of nucleotides on the structure and function of DnaK have been extensively studied. However, the key residues that govern the conformational motions that define the apo, ATP-bound, and ADP-bound states are not entirely clear. Here, we used molecular dynamics simulations, mutagenesis, and enzymatic assays to explore the molecular basis of this process. Simulations of DnaK's nucleotide-binding domain (NBD) in the apo, ATP-bound, and ADP/Pi-bound states suggested that each state has a distinct conformation, consistent with available biochemical and structural information. The simulations further suggested that large shearing motions between subdomains I-A and II-A dominated the conversion between these conformations. We found that several evolutionally conserved residues, especially G228 and G229, appeared to function as a hinge for these motions, because they predominantly populated two distinct states depending on whether ATP or ADP/Pi was bound. Consistent with the importance of these “hinge” residues, alanine point mutations caused DnaK to have reduced chaperone activities in vitro and in vivo. Together, these results clarify how sub-domain motions communicate allostery in DnaK. DnaK belongs to the highly conserved heat shock protein 70 (Hsp70) family, a group of ATP-dependent molecular chaperones that regulates proteostasis. Studies have suggested that global movements of the subdomains in the nucleotide-binding domain (NBD) of DnaK regulate its catalytic activity. However, there is less known about the key residues involved in these subdomain motions and whether these residues might also regulate inter-domain allostery with the substrate-binding domain (SBD). To examine the motions in the NBD, dynamics simulations of DnaK's NBD in the apo, ATP-bound, and ADP/Pi-bound states were performed. Through essential dynamics and torsion angle analyses, we identified motions and highly conserved hinge residues between subdomains IIA and IIB that are likely to be important for nucleotide cycling and for communicating the nucleotide state to the SBD. Supporting this model, mutating these conserved hinge residues affected ATPase activity and chaperone functions in vitro and in bacteria, suggesting their importance in the nucleotide-dependent motions in DnaK.
Collapse
|
49
|
Nicolaï A, Delarue P, Senet P. Conformational dynamics of full-length inducible human Hsp70 derived from microsecond molecular dynamics simulations in explicit solvent. J Biomol Struct Dyn 2013; 31:1111-26. [DOI: 10.1080/07391102.2012.726190] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
50
|
Mashaghi A, Kramer G, Lamb DC, Mayer MP, Tans SJ. Chaperone Action at the Single-Molecule Level. Chem Rev 2013; 114:660-76. [DOI: 10.1021/cr400326k] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Alireza Mashaghi
- AMOLF Institute, Science Park
104, 1098 XG Amsterdam, The Netherlands
| | - Günter Kramer
- Zentrum
für Molekulare Biologie der Universität Heidelberg (ZMBH), DKFZ-ZMBH Allianz, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany
| | - Don C. Lamb
- Physical
Chemistry, Department of Chemistry, Munich Center for Integrated Protein
Science (CiPSM) and Center for Nanoscience, Ludwig-Maximilians-Universität München, Butenandtstrasse 5-13, Gerhard-Ertl-Building, 81377 Munich, Germany
| | - Matthias P. Mayer
- Zentrum
für Molekulare Biologie der Universität Heidelberg (ZMBH), DKFZ-ZMBH Allianz, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany
| | - Sander J. Tans
- AMOLF Institute, Science Park
104, 1098 XG Amsterdam, The Netherlands
| |
Collapse
|