51
|
Jain A, Dokholyan NV, Lee AL. Allosteric inactivation of an engineered optogenetic GTPase. Proc Natl Acad Sci U S A 2023; 120:e2219254120. [PMID: 36972433 PMCID: PMC10083549 DOI: 10.1073/pnas.2219254120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 02/14/2023] [Indexed: 03/29/2023] Open
Abstract
Optogenetics is a technique for establishing direct spatiotemporal control over molecular function within living cells using light. Light application induces conformational changes within targeted proteins that produce changes in function. One of the applications of optogenetic tools is an allosteric control of proteins via light-sensing domain (LOV2), which allows direct and robust control of protein function. Computational studies supported by cellular imaging demonstrated that application of light allosterically inhibited signaling proteins Vav2, ITSN, and Rac1, but the structural and dynamic basis of such control has yet to be elucidated by experiment. Here, using NMR spectroscopy, we discover principles of action of allosteric control of cell division control protein 42 (CDC42), a small GTPase involved in cell signaling. Both LOV2 and Cdc42 employ flexibility in their function to switch between "dark"/"lit" or active/inactive states, respectively. By conjoining Cdc42 and phototropin1 LOV2 domains into the bi-switchable fusion Cdc42Lov, application of light-or alternatively, mutation in LOV2 to mimic light absorption-allosterically inhibits Cdc42 downstream signaling. The flow and patterning of allosteric transduction in this flexible system are well suited to observation by NMR. Close monitoring of the structural and dynamic properties of dark versus "lit" states of Cdc42Lov revealed lit-induced allosteric perturbations that extend to Cdc42's downstream effector binding site. Chemical shift perturbations for lit mimic, I539E, have distinct regions of sensitivity, and both the domains are coupled together, leading to bidirectional interdomain signaling. Insights gained from this optoallosteric design will increase our ability to control response sensitivity in future designs.
Collapse
Affiliation(s)
- Abha Jain
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC27599
| | - Nikolay V. Dokholyan
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA17033
- Department of Biochemistry and Molecular Biology, Penn State College of Medicine, Hershey, PA17033
| | - Andrew L. Lee
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC27599
| |
Collapse
|
52
|
Xiao S, Verkhivker GM, Tao P. Machine learning and protein allostery. Trends Biochem Sci 2023; 48:375-390. [PMID: 36564251 PMCID: PMC10023316 DOI: 10.1016/j.tibs.2022.12.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Revised: 11/16/2022] [Accepted: 12/02/2022] [Indexed: 12/23/2022]
Abstract
The fundamental biological importance and complexity of allosterically regulated proteins stem from their central role in signal transduction and cellular processes. Recently, machine-learning approaches have been developed and actively deployed to facilitate theoretical and experimental studies of protein dynamics and allosteric mechanisms. In this review, we survey recent developments in applications of machine-learning methods for studies of allosteric mechanisms, prediction of allosteric effects and allostery-related physicochemical properties, and allosteric protein engineering. We also review the applications of machine-learning strategies for characterization of allosteric mechanisms and drug design targeting SARS-CoV-2. Continuous development and task-specific adaptation of machine-learning methods for protein allosteric mechanisms will have an increasingly important role in bridging a wide spectrum of data-intensive experimental and theoretical technologies.
Collapse
Affiliation(s)
- Sian Xiao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX 75205, USA.
| | - Gennady M Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA; Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX 75205, USA.
| |
Collapse
|
53
|
Kish M, Subramanian S, Smith V, Lethbridge N, Cole L, Vollmer F, Bond NJ, Phillips JJ. Allosteric Regulation of Glycogen Phosphorylase by Order/Disorder Transition of the 250' and 280s Loops. Biochemistry 2023; 62:1360-1368. [PMID: 36989206 PMCID: PMC10116597 DOI: 10.1021/acs.biochem.2c00671] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Allostery is a fundamental mechanism of protein activation, yet the precise dynamic changes that underlie functional regulation of allosteric enzymes, such as glycogen phosphorylase (GlyP), remain poorly understood. Despite being the first allosteric enzyme described, its structural regulation is still a challenging problem: the key regulatory loops of the GlyP active site (250' and 280s) are weakly stable and often missing density or have large b-factors in structural models. This led to the longstanding hypothesis that GlyP regulation is achieved through gating of the active site by (dis)order transitions, as first proposed by Barford and Johnson. However, testing this requires a quantitative measurement of weakly stable local structure which, to date, has been technically challenging in such a large protein. Hydrogen-deuterium-exchange mass spectrometry (HDX-MS) is a powerful tool for studying protein dynamics, and millisecond HDX-MS has the ability to measure site-localized stability differences in weakly stable structures, making it particularly valuable for investigating allosteric regulation in GlyP. Here, we used millisecond HDX-MS to measure the local structural perturbations of glycogen phosphorylase b (GlyPb), the phosphorylated active form (GlyPa), and the inhibited glucose-6 phosphate complex (GlyPb:G6P) at near-amino acid resolution. Our results support the Barford and Johnson hypothesis for GlyP regulation by providing insight into the dynamic changes of the key regulatory loops.
Collapse
Affiliation(s)
- Monika Kish
- Living Systems Institute, Department of Biosciences, University of Exeter, Stocker Road, Exeter, EX4 4QD, U.K
| | - Sivaraman Subramanian
- Living Systems Institute, Department of Physics, University of Exeter, Stocker Road, Exeter, EX4 6QD, U.K
| | | | | | - Lindsay Cole
- Applied Photophysics Ltd, Leatherhead, KT227BA, U.K
| | - Frank Vollmer
- Living Systems Institute, Department of Physics, University of Exeter, Stocker Road, Exeter, EX4 6QD, U.K
| | - Nicholas J Bond
- Analytical Sciences, Biopharmaceutical Development, BioPharmaceuticals R&D, AstraZeneca, Milstein Building, Granta Park, Cambridge, CB21 6GH, U.K
| | - Jonathan J Phillips
- Living Systems Institute, Department of Biosciences, University of Exeter, Stocker Road, Exeter, EX4 4QD, U.K
- Alan Turing Institute, British Library, London, NW1 2DB, U.K
| |
Collapse
|
54
|
Hu Y, Li C, Hu M, Zhang Z, Fu R, Tang X, Wu T. Allosteric Nucleic Acid Enzyme: A Versatile Stimuli-Responsive Tool for Molecular Computing and Biosensing Nanodevices. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023:e2300207. [PMID: 36978231 DOI: 10.1002/smll.202300207] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 03/09/2023] [Indexed: 06/18/2023]
Abstract
Allostery is a naturally occurring mechanism in which effector binding induces the modulation and fine control of a related biomolecule function. Deoxyribozyme (DNAzyme) with catalytic activity and substrate recognition ability is ideal to be regulated by allosteric strategies. However, the current regulations frequently confront various obstacles, such as severe activity decay, signal leakage, and limited effectors. In this work, a rational regulation strategy for developing versatile effectors-responsive allosteric nucleic acid enzyme (ANAzyme) by introducing an allosteric domain in response to diverse effectors is established. The enzyme-like activity of this re-engineered ANAzyme can be modulated in a more predictable and fine way compared with the previous DNAzyme regulation strategies. Based on the allosteric strategy, the construction of allosterically coregulatory nanodevices and a series of basic logic gates and logic circuits are achieved, demonstrating that the proposed ANAzyme-regulated strategy showed great potential in molecular computing. Given these facts, the rational design of ANAzyme with the allosteric domain presented here can expand the available toolbox to develop a variety of stimuli-responsive allosteric DNA materials, including molecular machines, computing systems, biosensing platforms, and gene-silencing tools.
Collapse
Affiliation(s)
- Yuqiang Hu
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
- State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University, Beijing, 100191, P. R. China
| | - Changjiang Li
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
| | - Minghao Hu
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
| | - Zhen Zhang
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
| | - Ruolan Fu
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
| | - Xinjing Tang
- State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University, Beijing, 100191, P. R. China
| | - Tongbo Wu
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430030, P. R. China
| |
Collapse
|
55
|
Wand AJ. Deep mining of the protein energy landscape. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2023; 10:020901. [PMID: 37124940 PMCID: PMC10147411 DOI: 10.1063/4.0000180] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Accepted: 03/30/2023] [Indexed: 05/03/2023]
Abstract
For over half a century, it has been known that protein molecules naturally undergo extensive structural fluctuations, and that these internal motions are intimately related to their functional properties. The energy landscape view has provided a powerful framework for describing the various physical states that proteins visit during their lifetimes. This Perspective focuses on the commonly neglected and often disparaged axis of the protein energy landscape: entropy. Initially seen largely as a barrier to functionally relevant states of protein molecules, it has recently become clear that proteins retain considerable conformational entropy in the "native" state, and that this entropy can and often does contribute significantly to the free energy of fundamental protein properties, processes, and functions. NMR spectroscopy, molecular dynamics simulations, and emerging crystallographic views have matured in parallel to illuminate dynamic disorder of the "ground state" of proteins and their importance in not only transiting between biologically interesting structures but also greatly influencing their stability, cooperativity, and contribution to critical properties such as allostery.
Collapse
|
56
|
Gianni S, Jemth P. Allostery Frustrates the Experimentalist. J Mol Biol 2023; 435:167934. [PMID: 36586463 DOI: 10.1016/j.jmb.2022.167934] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 12/20/2022] [Accepted: 12/22/2022] [Indexed: 12/29/2022]
Abstract
Proteins interact with other proteins, with nucleic acids, lipids, carbohydrates and various small molecules in the living cell. These interactions have been quantified and structurally characterized in numerous studies such that we today have a comprehensive picture of protein structure and function. However, proteins are dynamic and even folded proteins are likely more heterogeneous than they appear in most descriptions. One property of proteins that relies on dynamics and heterogeneity is allostery, the ability of a protein to change structure and function upon ligand binding to an allosteric site. Over the last decades the concept of allostery was broadened to embrace all types of long-range interactions across a protein including purely entropic changes without a conformational change in single protein domains. But with this re-definition came a problem: How do we measure allostery? In this opinion, we discuss some caveats arising from the quantitative description of single-domain allostery from an experimental perspective and how the limitations cannot be separated from the definition of allostery per se. Furthermore, we attempt to tie together allostery with the concept of frustration in an effort to investigate the links between these two complex, and yet general, properties of proteins. We arrive at the conclusion that the sensitivity to perturbation of allosteric networks in single protein domains is too large for the networks to be of significant biological relevance.
Collapse
Affiliation(s)
- Stefano Gianni
- Istituto Pasteur-Fondazione Cenci Bolognetti and Istituto di Biologia e Patologia Molecolari del CNR, Dipartimento di Scienze Biochimiche "A. Rossi Fanelli," Sapienza Università di Roma, 00185 Rome, Italy.
| | - Per Jemth
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123 Uppsala, Sweden.
| |
Collapse
|
57
|
Madan LK, Welsh CL, Kornev AP, Taylor SS. The "violin model": Looking at community networks for dynamic allostery. J Chem Phys 2023; 158:081001. [PMID: 36859094 PMCID: PMC9957607 DOI: 10.1063/5.0138175] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2022] [Accepted: 02/03/2023] [Indexed: 02/09/2023] Open
Abstract
Allosteric regulation of proteins continues to be an engaging research topic for the scientific community. Models describing allosteric communication have evolved from focusing on conformation-based descriptors of protein structural changes to appreciating the role of internal protein dynamics as a mediator of allostery. Here, we explain a "violin model" for allostery as a contemporary method for approaching the Cooper-Dryden model based on redistribution of protein thermal fluctuations. Based on graph theory, the violin model makes use of community network analysis to functionally cluster correlated protein motions obtained from molecular dynamics simulations. This Review provides the theory and workflow of the methodology and explains the application of violin model to unravel the workings of protein kinase A.
Collapse
Affiliation(s)
- Lalima K. Madan
- Author to whom correspondence should be addressed: and . Telephone: 843.792.4525. Fax: 843.792.0481
| | - Colin L. Welsh
- Department of Cell and Molecular Pharmacology and Experimental Therapeutics, Medical University of South Carolina, 173 Ashley Ave., Charleston, South Carolina 29425, USA
| | - Alexandr P. Kornev
- Department of Pharmacology, University of California San Diego, 9500 Gilman Drive, San Diego, California, 92093, USA
| | | |
Collapse
|
58
|
Dynamics and mechanistic interpretations of nonribosomal peptide synthetase cyclization domains. Curr Opin Chem Biol 2023; 72:102228. [PMID: 36402006 DOI: 10.1016/j.cbpa.2022.102228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 10/10/2022] [Accepted: 10/12/2022] [Indexed: 11/17/2022]
Abstract
Ox-/thiazoline groups in nonribosomal peptides are formed by a variant of peptide-forming condensation domains called heterocyclization (Cy) domains and appear in a range of pharmaceutically important natural products and virulence factors. Recent cryo-EM, crystallographic, and NMR studies of Cy domains make it opportune to revisit outstanding questions regarding their molecular mechanisms. This review covers structural and dynamical findings about Cy domains that will inform future bioengineering efforts and our understanding of natural product synthesis.
Collapse
|
59
|
Liu J, Amaral LAN, Keten S. A new approach for extracting information from protein dynamics. Proteins 2023; 91:183-195. [PMID: 36094321 PMCID: PMC9844508 DOI: 10.1002/prot.26421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 08/25/2022] [Accepted: 09/06/2022] [Indexed: 01/19/2023]
Abstract
Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.
Collapse
Affiliation(s)
- Jenny Liu
- Department of Mechanical Engineering, Northwestern University
| | - Luís A. N. Amaral
- Department of Chemical and Biological Engineering, Northwestern University
| | - Sinan Keten
- Department of Mechanical Engineering, Northwestern University
| |
Collapse
|
60
|
Modi T, Campitelli P, Heyden M, Ozkan SB. Correlated Evolution of Low-Frequency Vibrations and Function in Enzymes. J Phys Chem B 2023; 127:616-622. [PMID: 36633931 DOI: 10.1021/acs.jpcb.2c05983] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Previous studies of the flexibility of ancestral proteins suggest that proteins evolve their function by altering their native state ensemble. Here, we propose a more direct method to analyze such changes during protein evolution by comparing thermally activated vibrations at frequencies below 6 THz, which report on the dynamics of collective protein modes. We analyzed the backbone vibrational density of states of ancestral and extant β-lactamases and thioredoxins and observed marked changes in the vibrational spectrum in response to evolution. Coupled with previously observed changes in protein flexibility, the observed shifts of vibrational mode densities suggest that protein dynamics and dynamical allostery are critical factors for the evolution of enzymes with specialized catalytic and biophysical properties.
Collapse
Affiliation(s)
- Tushar Modi
- Department of Physics, Arizona State University, Tempe, Arizona85287, United States
| | - Paul Campitelli
- Department of Physics, Arizona State University, Tempe, Arizona85287, United States
| | - Matthias Heyden
- School of Molecular Sciences, Arizona State University, Tempe, Arizona85287, United States
| | - S Banu Ozkan
- Department of Physics, Arizona State University, Tempe, Arizona85287, United States
| |
Collapse
|
61
|
Krishnan K, Tian H, Tao P, Verkhivker GM. Probing conformational landscapes and mechanisms of allosteric communication in the functional states of the ABL kinase domain using multiscale simulations and network-based mutational profiling of allosteric residue potentials. J Chem Phys 2022; 157:245101. [PMID: 36586979 PMCID: PMC11184971 DOI: 10.1063/5.0133826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Accepted: 12/05/2022] [Indexed: 12/12/2022] Open
Abstract
In the current study, multiscale simulation approaches and dynamic network methods are employed to examine the dynamic and energetic details of conformational landscapes and allosteric interactions in the ABL kinase domain that determine the kinase functions. Using a plethora of synergistic computational approaches, we elucidate how conformational transitions between the active and inactive ABL states can employ allosteric regulatory switches to modulate intramolecular communication networks between the ATP site, the substrate binding region, and the allosteric binding pocket. A perturbation-based network approach that implements mutational profiling of allosteric residue propensities and communications in the ABL states is proposed. Consistent with biophysical experiments, the results reveal functionally significant shifts of the allosteric interaction networks in which preferential communication paths between the ATP binding site and substrate regions in the active ABL state become suppressed in the closed inactive ABL form, which in turn features favorable allosteric coupling between the ATP site and the allosteric binding pocket. By integrating the results of atomistic simulations with dimensionality reduction methods and Markov state models, we analyze the mechanistic role of macrostates and characterize kinetic transitions between the ABL conformational states. Using network-based mutational scanning of allosteric residue propensities, this study provides a comprehensive computational analysis of long-range communications in the ABL kinase domain and identifies conserved regulatory hotspots that modulate kinase activity and allosteric crosstalk between the allosteric pocket, ATP binding site, and substrate binding regions.
Collapse
Affiliation(s)
| | - Hao Tian
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75205, USA
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75205, USA
| | - Gennady M. Verkhivker
- Author to whom correspondence should be addressed: . Telephone: 714-516-4586. Fax: 714-532-6048
| |
Collapse
|
62
|
Young BD, Cook ME, Costabile BK, Samanta R, Zhuang X, Sevdalis SE, Varney KM, Mancia F, Matysiak S, Lattman E, Weber DJ. Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery. J Mol Biol 2022; 434:167872. [PMID: 36354074 PMCID: PMC10871162 DOI: 10.1016/j.jmb.2022.167872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Revised: 09/20/2022] [Accepted: 10/24/2022] [Indexed: 11/06/2022]
Abstract
EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.
Collapse
Affiliation(s)
- Brianna D Young
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
| | - Mary E Cook
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
| | - Brianna K Costabile
- Department of Physiology and Cellular Biophysics, Columbia University, New York, NY 10032, USA
| | - Riya Samanta
- Biophysics Graduate Program, University of Maryland, College Park, MD 20742, USA; Fischell Department of Bioengineering, University of Maryland, College Park, MD 20742, USA
| | - Xinhao Zhuang
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
| | - Spiridon E Sevdalis
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
| | - Kristen M Varney
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
| | - Filippo Mancia
- Department of Physiology and Cellular Biophysics, Columbia University, New York, NY 10032, USA
| | - Silvina Matysiak
- Biophysics Graduate Program, University of Maryland, College Park, MD 20742, USA; Fischell Department of Bioengineering, University of Maryland, College Park, MD 20742, USA
| | - Eaton Lattman
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA; Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| | - David J Weber
- The Center for Biomolecular Therapeutics (CBT), Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA; The Institute of Bioscience and Biotechnology Research (IBBR), Rockville, MD 20850, USA.
| |
Collapse
|
63
|
Schay G, Fidy J, Herenyi L. Slow dynamics measured by phosphorescence lifetime reveals global conformational changes in human adult hemoglobin induced by allosteric effectors. PLoS One 2022; 17:e0278417. [PMID: 36454779 PMCID: PMC9714750 DOI: 10.1371/journal.pone.0278417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 11/16/2022] [Indexed: 12/03/2022] Open
Abstract
The mechanism underlying allostery in hemoglobin (Hb) is still not completely understood. Various models describing the action of allosteric effectors on Hb function have been published in the literature. It has also been reported that some allosteric effectors-such as chloride ions, inositol hexaphosphate, 2,3-diphospho-glycerate and bezafibrate-considerably lower the oxygen affinity of Hb. In this context, an important question is the extent to which these changes influence the conformational dynamics of the protein. Earlier, we elaborated a challenging method based on phosphorescence quenching, which makes characterizing protein-internal dynamics possible in the ms time range. The experimental technique involves phosphorescence lifetime measurements in thermal equilibrium at varied temperatures from 10 K up to 273 K, based on the signal of Zn-protoporphyrin substituted for the heme in the β-subunits of Hb. The thermal activation of protein dynamics was observed by the enhancement of phosphorescence quenching attributed to O2 diffusion. It was shown that the thermal activation of protein matrix dynamics was clearly distinguishable from the dynamic activation of the aqueous solvent, and was therefore highly specific for the protein. In the present work, the same method was used to study the changes in the parameters of the dynamic activation of human HbA induced by binding allosteric effectors. We interpreted the phenomenon as phase transition between two states. The fitting of this model to lifetime data yielded the change of energy and entropy in the activation process and the quenching rate in the dynamically activated state. The fitted parameters were particularly sensitive to the presence of allosteric effectors and could be interpreted in line with results from earlier experimental studies. The results suggest that allosteric effectors are tightly coupled to the dynamics of the whole protein, and thus underline the importance of global dynamics in the regulation of Hb function.
Collapse
Affiliation(s)
- Gusztáv Schay
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Judit Fidy
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Levente Herenyi
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
- * E-mail:
| |
Collapse
|
64
|
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
|
65
|
Ali AAAI, Gulzar A, Wolf S, Stock G. Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication. J Phys Chem Lett 2022; 13:9862-9868. [PMID: 36251493 DOI: 10.1021/acs.jpclett.2c02821] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α3-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α3-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns. Performing extensive nonequilibrium molecular dynamics simulations, the modeling of the experiment reproduces the measured time scales and reveals a detailed picture of the allosteric communication in PDZ3. In particular, a correlation analysis identifies a network of contacts connecting the α3-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α3-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.
Collapse
Affiliation(s)
- Ahmed A A I Ali
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Adnan Gulzar
- 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
|
66
|
Leander M, Liu Z, Cui Q, Raman S. Deep mutational scanning and machine learning reveal structural and molecular rules governing allosteric hotspots in homologous proteins. eLife 2022; 11:e79932. [PMID: 36226916 PMCID: PMC9662819 DOI: 10.7554/elife.79932] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 10/13/2022] [Indexed: 01/29/2023] Open
Abstract
A fundamental question in protein science is where allosteric hotspots - residues critical for allosteric signaling - are located, and what properties differentiate them. We carried out deep mutational scanning (DMS) of four homologous bacterial allosteric transcription factors (aTFs) to identify hotspots and built a machine learning model with this data to glean the structural and molecular properties of allosteric hotspots. We found hotspots to be distributed protein-wide rather than being restricted to 'pathways' linking allosteric and active sites as is commonly assumed. Despite structural homology, the location of hotspots was not superimposable across the aTFs. However, common signatures emerged when comparing hotspots coincident with long-range interactions, suggesting that the allosteric mechanism is conserved among the homologs despite differences in molecular details. Machine learning with our large DMS datasets revealed global structural and dynamic properties to be a strong predictor of whether a residue is a hotspot than local and physicochemical properties. Furthermore, a model trained on one protein can predict hotspots in a homolog. In summary, the overall allosteric mechanism is embedded in the structural fold of the aTF family, but the finer, molecular details are sequence-specific.
Collapse
Affiliation(s)
- Megan Leander
- Department of Biochemistry, University of Wisconsin-MadisonMadisonUnited States
| | - Zhuang Liu
- Department of Physics, Boston UniversityBostonUnited States
| | - Qiang Cui
- Department of Physics, Boston UniversityBostonUnited States
- Department of Chemistry, Boston UniversityBostonUnited States
| | - Srivatsan Raman
- Department of Biochemistry, University of Wisconsin-MadisonMadisonUnited States
- Department of Bacteriology, University of Wisconsin-MadisonMadisonUnited States
- Department of Chemical and Biological Engineering, University of Wisconsin-MadisonMadisonUnited States
| |
Collapse
|
67
|
Wen B, Zhang W, Zhang Y, Lei H, Cao Y, Li W, Wang W. Self-Effected Allosteric Coupling and Cooperativity in Hypoxic Response Regulation with Disordered Proteins. J Phys Chem Lett 2022; 13:9201-9209. [PMID: 36170455 DOI: 10.1021/acs.jpclett.2c02065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Hypersensitive regulation of cellular hypoxic response relies on cooperative displacement of one disordered protein (HIF-1α) by another disordered protein (CITED2) from the target in a negative feedback loop. Considering the weak intramolecule coupling in disordered proteins, the molecular mechanism of high cooperativity in the molecular displacement event remains elusive. Herein, we show that disordered proteins utilize a "self-effected allostery" mechanism to achieve high binding cooperativity. Different from the conventional allostery mechanisms shown by many structured or disordered proteins, this mechanism utilizes one part of the disordered protein as the effector to trigger the allosteric coupling and enhance the binding of the remaining part of the same disordered protein, contributing to high cooperativity of the displacement event. The conserved charge motif of CITED2 is the key determinant of the molecular displacement event by serving as the effector of allosteric coupling. Such self-effected allostery provides an efficient strategy to achieve high cooperativity in the molecular events involving disordered proteins.
Collapse
Affiliation(s)
- Bin Wen
- Wenzhou Key Laboratory of Biophysics, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Weiwei Zhang
- Wenzhou Key Laboratory of Biophysics, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
| | - Yangyang Zhang
- Wenzhou Key Laboratory of Biophysics, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Hai Lei
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Yi Cao
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Wenfei Li
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Wei Wang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| |
Collapse
|
68
|
Tan ZW, Tee WV, Guarnera E, Berezovsky IN. AlloMAPS 2: allosteric fingerprints of the AlphaFold and Pfam-trRosetta predicted structures for engineering and design. Nucleic Acids Res 2022; 51:D345-D351. [PMID: 36169226 PMCID: PMC9825619 DOI: 10.1093/nar/gkac828] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 08/26/2022] [Accepted: 09/15/2022] [Indexed: 01/29/2023] Open
Abstract
AlloMAPS 2 is an update of the Allosteric Mutation Analysis and Polymorphism of Signalling database, which contains data on allosteric communication obtained for predicted structures in the AlphaFold database (AFDB) and trRosetta-predicted Pfam domains. The data update contains Allosteric Signalling Maps (ASMs) and Allosteric Probing Maps (APMs) quantifying allosteric effects of mutations and of small probe binding, respectively. To ensure quality of the ASMs and APMs, we performed careful and accurate selection of protein sets containing high-quality predicted structures in both databases for each organism/structure, and the data is available for browsing and download. The data for remaining structures are available for download and should be used at user's discretion and responsibility. We believe these massive data can facilitate both diagnostics and drug design within the precision medicine paradigm. Specifically, it can be instrumental in the analysis of allosteric effects of pathological and rescue mutations, providing starting points for fragment-based design of allosteric effectors. The exhaustive character of allosteric signalling and probing fingerprints will be also useful in future developments of corresponding machine learning applications. The database is freely available at: http://allomaps.bii.a-star.edu.sg.
Collapse
Affiliation(s)
- Zhen Wah Tan
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01, Matrix, 138671, Singapore
| | - Wei-Ven Tee
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01, Matrix, 138671, Singapore
| | - Enrico Guarnera
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01, Matrix, 138671, Singapore
| | - Igor N Berezovsky
- To whom correspondence should be addressed. Tel: +65 6478 8269; Fax: +65 6478 9047;
| |
Collapse
|
69
|
Ghosh RK, Hilario E, Chang CEA, Mueller LJ, Dunn MF. Allosteric regulation of substrate channeling: Salmonella typhimurium tryptophan synthase. Front Mol Biosci 2022; 9:923042. [PMID: 36172042 PMCID: PMC9512447 DOI: 10.3389/fmolb.2022.923042] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Accepted: 08/11/2022] [Indexed: 11/13/2022] Open
Abstract
The regulation of the synthesis of L-tryptophan (L-Trp) in enteric bacteria begins at the level of gene expression where the cellular concentration of L-Trp tightly controls expression of the five enzymes of the Trp operon responsible for the synthesis of L-Trp. Two of these enzymes, trpA and trpB, form an αββα bienzyme complex, designated as tryptophan synthase (TS). TS carries out the last two enzymatic processes comprising the synthesis of L-Trp. The TS α-subunits catalyze the cleavage of 3-indole D-glyceraldehyde 3′-phosphate to indole and D-glyceraldehyde 3-phosphate; the pyridoxal phosphate-requiring β-subunits catalyze a nine-step reaction sequence to replace the L-Ser hydroxyl by indole giving L-Trp and a water molecule. Within αβ dimeric units of the αββα bienzyme complex, the common intermediate indole is channeled from the α site to the β site via an interconnecting 25 Å-long tunnel. The TS system provides an unusual example of allosteric control wherein the structures of the nine different covalent intermediates along the β-reaction catalytic path and substrate binding to the α-site provide the allosteric triggers for switching the αββα system between the open (T) and closed (R) allosteric states. This triggering provides a linkage that couples the allosteric conformational coordinate to the covalent chemical reaction coordinates at the α- and β-sites. This coupling drives the α- and β-sites between T and R conformations to achieve regulation of substrate binding and/or product release, modulation of the α- and β-site catalytic activities, prevention of indole escape from the confines of the active sites and the interconnecting tunnel, and synchronization of the α- and β-site catalytic activities. Here we review recent advances in the understanding of the relationships between structure, function, and allosteric regulation of the complex found in Salmonella typhimurium.
Collapse
Affiliation(s)
- Rittik K. Ghosh
- Department of Biochemistry, University of California, Riverside, Riverside, CA, United States
| | - Eduardo Hilario
- Department of Chemistry, University of California, Riverside, Riverside, CA, United States
| | - Chia-en A. Chang
- Department of Chemistry, University of California, Riverside, Riverside, CA, United States
| | - Leonard J. Mueller
- Department of Chemistry, University of California, Riverside, Riverside, CA, United States
- *Correspondence: Leonard J. Mueller, ; Michael F. Dunn,
| | - Michael F. Dunn
- Department of Biochemistry, University of California, Riverside, Riverside, CA, United States
- *Correspondence: Leonard J. Mueller, ; Michael F. Dunn,
| |
Collapse
|
70
|
Panecka-Hofman J, Poehner I, Wade R. Anti-trypanosomatid structure-based drug design - lessons learned from targeting the folate pathway. Expert Opin Drug Discov 2022; 17:1029-1045. [PMID: 36073204 DOI: 10.1080/17460441.2022.2113776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
INTRODUCTION Trypanosomatidic parasitic infections of humans and animals caused by Trypanosoma brucei, Trypanosoma cruzi, and Leishmania species pose a significant health and economic burden in developing countries. There are few effective and accessible treatments for these diseases, and the existing therapies suffer from problems such as parasite resistance and side effects. Structure-based drug design (SBDD) is one of the strategies that has been applied to discover new compounds targeting trypanosomatid-borne diseases. AREAS COVERED We review the current literature (mostly over the last 5 years, searched in PubMed database on Nov 11th 2021) on the application of structure-based drug design approaches to identify new anti-trypanosomatidic compounds that interfere with a validated target biochemical pathway, the trypanosomatid folate pathway. EXPERT OPINION The application of structure-based drug design approaches to perturb the trypanosomatid folate pathway has successfully provided many new inhibitors with good selectivity profiles, most of which are natural products or their derivatives or have scaffolds of known drugs. However, the inhibitory effect against the target protein(s) often does not translate to anti-parasitic activity. Further progress is hampered by our incomplete understanding of parasite biology and biochemistry, which is necessary to complement SBDD in a multiparameter optimization approach to discovering selective anti-parasitic drugs.
Collapse
Affiliation(s)
- Joanna Panecka-Hofman
- Division of Biophysics, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5a, 02-097 Warsaw, Poland
| | - Ina Poehner
- School of Pharmacy, University of Eastern Finland, Kuopio, Yliopistonranta 1C, PO Box 1627, FI-70211 Kuopio, Finland
| | - Rebecca Wade
- Center for Molecular Biology (ZMBH), Heidelberg University, Im Neuenheimer Feld 282, Heidelberg 69120, Germany.,Heidelberg Institute for Theoretical Studies (HITS), Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany.,DKFZ-ZMBH Alliance and Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| |
Collapse
|
71
|
Iorio A, Brochier-Armanet C, Mas C, Sterpone F, Madern D. Protein Conformational Space at the Edge of Allostery: Turning a Non-allosteric Malate Dehydrogenase into an "Allosterized" Enzyme using Evolution Guided Punctual Mutations. Mol Biol Evol 2022; 39:6691310. [PMID: 36056899 PMCID: PMC9486893 DOI: 10.1093/molbev/msac186] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.
Collapse
Affiliation(s)
- Antonio Iorio
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Céline Brochier-Armanet
- Univ Lyon, Université Lyon 1, CNRS, UMR5558, Laboratoire de Biométrie et Biologie Évolutive, 43 bd du 11 novembre 1918, F-69622, Villeurbanne, France
| | - Caroline Mas
- Univ. Grenoble Alpes, CEA, CNRS, IBS, 38000 Grenoble, France
| | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | | |
Collapse
|
72
|
Phenol sensing in nature is modulated via a conformational switch governed by dynamic allostery. J Biol Chem 2022; 298:102399. [PMID: 35988639 PMCID: PMC9556785 DOI: 10.1016/j.jbc.2022.102399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Revised: 08/11/2022] [Accepted: 08/12/2022] [Indexed: 11/21/2022] Open
Abstract
The NtrC family of proteins senses external stimuli and accordingly stimulates stress and virulence pathways via activation of associated σ54-dependent RNA polymerases. However, the structural determinants that mediate this activation are not well understood. Here, we establish using computational, structural, biochemical, and biophysical studies that MopR, an NtrC protein, harbors a dynamic bidirectional electrostatic network that connects the phenol pocket to two distal regions, namely the “G-hinge” and the “allosteric linker.” While the G-hinge influences the entry of phenol into the pocket, the allosteric linker passes the signal to the downstream ATPase domain. We show that phenol binding induces a rewiring of the electrostatic connections by eliciting dynamic allostery and demonstrates that perturbation of the core relay residues results in a complete loss of ATPase stimulation. Furthermore, we found a mutation of the G-hinge, ∼20 Å from the phenol pocket, promotes altered flexibility by shifting the pattern of conformational states accessed, leading to a protein with 7-fold enhanced phenol binding ability and enhanced transcriptional activation. Finally, we conducted a global analysis that illustrates that dynamic allostery-driven conserved community networks are universal and evolutionarily conserved across species. Taken together, these results provide insights into the mechanisms of dynamic allostery-mediated conformational changes in NtrC sensor proteins.
Collapse
|
73
|
Abstract
![]()
AlphaFold has burst into our lives. A powerful algorithm
that underscores
the strength of biological sequence data and artificial intelligence
(AI). AlphaFold has appended projects and research directions. The
database it has been creating promises an untold number of applications
with vast potential impacts that are still difficult to surmise. AI
approaches can revolutionize personalized treatments and usher in
better-informed clinical trials. They promise to make giant leaps
toward reshaping and revamping drug discovery strategies, selecting
and prioritizing combinations of drug targets. Here, we briefly overview
AI in structural biology, including in molecular dynamics simulations
and prediction of microbiota–human protein–protein interactions.
We highlight the advancements accomplished by the deep-learning-powered
AlphaFold in protein structure prediction and their powerful impact
on the life sciences. At the same time, AlphaFold does not resolve
the decades-long protein folding challenge, nor does it identify the
folding pathways. The models that AlphaFold provides do not capture
conformational mechanisms like frustration and allostery, which are
rooted in ensembles, and controlled by their dynamic distributions.
Allostery and signaling are properties of populations. AlphaFold also
does not generate ensembles of intrinsically disordered proteins and
regions, instead describing them by their low structural probabilities.
Since AlphaFold generates single ranked structures, rather than conformational
ensembles, it cannot elucidate the mechanisms of allosteric activating
driver hotspot mutations nor of allosteric drug resistance. However,
by capturing key features, deep learning techniques can use the single
predicted conformation as the basis for generating a diverse ensemble.
Collapse
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, United States.,Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Mingzhen Zhang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, United States
| | - Yonglan Liu
- Cancer Innovation Laboratory, National Cancer Institute, Frederick, Maryland 21702, United States
| | - Hyunbum Jang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, United States
| |
Collapse
|
74
|
Nam KM, Martinez-Corral R, Gunawardena J. The linear framework: using graph theory to reveal the algebra and thermodynamics of biomolecular systems. Interface Focus 2022; 12:20220013. [PMID: 35860006 PMCID: PMC9184966 DOI: 10.1098/rsfs.2022.0013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Accepted: 04/25/2022] [Indexed: 12/25/2022] Open
Abstract
The linear framework uses finite, directed graphs with labelled edges to model biomolecular systems. Graph vertices represent biochemical species or molecular states, edges represent reactions or transitions and labels represent rates. The graph yields a linear dynamics for molecular concentrations or state probabilities, with the graph Laplacian as the operator, and the labels encode the nonlinear interactions between system and environment. The labels can be specified by vertices of other graphs or by conservation laws or, when the environment consists of thermodynamic reservoirs, they may be constants. In the latter case, the graphs correspond to infinitesimal generators of Markov processes. The key advantage of the framework has been that steady states are determined as rational algebraic functions of the labels by the Matrix-Tree theorems of graph theory. When the system is at thermodynamic equilibrium, this prescription recovers equilibrium statistical mechanics but it continues to hold for non-equilibrium steady states. The framework goes beyond other graph-based approaches in treating the graph as a mathematical object, for which general theorems can be formulated that accommodate biomolecular complexity. It has been particularly effective at analysing enzyme-catalysed modification systems and input-output responses.
Collapse
Affiliation(s)
- Kee-Myoung Nam
- Department of Systems Biology, Harvard Medical School, Boston, MA 02115, USA
| | | | - Jeremy Gunawardena
- Department of Systems Biology, Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
75
|
Balasco N, Paladino A, Graziano G, D'Abramo M, Vitagliano L. Atomic-Level View of the Functional Transition in Vertebrate Hemoglobins: The Case of Antarctic Fish Hbs. J Chem Inf Model 2022; 62:3874-3884. [PMID: 35930673 PMCID: PMC9400108 DOI: 10.1021/acs.jcim.2c00727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Tetrameric hemoglobins (Hbs) are prototypal systems for studies aimed at unveiling basic structure-function relationships as well as investigating the molecular/structural basis of adaptation of living organisms to extreme conditions. However, a chronological analysis of decade-long studies conducted on Hbs is illuminating on the difficulties associated with the attempts of gaining functional insights from static structures. Here, we applied molecular dynamics (MD) simulations to explore the functional transition from the T to the R state of the hemoglobin of the Antarctic fish Trematomus bernacchii (HbTb). Our study clearly demonstrates the ability of the MD technique to accurately describe the transition of HbTb from the T to R-like states, as shown by a number of global and local structural indicators. A comparative analysis of the structural states that HbTb assumes in the simulations with those detected in previous MD analyses conducted on HbA (human Hb) highlights interesting analogies (similarity of the transition pathway) and differences (distinct population of intermediate states). In particular, the ability of HbTb to significantly populate intermediate states along the functional pathway explains the observed propensity of this protein to assume these structures in the crystalline state. It also explains some functional data reported on the protein that indicate the occurrence of other functional states in addition to the canonical R and T ones. These findings are in line with the emerging idea that the classical two-state view underlying tetrameric Hb functionality is probably an oversimplification and that other structural states play important roles in these proteins. The ability of MD simulations to accurately describe the functional pathway in tetrameric Hbs suggests that this approach may be effectively applied to unravel the molecular and structural basis of Hbs exhibiting peculiar functional properties as a consequence of the environmental adaptation of the host organism.
Collapse
Affiliation(s)
- Nicole Balasco
- Institute of Molecular Biology and Pathology, CNR c/o Dep. Chemistry, University of Rome, Sapienza, P.le A. Moro 5, 00185 Rome, Italy
| | - Antonella Paladino
- Institute of Biostructures and Bioimaging, CNR, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Giuseppe Graziano
- Department of Science and Technology, University of Sannio, via Francesco de Sanctis snc, Benevento 82100, Italy
| | - Marco D'Abramo
- Department of Chemistry, University of Rome Sapienza, P.le A.Moro 5, 00185 Rome, Italy
| | - Luigi Vitagliano
- Institute of Biostructures and Bioimaging, CNR, Via Pietro Castellino 111, 80131 Naples, Italy
| |
Collapse
|
76
|
Chen S, Giardina F, Choi GPT, Mahadevan L. Modular representation and control of floppy networks. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Geometric graph models of systems as diverse as proteins, DNA assemblies, architected materials and robot swarms are useful abstract representations of these objects that also unify ways to study their properties and control them in space and time. While much work has been done in the context of characterizing the behaviour of these networks close to critical points associated with bond and rigidity percolation, isostaticity, etc., much less is known about floppy, underconstrained networks that are far more common in nature and technology. Here, we combine geometric rigidity and algebraic sparsity to provide a framework for identifying the zero energy floppy modes via a representation that illuminates the underlying hierarchy and modularity of the network and thence the control of its nestedness and locality. Our framework allows us to demonstrate a range of applications of this approach that include robotic reaching tasks with motion primitives, and predicting the linear and nonlinear response of elastic networks based solely on infinitesimal rigidity and sparsity, which we test using physical experiments. Our approach is thus likely to be of use broadly in dissecting the geometrical properties of floppy networks using algebraic sparsity to optimize their function and performance.
Collapse
Affiliation(s)
- Siheng Chen
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
| | - Fabio Giardina
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
| | - Gary P. T. Choi
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge MA 02139, USA
| | - L. Mahadevan
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
- Department of Physics, Harvard University, Cambridge MA 02138, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge MA 02138, USA
| |
Collapse
|
77
|
Cai J, Zhao L, Li Y, He C, Wang C, Duan C. Binding of Dual-Function Hybridized Metal -Organic Capsules to Enzymes for Cascade Catalysis. JACS AU 2022; 2:1736-1746. [PMID: 35911460 PMCID: PMC9327082 DOI: 10.1021/jacsau.2c00322] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The combination of chemo- and biocatalysis for multistep syntheses provides attractive advantages in terms of evolvability, promiscuity, and sustainability striving for desirable catalytic performance. Through the encapsulation of flavin analogues by both NADH and heme mimics codecorated heteroleptic metal-organic capsules, herein, we report a progressive host-guest strategy to imitate cytochrome P450s catalysis for cascade oxidative coupling catalysis. Besides the construction of stable dual-function metal-organic capsules and the modification of cofactor-decorated capsules at the domain of enzymes, this supramolecular strategy involves multistage directional electron flow, affording reactive ferric peroxide species for inducing oxygenation. Under light irradiation, the metal-organic capsule selectively converts stilbene to oxidative coupling products (including 2-oxo-1,2-diphenylethyl formate, 2-alkoxy-1,2-diphenylethanone) in tandem with enzymatic reactions respectively, at the domain of natural enzymes. The ingenious combination of capsules and enzymes with the in situ-regenerated capsule-loaded NADH cofactor promises non-native coupling reactions by forming regional cooperation and division. This abiotic-biotic conjugated host-guest strategy is conducive to the de novo creation of multifunctional components approaching active enzymatic sites for reinforced matter and energy transporting, demonstrating a key role of multicomponent supramolecular catalysts for one-pot integrated catalytic conversions.
Collapse
Affiliation(s)
- Junkai Cai
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
- State
Key Laboratory of Coordination Chemistry, Nanjing University, Nanjing 210023, People’s Republic
of China
| | - Liang Zhao
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
| | - Yanan Li
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
| | - Cheng He
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
| | - Chong Wang
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
| | - Chunying Duan
- State
Key Laboratory of Fine Chemicals, Zhang Dayu School of Chemistry, Dalian University of Technology, Dalian 116024, People’s Republic of China
- State
Key Laboratory of Coordination Chemistry, Nanjing University, Nanjing 210023, People’s Republic
of China
| |
Collapse
|
78
|
Mishra SH, Kancherla AK, Marincin KA, Bouvignies G, Nerli S, Sgourakis N, Dowling DP, Frueh DP. Global protein dynamics as communication sensors in peptide synthetase domains. SCIENCE ADVANCES 2022; 8:eabn6549. [PMID: 35857508 PMCID: PMC9286511 DOI: 10.1126/sciadv.abn6549] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 05/28/2022] [Indexed: 05/04/2023]
Abstract
Biological activity is governed by the timely redistribution of molecular interactions, and static structural snapshots often appear insufficient to provide the molecular determinants that choreograph communication. This conundrum applies to multidomain enzymatic systems called nonribosomal peptide synthetases (NRPSs), which assemble simple substrates into complex metabolites, where a dynamic domain organization challenges rational design to produce new pharmaceuticals. Using a nuclear magnetic resonance (NMR) atomic-level readout of biochemical transformations, we demonstrate that global structural fluctuations help promote substrate-dependent communication and allosteric responses, and impeding these global dynamics by a point-site mutation hampers allostery and molecular recognition. Our results establish global structural dynamics as sensors of molecular events that can remodel domain interactions, and they provide new perspectives on mechanisms of allostery, protein communication, and NRPS synthesis.
Collapse
Affiliation(s)
- Subrata H. Mishra
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Aswani K. Kancherla
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Kenneth A. Marincin
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Guillaume Bouvignies
- Laboratoire des Biomolécules (LBM), Département de Chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, Paris, France
| | - Santrupti Nerli
- Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA, USA
| | - Nikolaos Sgourakis
- Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Daniel P. Dowling
- Department of Chemistry, University of Massachusetts Boston, Boston, MA, USA
| | - Dominique P. Frueh
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| |
Collapse
|
79
|
Wu N, Yaliraki SN, Barahona M. Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity. J Mol Biol 2022; 434:167749. [PMID: 35841931 DOI: 10.1016/j.jmb.2022.167749] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Revised: 07/09/2022] [Accepted: 07/11/2022] [Indexed: 11/28/2022]
Abstract
Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.
Collapse
Affiliation(s)
- Nan Wu
- Department of Chemistry Imperial College London, United Kingdom
| | | | | |
Collapse
|
80
|
Dubanevics I, McLeish TCB. Optimising Elastic Network Models for Protein Dynamics and Allostery: Spatial and Modal Cut-offs and Backbone Stiffness. J Mol Biol 2022; 434:167696. [PMID: 35810792 DOI: 10.1016/j.jmb.2022.167696] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 06/16/2022] [Accepted: 06/19/2022] [Indexed: 01/15/2023]
Abstract
The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.
Collapse
|
81
|
Diez G, Nagel D, Stock G. Correlation-Based Feature Selection to Identify Functional Dynamics in Proteins. J Chem Theory Comput 2022; 18:5079-5088. [PMID: 35793551 DOI: 10.1021/acs.jctc.2c00337] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To interpret molecular dynamics simulations of biomolecular systems, systematic dimensionality reduction methods are commonly employed. Among others, this includes principal component analysis (PCA) and time-lagged independent component analysis (TICA), which aim to maximize the variance and the time scale of the first components, respectively. A crucial first step of such an analysis is the identification of suitable and relevant input coordinates (the so-called features), such as backbone dihedral angles and interresidue distances. As typically only a small subset of those coordinates is involved in a specific biomolecular process, it is important to discard the remaining uncorrelated motions or weakly correlated noise coordinates. This is because they may exhibit large amplitudes or long time scales and therefore will be erroneously considered important by PCA and TICA, respectively. To discriminate collective motions underlying functional dynamics from uncorrelated motions, the correlation matrix of the input coordinates is block-diagonalized by a clustering method. This strategy avoids possible bias due to presumed functional observables and conformational states or variation principles that maximize variance or time scales. Considering several linear and nonlinear correlation measures and various clustering algorithms, it is shown that the combination of linear correlation and the Leiden community detection algorithm yields excellent results for all considered model systems. These include the functional motion of T4 lysozyme to demonstrate the successful identification of collective motion, as well as the folding of the villin headpiece to highlight the physical interpretation of the correlated motions in terms of a functional mechanism.
Collapse
Affiliation(s)
- Georg Diez
- Biomolecular Dynamics, Institute of Physics, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Daniel Nagel
- Biomolecular Dynamics, Institute of Physics, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| |
Collapse
|
82
|
MacKenzie DWS, Schaefer A, Steckner J, Leo CA, Naser D, Artikis E, Broom A, Ko T, Shah P, Ney MQ, Tran E, Smith MTJ, Fuglestad B, Wand AJ, Brooks CL, Meiering EM. A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein. Proc Natl Acad Sci U S A 2022; 119:e2119686119. [PMID: 35737838 PMCID: PMC9245636 DOI: 10.1073/pnas.2119686119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Accepted: 03/29/2022] [Indexed: 12/24/2022] Open
Abstract
Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.
Collapse
Affiliation(s)
| | - Anna Schaefer
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Julia Steckner
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Christopher A. Leo
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Dalia Naser
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Efrosini Artikis
- Department of Chemistry and Biophysics, University of Michigan, Ann Arbor, MI 48109
| | - Aron Broom
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Travis Ko
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Purnank Shah
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Mikaela Q. Ney
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Elisa Tran
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Martin T. J. Smith
- Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Brian Fuglestad
- Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104
| | - A. Joshua Wand
- Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104
| | - Charles L. Brooks
- Department of Chemistry and Biophysics, University of Michigan, Ann Arbor, MI 48109
| | | |
Collapse
|
83
|
Yuan Y, Deng J, Cui Q. Molecular Dynamics Simulations Establish the Molecular Basis for the Broad Allostery Hotspot Distributions in the Tetracycline Repressor. J Am Chem Soc 2022; 144:10870-10887. [PMID: 35675441 DOI: 10.1021/jacs.2c03275] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It is imperative to identify the network of residues essential to the allosteric coupling for the purpose of rationally engineering allostery in proteins. Deep mutational scanning analysis has emerged as a function-centric approach for identifying such allostery hotspots in a comprehensive and unbiased fashion, leading to observations that challenge our understanding of allostery at the molecular level. Specifically, a recent deep mutational scanning study of the tetracycline repressor (TetR) revealed an unexpectedly broad distribution of allostery hotspots throughout the protein structure. Using extensive molecular dynamics simulations (up to 50 μs) and free energy computations, we establish the molecular and energetic basis for the strong anticooperativity between the ligand and DNA binding sites. The computed free energy landscapes in different ligation states illustrate that allostery in TetR is well described by a conformational selection model, in which the apo state samples a broad set of conformations, and specific ones are selectively stabilized by either ligand or DNA binding. By examining a range of structural and dynamic properties of residues at both local and global scales, we observe that various analyses capture different subsets of experimentally identified hotspots, suggesting that these residues modulate allostery in distinct ways. These results motivate the development of a thermodynamic model that qualitatively explains the broad distribution of hotspot residues and their distinct features in molecular dynamics simulations. The multifaceted strategy that we establish here for hotspot evaluations and our insights into their mechanistic contributions are useful for modulating protein allostery in mechanistic and engineering studies.
Collapse
Affiliation(s)
- Yuchen Yuan
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
84
|
Post M, Lickert B, Diez G, Wolf S, Stock G. Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network. J Mol Biol 2022; 434:167679. [PMID: 35690098 DOI: 10.1016/j.jmb.2022.167679] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Revised: 06/02/2022] [Accepted: 06/06/2022] [Indexed: 12/13/2022]
Abstract
Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.
Collapse
Affiliation(s)
- Matthias Post
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany. https://twitter.com/@_posti
| | - Benjamin Lickert
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany. https://twitter.com/@BenjaminLickert
| | - Georg Diez
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany. https://twitter.com/@gegadiez
| | - 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
|
85
|
Huang J, Chu X, Luo Y, Wang Y, Zhang Y, Zhang Y, Li H. Insights into Phosphorylation-Induced Protein Allostery and Conformational Dynamics of Glycogen Phosphorylase via Integrative Structural Mass Spectrometry and In Silico Modeling. ACS Chem Biol 2022; 17:1951-1962. [PMID: 35675581 DOI: 10.1021/acschembio.2c00393] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Allosteric regulation plays a fundamental role in innumerable biological processes. Understanding its dynamic mechanism and impact at the molecular level is of great importance in disease diagnosis and drug discovery. Glycogen phosphorylase (GP) is a phosphoprotein responding to allosteric regulation and has significant biological importance to glycogen metabolism. Although the atomic structures of GP have been previously solved, the conformational dynamics of GP related to allostery regulation remain largely elusive due to its macromolecular size (∼196 kDa). Here, we integrated native top-down mass spectrometry (nTD-MS), hydrogen-deuterium exchange MS (HDX-MS), protection factor (PF) analysis, molecular dynamics (MD) simulations, and allostery signaling analysis to examine the structural basis and dynamics for the allosteric regulation of GP by phosphorylation. nTD-MS reveals differences in structural stability as well as oligomeric state between the unphosphorylated (GPb) and phosphorylated (GPa) forms. HDX-MS, PF analysis, and MD simulations further pinpoint the structural differences between GPb and GPa involving the binding interfaces (the N-terminal and tower-tower helices), catalytic site, and PLP-binding region. More importantly, it also allowed us to complete the missing link of the long-range communication process from the N-terminal tail to the catalytic site caused by phosphorylation. This integrative MS and in silico-based platform is highly complementary to biophysical approaches and yields valuable insights into protein structures and dynamic regulation.
Collapse
Affiliation(s)
- Jing Huang
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou Higher Education Mega Center, No. 132 Wai Huan Dong Lu, Guangzhou 510006, China
| | - Xiakun Chu
- Advanced Materials Thrust, The Hong Kong University of Science and Technology (Guangzhou), Nansha, Guangzhou, Guangdong 511400, China
| | - Yuxiang Luo
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou Higher Education Mega Center, No. 132 Wai Huan Dong Lu, Guangzhou 510006, China
| | - Yong Wang
- The Provincial International Science and Technology Cooperation Base on Engineering Biology, International Campus of Zhejiang University, College of Life Sciences, Shanghai Institute for Advanced Study, Institute of Quantitative Biology, Zhejiang University, Haining 314400, China
| | - Ying Zhang
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou Higher Education Mega Center, No. 132 Wai Huan Dong Lu, Guangzhou 510006, China
| | - Yu Zhang
- State Key Laboratory of Crop Stress Adaptation and Improvement, Key Laboratory of Plant Stress Biology, School of Life Sciences, Henan University, Kaifeng 475004, China
| | - Huilin Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou Higher Education Mega Center, No. 132 Wai Huan Dong Lu, Guangzhou 510006, China.,Guangdong Key Laboratory of Chiral Molecule and Drug Discovery, School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, Guangdong 510006, China
| |
Collapse
|
86
|
Bern D, Tobi D. The effect of dimerization and ligand binding on the dynamics of Kaposi's sarcoma‐associated herpesvirus protease. Proteins 2022; 90:1267-1277. [PMID: 35084062 PMCID: PMC9305915 DOI: 10.1002/prot.26307] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2021] [Revised: 11/18/2021] [Accepted: 01/24/2022] [Indexed: 11/10/2022]
Abstract
The Kaposi's sarcoma‐associated herpesvirus protease is essential for virus maturation. This protease functions under allosteric regulation that establishes its enzymatic activity upon dimerization. It exists in equilibrium between an inactive monomeric state and an active, weakly associating, dimeric state that is stabilized upon ligand binding. The dynamics of the protease dimer and its monomer were studied using the Gaussian network model and the anisotropic network model , and its role in mediating the allosteric regulation is demonstrated. We show that the dimer is composed of five dynamical domains. The central domain is formed upon dimerization and composed of helix five of each monomer, in addition to proximal and distal domains of each monomer. Dimerization reduces the mobility of the central domains and increases the mobility of the distal domains, in particular the binding site within them. The three slowest ANM modes of the dimer assist the protease in ligand binding, motion of the conserved Arg142 and Arg143 toward the oxyanion, and reducing the activation barrier for the tetrahedral transition state by stretching the bond that is cleaved by the protease. In addition, we show that ligand binding reduces the motion of helices α1 and α5 at the interface and explain how ligand binding can stabilize the dimer.
Collapse
Affiliation(s)
- David Bern
- Department of Molecular Biology Ariel University Ariel Israel
| | - Dror Tobi
- Department of Molecular Biology Ariel University Ariel Israel
- Department of Computer Sciences Ariel University Ariel Israel
| |
Collapse
|
87
|
Guan X, Tan C, Li W, Wang W, Thirumalai D. Role of water-bridged interactions in metal ion coupled protein allostery. PLoS Comput Biol 2022; 18:e1010195. [PMID: 35653400 PMCID: PMC9197054 DOI: 10.1371/journal.pcbi.1010195] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 06/14/2022] [Accepted: 05/10/2022] [Indexed: 11/18/2022] Open
Abstract
Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca2+ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca2+ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca2+ binding and protein allostery. The bridging water molecules between Ca2+ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca2+ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca2+ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery. Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca2+ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca2+ to the binding site involves an intermediate in which the water molecule bridges the Ca2+ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca2+ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca2+ ions and is needed for efficient allosteric communication.
Collapse
Affiliation(s)
- Xingyue Guan
- Department of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing, China
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang, China
| | - Cheng Tan
- Department of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing, China
| | - Wenfei Li
- Department of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing, China
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang, China
- * E-mail: (WL); (WW); (DT)
| | - Wei Wang
- Department of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing, China
- * E-mail: (WL); (WW); (DT)
| | - D. Thirumalai
- Department of Chemistry, University of Texas, Texas, United States of America
- * E-mail: (WL); (WW); (DT)
| |
Collapse
|
88
|
Casadó-Anguera V, Casadó V. Unmasking allosteric binding sites: Novel targets for GPCR drug discovery. Expert Opin Drug Discov 2022; 17:897-923. [PMID: 35649692 DOI: 10.1080/17460441.2022.2085684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
INTRODUCTION Unexpected non-apparent and hidden allosteric binding sites are non-classical and non-apparent allosteric centers in 3-D X-ray protein structures until orthosteric or allosteric ligands bind to them. The orthosteric center of one protomer that modulates binding centers of the other protomers within an oligomer is also an unexpected allosteric site. Furthermore, another partner protein can also produce these effects, acting as an unexpected allosteric modulator. AREAS COVERED This review summarizes both classical and non-classical allosterism. The authors focus on G protein-coupled receptor (GPCR) oligomers as a paradigm of allosteric molecules. Moreover, they show several examples of unexpected allosteric sites such as hidden allosteric sites in a protomer that appear after the interaction with other molecules and the allosterism exerted between orthosteric sites within GPCR oligomer, emphasizing on the allosteric modulations that can occur between binding sites. EXPERT OPINION The study of these new non-classical allosteric sites will expand the diversity of allosteric control on the function of orthosteric sites within proteins, whether GPCRs or other receptors, enzymes or transporters. Moreover, the design of new drugs targeting these hidden allosteric sites or already known orthosteric sites acting as allosteric sites in protein homo- or hetero-oligomers will increase the therapeutic potential of allosterism.
Collapse
Affiliation(s)
- Verònica Casadó-Anguera
- Laboratory of Molecular Neuropharmacology, Department of Biochemistry and Molecular Biomedicine, Faculty of Biology, Universitat de Barcelona, and Institute of Biomedicine of the Universitat de Barcelona, Barcelona, Spain.,Laboratory of Neuropharmacology-Neurophar, Department of Experimental and Health Sciences, Universitat Pompeu Fabra, Barcelona, Spain
| | - Vicent Casadó
- Laboratory of Molecular Neuropharmacology, Department of Biochemistry and Molecular Biomedicine, Faculty of Biology, Universitat de Barcelona, and Institute of Biomedicine of the Universitat de Barcelona, Barcelona, Spain
| |
Collapse
|
89
|
Mingione VR, Foda Z, Paung Y, Philipose H, Rangwala AM, Shan Y, Seeliger MA. Validation of an allosteric binding site of Src kinase identified by unbiased ligand binding simulations. J Mol Biol 2022; 434:167628. [PMID: 35595169 DOI: 10.1016/j.jmb.2022.167628] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 04/13/2022] [Accepted: 05/04/2022] [Indexed: 10/18/2022]
Abstract
Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a v-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.
Collapse
Affiliation(s)
- Victoria R Mingione
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA
| | - Zachariah Foda
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA
| | - YiTing Paung
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA
| | - Hannah Philipose
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA
| | - Aziz M Rangwala
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA
| | - Yibing Shan
- Antidote Health Foundation for Cure of Cancer, Cambridge, MA 02139, USA.
| | - Markus A Seeliger
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794, USA.
| |
Collapse
|
90
|
Domino effect in allosteric signaling of peptide binding. J Mol Biol 2022; 434:167661. [DOI: 10.1016/j.jmb.2022.167661] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 05/24/2022] [Accepted: 05/24/2022] [Indexed: 11/22/2022]
|
91
|
Subsets of Slow Dynamic Modes Reveal Global Information Sources as Allosteric Sites. J Mol Biol 2022; 434:167644. [DOI: 10.1016/j.jmb.2022.167644] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 05/17/2022] [Accepted: 05/17/2022] [Indexed: 02/06/2023]
|
92
|
Kremer DM, Lyssiotis CA. Targeting allosteric regulation of cancer metabolism. Nat Chem Biol 2022; 18:441-450. [PMID: 35484254 DOI: 10.1038/s41589-022-00997-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Accepted: 02/14/2022] [Indexed: 12/13/2022]
Abstract
Metabolic reprogramming is observed across all cancer types. Indeed, the success of many classic chemotherapies stems from their targeting of cancer metabolism. Contemporary research in this area has refined our understanding of tumor-specific metabolic mechanisms and has revealed strategies for exploiting these vulnerabilities selectively. Based on this growing understanding, new small-molecule tools and drugs have been developed to study and target tumor metabolism. Here, we highlight allosteric modulation of metabolic enzymes as an attractive mechanism of action for small molecules that target metabolic enzymes. We then discuss the mechanistic insights garnered from their application in cancer studies and highlight the achievements of this approach in targeting cancer metabolism. Finally, we discuss technological advances in drug discovery for allosteric modulators of enzyme activity.
Collapse
Affiliation(s)
- Daniel M Kremer
- Department of Molecular & Integrative Physiology, University of Michigan, Ann Arbor, MI, USA.,Graduate Program in Chemical Biology, University of Michigan, Ann Arbor, MI, USA.,Department of Chemistry, the Scripps Research Institute, La Jolla, CA, USA
| | - Costas A Lyssiotis
- Department of Molecular & Integrative Physiology, University of Michigan, Ann Arbor, MI, USA. .,Department of Internal Medicine, Division of Gastroenterology and Hepatology, University of Michigan, Ann Arbor, MI, USA. .,Rogel Cancer Center, University of Michigan, Ann Arbor, MI, USA.
| |
Collapse
|
93
|
Abstract
A substantial portion of molecules in an organism are involved in regulation of a wide spectrum of biological processes. Several models have been presented for various forms of biological regulation, including gene expression regulation and physiological regulation; however, a generic model is missing. Recently a new unifying theory in biology, poikilosis, was presented. Poikilosis indicates that all systems display intrinsic heterogeneity, which is a normal state. The concept of poikilosis allowed development of a model for biological regulation applicable to all types of regulated systems. The perturbation-lagom-TATAR countermeasures-regulator (PLTR) model combines the effects of perturbation and lagom (allowed and sufficient extent of heterogeneity) in a system with tolerance, avoidance, repair, attenuation and resistance (TARAR) countermeasures, and possible regulators. There are three modes of regulation, two of which are lagom-related. In the first scenario, lagom is maintained, both intrinsic (passive) and active TARAR countermeasures can be involved. In the second mode, there is a shift from one lagom to another. In the third mode, reguland regulation, the regulated entity is the target of a regulatory shift, which is often irreversible or requires action of another regulator to return to original state. After the shift, the system enters to lagom maintenance mode, but at new lagom extent. The model is described and elaborated with examples and applications, including medicine and systems biology. Consequences of non-lagom extent of heterogeneity are introduced, along with a novel idea for therapy by reconstituting biological processes to lagom extent, even when the primary effect cannot be treated.
Collapse
Affiliation(s)
- Mauno Vihinen
- Department of Experimental Medical Science, Lund University, BMC B13, Lund, SE-221 84, Sweden
| |
Collapse
|
94
|
Abstract
A substantial portion of molecules in an organism are involved in regulation of a wide spectrum of biological processes. Several models have been presented for various forms of biological regulation, including gene expression regulation and physiological regulation; however, a generic model is missing. Recently a new unifying theory in biology, poikilosis, was presented. Poikilosis indicates that all systems display intrinsic heterogeneity. The concept of poikilosis allowed development of a model for biological regulation applicable to all types of regulated systems. The perturbation-lagom-TATAR countermeasures-regulator (PLTR) model combines the effects of perturbation and lagom (allowed and sufficient extent of heterogeneity) in a system with tolerance, avoidance, repair, attenuation and resistance (TARAR) countermeasures, and possible regulators. There are three modes of regulation, two of which are lagom-related. In the first scenario, lagom is maintained, both intrinsic (passive) and active TARAR countermeasures can be involved. In the second mode, there is a shift from one lagom to another. In the third mode, reguland regulation, the regulated entity is the target of a regulatory shift, which is often irreversible or requires action of another regulator to return to original state. After the shift, the system enters to lagom maintenance mode, but at new lagom extent. The model is described and elaborated with examples and applications, including medicine and systems biology. Consequences of non-lagom extent of heterogeneity are introduced, along with a novel idea for therapy by reconstituting biological processes to lagom extent, even when the primary effect cannot be treated.
Collapse
Affiliation(s)
- Mauno Vihinen
- Department of Experimental Medical Science, Lund University, BMC B13, Lund, SE-221 84, Sweden
| |
Collapse
|
95
|
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
|
96
|
Verkhivker GM, Agajanian S, Kassab R, Krishnan K. Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants. J Chem Inf Model 2022; 62:1956-1978. [PMID: 35377633 DOI: 10.1021/acs.jcim.2c00124] [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/22/2022]
Abstract
The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Steve Agajanian
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Ryan Kassab
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Keerthi Krishnan
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| |
Collapse
|
97
|
Tastan Bishop Ö, Mutemi Musyoka T, Barozi V. Allostery and missense mutations as intermittently linked promising aspects of modern computational drug discovery. J Mol Biol 2022; 434:167610. [DOI: 10.1016/j.jmb.2022.167610] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2022] [Revised: 04/21/2022] [Accepted: 04/22/2022] [Indexed: 12/15/2022]
|
98
|
Nussinov R, Tsai CJ, Jang H. Allostery, and how to define and measure signal transduction. Biophys Chem 2022; 283:106766. [PMID: 35121384 PMCID: PMC8898294 DOI: 10.1016/j.bpc.2022.106766] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Revised: 01/21/2022] [Accepted: 01/24/2022] [Indexed: 12/15/2022]
Abstract
Here we ask: What is productive signaling? How to define it, how to measure it, and most of all, what are the parameters that determine it? Further, what determines the strength of signaling from an upstream to a downstream node in a specific cell? These questions have either not been considered or not entirely resolved. The requirements for the signal to propagate downstream to activate (repress) transcription have not been considered either. Yet, the questions are pivotal to clarify, especially in diseases such as cancer where determination of signal propagation can point to cell proliferation and to emerging drug resistance, and to neurodevelopmental disorders, such as RASopathy, autism, attention-deficit/hyperactivity disorder (ADHD), and cerebral palsy. Here we propose a framework for signal transduction from an upstream to a downstream node addressing these questions. Defining cellular processes, experimentally measuring them, and devising powerful computational AI-powered algorithms that exploit the measurements, are essential for quantitative science.
Collapse
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA; Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel.
| | - Chung-Jung Tsai
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| | - Hyunbum Jang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| |
Collapse
|
99
|
Karlsson E, Sorgenfrei FA, Andersson E, Dogan J, Jemth P, Chi CN. The dynamic properties of a nuclear coactivator binding domain are evolutionarily conserved. Commun Biol 2022; 5:286. [PMID: 35354917 PMCID: PMC8967867 DOI: 10.1038/s42003-022-03217-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 03/02/2022] [Indexed: 12/21/2022] Open
Abstract
Evolution of proteins is constrained by their structure and function. While there is a consensus that the plasticity of intrinsically disordered proteins relaxes the structural constraints on evolution there is a paucity of data on the molecular details of these processes. The Nuclear Coactivator Binding Domain (NCBD) from CREB-binding protein is a protein interaction domain, which contains a hydrophobic core but is not behaving as a typical globular domain, and has been described as 'molten-globule like'. The highly dynamic properties of NCBD makes it an interesting model system for evolutionary structure-function investigation of intrinsically disordered proteins. We have here compared the structure and biophysical properties of an ancient version of NCBD present in a bilaterian animal ancestor living around 600 million years ago with extant human NCBD. Using a combination of NMR spectroscopy, circular dichroism and kinetics we show that although NCBD has increased its thermodynamic stability, it has retained its dynamic biophysical properties in the ligand-free state in the evolutionary lineage leading from the last common bilaterian ancestor to humans. Our findings suggest that the dynamic properties of NCBD have been maintained by purifying selection and thus are important for its function, which includes mediating several distinct protein-protein interactions.
Collapse
Affiliation(s)
- Elin Karlsson
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden
| | - Frieda A Sorgenfrei
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden.,acib GmbH, Krenngasse 37, 8010 Graz c/o University of Graz, Institute of Chemistry, NAWI Graz, BioTechMed Graz, Heinrichstrasse 28, 8010, Graz, Austria
| | - Eva Andersson
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden
| | - Jakob Dogan
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden
| | - Per Jemth
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden.
| | - Celestine N Chi
- Department of Medical Biochemistry and Microbiology, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden. .,Department of Pharmaceutical Biosciences, Uppsala University, BMC Box 582, SE-75123, Uppsala, Sweden.
| |
Collapse
|
100
|
Yuce M, Sarica Z, Ates B, Kurkcuoglu O. Exploring species-specific inhibitors with multiple target sites on S. aureus pyruvate kinase using a computational workflow. J Biomol Struct Dyn 2022; 41:3496-3510. [PMID: 35302925 DOI: 10.1080/07391102.2022.2051743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Experimental evidence indicated that bacterial pyruvate kinase of glycolysis can be evaluated as an alternative target to eliminate infections, while antibiotic resistance poses a global threat. Here, we use a computational workflow to reveal and investigate the potential allosteric sites of methicillin-resistant S. aureus PK, which can help in designing species-specific drugs to inhibit activity of this organism. Residue interaction networks point to a known allosteric site at the small C-C interface, a potential allosteric site near the small interface (site #1), and a second potential allosteric site at the large interface (site #2). 2 µs-long molecular dynamics (MD) simulations with AMBER16 generate different conformations of one narrow target site. Known and potential allosteric sites on the selected conformers are investigated using ensemble docking with AutoDock Vina and a library of 2447 FDA-approved drugs. We determine 18 hits, comprising ergot-alkaloids, anti-cancer-agents, antivirals, analgesics, cardiac glycosides, all with a high docking z-score for three sites. 5 selected compounds with high, average and low z-scores are subjected to 50 ns-long MD simulations for MM-GBSA calculations. ΔGbind values up to -49.3 kcal/mol at the C-C interface, up to -32.7 kcal/mol at site #1, and up to -53.3 kcal/mol at site #2 support the docking calculations. We investigate mitapivat and TT-232 as reference compounds under clinical trial, targeting human PK isomers. We suggest 18 FDA-approved hits from the docking calculations and TT-232 as potential inhibitors with multiple target sites on S. aureus PK. This study also proposes pharmacophores models for de novo drug design.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Merve Yuce
- Department of Chemical Engineering, Istanbul Technical University, Istanbul, Turkey
| | - Zehra Sarica
- Computational Science and Engineering Division, Informatics Institute, Istanbul Technical University, Istanbul, Turkey
| | - Beril Ates
- Department of Chemical Engineering, Istanbul Technical University, Istanbul, Turkey
| | - Ozge Kurkcuoglu
- Department of Chemical Engineering, Istanbul Technical University, Istanbul, Turkey
| |
Collapse
|