1
|
Erkip A, Erman B. Dynamically driven correlations in elastic net models reveal sequence of events and causality in proteins. Proteins 2024; 92:1113-1126. [PMID: 38687146 DOI: 10.1002/prot.26697] [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: 01/18/2024] [Revised: 04/07/2024] [Accepted: 04/16/2024] [Indexed: 05/02/2024]
Abstract
An explicit analytic solution is given for the Langevin equation applied to the Gaussian Network Model of a protein subjected to both a random and a deterministic periodic force. Synchronous and asynchronous components of time correlation functions are derived and an expression for phase differences in the time correlations of residue pairs is obtained. The synchronous component enables the determination of dynamic communities within the protein structure. The asynchronous component reveals causality, where the time correlation function between residues i and j differs depending on whether i is observed before j or vice versa, resulting in directional information flow. Driver and driven residues in the allosteric process of cyclophilin A and human NAD-dependent isocitrate dehydrogenase are determined by a perturbation-scanning technique. Factors affecting phase differences between fluctuations of residues, such as network topology, connectivity, and residue centrality, are identified. Within the constraints of the isotropic Gaussian Network Model, our results show that asynchronicity increases with viscosity and distance between residues, decreases with increasing connectivity, and decreases with increasing levels of eigenvector centrality.
Collapse
Affiliation(s)
- Albert Erkip
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Burak Erman
- Department of Chemical and Biological Engineering, Koc University, Rumelifeneri Yolu, Istanbul, Turkey
| |
Collapse
|
2
|
Böttcher L, Porter MA. Complex networks with complex weights. Phys Rev E 2024; 109:024314. [PMID: 38491610 DOI: 10.1103/physreve.109.024314] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 12/20/2023] [Indexed: 03/18/2024]
Abstract
In many studies, it is common to use binary (i.e., unweighted) edges to examine networks of entities that are either adjacent or not adjacent. Researchers have generalized such binary networks to incorporate edge weights, which allow one to encode node-node interactions with heterogeneous intensities or frequencies (e.g., in transportation networks, supply chains, and social networks). Most such studies have considered real-valued weights, despite the fact that networks with complex weights arise in fields as diverse as quantum information, quantum chemistry, electrodynamics, rheology, and machine learning. Many of the standard network-science approaches in the study of classical systems rely on the real-valued nature of edge weights, so it is necessary to generalize them if one seeks to use them to analyze networks with complex edge weights. In this paper, we examine how standard network-analysis methods fail to capture structural features of networks with complex edge weights. We then generalize several network measures to the complex domain and show that random-walk centralities provide a useful approach to examine node importances in networks with complex weights.
Collapse
Affiliation(s)
- Lucas Böttcher
- Department of Computational Science and Philosophy, Frankfurt School of Finance and Management, 60322 Frankfurt am Main, Germany
- Department of Medicine, University of Florida, Gainesville, Florida, 32610, USA
| | - Mason A Porter
- Department of Mathematics, University of California, Los Angeles, California 90095, USA
- Department of Sociology, University of California, Los Angeles, California 90095, USA
- Santa Fe Institute, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
3
|
Klem H, Alegre-Requena JV, Paton RS. Catalytic Effects of Active Site Conformational Change in the Allosteric Activation of Imidazole Glycerol Phosphate Synthase. ACS Catal 2023; 13:16249-16257. [PMID: 38125975 PMCID: PMC10729027 DOI: 10.1021/acscatal.3c04176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Revised: 11/10/2023] [Accepted: 11/13/2023] [Indexed: 12/23/2023]
Abstract
Imidazole glycerol phosphate synthase (IGPS) is a class-I glutamine amidotransferase (GAT) that hydrolyzes glutamine. Ammonia is produced and transferred to a second active site, where it reacts with N1-(5'-phosphoribosyl)-formimino-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) to form precursors to purine and histidine biosynthesis. Binding of PrFAR over 25 Å away from the active site increases glutaminase efficiency by ∼4500-fold, primarily altering the glutamine turnover number. IGPS has been the focus of many studies on allosteric communication; however, atomic details for how the glutamine hydrolysis rate increases in the presence of PrFAR are lacking. We present a density functional theory study on 237-atom active site cluster models of IGPS based on crystallized structures representing the inactive and allosterically active conformations and investigate the multistep reaction leading to thioester formation and ammonia production. The proposed mechanism is supported by similar, well-studied enzyme mechanisms, and the corresponding energy profile is consistent with steady-state kinetic studies of PrFAR + IGPS. Additional active site models are constructed to examine the relationship between active site structural change and transition-state stabilization via energy decomposition schemes. The results reveal that the inactive IGPS conformation does not provide an adequately formed oxyanion hole structure and that repositioning of the oxyanion strand relative to the substrate is vital for a catalysis-competent oxyanion hole, with or without the hVal51 dihedral flip. These findings are valuable for future endeavors in modeling the IGPS allosteric mechanism by providing insight into the atomistic changes required for rate enhancement that can inform suitable reaction coordinates for subsequent investigations.
Collapse
Affiliation(s)
- Heidi Klem
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Juan V Alegre-Requena
- Dpto.de Química Inorgánica, Instituto de Síntesis Química y Catálisis Homogénea (ISQCH), CSIC, Universidad de Zaragoza, Zaragoza 50009, Spain
| | - Robert S Paton
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| |
Collapse
|
4
|
La Sala G, Pfleger C, Käck H, Wissler L, Nevin P, Böhm K, Janet JP, Schimpl M, Stubbs CJ, De Vivo M, Tyrchan C, Hogner A, Gohlke H, Frolov AI. Combining structural and coevolution information to unveil allosteric sites. Chem Sci 2023; 14:7057-7067. [PMID: 37389247 PMCID: PMC10306073 DOI: 10.1039/d2sc06272k] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Accepted: 06/02/2023] [Indexed: 07/01/2023] Open
Abstract
Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.
Collapse
Affiliation(s)
- Giuseppina La Sala
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
| | - Helena Käck
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Lisa Wissler
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Philip Nevin
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Kerstin Böhm
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Jon Paul Janet
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Marianne Schimpl
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Christopher J Stubbs
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Design, Istituto Italiano di Tecnologia Via Morego 30 16163 Genoa Italy
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory & Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Anders Hogner
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Institute of Bio- and Geosciences (IBG-4: Bioinformatics) Forschungszentrum Jülich GmbH 52425 Jülich Germany
| | - Andrey I Frolov
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| |
Collapse
|
5
|
Rosignoli S, Paiardini A. Boosting the Full Potential of PyMOL with Structural Biology Plugins. Biomolecules 2022; 12:biom12121764. [PMID: 36551192 PMCID: PMC9775141 DOI: 10.3390/biom12121764] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Revised: 11/23/2022] [Accepted: 11/24/2022] [Indexed: 11/29/2022] Open
Abstract
Over the past few decades, the number of available structural bioinformatics pipelines, libraries, plugins, web resources and software has increased exponentially and become accessible to the broad realm of life scientists. This expansion has shaped the field as a tangled network of methods, algorithms and user interfaces. In recent years PyMOL, widely used software for biomolecules visualization and analysis, has started to play a key role in providing an open platform for the successful implementation of expert knowledge into an easy-to-use molecular graphics tool. This review outlines the plugins and features that make PyMOL an eligible environment for supporting structural bioinformatics analyses.
Collapse
|
6
|
Helmer N, Wolf S, Stock G. Energy Transport and Its Function in Heptahelical Transmembrane Proteins. J Phys Chem B 2022; 126:8735-8746. [PMID: 36261792 DOI: 10.1021/acs.jpcb.2c05892] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β2 adrenergic receptor (β2AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β2AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β2AR, it seems not particularly suited to explain allosteric phenomena.
Collapse
Affiliation(s)
- Nadja Helmer
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| |
Collapse
|
7
|
Calvó-Tusell C, Maria-Solano MA, Osuna S, Feixas F. Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase. J Am Chem Soc 2022; 144:7146-7159. [PMID: 35412310 PMCID: PMC9052757 DOI: 10.1021/jacs.1c12629] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Deciphering the molecular
mechanisms of enzymatic allosteric regulation
requires the structural characterization of functional states and
also their time evolution toward the formation of the allosterically
activated ternary complex. The transient nature and usually slow millisecond
time scale interconversion between these functional states hamper
their experimental and computational characterization. Here, we combine
extensive molecular dynamics simulations, enhanced sampling techniques,
and dynamical networks to describe the allosteric activation of imidazole
glycerol phosphate synthase (IGPS) from the substrate-free form to
the active ternary complex. IGPS is a heterodimeric bienzyme complex
whose HisH subunit is responsible for hydrolyzing glutamine and delivering
ammonia for the cyclase activity in HisF. Despite significant advances
in understanding the underlying allosteric mechanism, essential molecular
details of the long-range millisecond allosteric activation of IGPS
remain hidden. Without using a priori information
of the active state, our simulations uncover how IGPS, with the allosteric
effector bound in HisF, spontaneously captures glutamine in a catalytically
inactive HisH conformation, subsequently attains a closed HisF:HisH
interface, and finally forms the oxyanion hole in HisH for efficient
glutamine hydrolysis. We show that the combined effector and substrate
binding dramatically decreases the conformational barrier associated
with oxyanion hole formation, in line with the experimentally observed
4500-fold activity increase in glutamine hydrolysis. The allosteric
activation is controlled by correlated time-evolving dynamic networks
connecting the effector and substrate binding sites. This computational
strategy tailored to describe millisecond events can be used to rationalize
the effect of mutations on the allosteric regulation and guide IGPS
engineering efforts.
Collapse
Affiliation(s)
- Carla Calvó-Tusell
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| | - Miguel A Maria-Solano
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain.,Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Science, Ewha Womans University, 03760 Seoul, Republic of Korea
| | - Sílvia Osuna
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain.,Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Catalonia, Spain
| | - Ferran Feixas
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| |
Collapse
|
8
|
Yao XQ, Hamelberg D. From Distinct to Differential Conformational Dynamics to Map Allosteric Communication Pathways in Proteins. J Phys Chem B 2022; 126:2612-2620. [PMID: 35319195 DOI: 10.1021/acs.jpcb.2c00199] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Initiation of biological processes involving protein-ligand binding, transient protein-protein interactions, or amino acid modifications alters the conformational dynamics of proteins. Accompanying these biological processes are ensuing coupled atomic level conformational changes within the proteins. These conformational changes collectively connect multiple amino acid residues at distal allosteric, binding, and/or active sites. Local changes due to, for example, binding of a regulatory ligand at an allosteric site initiate the allosteric regulation. The allosteric signal propagates throughout the protein structure, causing changes at distal sites, activating, deactivating, or modifying the function of the protein. Hence, dynamical responses within protein structures to stimuli contain critical information on protein function. In this Perspective, we examine the description of allosteric regulation from protein dynamical responses and associated alternative and emerging computational approaches to map allosteric communication pathways between distal sites in proteins at the atomic level.
Collapse
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
9
|
Yao XQ, Hamelberg D. Residue–Residue Contact Changes during Functional Processes Define Allosteric Communication Pathways. J Chem Theory Comput 2022; 18:1173-1187. [DOI: 10.1021/acs.jctc.1c00669] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
10
|
Sladek V, Yamamoto Y, Harada R, Shoji M, Shigeta Y, Sladek V. pyProGA-A PyMOL plugin for protein residue network analysis. PLoS One 2021; 16:e0255167. [PMID: 34329304 PMCID: PMC8323899 DOI: 10.1371/journal.pone.0255167] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2021] [Accepted: 07/11/2021] [Indexed: 11/18/2022] Open
Abstract
The field of protein residue network (PRN) research has brought several useful methods and techniques for structural analysis of proteins and protein complexes. Many of these are ripe and ready to be used by the proteomics community outside of the PRN specialists. In this paper we present software which collects an ensemble of (network) methods tailored towards the analysis of protein-protein interactions (PPI) and/or interactions of proteins with ligands of other type, e.g. nucleic acids, oligosaccharides etc. In parallel, we propose the use of the network differential analysis as a method to identify residues mediating key interactions between proteins. We use a model system, to show that in combination with other, already published methods, also included in pyProGA, it can be used to make such predictions. Such extended repertoire of methods allows to cross-check predictions with other methods as well, as we show here. In addition, the possibility to construct PRN models from various kinds of input is so far a unique asset of our code. One can use structural data as defined in PDB files and/or from data on residue pair interaction energies, either from force-field parameters or fragment molecular orbital (FMO) calculations. pyProGA is a free open-source software available from https://gitlab.com/Vlado_S/pyproga.
Collapse
Affiliation(s)
- Vladimir Sladek
- Institute of Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Yuta Yamamoto
- Department of Chemistry, Rikkyo University, Nishi-Ikebukuro, Tokyo, Japan
| | - Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki, Japan
| | - Mitsuo Shoji
- Center for Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki, Japan
| | - Vladimir Sladek
- Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia
| |
Collapse
|
11
|
Foster CA, Silversmith RE, Immormino RM, Vass LR, Kennedy EN, Pazy Y, Collins EJ, Bourret RB. Role of Position K+4 in the Phosphorylation and Dephosphorylation Reaction Kinetics of the CheY Response Regulator. Biochemistry 2021; 60:2130-2151. [PMID: 34167303 DOI: 10.1021/acs.biochem.1c00246] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Two-component signaling is a primary method by which microorganisms interact with their environments. A kinase detects stimuli and modulates autophosphorylation activity. The signal propagates by phosphotransfer from the kinase to a response regulator, eliciting a response. Response regulators operate over a range of time scales, corresponding to their related biological processes. Response regulator active site chemistry is highly conserved, but certain variable residues can influence phosphorylation kinetics. An Ala-to-Pro substitution (K+4, residue 113) in the Escherichia coli response regulator CheY triggers a constitutively active phenotype; however, the A113P substitution is too far from the active site to directly affect phosphochemistry. To better understand the activating mechanism(s) of the substitution, we analyzed receiver domain sequences to characterize the evolutionary role of the K+4 position. Although most featured Pro, Leu, Ile, and Val residues, chemotaxis-related proteins exhibited atypical Ala, Gly, Asp, and Glu residues at K+4. Structural and in silico analyses revealed that CheY A113P adopted a partially active configuration. Biochemical data showed that A113P shifted CheY toward a more activated state, enhancing autophosphorylation. By characterizing CheY variants, we determined that this functionality was transmitted through a hydrophobic network bounded by the β5α5 loop and the α1 helix of CheY. This region also interacts with the phosphodonor CheAP1, suggesting that binding generates an activating perturbation similar to the A113P substitution. Atypical residues like Ala at the K+4 position likely serve two purposes. First, restricting autophosphorylation may minimize background noise generated by intracellular phosphodonors such as acetyl phosphate. Second, optimizing interactions with upstream partners may help prime the receiver domain for phosphorylation.
Collapse
Affiliation(s)
- Clay A Foster
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Ruth E Silversmith
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Robert M Immormino
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Luke R Vass
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Emily N Kennedy
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yael Pazy
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Edward J Collins
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Robert B Bourret
- Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
12
|
Sobieraj M, Setny P. Entropy-based distance cutoff for protein internal contact networks. Proteins 2021; 89:1333-1339. [PMID: 34053102 DOI: 10.1002/prot.26154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Revised: 05/17/2021] [Accepted: 05/19/2021] [Indexed: 11/05/2022]
Abstract
Protein structure networks (PSNs) have long been used to provide a coarse yet meaningful representation of protein structure, dynamics, and internal communication pathways. An important question is what criteria should be applied to construct the network so that to include relevant interresidue contacts while avoiding unnecessary connections. To address this issue, we systematically considered varying residue distance cutoff length and the probability threshold for contact formation to construct PSNs based on atomistic molecular dynamics in order to assess the amount of mutual information within the resulting representations. We found that the minimum in mutual information is universally achieved at the cutoff length of 5 Å, irrespective of the applied contact formation probability threshold in all considered, distinct proteins. Assuming that the optimal PSNs should be characterized by the least amount of redundancy, which corresponds to the minimum in mutual information, this finding suggests an objective criterion for cutoff distance and supports the existing preference towards its customary selection around 5 Å length, typically based to date on heuristic criteria.
Collapse
Affiliation(s)
- Marcin Sobieraj
- Center of New Technologies, University of Warsaw, Warsaw, Poland.,Department of Physics, University of Warsaw, Warsaw, Poland
| | - Piotr Setny
- Center of New Technologies, University of Warsaw, Warsaw, Poland
| |
Collapse
|
13
|
Dodd T, Yao XQ, Hamelberg D, Ivanov I. Subsets of adjacent nodes (SOAN): a fast method for computing suboptimal paths in protein dynamic networks. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1893847] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Thomas Dodd
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Ivaylo Ivanov
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| |
Collapse
|
14
|
Botello-Smith WM, Luo YL. Concepts, Practices, and Interactive Tutorial for Allosteric Network Analysis of Molecular Dynamics Simulations. Methods Mol Biol 2021; 2302:311-334. [PMID: 33877635 DOI: 10.1007/978-1-0716-1394-8_17] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Over the past decade, concepts of network theory in combination with dynamical information from conformational ensembles have been widely applied to gain insights in understanding allosteric regulation in biomolecules. In this chapter, we introduce the basic theories and protocols used in protein dynamics network analysis through a series of interactive python Jupyter notebook scripts. While various network analysis methods exist in the literature, here we focus on the two popular methods based on correlated atomic motions and pairwise interaction energies. While the tutorial is based on a small prototypic protein, the workflow and protocol introduced here are optimized to handle large membrane proteins.
Collapse
Affiliation(s)
- Wesley M Botello-Smith
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, CA, USA.
| | - Yun Lyna Luo
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, CA, USA
| |
Collapse
|
15
|
Rocha GV, Bastos LS, Costa MGS. Identification of potential allosteric binding sites in cathepsin K based on intramolecular communication. Proteins 2020; 88:1675-1687. [PMID: 32683717 DOI: 10.1002/prot.25985] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Revised: 06/02/2020] [Accepted: 07/12/2020] [Indexed: 12/18/2022]
Abstract
Network theory methods and molecular dynamics (MD) simulations are accepted tools to study allosteric regulation. Indeed, dynamic networks built upon correlation analysis of MD trajectories provide detailed information about communication paths between distant sites. In this context, we aimed to understand whether the efficiency of intramolecular communication could be used to predict the allosteric potential of a given site. To this end, we performed MD simulations and network theory analyses in cathepsin K (catK), whose allosteric sites are well defined. To obtain a quantitative measure of the efficiency of communication, we designed a new protocol that enables the comparison between properties related to ensembles of communication paths obtained from different sites. Further, we applied our strategy to evaluate the allosteric potential of different catK cavities not yet considered for drug design. Our predictions of the allosteric potential based on intramolecular communication correlate well with previous catK experimental and theoretical data. We also discuss the possibility of applying our approach to other proteins from the same family.
Collapse
Affiliation(s)
- Gisele V Rocha
- Programa de Computação Científica, Vice-Presidência de Educação, Informação e Comunicação, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratoire de Biologie et de Pharmacologie Appliquée, Ecole Normale Supérieure Paris Saclay, Centre National de la Recherche Scientifique, Cachan, France
| | - Leonardo S Bastos
- Programa de Computação Científica, Vice-Presidência de Educação, Informação e Comunicação, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Department of Infectious Diseases Epidemiology, London School of Hygiene and Tropical Medicine, London, UK
| | - Mauricio G S Costa
- Programa de Computação Científica, Vice-Presidência de Educação, Informação e Comunicação, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratoire de Biologie et de Pharmacologie Appliquée, Ecole Normale Supérieure Paris Saclay, Centre National de la Recherche Scientifique, Cachan, France
| |
Collapse
|
16
|
Lake PT, Davidson RB, Klem H, Hocky GM, McCullagh M. Residue-Level Allostery Propagates through the Effective Coarse-Grained Hessian. J Chem Theory Comput 2020; 16:3385-3395. [PMID: 32251581 DOI: 10.1021/acs.jctc.9b01149] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The long-ranged coupling between residues that gives rise to allostery in a protein is built up from short-ranged physical interactions. Computational tools used to predict this coupling and its functional relevance have relied on the application of graph theoretical metrics to residue-level correlations measured from all-atom molecular dynamics simulations. The short-ranged interactions that yield these long-ranged residue-level correlations are quantified by the effective coarse-grained Hessian. Here we compute an effective harmonic coarse-grained Hessian from simulations of a benchmark allosteric protein, IGPS, and demonstrate the improved locality of this graph Laplacian over two other connectivity matrices. Additionally, two centrality metrics are developed that indicate the direct and indirect importance of each residue at producing the covariance between the effector binding pocket and the active site. The residue importance indicated by these two metrics is corroborated by previous mutagenesis experiments and leads to unique functional insights; in contrast to previous computational analyses, our results suggest that fP76-hK181 is the most important contact for conveying direct allosteric paths across the HisF-HisH interface. The connectivity around fD98 is found to be important at affecting allostery through indirect means.
Collapse
Affiliation(s)
- Peter T Lake
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Russell B Davidson
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Heidi Klem
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Glen M Hocky
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Martin McCullagh
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| |
Collapse
|
17
|
Wang J, Samanta R, Custer G, Look C, Matysiak S, Beckett D. Tuning Allostery through Integration of Disorder to Order with a Residue Network. Biochemistry 2020; 59:790-801. [PMID: 31899864 DOI: 10.1021/acs.biochem.9b01006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In allostery, a signal from one site in a protein is transmitted to a second site to alter its function. Due to its ubiquity in biology and the potential for its exploitation in drug and protein design, the molecular basis of allosteric communication continues to be the subject of intense research. Although allosterically coupled sites are frequently characterized by disorder, how communication between disordered segments occurs remains obscure. Allosteric activation of Escherichia coli BirA dimerization occurs via coupled distant disorder-to-order transitions. In this work, combined structural and computational studies reveal an extensive residue network in BirA. Substitution of several network residues yields large perturbations to allostery. Force distribution analysis reveals that disruptions to the disorder-to-order transitions through amino acid substitution are manifested in shifts in the energy experienced by network residues as well as alterations in packing of an α-helix that plays a critical role in allostery. The combined results reveal a highly distributed allosteric mechanism that is robust to sequence change.
Collapse
Affiliation(s)
- Jingheng Wang
- Department of Chemistry & Biochemistry , University of Maryland , College Park , Maryland 20742 , United States
| | - Riya Samanta
- Biophysics Graduate Program , University of Maryland , College Park , Maryland 20742 , United States
| | - Gregory Custer
- Fischell Department of Bioengineering , University of Maryland , College Park , Maryland 20742 , United States
| | - Christopher Look
- Fischell Department of Bioengineering , University of Maryland , College Park , Maryland 20742 , United States
| | - Silvina Matysiak
- Fischell Department of Bioengineering , University of Maryland , College Park , Maryland 20742 , United States
| | - Dorothy Beckett
- Department of Chemistry & Biochemistry , University of Maryland , College Park , Maryland 20742 , United States
| |
Collapse
|
18
|
Valiño Borau L, Gulzar A, Stock G. Master equation model to predict energy transport pathways in proteins. J Chem Phys 2020; 152:045103. [DOI: 10.1063/1.5140070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Luis Valiño Borau
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| |
Collapse
|
19
|
Leitner DM, Pandey HD, Reid KM. Energy Transport across Interfaces in Biomolecular Systems. J Phys Chem B 2019; 123:9507-9524. [DOI: 10.1021/acs.jpcb.9b07086] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- David M. Leitner
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Hari Datt Pandey
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Korey M. Reid
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| |
Collapse
|
20
|
Gulzar A, Valiño Borau L, Buchenberg S, Wolf S, Stock G. Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study. J Chem Theory Comput 2019; 15:5750-5757. [PMID: 31433644 DOI: 10.1021/acs.jctc.9b00598] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.
Collapse
Affiliation(s)
- Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Luis Valiño Borau
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Sebastian Buchenberg
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| |
Collapse
|
21
|
Serçinoglu O, Ozbek P. gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations. Nucleic Acids Res 2019; 46:W554-W562. [PMID: 29800260 PMCID: PMC6030995 DOI: 10.1093/nar/gky381] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 05/22/2018] [Indexed: 11/12/2022] Open
Abstract
Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.
Collapse
Affiliation(s)
- Onur Serçinoglu
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| |
Collapse
|
22
|
Yao XQ, Momin M, Hamelberg D. Establishing a Framework of Using Residue–Residue Interactions in Protein Difference Network Analysis. J Chem Inf Model 2019; 59:3222-3228. [DOI: 10.1021/acs.jcim.9b00320] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Mohamed Momin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
23
|
Xia Q, Ding Y. Thermostability of Lipase A and Dynamic Communication Based on Residue Interaction Network. Protein Pept Lett 2019; 26:702-716. [PMID: 31215367 DOI: 10.2174/0929866526666190617091812] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Revised: 04/10/2019] [Accepted: 04/25/2019] [Indexed: 11/22/2022]
Abstract
OBJECTIVE Dynamic communication caused by mutation affects protein stability. The main objective of this study is to explore how mutations affect communication and to provide further insight into the relationship between heat resistance and signal propagation of Bacillus subtilis lipase (Lip A). METHODS The relationship between dynamic communication and Lip A thermostability is studied by long-time MD simulation and residue interaction network. The Dijkstra algorithm is used to get the shortest path of each residue pair. Subsequently, time-series frequent paths and spatio-temporal frequent paths are mined through an Apriori-like algorithm. RESULTS Time-series frequent paths show that the communication between residue pairs, both in wild-type lipase (WTL) and mutant 6B, becomes chaotic with an increase in temperature; however, more residues in 6B can maintain stable communication at high temperature, which may be associated with the structural rigidity. Furthermore, spatio-temporal frequent paths reflect the interactions among secondary structures. For WTL at 300K, β7, αC, αB, the longest loop, αA and αF contact frequently. The 310-helix between β3 and αA is penetrated by spatio-temporal frequent paths. At 400K, only αC can be frequently transmitted. For 6B, when at 300K, αA and αF are in more tight contact by spatio-temporal frequent paths though I157M and N166Y. Moreover, the rigidity of the active site His156 and the C-terminal of Lip A are increased, as reflected by the spatio-temporal frequent paths. At 400K, αA and αF, 310-helix between β3 and αA, the longest loop, and the loop where the active site Asp133 is located can still maintain stable communication. CONCLUSION From the perspective of residue dynamic communication, it is obviously found that mutations cause changes in interactions between secondary structures and enhance the rigidity of the structure, contributing to the thermal stability and functional activity of 6B.
Collapse
Affiliation(s)
- Qian Xia
- Laboratory of Media Design and Software Technology, Jiangnan University, Wuxi, Jiangsu, 214122, China
| | - Yanrui Ding
- Laboratory of Media Design and Software Technology, Jiangnan University, Wuxi, Jiangsu, 214122, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Wuxi, Jiangsu, 214122, China
| |
Collapse
|
24
|
Liang Z, Verkhivker GM, Hu G. Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications. Brief Bioinform 2019; 21:815-835. [DOI: 10.1093/bib/bbz029] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Revised: 02/04/2019] [Accepted: 02/21/2019] [Indexed: 12/24/2022] Open
Abstract
Abstract
Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.
Collapse
Affiliation(s)
- Zhongjie Liang
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, USA
| | - Guang Hu
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| |
Collapse
|
25
|
Botello-Smith WM, Luo Y. Robust Determination of Protein Allosteric Signaling Pathways. J Chem Theory Comput 2019; 15:2116-2126. [DOI: 10.1021/acs.jctc.8b01197] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Wesley M. Botello-Smith
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Yun Luo
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| |
Collapse
|
26
|
Sladek V, Tokiwa H, Shimano H, Shigeta Y. Protein Residue Networks from Energetic and Geometric Data: Are They Identical? J Chem Theory Comput 2018; 14:6623-6631. [PMID: 30500196 DOI: 10.1021/acs.jctc.8b00733] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein residue networks (PRN) from energetic and geometric data are probably not identical. PRNs constructed from ab initio pair interaction energies are analyzed for the first time and compared to PRN based on center of mass separation. We use modern, previously unused algorithms such as global and local efficiencies to quantitatively confirm that both types of PRNs do exhibit small-world character. The main novelty finding is that interaction energy-based PRNs preserve small-world character even when clustered. A node hierarchy independent of the cutoff energy used for the edge creation is characteristic for them. Efficiency centrality identifies hubs responsible for such behavior. The interaction energy-based PRNs seem to comply with the scale-free network model with respect to efficiency centrality distribution as opposed to distance based PRNs. Community detection is introduced into protein network research as an extension beyond cluster analysis to study tertiary and quaternary structures.
Collapse
Affiliation(s)
- Vladimir Sladek
- Institute of Chemistry - Centre for Glycomics , Dubravska cesta 9 , 84538 Bratislava , Slovakia.,Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan
| | - Hiroaki Tokiwa
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Chemistry , Rikkyo University , Nishi-Ikebukuro , Toshima, Tokyo 171-8501 , Japan
| | - Hitoshi Shimano
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Internal Medicine, Faculty of Medicine , University of Tsukuba , 1-1-1 Tennodai , Tsukuba, Ibaraki 305-8575 , Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences , University of Tsukuba , Tennodai 1-1-1 , Tsukuba, Ibaraki 305-8577 , Japan
| |
Collapse
|
27
|
Bingöl EN, Serçinoğlu O, Ozbek P. How do mutations and allosteric inhibitors modulate caspase-7 activity? A molecular dynamics study. J Biomol Struct Dyn 2018; 37:3456-3466. [DOI: 10.1080/07391102.2018.1517611] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Elif Naz Bingöl
- Department of Bioengineering, Institute of Pure and Applied Sciences, Marmara University, Istanbul, Turkey
| | - Onur Serçinoğlu
- Department of Bioengineering, Institute of Pure and Applied Sciences, Marmara University, Istanbul, Turkey
- Faculty of Engineering, Department of Bioengineering, Marmara University, Istanbul, Turkey
| | - Pemra Ozbek
- Faculty of Engineering, Department of Bioengineering, Marmara University, Istanbul, Turkey
| |
Collapse
|
28
|
Leitner DM, Yamato T. MAPPING ENERGY TRANSPORT NETWORKS IN PROTEINS. REVIEWS IN COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1002/9781119518068.ch2] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
29
|
Hodges M, Barahona M, Yaliraki SN. Allostery and cooperativity in multimeric proteins: bond-to-bond propensities in ATCase. Sci Rep 2018; 8:11079. [PMID: 30038211 PMCID: PMC6056424 DOI: 10.1038/s41598-018-27992-z] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2018] [Accepted: 06/13/2018] [Indexed: 11/08/2022] Open
Abstract
Aspartate carbamoyltransferase (ATCase) is a large dodecameric enzyme with six active sites that exhibits allostery: its catalytic rate is modulated by the binding of various substrates at distal points from the active sites. A recently developed method, bond-to-bond propensity analysis, has proven capable of predicting allosteric sites in a wide range of proteins using an energy-weighted atomistic graph obtained from the protein structure and given knowledge only of the location of the active site. Bond-to-bond propensity establishes if energy fluctuations at given bonds have significant effects on any other bond in the protein, by considering their propagation through the protein graph. In this work, we use bond-to-bond propensity analysis to study different aspects of ATCase activity using three different protein structures and sources of fluctuations. First, we predict key residues and bonds involved in the transition between inactive (T) and active (R) states of ATCase by analysing allosteric substrate binding as a source of energy perturbations in the protein graph. Our computational results also indicate that the effect of multiple allosteric binding is non linear: a switching effect is observed after a particular number and arrangement of substrates is bound suggesting a form of long range communication between the distantly arranged allosteric sites. Second, cooperativity is explored by considering a bisubstrate analogue as the source of energy fluctuations at the active site, also leading to the identification of highly significant residues to the T ↔ R transition that enhance cooperativity across active sites. Finally, the inactive (T) structure is shown to exhibit a strong, non linear communication between the allosteric sites and the interface between catalytic subunits, rather than the active site. Bond-to-bond propensity thus offers an alternative route to explain allosteric and cooperative effects in terms of detailed atomistic changes to individual bonds within the protein, rather than through phenomenological, global thermodynamic arguments.
Collapse
Affiliation(s)
- Maxwell Hodges
- Department of Chemistry, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
| | - Mauricio Barahona
- Department of Mathematics, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
| | - Sophia N Yaliraki
- Department of Chemistry, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom.
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom.
| |
Collapse
|
30
|
Yao XQ, Momin M, Hamelberg D. Elucidating Allosteric Communications in Proteins with Difference Contact Network Analysis. J Chem Inf Model 2018; 58:1325-1330. [DOI: 10.1021/acs.jcim.8b00250] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Mohamed Momin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
31
|
Álvarez LD, Presman DM, Pecci A. Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain suggest a role of the lever-arm mobility in transcriptional output. PLoS One 2017; 12:e0189588. [PMID: 29244866 PMCID: PMC5731742 DOI: 10.1371/journal.pone.0189588] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 11/29/2017] [Indexed: 11/19/2022] Open
Abstract
One of the first and essential steps in gene expression regulation involves the recruitment of transcription factors (TFs) to specific response elements located at enhancers and/or promoters of targeted genes. These DNA elements have a certain variability in both sequence and length, which may affect the final transcriptional output. The molecular mechanisms in which TFs integrate the subtle differences within specific recognition sequences to offer different transcriptional responses is still largely unknown. Here we used molecular dynamics simulations to study the DNA binding behavior of the glucocorticoid receptor (GR), a ligand-regulated TF with pleiotropic effects in almost all cells. By comparing the behavior of the wild type receptor and a well characterized Ala477Thr substitution within the rat GR DNA binding domain, we found that the region that connects the two-zinc fingers (i.e. the lever arm) would likely play a key role in GR transcriptional output.
Collapse
Affiliation(s)
- Lautaro Damián Álvarez
- Universidad de Buenos Aires, CONICET, UMYMFOR and Departamento de Química Orgánica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina
- * E-mail:
| | - Diego Martín Presman
- Laboratory of Receptor Biology and Gene Expression, Building 41, 41 Library Drive, National Cancer Institute, National Institutes of Health, Bethesda, MD, Unitec States of America
| | - Adalí Pecci
- Universidad de Buenos Aires, CONICET, IFIBYNE and Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Buenos Aires, Argentina
| |
Collapse
|
32
|
Gnanasekaran R. Computational study to understand the energy transfer pathways within amicyanin. J Mol Graph Model 2017; 78:88-95. [PMID: 29054098 DOI: 10.1016/j.jmgm.2017.09.023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2017] [Revised: 09/28/2017] [Accepted: 09/30/2017] [Indexed: 11/27/2022]
Abstract
Vibrational energy diffusivities between the residues present in Amicyanin copper protein are calculated and presented in form of communication map. From those results energy flow pathways from the copper metal ion to the inter protein residue Glu31 are identified. Our finding suggests many different pathways are possible and copper metal ion in oxidized and reduced state switches the pathways. Our finding also suggests the cooperative nature of surrounding residues and water molecules towards selecting the pathways. The major transport channels in the oxidised state are, Cu2+---> MET28---> LYS29---> TYR30---> GLU31 and Cu2+---> MET98---> TYR30--- > GLU31. And in the reduced state Cu+---> CYS9---> TYR30---> GLU31 and Cu+---> MET28---> LYS2---> TYR30---> GLU31. We studied further the interaction energies between the copper ion and neighbouring residues using B3LYP/QZVP method. Both the methods complement each other in predicting the energy flow pathways and the cooperative nature of residues.
Collapse
|
33
|
Abstract
Allostery represents a fundamental mechanism of biological regulation that is mediated via long-range communication between distant protein sites. Although little is known about the underlying dynamical process, recent time-resolved infrared spectroscopy experiments on a photoswitchable PDZ domain (PDZ2S) have indicated that the allosteric transition occurs on multiple timescales. Here, using extensive nonequilibrium molecular dynamics simulations, a time-dependent picture of the allosteric communication in PDZ2S is developed. The simulations reveal that allostery amounts to the propagation of structural and dynamical changes that are genuinely nonlinear and can occur in a nonlocal fashion. A dynamic network model is constructed that illustrates the hierarchy and exceeding structural heterogeneity of the process. In compelling agreement with experiment, three physically distinct phases of the time evolution are identified, describing elastic response ([Formula: see text] ns), inelastic reorganization ([Formula: see text] ns), and structural relaxation ([Formula: see text]s). Issues such as the similarity to downhill folding as well as the interpretation of allosteric pathways are discussed.
Collapse
|
34
|
Lee Y, Kim S, Choi S, Hyeon C. Ultraslow Water-Mediated Transmembrane Interactions Regulate the Activation of A2A Adenosine Receptor. Biophys J 2017; 111:1180-1191. [PMID: 27653477 DOI: 10.1016/j.bpj.2016.08.002] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Revised: 07/09/2016] [Accepted: 08/02/2016] [Indexed: 01/04/2023] Open
Abstract
Water molecules inside a G-protein coupled receptor (GPCR) have recently been spotlighted in a series of crystal structures. To decipher the dynamics and functional roles of internal water molecules in GPCR activity, we studied the A2A adenosine receptor using microsecond molecular-dynamics simulations. Our study finds that the amount of water flux across the transmembrane (TM) domain varies depending on the receptor state, and that the water molecules of the TM channel in the active state flow three times more slowly than those in the inactive state. Depending on the location in solvent-protein interface as well as the receptor state, the average residence time of water in each residue varies from ∼O(10(2)) ps to ∼O(10(2)) ns. Especially, water molecules, exhibiting ultraslow relaxation (∼O(10(2)) ns) in the active state, are found around the microswitch residues that are considered activity hotspots for GPCR function. A continuous allosteric network spanning the TM domain, arising from water-mediated contacts, is unique in the active state, underscoring the importance of slow water molecules in the activation of GPCRs.
Collapse
Affiliation(s)
- Yoonji Lee
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul, Korea
| | - Songmi Kim
- Korea Institute for Advanced Study, Seoul, Korea
| | - Sun Choi
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul, Korea.
| | | |
Collapse
|
35
|
Salamanca Viloria J, Allega MF, Lambrughi M, Papaleo E. An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass. Sci Rep 2017; 7:2838. [PMID: 28588190 PMCID: PMC5460117 DOI: 10.1038/s41598-017-01498-6] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 03/28/2017] [Indexed: 02/05/2023] Open
Abstract
Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.
Collapse
Affiliation(s)
- Juan Salamanca Viloria
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Maria Francesca Allega
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark.
| |
Collapse
|
36
|
Ahalawat N, Murarka RK. Molecular Mechanism of Nucleotide-Dependent Allosteric Regulation in AMP-Activated Protein Kinase. J Phys Chem B 2017; 121:2919-2930. [PMID: 28345916 DOI: 10.1021/acs.jpcb.6b11223] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The AMP-activated protein kinase (AMPK), a central enzyme in the regulation of energy homeostasis, is an important drug target for type 2 diabetes, obesity, and cancer. Binding of adenosine nucleotides to the regulatory γ-subunit tightly regulates the activity of this enzyme. Though recent crystal structures of AMPK have provided important insights into the allosteric activation of AMPK, molecular details of the regulatory mechanism of AMPK activation is still elusive. Here, we have performed extensive all-atom molecular dynamics (MD) simulations and shown that the kinase domain (KD) and γ-subunit come closer resulting in a more compact heterotrimeric AMPK complex in AMP-bound state compared to the ATP-bound state. The binding of ATP at site 3 of regulatory γ-subunit allosterically inhibits AMPK by destabilizing different regulatory regions of α-subunit: the autoinhibitory domain, the linker region, and the activation loop of the kinase core. The catalytically important residues experience a change in mechanical stress, and major rearrangements in community structure derived from residue-residue interaction energy-based network are observed in KD and α-linker region upon binding of different nucleotides. Our results also highlight the role of conserved charged residues forming an ionic network near the site 3 of γ-subunit in allosteric communications.
Collapse
Affiliation(s)
- Navjeet Ahalawat
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal , Bhopal By-pass Road, Bhauri, Bhopal 462 066, MP, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal , Bhopal By-pass Road, Bhauri, Bhopal 462 066, MP, India
| |
Collapse
|
37
|
Karain WI, Qaraeen NI. The adaptive nature of protein residue networks. Proteins 2017; 85:917-923. [PMID: 28168745 DOI: 10.1002/prot.25261] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Revised: 01/19/2017] [Accepted: 01/30/2017] [Indexed: 11/09/2022]
Abstract
Protein residue networks PRNs are used to describe proteins. These networks are usually based on an average structure for the protein. However, proteins are dynamic entities that are affected by their surroundings. In this work, we study the effect of temperatures above and below the protein dynamical transition temperature(≈200 K), on three important network parameters gleaned from weighted PRNs for the solvated β-lactamase inhibitory protein BLIP: the betweenness centrality B, the closeness centrality C, and the clustering coefficient CC. The B and C values will be extracted for each node from PRNs at six different temperatures: 150 K, 180 K, 200 K, 220 K, 250 K, and 310 K respectively. The average value for the CC for each PRN will also be calculated at each temperature, respectively. We find that at temperatures ≤200 K, the network nodes with the most significant B and C values tend to have lower relative solvent accessibility RSA values, and to fall within the protein secondary structure elements (α helices and β sheets). At temperatures >200 K, the significant nodes in terms of B and C tend to have larger RSA values, and to fall on the connecting loops in the protein. The average CC decreases in value for the PRNs up to 200 K, and then remains basically constant above 200 K. This clearly shows that any conclusions based on static PRNs should be handled with care. The dynamic nature of proteins and its coupling to the surrounding environment should be taken into consideration when using the PRN paradigm. Proteins 2017; 85:917-923. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Wael I Karain
- Department of Physics, Birzeit University, Birzeit, Palestine
| | - Nael I Qaraeen
- Department of Computer Science, Birzeit University, Birzeit, Palestine
| |
Collapse
|
38
|
Long S, Tian P. Nonlinear backbone torsional pair correlations in proteins. Sci Rep 2016; 6:34481. [PMID: 27708342 PMCID: PMC5052647 DOI: 10.1038/srep34481] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2016] [Accepted: 09/14/2016] [Indexed: 12/27/2022] Open
Abstract
Protein allostery requires dynamical structural correlations. Physical origin of which, however, remain elusive despite intensive studies during last two and half decades. Based on analysis of molecular dynamics (MD) simulation trajectories for ten proteins with different sizes and folds, we found that nonlinear backbone torsional pair (BTP) correlations, which are mainly spatially long-ranged and are dominantly executed by loop residues, exist extensively in most analyzed proteins. Examination of torsional motion for correlated BTPs suggested that such nonlinear correlations are mainly associated aharmonic torsional state transitions and in some cases strongly anisotropic local torsional motion of participating torsions, and occur on widely different and relatively longer time scales. In contrast, correlations between backbone torsions in stable α helices and β strands are mainly linear and spatially short-ranged, and are more likely to associate with harmonic local torsional motion. Further analysis revealed that the direct cause of nonlinear contributions are heterogeneous linear correlations. These findings implicate a general search strategy for novel allosteric modulation sites of protein activities.
Collapse
Affiliation(s)
- Shiyang Long
- School of Life Sciences, Jilin University, Changchun, 130012 China
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun, 130012 China.,MOE Key Laboratory of Molecular Enzymology and Engineering, Jilin University, Changchun, 130012 China
| |
Collapse
|
39
|
Amor BRC, Schaub MT, Yaliraki SN, Barahona M. Prediction of allosteric sites and mediating interactions through bond-to-bond propensities. Nat Commun 2016; 7:12477. [PMID: 27561351 PMCID: PMC5007447 DOI: 10.1038/ncomms12477] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2015] [Accepted: 07/05/2016] [Indexed: 11/09/2022] Open
Abstract
Allostery is a fundamental mechanism of biological regulation, in which binding of a molecule at a distant location affects the active site of a protein. Allosteric sites provide targets to fine-tune protein activity, yet we lack computational methodologies to predict them. Here we present an efficient graph-theoretical framework to reveal allosteric interactions (atoms and communication pathways strongly coupled to the active site) without a priori information of their location. Using an atomistic graph with energy-weighted covalent and weak bonds, we define a bond-to-bond propensity quantifying the non-local effect of instantaneous bond fluctuations propagating through the protein. Significant interactions are then identified using quantile regression. We exemplify our method with three biologically important proteins: caspase-1, CheY, and h-Ras, correctly predicting key allosteric interactions, whose significance is additionally confirmed against a reference set of 100 proteins. The almost-linear scaling of our method renders it suitable for high-throughput searches for candidate allosteric sites.
Collapse
Affiliation(s)
- B R C Amor
- Department of Chemistry, Imperial College London, London SW7 2AZ, UK.,Institute of Chemical Biology, Imperial College London, London SW7 2AZ, UK
| | - M T Schaub
- Department of Mathematics, Imperial College London, London SW7 2AZ, UK
| | - S N Yaliraki
- Department of Chemistry, Imperial College London, London SW7 2AZ, UK.,Institute of Chemical Biology, Imperial College London, London SW7 2AZ, UK
| | - M Barahona
- Institute of Chemical Biology, Imperial College London, London SW7 2AZ, UK.,Department of Mathematics, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
40
|
Abstract
We examine energy dynamics in the unliganded and liganded states of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI), which exhibits cooperativity mediated by the cluster of water molecules at the interface upon ligand binding and dissociation. We construct and analyze a dynamic network in which nodes representing the residues, hemes, and water cluster are connected by edges that represent energy transport times, as well as a nonbonded network (NBN) indicating regions that respond rapidly to local strain within the protein via nonbonded interactions. One of the two largest NBNs includes the Lys30-Asp89 salt bridge critical for stabilizing the dimer. The other includes the hemes and surrounding residues, as well as, in the unliganded state, the cluster of water molecules between the globules. Energy transport in the protein appears to be controlled by the Lys30-Asp89 salt bridge critical for stabilizing the dimer, as well as the interface water cluster in the unliganded state. Possible connections between energy transport dynamics in response to local strain identified here and allosteric transitions in HbI are discussed.
Collapse
Affiliation(s)
- David M Leitner
- Department of Chemistry and Chemical Physics Program, University of Nevada , Reno, Nevada 89557, United States
| |
Collapse
|
41
|
Ricci CG, Silveira RL, Rivalta I, Batista VS, Skaf MS. Allosteric Pathways in the PPARγ-RXRα nuclear receptor complex. Sci Rep 2016; 6:19940. [PMID: 26823026 PMCID: PMC4731802 DOI: 10.1038/srep19940] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Accepted: 12/07/2015] [Indexed: 02/08/2023] Open
Abstract
Understanding the nature of allostery in DNA-nuclear receptor (NR) complexes is of fundamental importance for drug development since NRs regulate the transcription of a myriad of genes in humans and other metazoans. Here, we investigate allostery in the peroxisome proliferator-activated/retinoid X receptor heterodimer. This important NR complex is a target for antidiabetic drugs since it binds to DNA and functions as a transcription factor essential for insulin sensitization and lipid metabolism. We find evidence of interdependent motions of Ω-loops and PPARγ-DNA binding domain with contacts susceptible to conformational changes and mutations, critical for regulating transcriptional functions in response to sequence-dependent DNA dynamics. Statistical network analysis of the correlated motions, observed in molecular dynamics simulations, shows preferential allosteric pathways with convergence centers comprised of polar amino acid residues. These findings are particularly relevant for the design of allosteric modulators of ligand-dependent transcription factors.
Collapse
Affiliation(s)
- Clarisse G Ricci
- Institute of Chemistry, University of Campinas-UNICAMP, Cx. P. 6154, Campinas SP 13084-862, Brazil
| | - Rodrigo L Silveira
- Institute of Chemistry, University of Campinas-UNICAMP, Cx. P. 6154, Campinas SP 13084-862, Brazil
| | - Ivan Rivalta
- Université de Lyon, CNRS, Laboratoire de Chimie, École Normale Supérieure de Lyon, UMR 5182, 46 Allée d'Italie, 69364 Lyon, Cedex 07, France.,Department of Chemistry, Yale University, P.O. Box 208107, New Haven, CT 06520-8167, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, CT 06520-8167, United States
| | - Munir S Skaf
- Institute of Chemistry, University of Campinas-UNICAMP, Cx. P. 6154, Campinas SP 13084-862, Brazil
| |
Collapse
|
42
|
Affiliation(s)
- Andre A. S. T. Ribeiro
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Vanessa Ortiz
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|
43
|
Dynamic Transmission of Protein Allostery without Structural Change: Spatial Pathways or Global Modes? Biophys J 2015; 109:1240-50. [PMID: 26338443 DOI: 10.1016/j.bpj.2015.08.009] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2015] [Revised: 07/02/2015] [Accepted: 08/07/2015] [Indexed: 11/23/2022] Open
Abstract
We examine the contrast between mechanisms for allosteric signaling that involve structural change, and those that do not, from the perspective of allosteric pathways. In particular we treat in detail the case of fluctuation-allostery by which amplitude modulation of the thermal fluctuations of the elastic normal modes conveys the allosteric signal, and address the question of what an allosteric pathway means in this case. We find that a perturbation theory of thermal elastic solids and nonperturbative approach (by super-coarse-graining elasticity into internal bending modes) have opposite signatures in their structure of correlated pathways. We illustrate the effect from analysis of previous results from GlxR of Corynebacterium glutamicum, an example of the CRP/FNR transcription family of allosteric homodimers. We find that the visibility of both correlated pathways and disconnected sites of correlated motion in this protein suggests that mechanisms of local elastic stretch and bend are recruited for the purpose of creating and controlling allosteric cooperativity.
Collapse
|
44
|
Tiberti M, Invernizzi G, Papaleo E. (Dis)similarity Index To Compare Correlated Motions in Molecular Simulations. J Chem Theory Comput 2015; 11:4404-14. [PMID: 26575932 DOI: 10.1021/acs.jctc.5b00512] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Molecular dynamics (MD) simulations are widely used to complement or guide experimental studies in the characterization of protein dynamics, thanks to improvements in force-field accuracy, along with in the software and hardware to sample the conformational landscape of proteins. Among the different applications of MD simulations, the study of correlated motions is largely employed for different purposes. Several metrics have been developed to describe correlated motions in the MD ensemble, such as methods based on Pearson Correlation or Mutual Information. Cross-correlation analysis of MD trajectories is indeed appealing not only to identify residues characterized by coupled fluctuations in protein structures but also since it can be used to extrapolate motions along directions in which major conformational changes should occur, for example on longer time scales than the ones that are actually simulated. Nevertheless, most of the MD studies employ average correlation maps and mostly in a qualitative way, even when different systems or different replicates of the same system are compared. The broad application of correlation metrics in the analysis of MD simulations, especially for comparative purposes, requires a step forward toward more quantitative and accurate comparisons. We thus here employed a simple but effective index, which is based on a normalized Frobenius norm of the differences between protein correlation maps, to compare correlated motions. We applied this index for a quantitative comparison of correlated motions from MD simulations of seven proteins of different size and fold. We also employed the index to assess the robustness of correlation description when multi-replicate MD simulations of a same system are used, and we compared our index to metrics for comparison of structural ensembles such as Root Mean Square Inner Product and the Bhattacharyya Coefficient.
Collapse
Affiliation(s)
- Matteo Tiberti
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Gaetano Invernizzi
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Elena Papaleo
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| |
Collapse
|
45
|
Ribeiro AAST, Ortiz V. MDN: A Web Portal for Network Analysis of Molecular Dynamics Simulations. Biophys J 2015; 109:1110-6. [PMID: 26143656 DOI: 10.1016/j.bpj.2015.06.013] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2015] [Revised: 05/08/2015] [Accepted: 06/03/2015] [Indexed: 01/01/2023] Open
Abstract
We introduce a web portal that employs network theory for the analysis of trajectories from molecular dynamics simulations. Users can create protein energy networks following methodology previously introduced by our group, and can identify residues that are important for signal propagation, as well as measure the efficiency of signal propagation by calculating the network coupling. This tool, called MDN, was used to characterize signal propagation in Escherichia coli heat-shock protein 70-kDa. Two variants of this protein experimentally shown to be allosterically active exhibit higher network coupling relative to that of two inactive variants. In addition, calculations of partial coupling suggest that this quantity could be used as part of the criteria to determine pocket druggability in drug discovery studies.
Collapse
Affiliation(s)
| | - Vanessa Ortiz
- Department of Chemical Engineering, Columbia University, New York, New York.
| |
Collapse
|
46
|
Papaleo E. Integrating atomistic molecular dynamics simulations, experiments, and network analysis to study protein dynamics: strength in unity. Front Mol Biosci 2015; 2:28. [PMID: 26075210 PMCID: PMC4445042 DOI: 10.3389/fmolb.2015.00028] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2015] [Accepted: 05/08/2015] [Indexed: 12/11/2022] Open
Abstract
In the last years, we have been observing remarkable improvements in the field of protein dynamics. Indeed, we can now study protein dynamics in atomistic details over several timescales with a rich portfolio of experimental and computational techniques. On one side, this provides us with the possibility to validate simulation methods and physical models against a broad range of experimental observables. On the other side, it also allows a complementary and comprehensive view on protein structure and dynamics. What is needed now is a better understanding of the link between the dynamic properties that we observe and the functional properties of these important cellular machines. To make progresses in this direction, we need to improve the physical models used to describe proteins and solvent in molecular dynamics, as well as to strengthen the integration of experiments and simulations to overcome their own limitations. Moreover, now that we have the means to study protein dynamics in great details, we need new tools to understand the information embedded in the protein ensembles and in their dynamic signature. With this aim in mind, we should enrich the current tools for analysis of biomolecular simulations with attention to the effects that can be propagated over long distances and are often associated to important biological functions. In this context, approaches inspired by network analysis can make an important contribution to the analysis of molecular dynamics simulations.
Collapse
Affiliation(s)
- Elena Papaleo
- Structural Biology and Nuclear Magnetic Resonance Laboratory, Department of Biology, University of Copenhagen Copenhagen, Denmark
| |
Collapse
|
47
|
Karain WI, Qaraeen NI. Weighted protein residue networks based on joint recurrences between residues. BMC Bioinformatics 2015; 16:173. [PMID: 26003989 PMCID: PMC4491895 DOI: 10.1186/s12859-015-0621-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2015] [Accepted: 05/18/2015] [Indexed: 01/07/2023] Open
Abstract
BACKGROUND Weighted and un-weighted protein residue networks can predict key functional residues in proteins based on the closeness centrality C and betweenness centrality B values for each residue. A static snapshot of the protein structure, and a cutoff distance, are used to define edges between the network nodes. In this work we apply the weighted network approach to study the β-Lactamase Inhibitory Protein (BLIP). Joint recurrences extracted from molecular dynamics MD trajectory positions of the protein residue carbon alpha atoms are used to define edge weights between nodes, and no cutoff distance is used. The results for B and C from our approach are compared with those extracted from an un-weighted network, and a weighted network that uses interatomic contacts to define edge weights between nodes, respectively. RESULTS The joint recurrence weighted network approach performs well in pointing out key protein residues. Furthermore, it seems to emphasize residues with medium to high relative solvent accessibility that lie in loop regions between secondary structure elements of the protein. CONCLUSIONS Protein residue networks that use joint recurrences extracted from molecular dynamics simulations of a solvated protein perform well in pointing to hotspot residues and hotspot clusters. This approach uses no distance cutoff threshold, and does not exclude any interactions between the residues, including water-mediated interactions.
Collapse
Affiliation(s)
- Wael I Karain
- Department of Physics, Birzeit University, Birzeit, Palestine.
| | - Nael I Qaraeen
- Department of Computer Science, Birzeit University, Birzeit, Palestine.
| |
Collapse
|
48
|
Ribeiro AA, Ortiz V. Local elastic constants of LacI and implications for allostery. J Mol Graph Model 2015; 57:106-13. [DOI: 10.1016/j.jmgm.2015.01.013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2014] [Revised: 01/28/2015] [Accepted: 01/29/2015] [Indexed: 11/28/2022]
|
49
|
Ribeiro AAST, Ortiz V. Energy Propagation and Network Energetic Coupling in Proteins. J Phys Chem B 2015; 119:1835-46. [DOI: 10.1021/jp509906m] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andre A. S. T. Ribeiro
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Vanessa Ortiz
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|