1
|
Sandoval-Velasco M, Dudchenko O, Rodríguez JA, Pérez Estrada C, Dehasque M, Fontsere C, Mak SST, Khan R, Contessoto VG, Oliveira Junior AB, Kalluchi A, Zubillaga Herrera BJ, Jeong J, Roy RP, Christopher I, Weisz D, Omer AD, Batra SS, Shamim MS, Durand NC, O'Connell B, Roca AL, Plikus MV, Kusliy MA, Romanenko SA, Lemskaya NA, Serdyukova NA, Modina SA, Perelman PL, Kizilova EA, Baiborodin SI, Rubtsov NB, Machol G, Rath K, Mahajan R, Kaur P, Gnirke A, Garcia-Treviño I, Coke R, Flanagan JP, Pletch K, Ruiz-Herrera A, Plotnikov V, Pavlov IS, Pavlova NI, Protopopov AV, Di Pierro M, Graphodatsky AS, Lander ES, Rowley MJ, Wolynes PG, Onuchic JN, Dalén L, Marti-Renom MA, Gilbert MTP, Aiden EL. Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample. Cell 2024; 187:3541-3562.e51. [PMID: 38996487 DOI: 10.1016/j.cell.2024.06.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 03/07/2024] [Accepted: 06/03/2024] [Indexed: 07/14/2024]
Abstract
Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.
Collapse
Affiliation(s)
| | - Olga Dudchenko
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA.
| | - Juan Antonio Rodríguez
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark; Centre Nacional d'Anàlisi Genòmica, CNAG, 08028 Barcelona, Spain
| | - Cynthia Pérez Estrada
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA
| | - Marianne Dehasque
- Centre for Palaeogenetics, SE-106 91 Stockholm, Sweden; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 10405 Stockholm, Sweden; Department of Zoology, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Claudia Fontsere
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark
| | - Sarah S T Mak
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark
| | - Ruqayya Khan
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | | | | | - Achyuth Kalluchi
- Department of Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE 68198, USA
| | - Bernardo J Zubillaga Herrera
- Department of Physics, Northeastern University, Boston, MA 02115, USA; Center for Theoretical Biological Physics, Northeastern University, Boston, MA 02215, USA
| | - Jiyun Jeong
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Renata P Roy
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Departments of Biology and Physics, Texas Southern University, Houston, TX 77004, USA
| | - Ishawnia Christopher
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - David Weisz
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Arina D Omer
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Sanjit S Batra
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Muhammad S Shamim
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Neva C Durand
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Brendan O'Connell
- Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95064, USA; Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR 97239, USA
| | - Alfred L Roca
- Department of Animal Sciences and Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Maksim V Plikus
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA 92697, USA
| | - Mariya A Kusliy
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | | | - Natalya A Lemskaya
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | | | - Svetlana A Modina
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | - Polina L Perelman
- Institute of Molecular and Cellular Biology SB RAS, Novosibirsk 630090, Russia
| | - Elena A Kizilova
- Institute of Cytology and Genetics SB RAS, Novosibirsk 630090, Russia
| | | | - Nikolai B Rubtsov
- Institute of Cytology and Genetics SB RAS, Novosibirsk 630090, Russia
| | - Gur Machol
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Krisha Rath
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Ragini Mahajan
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA
| | - Parwinder Kaur
- UWA School of Agriculture and Environment, University of Western Australia, Perth, WA 6009, Australia
| | - Andreas Gnirke
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | | | - Rob Coke
- San Antonio Zoo, San Antonio, TX 78212, USA
| | | | | | - Aurora Ruiz-Herrera
- Departament de Biologia Cel·lular, Fisiologia i Immunologia and Genome Integrity and Instability Group, Institut de Biotecnologia i Biomedicina, Universitat Autònoma de Barcelona, 08193 Cerdanyola del Vallès, Spain
| | | | | | - Naryya I Pavlova
- Institute of Biological Problems of Cryolitezone SB RAS, Yakutsk 677000, Russia
| | - Albert V Protopopov
- Academy of Sciences of Sakha Republic, Yakutsk 677000, Russia; North-Eastern Federal University, Yakutsk 677027, Russia
| | - Michele Di Pierro
- Department of Physics, Northeastern University, Boston, MA 02115, USA; Center for Theoretical Biological Physics, Northeastern University, Boston, MA 02215, USA
| | | | - Eric S Lander
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Department of Systems Biology, Harvard Medical School, Boston, MA 02115, USA; Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - M Jordan Rowley
- Department of Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE 68198, USA
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA; Departments of Physics, Astronomy, & Chemistry, Rice University, Houston, TX 77005, USA
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Department of Biosciences, Rice University, Houston, TX 77005, USA; Departments of Physics, Astronomy, & Chemistry, Rice University, Houston, TX 77005, USA
| | - Love Dalén
- Centre for Palaeogenetics, SE-106 91 Stockholm, Sweden; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 10405 Stockholm, Sweden; Department of Zoology, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Marc A Marti-Renom
- Centre Nacional d'Anàlisi Genòmica, CNAG, 08028 Barcelona, Spain; Centre for Genomic Regulation, The Barcelona Institute for Science and Technology, 08003 Barcelona, Spain; ICREA, 08010 Barcelona, Spain; Universitat Pompeu Fabra, 08002 Barcelona, Spain.
| | - M Thomas P Gilbert
- Center for Evolutionary Hologenomics, University of Copenhagen, DK-1353 Copenhagen, Denmark; University Museum NTNU, 7012 Trondheim, Norway.
| | - Erez Lieberman Aiden
- The Center for Genome Architecture and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA; Center for Theoretical Biological Physics, Rice University, Houston, TX 77030, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA.
| |
Collapse
|
2
|
Viegas RG, Martins IBS, Leite VBP. Understanding the Energy Landscape of Intrinsically Disordered Protein Ensembles. J Chem Inf Model 2024; 64:4149-4157. [PMID: 38713459 DOI: 10.1021/acs.jcim.4c00080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
A substantial portion of various organisms' proteomes comprises intrinsically disordered proteins (IDPs) that lack a defined three-dimensional structure. These IDPs exhibit a diverse array of conformations, displaying remarkable spatiotemporal heterogeneity and exceptional conformational flexibility. Characterizing the structure or structural ensemble of IDPs presents significant conceptual and methodological challenges owing to the absence of a well-defined native structure. While databases such as the Protein Ensemble Database (PED) provide IDP ensembles obtained through a combination of experimental data and molecular modeling, the absence of reaction coordinates poses challenges in comprehensively understanding pertinent aspects of the system. In this study, we leverage the energy landscape visualization method (JCTC, 6482, 2019) to scrutinize four IDP ensembles sourced from PED. ELViM, a methodology that circumvents the need for a priori reaction coordinates, aids in analyzing the ensembles. The specific IDP ensembles investigated are as follows: two fragments of nucleoporin (NUL: 884-993 and NUS: 1313-1390), yeast sic 1 N-terminal (1-90), and the N-terminal SH3 domain of Drk (1-59). Utilizing ELViM enables the comprehensive validation of ensembles, facilitating the detection of potential inconsistencies in the sampling process. Additionally, it allows for identifying and characterizing the most prevalent conformations within an ensemble. Moreover, ELViM facilitates the comparative analysis of ensembles obtained under diverse conditions, thereby providing a powerful tool for investigating the functional mechanisms of IDPs.
Collapse
Affiliation(s)
- Rafael G Viegas
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, São Paulo 15.808-305, Brazil
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Ingrid B S Martins
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Vitor B P Leite
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
3
|
Pereira AF, Martínez L. Helical Content Correlations and Hydration Structures of the Folding Ensemble of the B Domain of Protein A. J Chem Inf Model 2024; 64:3350-3359. [PMID: 38566451 DOI: 10.1021/acs.jcim.3c01822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
The B domain of protein A (BdpA), a small three-helix bundle, folds on a time scale of a few microseconds with heterogeneous native and unfolded states. It is widely used as a model for understanding protein folding mechanisms. In this work, we use structure-based models (SBMs) and atomistic simulations to comprehensively investigate how BdpA folding is associated with the formation of its secondary structure. The energy landscape visualization method (ELViM) was used to characterize the pathways that connect the folded and unfolded states of BdpA as well as the sets of structures displaying specific ellipticity patterns. We show that the native state conformational diversity is due mainly to the conformational variability of helix I. Helices I, II, and III occur in a weakly correlated manner, with Spearman's rank correlation coefficients of 0.1539 (I and II), 0.1259 (I and III), and 0.2561 (II and III). These results, therefore, suggest the highest cooperativity between helices II and III. Our results allow the clustering of partially folded structures of folding of the B domain of protein A on the basis of its secondary structure, paving the way to an understanding of environmental factors in the relative stability of the basins of the folding ensemble, which are illustrated by the structural dependency of the protein hydration structures, as computed with minimum-distance distribution functions.
Collapse
Affiliation(s)
- Ander Francisco Pereira
- Institute of Chemistry and Center for Computing in Engineering & Science, Universidade Estadual de Campinas (UNICAMP), 13083-861 Campinas, SP, Brazil
| | - Leandro Martínez
- Institute of Chemistry and Center for Computing in Engineering & Science, Universidade Estadual de Campinas (UNICAMP), 13083-861 Campinas, SP, Brazil
| |
Collapse
|
4
|
Viegas RG, Martins IBS, Sanches MN, Oliveira Junior AB, Camargo JBD, Paulovich FV, Leite VBP. ELViM: Exploring Biomolecular Energy Landscapes through Multidimensional Visualization. J Chem Inf Model 2024; 64:3443-3450. [PMID: 38506664 DOI: 10.1021/acs.jcim.4c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024]
Abstract
Molecular dynamics (MD) simulations provide a powerful means of exploring the dynamic behavior of biomolecular systems at the atomic level. However, analyzing the vast data sets generated by MD simulations poses significant challenges. This article discusses the energy landscape visualization method (ELViM), a multidimensional reduction technique inspired by the energy landscape theory. ELViM transcends one-dimensional representations, offering a comprehensive analysis of the effective conformational phase space without the need for predefined reaction coordinates. We apply the ELViM to study the folding landscape of the antimicrobial peptide Polybia-MP1, showcasing its versatility in capturing complex biomolecular dynamics. Using dissimilarity matrices and a force-scheme approach, the ELViM provides intuitive visualizations, revealing structural correlations and local conformational signatures. The method is demonstrated to be adaptable, robust, and applicable to various biomolecular systems.
Collapse
Affiliation(s)
- Rafael Giordano Viegas
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, São Paulo 15.808-305, Brazil
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Ingrid B S Martins
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Murilo Nogueira Sanches
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | | | - Juliana B de Camargo
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Fernando V Paulovich
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, Eindhoven 5600 MB, The Netherlands
| | - Vitor B P Leite
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
5
|
da Silva FB, Simien JM, Viegas RG, Haglund E, Leite VBP. Exploring the folding landscape of leptin: Insights into threading pathways. J Struct Biol 2024; 216:108054. [PMID: 38065428 DOI: 10.1016/j.jsb.2023.108054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Revised: 11/27/2023] [Accepted: 12/04/2023] [Indexed: 12/17/2023]
Abstract
The discovery of new protein topologies with entanglements and loop-crossings have shown the impact of local amino acid arrangement and global three-dimensional structures. This phenomenon plays a crucial role in understanding how protein structure relates to folding and function, affecting the global stability, and biological activity. Protein entanglements encompassing knots and non-trivial topologies add complexity to their folding free energy landscapes. However, the initial native contacts driving the threading event for entangled proteins remains elusive. The Pierced Lasso Topology (PLT) represents an entangled topology where a covalent linker creates a loop in which the polypeptide backbone is threaded through. Compared to true knotted topologies, PLTs are simpler topologies where the covalent-loop persists in all conformations. In this work, the PLT protein leptin, is used to visualize and differentiate the preference for slipknotting over plugging transition pathways along the folding route. We utilize the Energy Landscape Visualization Method (ELViM), a multidimensional projection technique, to visualize and distinguish early threaded conformations that cannot be observed in an in vitro experiment. Critical contacts for the leptin threading mechanisms were identified where the competing pathways are determined by the formation of a hairpin loop in the unfolded basin. Thus, prohibiting the dominant slipknotting pathway. Furthermore, ELViM offers insights into distinct folding pathways associated with slipknotting and plugging providing a novel tool for de novo design and in vitro experiments with residue specific information of threading events in silico.
Collapse
Affiliation(s)
- Fernando Bruno da Silva
- Centre of New Technologies, University of Warsaw, Banacha 2c, Warsaw, Poland; Institute of Biosciences, Humanities and Exact Sciences (IBILCE), São Paulo State University (UNESP), São José do Rio Preto, SP, Brazil
| | - Jennifer M Simien
- Department of Chemistry, University of Hawaii at Manoa, Honolulu, Hawaii 96822, United States
| | - Rafael G Viegas
- Institute of Biosciences, Humanities and Exact Sciences (IBILCE), São Paulo State University (UNESP), São José do Rio Preto, SP, Brazil; Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, SP 15.808-305, Brazil
| | - Ellinor Haglund
- Department of Chemistry, University of Hawaii at Manoa, Honolulu, Hawaii 96822, United States.
| | - Vitor Barbanti Pereira Leite
- Institute of Biosciences, Humanities and Exact Sciences (IBILCE), São Paulo State University (UNESP), São José do Rio Preto, SP, Brazil.
| |
Collapse
|
6
|
Martins IBS, Viegas RG, Sanches MN, de Araujo AS, Leite VBP. Probing Mastoparan-like Antimicrobial Peptides Interaction with Model Membrane Through Energy Landscape Analysis. J Phys Chem B 2024; 128:163-171. [PMID: 38159056 DOI: 10.1021/acs.jpcb.3c05852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
Antimicrobial Peptides (AMPs) have emerged as promising alternatives to conventional antibiotics due to their capacity to disrupt the lipid packing of bacterial cell membranes. This mechanism of action may prevent the development of resistance by bacteria. Understanding their role in lipid packing disruption and their structural properties upon interaction with bacterial membranes is highly desirable. In this study, we employed Molecular Dynamics simulations and the Energy Landscape Visualization Method (ELViM) to characterize and compare the conformational ensembles of mastoparan-like Polybia-MP1 and its analogous H-MP1, in which histidines replace lysine residues. Two situations were analyzed: (i) the peptides in their free state in an aqueous solution containing water and ions and (ii) the peptides spontaneously adsorbing onto an anionic lipid bilayer, used as a bacteria membrane mimetic. ELViM was used to project a single effective conformational phase space for both peptides, providing a comparative analysis. This projection enabled us to map the conformational ensembles of each peptide in an aqueous solution and assess the structural effects of substituting lysines with histidines in H-MP1. Furthermore, a single conformational phase space analysis was employed to describe structural changes during the adsorption process using the same framework. We show that ELViM provides a comprehensive analysis, able to identify discrepancies in the conformational ensembles of these peptides that may affect their affinity to the membrane and adsorption kinetics.
Collapse
Affiliation(s)
- Ingrid B S Martins
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
- Biophysics Institute Carlos Chagas Filho, Federal University of Rio de Janeiro, Rio de Janeiro, RJ 21941-902, Brazil
| | - Rafael G Viegas
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, SP 15.808-305, Brazil
| | - Murilo N Sanches
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Alexandre S de Araujo
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| |
Collapse
|
7
|
Dias RVR, Pedro RP, Sanches MN, Moreira GC, Leite VBP, Caruso IP, de Melo FA, de Oliveira LC. Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding. J Chem Inf Model 2023; 63:6344-6353. [PMID: 37824286 DOI: 10.1021/acs.jcim.3c00955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.
Collapse
Affiliation(s)
- Raphael V R Dias
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
- Multiuser Center for Biomolecular Innovation (CMIB), São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Renan P Pedro
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
- Multiuser Center for Biomolecular Innovation (CMIB), São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Murilo N Sanches
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
| | - Giovana C Moreira
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
- Multiuser Center for Biomolecular Innovation (CMIB), São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Vitor B P Leite
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
| | - Icaro P Caruso
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
- Multiuser Center for Biomolecular Innovation (CMIB), São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Fernando A de Melo
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
- Multiuser Center for Biomolecular Innovation (CMIB), São Paulo State University (UNESP), São José do Rio Preto, SP 15054-000, Brazil
| | - Leandro C de Oliveira
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities, and Exact Sciences, São José do Rio Preto, SP 15054-000, Brazil
| |
Collapse
|
8
|
Viegas RG, Sanches MN, Chen AA, Paulovich FV, Garcia AE, Leite VBP. Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop. J Chem Inf Model 2023; 63:5641-5649. [PMID: 37606640 DOI: 10.1021/acs.jcim.3c00426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.
Collapse
Affiliation(s)
- Rafael G Viegas
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, São Paulo 15.808-305, Brazil
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Murilo N Sanches
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Alan A Chen
- Department of Chemistry and the RNA Institute, University at Albany, Albany, New York 12222, United States
| | - Fernando V Paulovich
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, Eindhoven 5600 MB, the Netherlands
| | - Angel E Garcia
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Physics Division, National Science Foundation, 2415 Eisenhower Ave, Alexandria, Virginia 22314, United States
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
9
|
da Silva FB, Martins de Oliveira V, de Oliveira Junior AB, Contessoto VDG, Leite VBP. Probing the Energy Landscape of Spectrin R15 and R16 and the Effects of Non-native Interactions. J Phys Chem B 2023; 127:1291-1300. [PMID: 36723393 DOI: 10.1021/acs.jpcb.2c06178] [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/02/2023]
Abstract
Understanding the details of a protein folding mechanism can be a challenging and complex task. One system with an interesting folding behavior is the α-spectrin domain, where the R15 folds three-orders of magnitude faster than its homologues R16 and R17, despite having similar structures. The molecular origins that explain these folding rate differences remain unclear, but our previous work revealed that a combined effect produced by non-native interactions could be a reasonable cause for these differences. In this study, we explore further the folding process by identifying the molecular paths, metastable states, and the collective motions that lead these unfolded proteins to their native state conformation. Our results uncovered the differences between the folding pathways for the wild-type R15 and R16 and an R16 mutant. The metastable ensembles that speed down the folding were identified using an energy landscape visualization method (ELViM). These ensembles correspond to similar experimentally reported configurations. Our observations indicate that the non-native interactions are also associated with secondary structure misdocking. This computational methodology can be used as a fast, straightforward protocol for shedding light on systems with unclear folding or conformational traps.
Collapse
Affiliation(s)
- Fernando Bruno da Silva
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo15054-000, Brazil
| | - Vinícius Martins de Oliveira
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| | | | | | - Vitor B P Leite
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo15054-000, Brazil
| |
Collapse
|
10
|
Mouro PR, Sanches MN, Leite VBP, Chahine J. Exploring the Folding Mechanism of Dimeric Superoxide Dismutase. J Phys Chem B 2023; 127:1338-1349. [PMID: 36716437 DOI: 10.1021/acs.jpcb.2c08877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The Cu/Zn Human Superoxide Dismutase (SOD1) is a dimeric metalloenzyme whose genetic mutations are directly related to amyotrophic lateral sclerosis (ALS), so understanding its folding mechanism is of fundamental importance. Currently, the SOD1 dimer formation is studied via molecular dynamics simulations using a simplified structure-based model and an all-atom model. Results from the simplified model reveal a mechanism dependent on distances between monomers, which are limited by constraints to mimic concentration dependence. The stability of intermediates (during the int state) is significantly affected by this distance, as well as by the presence of two folded monomers prior to dimer formation. The kinetics of interface formation are also highly dependent on the separation distance. The folding temperature of the dimer is about 4.2% higher than that of the monomer, a value not too different from experimental data. All-atom simulations on the apo dimer give binding free energy between monomers similar to experimental values. An intermediate state is evident for the apo form at a separation distance between monomers slightly larger than the native distance which has little formed interface between monomers. We have shown that this intermediate is stabilized by non-native intra- and intercontacts.
Collapse
Affiliation(s)
- Paulo R Mouro
- São Paulo State University (UNESP), IBILCE, São José do Rio Preto15054-000, Brazil
| | - Murilo N Sanches
- São Paulo State University (UNESP), IBILCE, São José do Rio Preto15054-000, Brazil
| | - Vitor B P Leite
- São Paulo State University (UNESP), IBILCE, São José do Rio Preto15054-000, Brazil
| | - Jorge Chahine
- São Paulo State University (UNESP), IBILCE, São José do Rio Preto15054-000, Brazil
| |
Collapse
|
11
|
Tsutsumi T, Ono Y, Taketsugu T. Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene. J Chem Theory Comput 2022; 18:7483-7495. [PMID: 36351076 DOI: 10.1021/acs.jctc.2c00560] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.
Collapse
Affiliation(s)
- Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo060-0810, Japan.,L-Station, Creative Research Institution (CRI), Hokkaido University, Sapporo060-0812, Japan
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo001-0021, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo001-0021, Japan
| |
Collapse
|
12
|
Tsutsumi T, Ono Y, Taketsugu T. Reaction Space Projector (ReSPer) for Visualizing Dynamic Reaction Routes Based on Reduced-Dimension Space. Top Curr Chem (Cham) 2022; 380:19. [PMID: 35266073 DOI: 10.1007/s41061-022-00377-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 02/21/2022] [Indexed: 11/26/2022]
Abstract
To analyze chemical reaction dynamics based on a reaction path network, we have developed the "Reaction Space Projector" (ReSPer) method with the aid of the dimensionality reduction method. This program has two functions: the construction of a reduced-dimensionality reaction space from a molecular structure dataset, and the projection of dynamic trajectories into the low-dimensional reaction space. In this paper, we apply ReSPer to isomerization and bifurcation reactions of the Au5 cluster and succeed in analyzing dynamic reaction routes involved in multiple elementary reaction processes, constructing complicated networks (called "closed islands") of nuclear permutation-inversion (NPI) isomerization reactions, and elucidating dynamic behaviors in bifurcation reactions with reference to bundles of trajectories. Interestingly, in the second application, we find a correspondence between the contribution ratios in the ability to visualize and the symmetry of the morphology of closed islands. In addition, the third application suggests the existence of boundaries that determine the selectivity in bifurcation reactions, which was discussed in the phase space. The ReSPer program is a versatile and robust tool to clarify dynamic reaction mechanisms based on the reduced-dimensionality reaction space without prior knowledge of target reactions.
Collapse
Affiliation(s)
- Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, 001-0021, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan.
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, 001-0021, Japan.
| |
Collapse
|
13
|
Kulkarni P, Leite VBP, Roy S, Bhattacharyya S, Mohanty A, Achuthan S, Singh D, Appadurai R, Rangarajan G, Weninger K, Orban J, Srivastava A, Jolly MK, Onuchic JN, Uversky VN, Salgia R. Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma. BIOPHYSICS REVIEWS 2022; 3:011306. [PMID: 38505224 PMCID: PMC10903413 DOI: 10.1063/5.0080512] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Accepted: 02/17/2022] [Indexed: 03/21/2024]
Abstract
Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.
Collapse
Affiliation(s)
- Prakash Kulkarni
- Department of Medical Oncology and Therapeutics Research, City of Hope National Medical Center, Duarte, California 91010, USA
| | - Vitor B. P. Leite
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal 741246, India
| | - Supriyo Bhattacharyya
- Translational Bioinformatics, Center for Informatics, Department of Computational and Quantitative Medicine, City of Hope National Medical Center, Duarte, California 91010, USA
| | - Atish Mohanty
- Department of Medical Oncology and Therapeutics Research, City of Hope National Medical Center, Duarte, California 91010, USA
| | - Srisairam Achuthan
- Center for Informatics, Division of Research Informatics, City of Hope National Medical Center, Duarte, California 91010, USA
| | - Divyoj Singh
- Center for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India
| | - Rajeswari Appadurai
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka, India
| | - Govindan Rangarajan
- Department of Mathematics, Indian Institute of Science, Bangalore 560012, India
| | - Keith Weninger
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695, USA
| | | | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka, India
| | - Mohit Kumar Jolly
- Center for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India
| | - Jose N. Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005-1892, USA
| | | | - Ravi Salgia
- Department of Medical Oncology and Therapeutics Research, City of Hope National Medical Center, Duarte, California 91010, USA
| |
Collapse
|
14
|
Abstract
Abstract
Ankyrin (ANK) repeat proteins are coded by tandem occurrences of patterns with around 33 amino acids. They often mediate protein–protein interactions in a diversity of biological systems. These proteins have an elongated non-globular shape and often display complex folding mechanisms. This work investigates the energy landscape of representative proteins of this class made up of 3, 4 and 6 ANK repeats using the energy-landscape visualisation method (ELViM). By combining biased and unbiased coarse-grained molecular dynamics AWSEM simulations that sample conformations along the folding trajectories with the ELViM structure-based phase space, one finds a three-dimensional representation of the globally funnelled energy surface. In this representation, it is possible to delineate distinct folding pathways. We show that ELViMs can project, in a natural way, the intricacies of the highly dimensional energy landscapes encoded by the highly symmetric ankyrin repeat proteins into useful low-dimensional representations. These projections can discriminate between multiplicities of specific parallel folding mechanisms that otherwise can be hidden in oversimplified depictions.
Collapse
|
15
|
Sanches MN, Knapp K, Oliveira AB, Wolynes PG, Onuchic JN, Leite VBP. Examining the Ensembles of Amyloid-β Monomer Variants and Their Propensities to Form Fibers Using an Energy Landscape Visualization Method. J Phys Chem B 2021; 126:93-99. [PMID: 34968059 DOI: 10.1021/acs.jpcb.1c08525] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
The amyloid-β (Aβ) monomer, an intrinsically disordered peptide, is produced by the cleavage of the amyloid precursor protein, leading to Aβ-40 and Aβ-42 as major products. These two isoforms generate pathological aggregates, whose accumulation correlates with Alzheimer's disease (AD). Experiments have shown that even though the natural abundance of Aβ-42 is smaller than that for Aβ-40, the Aβ-42 is more aggregation-prone compared to Aβ-40. Moreover, several single-point mutations are associated with early onset forms of AD. This work analyzes coarse-grained associative-memory, water-mediated, structure and energy model (AWSEM) simulations of normal Aβ-40 and Aβ-42 monomers, along with six single-point mutations associated with early onset disease. We analyzed the simulations using the energy landscape visualization method (ELViM), a reaction-coordinate-free approach suited to explore the frustrated energy landscapes of intrinsically disordered proteins. ELViM is shown to distinguish the monomer ensembles of variants that rapidly form fibers from those that do not form fibers as readily. It also delineates the amino acid contacts characterizing each ensemble. The results shed light on the potential of ELViM to probe intrinsically disordered proteins.
Collapse
Affiliation(s)
- Murilo N Sanches
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Kaitlin Knapp
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Antonio B Oliveira
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States.,Departments of Physics and Astronomy, Chemistry, and Biosciences, Rice University, Houston, Texas 77005, United States
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
16
|
Oliveira Junior AB, Lin X, Kulkarni P, Onuchic JN, Roy S, Leite VBP. Exploring Energy Landscapes of Intrinsically Disordered Proteins: Insights into Functional Mechanisms. J Chem Theory Comput 2021; 17:3178-3187. [PMID: 33871257 DOI: 10.1021/acs.jctc.1c00027] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Intrinsically disordered proteins (IDPs) lack a rigid three-dimensional structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed energy landscape visualization method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, prostate-associated gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely, HIPK1 and CLK2, phosphorylate PAGE4, generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyperphosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conformational diversity among different phosphoforms, we have analyzed their atomistic trajectories simulated using AWSEM forcefield, and the energy landscapes were elucidated using ELViM. This method allows us to identify and compare the population distributions of different conformational ensembles of PAGE4 phosphoforms using the same effective phase space. The results reveal a predominant conformational ensemble with an extended C-terminal segment of WT PAGE4, which exposes a functional residue Thr51, implying its potential of undertaking a fly-casting mechanism while binding to its cognate partner. In contrast, for HIPK1-PAGE4, a compact conformational ensemble enhances its population sequestering phosphorylated-Thr51. This clearly explains the experimentally observed weaker affinity of HIPK1-PAGE4 for c-Jun. ELViM appears as a powerful tool, especially to analyze the highly frustrated energy landscape representation of IDPs where appropriate reaction coordinates are hard to apprehend.
Collapse
Affiliation(s)
- Antonio B Oliveira Junior
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005-1892, United States
| | - Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Prakash Kulkarni
- Department of Medical Oncology and Therapeutics Research, City of Hope National Medical Center, Duarte, California 91010, United States
| | - José N Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005-1892, United States
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal 741246, India
| | - Vitor B P Leite
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
17
|
de Oliveira VM, Caetano DLZ, da Silva FB, Mouro PR, de Oliveira AB, de Carvalho SJ, Leite VBP. pH and Charged Mutations Modulate Cold Shock Protein Folding and Stability: A Constant pH Monte Carlo Study. J Chem Theory Comput 2020; 16:765-772. [PMID: 31756296 DOI: 10.1021/acs.jctc.9b00894] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The folding and stability of proteins is a fundamental problem in several research fields. In the present paper, we have used different computational approaches to study the effects caused by changes in pH and for charged mutations in cold shock proteins from Bacillus subtilis (Bs-CspB). First, we have investigated the contribution of each ionizable residue for these proteins to their thermal stability using the TKSA-MC, a Web server for rational mutation via optimizing the protein charge interactions. Based on these results, we have proposed a new mutation in an already optimized Bs-CspB variant. We have evaluated the effects of this new mutation in the folding energy landscape using structure-based models in Monte Carlo simulation at constant pH, SBM-CpHMC. Our results using this approach have indicated that the charge rearrangements already in the unfolded state are critical to the thermal stability of Bs-CspB. Furthermore, the conjunction of these simplified methods was able not only to predict stabilizing mutations in different pHs but also to provide essential information about their effects in each stage of protein folding.
Collapse
Affiliation(s)
- Vinícius M de Oliveira
- Brazilian Biosciences National Laboratory, National Center for Research in Energy and Materials, LNBio/CNPEM , Campinas , São Paulo , 13083-970 , Brazil
| | - Daniel L Z Caetano
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil
| | - Fernando B da Silva
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil
| | - Paulo R Mouro
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil
| | - Antonio B de Oliveira
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil
| | - Sidney J de Carvalho
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil
| | - Vitor B P Leite
- Department of Physics , São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences , São José do Rio Preto , São Paulo , 15054-000 , Brazil.,Center for Theoretical Biological Physics , Rice University , Houston , Texas 77005 , United States
| |
Collapse
|