1
|
Pražnikar J. Using graphlet degree vectors to predict atomic displacement parameters in protein structures. Acta Crystallogr D Struct Biol 2023; 79:1109-1119. [PMID: 37987168 PMCID: PMC10833351 DOI: 10.1107/s2059798323009142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 10/17/2023] [Indexed: 11/22/2023] Open
Abstract
In structural biology, atomic displacement parameters, commonly used in the form of B values, describe uncertainties in atomic positions. Their distribution over the structure can provide hints on local structural reliability and mobility. A spatial macromolecular model can be represented by a graph whose nodes are atoms and whose edges correspond to all interatomic contacts within a certain distance. Small connected subgraphs, called graphlets, provide information about the wiring of a particular atom. The multiple linear regression approach based on this information aims to predict a distribution of values of isotropic atomic displacement parameters (B values) within a protein structure, given the atomic coordinates and molecular packing. By modeling the dynamic component of atomic uncertainties, this method allows the B values obtained from experimental crystallographic or cryo-electron microscopy studies to be reproduced relatively well.
Collapse
Affiliation(s)
- Jure Pražnikar
- Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, Koper, Slovenia
- Department of Biochemistry, Molecular and Structural Biology, Institute Jožef Stefan, Jamova 39, Ljubljana, Slovenia
| |
Collapse
|
2
|
La Sala G, Pfleger C, Käck H, Wissler L, Nevin P, Böhm K, Janet JP, Schimpl M, Stubbs CJ, De Vivo M, Tyrchan C, Hogner A, Gohlke H, Frolov AI. Combining structural and coevolution information to unveil allosteric sites. Chem Sci 2023; 14:7057-7067. [PMID: 37389247 PMCID: PMC10306073 DOI: 10.1039/d2sc06272k] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Accepted: 06/02/2023] [Indexed: 07/01/2023] Open
Abstract
Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.
Collapse
Affiliation(s)
- Giuseppina La Sala
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
| | - Helena Käck
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Lisa Wissler
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Philip Nevin
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Kerstin Böhm
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Jon Paul Janet
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Marianne Schimpl
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Christopher J Stubbs
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Design, Istituto Italiano di Tecnologia Via Morego 30 16163 Genoa Italy
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory & Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Anders Hogner
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Institute of Bio- and Geosciences (IBG-4: Bioinformatics) Forschungszentrum Jülich GmbH 52425 Jülich Germany
| | - Andrey I Frolov
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| |
Collapse
|
3
|
Muhammad N, Hussain SI, Rehman ZU, Khan SA, Jan S, Khan N, Muzammal M, Abbasi SW, Kakar N, Rehman ZU, Khan MA, Mirza MU, Muhammad N, Khan S, Wasif N. Autosomal recessive variants c.953A>C and c.97-1G>C in NSUN2 causing intellectual disability: a molecular dynamics simulation study of loss-of-function mechanisms. Front Neurol 2023; 14:1168307. [PMID: 37305761 PMCID: PMC10249782 DOI: 10.3389/fneur.2023.1168307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 04/28/2023] [Indexed: 06/13/2023] Open
Abstract
Introduction Intellectual disability (ID) is a clinically and genetically heterogeneous disorder. It drastically affects the learning capabilities of patients and eventually reduces their IQ level below 70. Methods The current genetic study ascertained two consanguineous Pakistani families suffering from autosomal recessive intellectual developmental disorder-5 (MRT5). We have used exome sequencing followed by Sanger sequencing to identify the disease-causing variants. Results and discussion Genetic analysis using whole exome sequencing in these families identified two novel mutations in the NSUN2 (NM_017755.5). Family-A segregated a novel missense variant c.953A>C; p.Tyr318Ser in exon-9 of the NSUN2. The variant substituted an amino acid Tyr318, highly conserved among different animal species and located in the functional domain of NSUN2 known as "SAM-dependent methyltransferase RsmB/NOP2-type". Whereas in family B, we identified a novel splice site variant c.97-1G>C that affects the splice acceptor site of NSUN2. The identified splice variant (c.97-1G>C) was predicted to result in the skipping of exon-2, which would lead to a frameshift followed by a premature stop codon (p. His86Profs*16). Furthermore, it could result in the termination of translation and synthesis of dysfunctional protein, most likely leading to nonsense-mediated decay. The dynamic consequences of NSUN2 missense variant was further explored together with wildtype through molecular dynamic simulations, which uncovered the disruption of NSUN2 function due to a gain in structural flexibility. The present molecular genetic study further extends the mutational spectrum of NSUN2 to be involved in ID and its genetic heterogeneity in the Pakistani population.
Collapse
Affiliation(s)
- Nazif Muhammad
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Syeda Iqra Hussain
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Zia Ur Rehman
- Department of General Medicine, Northwest General Hospital & Research Center, Peshawar, Khyber Pakhtunkhwa, Pakistan
| | - Sher Alam Khan
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Samin Jan
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Niamatullah Khan
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Muhammad Muzammal
- Gomal Center of Biochemistry and Biotechnology, Gomal University, D.I.Khan, Khyber Pakhtunkhwa, Pakistan
| | - Sumra Wajid Abbasi
- NUMS Department of Biological Sciences, National University of Medical Sciences, The Mall, Rawalpindi, Punjab, Pakistan
| | - Naseebullah Kakar
- Department of Biotechnology, Faculty of Life Sciences and Informatics, BUITEMS, Quetta, Pakistan
- Institute of Human Genetics, Universitätsklinikum Schleswig-Holstein, Lübeck, Germany
| | - Zia Ur Rehman
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Muzammil Ahmad Khan
- Gomal Center of Biochemistry and Biotechnology, Gomal University, D.I.Khan, Khyber Pakhtunkhwa, Pakistan
| | - Muhammad Usman Mirza
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON, Canada
| | - Noor Muhammad
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Saadullah Khan
- Department of Biotechnology and Genetic Engineering, Kohat University of Science and Technology (KUST), Kohat, Khyber Pakhtunkhwa, Pakistan
| | - Naveed Wasif
- Institute of Human Genetics, Ulm University, and Ulm University Medical Center, Ulm, Germany
- Institute of Human Genetics, University Hospital Schleswig-Holstein, Kiel, Germany
| |
Collapse
|
4
|
Wankowicz SA, de Oliveira SH, Hogan DW, van den Bedem H, Fraser JS. Ligand binding remodels protein side-chain conformational heterogeneity. eLife 2022; 11:e74114. [PMID: 35312477 PMCID: PMC9084896 DOI: 10.7554/elife.74114] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Accepted: 03/18/2022] [Indexed: 12/15/2022] Open
Abstract
While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent-exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decreases and relative hydrophobicity increases. Across a series of 13 inhibitor-bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from nuclear magnetic resonance studies suggesting that residual side-chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.
Collapse
Affiliation(s)
- Stephanie A Wankowicz
- Department of Bioengineering and Therapeutic Sciences, University of California, San FranciscoSan FranciscoUnited States
- Biophysics Graduate Program, University of California San FranciscoSan FranciscoUnited States
| | | | - Daniel W Hogan
- Department of Bioengineering and Therapeutic Sciences, University of California, San FranciscoSan FranciscoUnited States
| | - Henry van den Bedem
- Department of Bioengineering and Therapeutic Sciences, University of California, San FranciscoSan FranciscoUnited States
- Atomwise Inc.San FranciscoUnited States
| | - James S Fraser
- Department of Bioengineering and Therapeutic Sciences, University of California, San FranciscoSan FranciscoUnited States
| |
Collapse
|
5
|
Muneer I, Ahmad S, Naz A, Abbasi SW, Alblihy A, Aloliqi AA, Aba Alkhayl FF, Alrumaihi F, Ahmad S, El Bakri Y, Tahir Ul Qamar M. Discovery of Novel Inhibitors From Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation. Front Mol Biosci 2021; 8:716735. [PMID: 34765641 PMCID: PMC8576517 DOI: 10.3389/fmolb.2021.716735] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2021] [Accepted: 09/20/2021] [Indexed: 12/15/2022] Open
Abstract
V-domain Ig suppressor of T cell activation (VISTA) is an immune checkpoint and is a type I transmembrane protein. VISTA is linked to immunotherapy resistance, and it is a potential immune therapeutic target, especially for triple-negative breast cancer. It expresses at a high concentration in regulatory T cells and myeloid-derived suppressor cells, and its functional blockade is found to delay tumor growth. A useful medicinal plant database for drug designing (MPD3), which is a collection of phytochemicals from diverse plant families, was employed in virtual screening against VISTA to prioritize natural inhibitors against VISTA. Three compounds, Paratocarpin K (PubChem ID: 14187087), 3-(1H-Indol-3-yl)-2-(trimethylazaniumyl)propanoate (PubChem ID: 3861164), and 2-[(5-Benzyl-4-ethyl-1,2,4-triazol-3-yl)sulfanylmethyl]-5-methyl-1,3,4-oxadiazole (PubChem ID: 6494266), having binding energies stronger than -6 kcal/mol were found to have two common hydrogen bond interactions with VISTA active site residues: Arg54 and Arg127. The dynamics of the compound-VISTA complexes were further explored to infer binding stability of the systems. Results revealed that the compound 14187087 and 6494266 systems are highly stable with an average RMSD of 1.31 Å. Further affirmation on the results was achieved by running MM-GBSA on the MD simulation trajectories, which re-ranked 14187087 as the top-binder with a net binding energy value of -33.33 kcal/mol. In conclusion, the present study successfully predicted natural compounds that have the potential to block the function of VISTA and therefore can be utilized further in experimental studies to validate their real anti-VISTA activity.
Collapse
Affiliation(s)
- Iqra Muneer
- School of Life Sciences, University of Science and Technology of China, Hefei, China
| | - Sajjad Ahmad
- Department of Health and Biological Sciences, Abasyn University, Peshawar, Pakistan
| | - Anam Naz
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Sumra Wajid Abbasi
- NUMS Department of Biological Sciences, National University of Medical Sciences, Rawalpindi, Pakistan
| | - Adel Alblihy
- Medical Center, King Fahad Security College (KFSC), Riyadh, Saudi Arabia
| | - Abdulaziz A Aloliqi
- Department of Medical Biotechnology, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Faris F Aba Alkhayl
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Faris Alrumaihi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Sarfraz Ahmad
- Department of Chemistry, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia
| | - Youness El Bakri
- Department of Theoretical and Applied Chemistry, South Ural State University, Chelyabinsk, Russia
| | | |
Collapse
|
6
|
Hauptstein N, Pouyan P, Kehrein J, Dirauf M, Driessen MD, Raschig M, Licha K, Gottschaldt M, Schubert US, Haag R, Meinel L, Sotriffer C, Lühmann T. Molecular Insights into Site-Specific Interferon-α2a Bioconjugates Originated from PEG, LPG, and PEtOx. Biomacromolecules 2021; 22:4521-4534. [PMID: 34643378 DOI: 10.1021/acs.biomac.1c00775] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Conjugation of biologics with polymers modulates their pharmacokinetics, with polyethylene glycol (PEG) as the gold standard. We compared alternative polymers and two types of cyclooctyne linkers (BCN/DBCO) for bioconjugation of interferon-α2a (IFN-α2a) using 10 kDa polymers including linear mPEG, poly(2-ethyl-2-oxazoline) (PEtOx), and linear polyglycerol (LPG). IFN-α2a was azide functionalized via amber codon expansion and bioorthogonally conjugated to all cyclooctyne linked polymers. Polymer conjugation did not impact IFN-α2a's secondary structure and only marginally reduced IFN-α2a's bioactivity. In comparison to PEtOx, the LPG polymer attached via the less rigid cyclooctyne linker BCN was found to stabilize IFN-α2a against thermal stress. These findings were further detailed by molecular modeling studies which showed a modulation of protein flexibility upon PEtOx conjugation and a reduced amount of protein native contacts as compared to PEG and LPG originated bioconjugates. Polymer interactions with IFN-α2a were further assessed via a limited proteolysis (LIP) assay, which resulted in comparable proteolytic cleavage patterns suggesting weak interactions with the protein's surface. In conclusion, both PEtOx and LPG bioconjugates resulted in a similar biological outcome and may become promising PEG alternatives for bioconjugation.
Collapse
Affiliation(s)
- Niklas Hauptstein
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Paria Pouyan
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Josef Kehrein
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Michael Dirauf
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Marc D Driessen
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Martina Raschig
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Kai Licha
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Michael Gottschaldt
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Ulrich S Schubert
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Rainer Haag
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Lorenz Meinel
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany.,Helmholtz Institute for RNA-Based Infection Research (HIRI), Helmholtz Center for Infection Research (HZI), 97080 Würzburg, Germany
| | - Christoph Sotriffer
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Tessa Lühmann
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| |
Collapse
|
7
|
Bentham AR, Petit-Houdenot Y, Win J, Chuma I, Terauchi R, Banfield MJ, Kamoun S, Langner T. A single amino acid polymorphism in a conserved effector of the multihost blast fungus pathogen expands host-target binding spectrum. PLoS Pathog 2021; 17:e1009957. [PMID: 34758051 PMCID: PMC8608293 DOI: 10.1371/journal.ppat.1009957] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Revised: 11/22/2021] [Accepted: 09/14/2021] [Indexed: 12/14/2022] Open
Abstract
Accelerated gene evolution is a hallmark of pathogen adaptation and specialization following host-jumps. However, the molecular processes associated with adaptive evolution between host-specific lineages of a multihost plant pathogen remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), host specialization on different grass hosts is generally associated with dynamic patterns of gain and loss of virulence effector genes that tend to define the distinct genetic lineages of this pathogen. Here, we unravelled the biochemical and structural basis of adaptive evolution of APikL2, an exceptionally conserved paralog of the well-studied rice-lineage specific effector AVR-Pik. Whereas AVR-Pik and other members of the six-gene AVR-Pik family show specific patterns of presence/absence polymorphisms between grass-specific lineages of M. oryzae, APikL2 stands out by being ubiquitously present in all blast fungus lineages from 13 different host species. Using biochemical, biophysical and structural biology methods, we show that a single aspartate to asparagine polymorphism expands the binding spectrum of APikL2 to host proteins of the heavy-metal associated (HMA) domain family. This mutation maps to one of the APikL2-HMA binding interfaces and contributes to an altered hydrogen-bonding network. By combining phylogenetic ancestral reconstruction with an analysis of the structural consequences of allelic diversification, we revealed a common mechanism of effector specialization in the AVR-Pik/APikL2 family that involves two major HMA-binding interfaces. Together, our findings provide a detailed molecular evolution and structural biology framework for diversification and adaptation of a fungal pathogen effector family following host-jumps.
Collapse
Affiliation(s)
- Adam R. Bentham
- Department of Biological Chemistry, John Innes Centre, Norwich Research Park, Norwich, United Kingdom
| | - Yohann Petit-Houdenot
- The Sainsbury Laboratory, University of East Anglia, Norwich Research Park, Norwich, United Kingdom
- Université Paris-Saclay, INRAE, AgroParisTech, UMR BIOGER, Thiverval-Grignon, France
| | - Joe Win
- The Sainsbury Laboratory, University of East Anglia, Norwich Research Park, Norwich, United Kingdom
| | - Izumi Chuma
- Obihiro University of Agriculture and Veterinary Medicine, Obihiro, Japan
| | - Ryohei Terauchi
- Kyoto University, Kyoto, Japan
- Iwate Biotechnology Research Center, Kitakami, Japan
| | - Mark J. Banfield
- Department of Biological Chemistry, John Innes Centre, Norwich Research Park, Norwich, United Kingdom
| | - Sophien Kamoun
- The Sainsbury Laboratory, University of East Anglia, Norwich Research Park, Norwich, United Kingdom
| | - Thorsten Langner
- The Sainsbury Laboratory, University of East Anglia, Norwich Research Park, Norwich, United Kingdom
| |
Collapse
|
8
|
Nutschel C, Coscolín C, David B, Mulnaes D, Ferrer M, Jaeger KE, Gohlke H. Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones. J Chem Inf Model 2021; 61:2383-2395. [PMID: 33949194 DOI: 10.1021/acs.jcim.1c00152] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.
Collapse
Affiliation(s)
- Christina Nutschel
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Cristina Coscolín
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Benoit David
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Daniel Mulnaes
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Manuel Ferrer
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, 52425 Jülich, Germany.,Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Holger Gohlke
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
9
|
Chen J, Wang W, Pang L, Zhu W. Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics. Phys Chem Chem Phys 2021; 22:21238-21250. [PMID: 32930679 DOI: 10.1039/d0cp03766d] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China.
| |
Collapse
|
10
|
Temperature dependent aggregation mechanism and pathway of lysozyme: By all atom and coarse grained molecular dynamics simulation. J Mol Graph Model 2020; 103:107816. [PMID: 33291026 DOI: 10.1016/j.jmgm.2020.107816] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2020] [Revised: 11/18/2020] [Accepted: 11/18/2020] [Indexed: 11/21/2022]
Abstract
Aggregation of protein causes various diseases including Alzheimer's disease, Parkinson's disease, and type II diabetes. It was found that aggregation of protein depends on many factors like temperature, pH, salt type, salt concentration, ionic strength, protein concentration, co solutes. Here we have tried to capture the aggregation mechanism and pathway of hen egg white lysozyme using molecular dynamics simulations at two different temperatures; 300 K and 340 K. Along with the all atom simulations to get the atomistic details of aggregation mechanism, we have used coarse grained simulation with MARTINI force field to monitor the aggregation for longer duration. Our results suggest that due to the aggregation, changes in the conformation of lysozyme are more at 340 K than at 300 K. The change in the conformation of the lysozyme at 300 K is mainly due to aggregation where at 340 K change in conformation of lysozyme is due to both aggregation and temperature. Also, a more compact aggregated system is formed at 340 K.
Collapse
|
11
|
Grahame DSA, Dupuis JH, Bryksa BC, Tanaka T, Yada RY. Comparative bioinformatic and structural analyses of pepsin and renin. Enzyme Microb Technol 2020; 141:109632. [PMID: 33051007 DOI: 10.1016/j.enzmictec.2020.109632] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2020] [Revised: 06/25/2020] [Accepted: 07/08/2020] [Indexed: 11/16/2022]
Abstract
Pepsin, the archetypal pepsin-like aspartic protease, is irreversibly denatured when exposed to neutral pH conditions whereas renin, a structural homologue of pepsin, is fully stable and optimally active in the same conditions despite sharing highly similar enzyme architecture. To gain insight into the structural determinants of differential aspartic protease pH stability, the present study used comparative bioinformatic and structural analyses. In pepsin, an abundance of polar and aspartic acid residues were identified, a common trait with other acid-stable enzymes. Conversely, renin was shown to have increased levels of basic amino acids. In both pepsin and renin, the solvent exposure of these charged groups was high. Having similar overall acidic residue content, the solvent-exposed basic residues may allow for extensive salt bridge formation in renin, whereas in pepsin, these residues are protonated and serve to form stabilizing hydrogen bonds at low pH. Relative differences in structure and sequence in the turn and joint regions of the β-barrel and ψ-loop in both the N- and C-terminal lobes were identified as regions of interest in defining divergent pH stability. Compared to the structural rigidity of renin, pepsin has more instability associated with the N-terminus, specifically the B/C connector. By contrast, renin exhibits greater C-terminal instability in turn and connector regions. Overall, flexibility differences in connector regions, and amino acid composition, particularly in turn and joint regions of the β-barrel and ψ-loops, likely play defining roles in determining pH stability for renin and pepsin.
Collapse
Affiliation(s)
- Douglas S A Grahame
- Department of Food Science, Ontario Agricultural College, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - John H Dupuis
- Food, Nutrition, and Health Program, Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC, V6T 1Z4 Canada
| | - Brian C Bryksa
- Department of Food Science, Ontario Agricultural College, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Takuji Tanaka
- Department of Food and Bioproduct Sciences, College of Agriculture and Bioresources, University of Saskatchewan, Saskatoon, SK, S7N 5A8 Canada
| | - Rickey Y Yada
- Department of Food Science, Ontario Agricultural College, University of Guelph, Guelph, ON, N1G 2W1, Canada; Food, Nutrition, and Health Program, Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC, V6T 1Z4 Canada.
| |
Collapse
|
12
|
Deleterious Variants in WNT10A, EDAR, and EDA Causing Isolated and Syndromic Tooth Agenesis: A Structural Perspective from Molecular Dynamics Simulations. Int J Mol Sci 2019; 20:ijms20215282. [PMID: 31652981 PMCID: PMC6862269 DOI: 10.3390/ijms20215282] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Revised: 10/18/2019] [Accepted: 10/20/2019] [Indexed: 02/07/2023] Open
Abstract
The dental abnormalities are the typical features of many ectodermal dysplasias along with congenital malformations of nails, skin, hair, and sweat glands. However, several reports of non-syndromic/isolated tooth agenesis have also been found in the literature. The characteristic features of hypohidrotic ectodermal dysplasia (HED) comprise of hypodontia/oligodontia, along with hypohidrosis/anhidrosis, and hypotrichosis. Pathogenic variants in EDA, EDAR, EDARADD, and TRAF6, cause the phenotypic expression of HED. Genetic alterations in EDA and WNT10A cause particularly non-syndromic/isolated oligodontia. In the current project, we recruited 57 patients of 17 genetic pedigrees (A-Q) from different geographic regions of the world, including Pakistan, Egypt, Saudi Arabia, and Syria. The molecular investigation of different syndromic and non-syndromic dental conditions, including hypodontia, oligodontia, generalized odontodysplasia, and dental crowding was carried out by using exome and Sanger sequencing. We have identified a novel missense variant (c.311G>A; p.Arg104His) in WNT10A in three oligodontia patients of family A, two novel sequence variants (c.207delinsTT, p.Gly70Trpfs*25 and c.1300T>G; p.Try434Gly) in EDAR in three patients of family B and four patients of family C, respectively. To better understand the structural and functional consequences of missense variants in WNT10A and EDAR on the stability of the proteins, we have performed extensive molecular dynamic (MD) simulations. We have also identified three previously reported pathogenic variants (c.1076T>C; p.Met359Thr), (c.1133C>T; p.Thr378Met) and (c.594_595insC; Gly201Argfs*39) in EDA in family D (four patients), E (two patients) and F (one patient), correspondingly. Presently, our data explain the genetic cause of 18 syndromic and non-syndromic tooth agenesis patients in six autosomal recessive and X-linked pedigrees (A-F), which expand the mutational spectrum of these unique clinical manifestations.
Collapse
|
13
|
Erman B. A computational model for controlling conformational cooperativity and function in proteins. Proteins 2018; 86:1001-1009. [PMID: 30051502 DOI: 10.1002/prot.25535] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2018] [Revised: 05/08/2018] [Accepted: 06/02/2018] [Indexed: 01/02/2023]
Abstract
We present a computational model that allows for rapid prediction of correlations among a set of residue pairs when the fluctuations of another set of residues are perturbed. The simple theory presented here is based on the knowledge of the fluctuation covariance matrix only. In this sense, the theory is model independent and therefore universal. Perturbation of any set of fluctuations and the resulting response of the remaining set are calculated using conditional probabilities of a multivariate normal distribution. The model is expected to rapidly and accurately map the consequences of mutations in proteins, as well as allosteric activity and ligand binding. Knowledge of triple correlations of fluctuations of residues i, j, and k, 〈 Δ R i Δ R j Δ R k 〉 emerges as the necessary source of information for controlling residue pairs by perturbing a distant residue. Triple correlations have not received wide attention in literature. Perturbation-response-function relations for ubiquitin (UBQ) are discussed as an example. Covariance matrix for UBQ obtained from the Gaussian Network Model combined with the present computational algorithm is able to reflect the millisecond molecular dynamics correlations and observed NMR results. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Burak Erman
- Department of Chemical and Biological Engineering; Koc University; Sariyer Istanbul Turkey
| |
Collapse
|
14
|
Carugo O. Atomic displacement parameters in structural biology. Amino Acids 2018; 50:775-786. [DOI: 10.1007/s00726-018-2574-y] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Accepted: 04/19/2018] [Indexed: 01/14/2023]
|
15
|
Pfleger C, Minges A, Boehm M, McClendon CL, Torella R, Gohlke H. Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery. J Chem Theory Comput 2017; 13:6343-6357. [PMID: 29112408 DOI: 10.1021/acs.jctc.7b00529] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Allostery describes the functional coupling between sites in biomolecules. Recently, the role of changes in protein dynamics for allosteric communication has been highlighted. A quantitative and predictive description of allostery is fundamental for understanding biological processes. Here, we integrate an ensemble-based perturbation approach with the analysis of biomolecular rigidity and flexibility to construct a model of dynamic allostery. Our model, by definition, excludes the possibility of conformational changes, evaluates static, not dynamic, properties of molecular systems, and describes allosteric effects due to ligand binding in terms of a novel free-energy measure. We validated our model on three distinct biomolecular systems: eglin c, protein tyrosine phosphatase 1B, and the lymphocyte function-associated antigen 1 domain. In all cases, it successfully identified key residues for signal transmission in very good agreement with the experiment. It correctly and quantitatively discriminated between positively or negatively cooperative effects for one of the systems. Our model should be a promising tool for the rational discovery of novel allosteric drugs.
Collapse
Affiliation(s)
- Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf , Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Alexander Minges
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf , Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Markus Boehm
- Medicinal Sciences, Pfizer, Inc. , 1 Portland Street, Cambridge, Massachusetts 02139, United States
| | - Christopher L McClendon
- Medicinal Sciences, Pfizer, Inc. , 1 Portland Street, Cambridge, Massachusetts 02139, United States
| | - Rubben Torella
- Medicinal Sciences, Pfizer, Inc. , 1 Portland Street, Cambridge, Massachusetts 02139, United States
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf , Universitätsstr. 1, 40225 Düsseldorf, Germany
| |
Collapse
|
16
|
Hanke CA, Gohlke H. Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch. PLoS One 2017; 12:e0179271. [PMID: 28640851 PMCID: PMC5480868 DOI: 10.1371/journal.pone.0179271] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2016] [Accepted: 05/27/2017] [Indexed: 12/18/2022] Open
Abstract
Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.
Collapse
Affiliation(s)
- Christian A. Hanke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
17
|
Hermans SM, Pfleger C, Nutschel C, Hanke CA, Gohlke H. Rigidity theory for biomolecules: concepts, software, and applications. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1311] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Susanne M.A. Hermans
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christopher Pfleger
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christina Nutschel
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christian A. Hanke
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| |
Collapse
|
18
|
Gohlke H, Ben-Shalom IY, Kopitz H, Pfeiffer-Marek S, Baringhaus KH. Rigidity Theory-Based Approximation of Vibrational Entropy Changes upon Binding to Biomolecules. J Chem Theory Comput 2017; 13:1495-1502. [DOI: 10.1021/acs.jctc.7b00014] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Holger Gohlke
- Institute
for Pharmaceutical and Medicinal Chemistry, Department of Mathematics
and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Ido Y. Ben-Shalom
- Institute
for Pharmaceutical and Medicinal Chemistry, Department of Mathematics
and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Hannes Kopitz
- Institute
for Pharmaceutical and Medicinal Chemistry, Department of Mathematics
and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Stefania Pfeiffer-Marek
- R&D/Pre-Development Sciences, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| | - Karl-Heinz Baringhaus
- R&D Resources/Site Direction, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| |
Collapse
|
19
|
Raschka S, Bemister-Buffington J, Kuhn LA. Detecting the native ligand orientation by interfacial rigidity: SiteInterlock. Proteins 2016; 84:1888-1901. [DOI: 10.1002/prot.25172] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2016] [Revised: 09/19/2016] [Accepted: 09/27/2016] [Indexed: 01/27/2023]
Affiliation(s)
- Sebastian Raschka
- Department of Biochemistry and Molecular Biology; Michigan State University; East Lansing Michigan 48824 USA
| | - Joseph Bemister-Buffington
- Department of Biochemistry and Molecular Biology; Michigan State University; East Lansing Michigan 48824 USA
| | - Leslie A. Kuhn
- Department of Biochemistry and Molecular Biology; Michigan State University; East Lansing Michigan 48824 USA
- Department of Computer Science and Engineering; Michigan State University; East Lansing Michigan 48824 USA
| |
Collapse
|
20
|
Li F, Liu J, Liu N, Kuhn LA, Garavito RM, Ferguson-Miller S. Translocator Protein 18 kDa (TSPO): An Old Protein with New Functions? Biochemistry 2016; 55:2821-31. [PMID: 27074410 DOI: 10.1021/acs.biochem.6b00142] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Translocator protein 18 kDa (TSPO) was previously known as the peripheral benzodiazepine receptor (PBR) in eukaryotes, where it is mainly localized to the mitochondrial outer membrane. Considerable evidence indicates that it plays regulatory roles in steroidogenesis and apoptosis and is involved in various human diseases, such as metastatic cancer, Alzheimer's and Parkinson's disease, inflammation, and anxiety disorders. Ligands of TSPO are widely used as diagnostic tools and treatment options, despite there being no clear understanding of the function of TSPO. An ortholog in the photosynthetic bacterium Rhodobacter was independently discovered as the tryptophan-rich sensory protein (TspO) and found to play a role in the response to changes in oxygen and light conditions that regulate photosynthesis and respiration. As part of this highly conserved protein family found in all three kingdoms, the rat TSPO is able to rescue the knockout phenotype in Rhodobacter, indicating functional as well as structural conservation. Recently, a major breakthrough in the field was achieved: the determination of atomic-resolution structures of TSPO from different species by several independent groups. This now allows us to reexamine the function of TSPO with a molecular perspective. In this review, we focus on recently determined structures of TSPO and their implications for potential functions of this ubiquitous multifaceted protein. We suggest that TSPO is an ancient bacterial receptor/stress sensor that has developed additional interactions, partners, and roles in its mitochondrial outer membrane environment in eukaryotes.
Collapse
Affiliation(s)
- Fei Li
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Jian Liu
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Nan Liu
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States.,Department of Computer Science and Engineering, Michigan State University , East Lansing, Michigan 48824-1319, United States.,Department of Chemistry, Michigan State University , East Lansing, Michigan 48824-1319, United States
| | - Leslie A Kuhn
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States.,Department of Computer Science and Engineering, Michigan State University , East Lansing, Michigan 48824-1319, United States
| | - R Michael Garavito
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Shelagh Ferguson-Miller
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|
21
|
Rathi PC, Fulton A, Jaeger KE, Gohlke H. Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis. PLoS Comput Biol 2016; 12:e1004754. [PMID: 27003415 PMCID: PMC4803202 DOI: 10.1371/journal.pcbi.1004754] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2015] [Accepted: 01/14/2016] [Indexed: 11/29/2022] Open
Abstract
Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability. Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.
Collapse
Affiliation(s)
- Prakash Chandra Rathi
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University, Düsseldorf, Germany
| | - Alexander Fulton
- Institute of Molecular Enzyme Technology, Heinrich-Heine-University, Düsseldorf, Germany
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich-Heine-University, Düsseldorf, Germany
- Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany
- * E-mail: (KEJ); (HG)
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University, Düsseldorf, Germany
- * E-mail: (KEJ); (HG)
| |
Collapse
|
22
|
Abstract
Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.
Collapse
Affiliation(s)
- H. Pezeshgi Modarres
- Molecular Cell Biomechanics Laboratory
- Departments of Bioengineering and Mechanical Engineering
- University of California Berkeley
- Berkeley
- USA
| | - M. R. Mofrad
- Molecular Cell Biomechanics Laboratory
- Departments of Bioengineering and Mechanical Engineering
- University of California Berkeley
- Berkeley
- USA
| | - A. Sanati-Nezhad
- BioMEMS and Bioinspired Microfluidic Laboratory
- Department of Mechanical and Manufacturing Engineering
- University of Calgary
- Calgary
- Canada
| |
Collapse
|
23
|
P C, R R. A systematic molecular dynamics approach to the structural characterization of amyloid aggregation propensity of β2-microglobulin mutant D76N. MOLECULAR BIOSYSTEMS 2016; 12:850-9. [DOI: 10.1039/c5mb00759c] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Beta-2 microglobulin (β2m) is an amyloidogenic protein belonging to the immunoglobulin superfamily, responsible for the dialysis-related amyloidosis (DRA).
Collapse
Affiliation(s)
- Chandrasekaran P
- Bioinformatics Division
- School of Biosciences and Technology
- VIT University
- Vellore 632 014
- India
| | - Rajasekaran R
- Bioinformatics Division
- School of Biosciences and Technology
- VIT University
- Vellore 632 014
- India
| |
Collapse
|
24
|
Xie H, Li Y, Yu F, Xie X, Qiu K, Fu J. An Investigation of Molecular Docking and Molecular Dynamic Simulation on Imidazopyridines as B-Raf Kinase Inhibitors. Int J Mol Sci 2015; 16:27350-61. [PMID: 26580609 PMCID: PMC4661884 DOI: 10.3390/ijms161126026] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2015] [Revised: 10/31/2015] [Accepted: 11/06/2015] [Indexed: 11/16/2022] Open
Abstract
In the recent cancer treatment, B-Raf kinase is one of key targets. Nowadays, a group of imidazopyridines as B-Raf kinase inhibitors have been reported. In order to investigate the interaction between this group of inhibitors and B-Raf kinase, molecular docking, molecular dynamic (MD) simulation and binding free energy (ΔGbind) calculation were performed in this work. Molecular docking was carried out to identify the key residues in the binding site, and MD simulations were performed to determine the detail binding mode. The results obtained from MD simulation reveal that the binding site is stable during the MD simulations, and some hydrogen bonds (H-bonds) in MD simulations are different from H-bonds in the docking mode. Based on the obtained MD trajectories, ΔGbind was computed by using Molecular Mechanics Generalized Born Surface Area (MM-GBSA), and the obtained energies are consistent with the activities. An energetic analysis reveals that both electrostatic and van der Waals contributions are important to ΔGbind, and the unfavorable polar solvation contribution results in the instability of the inhibitor with the lowest activity. These results are expected to understand the binding between B-Raf and imidazopyridines and provide some useful information to design potential B-Raf inhibitors.
Collapse
Affiliation(s)
- Huiding Xie
- Department of Chemistry, Yunnan University, Kunming 650091, China.
- Department of Chemistry, School of Pharmaceutical Science & Yunnan Key Laboratory of Pharmacology for Natural Products, Kunming Medical University, Kunming 650500, China.
| | - Yupeng Li
- Department of Chemistry, School of Pharmaceutical Science & Yunnan Key Laboratory of Pharmacology for Natural Products, Kunming Medical University, Kunming 650500, China.
| | - Fang Yu
- Department of Chemistry, School of Pharmaceutical Science & Yunnan Key Laboratory of Pharmacology for Natural Products, Kunming Medical University, Kunming 650500, China.
| | - Xiaoguang Xie
- Department of Chemistry, Yunnan University, Kunming 650091, China.
| | - Kaixiong Qiu
- Department of Chemistry, School of Pharmaceutical Science & Yunnan Key Laboratory of Pharmacology for Natural Products, Kunming Medical University, Kunming 650500, China.
| | - Jijun Fu
- Department of Chemistry, School of Pharmaceutical Science & Yunnan Key Laboratory of Pharmacology for Natural Products, Kunming Medical University, Kunming 650500, China.
| |
Collapse
|
25
|
Aβ1-25-Derived Sphingolipid-Domain Tracer Peptide SBD Interacts with Membrane Ganglioside Clusters via a Coil-Helix-Coil Motif. Int J Mol Sci 2015; 16:26318-32. [PMID: 26540054 PMCID: PMC4661814 DOI: 10.3390/ijms161125955] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2015] [Revised: 10/19/2015] [Accepted: 10/22/2015] [Indexed: 12/30/2022] Open
Abstract
The Amyloid-β (Aβ)-derived, sphingolipid binding domain (SBD) peptide is a fluorescently tagged probe used to trace the diffusion behavior of sphingolipid-containing microdomains in cell membranes through binding to a constellation of glycosphingolipids, sphingomyelin, and cholesterol. However, the molecular details of the binding mechanism between SBD and plasma membrane domains remain unclear. Here, to investigate how the peptide recognizes the lipid surface at an atomically detailed level, SBD peptides in the environment of raft-like bilayers were examined in micro-seconds-long molecular dynamics simulations. We found that SBD adopted a coil-helix-coil structural motif, which binds to multiple GT1b gangliosides via salt bridges and CH–π interactions. Our simulation results demonstrate that the CH–π and electrostatic forces between SBD monomers and GT1b gangliosides clusters are the main driving forces in the binding process. The presence of the fluorescent dye and linker molecules do not change the binding mechanism of SBD probes with gangliosides, which involves the helix-turn-helix structural motif that was suggested to constitute a glycolipid binding domain common to some sphingolipid interacting proteins, including HIV gp120, prion, and Aβ.
Collapse
|
26
|
Tiberti M, Invernizzi G, Papaleo E. (Dis)similarity Index To Compare Correlated Motions in Molecular Simulations. J Chem Theory Comput 2015; 11:4404-14. [PMID: 26575932 DOI: 10.1021/acs.jctc.5b00512] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Molecular dynamics (MD) simulations are widely used to complement or guide experimental studies in the characterization of protein dynamics, thanks to improvements in force-field accuracy, along with in the software and hardware to sample the conformational landscape of proteins. Among the different applications of MD simulations, the study of correlated motions is largely employed for different purposes. Several metrics have been developed to describe correlated motions in the MD ensemble, such as methods based on Pearson Correlation or Mutual Information. Cross-correlation analysis of MD trajectories is indeed appealing not only to identify residues characterized by coupled fluctuations in protein structures but also since it can be used to extrapolate motions along directions in which major conformational changes should occur, for example on longer time scales than the ones that are actually simulated. Nevertheless, most of the MD studies employ average correlation maps and mostly in a qualitative way, even when different systems or different replicates of the same system are compared. The broad application of correlation metrics in the analysis of MD simulations, especially for comparative purposes, requires a step forward toward more quantitative and accurate comparisons. We thus here employed a simple but effective index, which is based on a normalized Frobenius norm of the differences between protein correlation maps, to compare correlated motions. We applied this index for a quantitative comparison of correlated motions from MD simulations of seven proteins of different size and fold. We also employed the index to assess the robustness of correlation description when multi-replicate MD simulations of a same system are used, and we compared our index to metrics for comparison of structural ensembles such as Root Mean Square Inner Product and the Bhattacharyya Coefficient.
Collapse
Affiliation(s)
- Matteo Tiberti
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Gaetano Invernizzi
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Elena Papaleo
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| |
Collapse
|
27
|
Karshikoff A, Nilsson L, Ladenstein R. Rigidity versus flexibility: the dilemma of understanding protein thermal stability. FEBS J 2015; 282:3899-917. [PMID: 26074325 DOI: 10.1111/febs.13343] [Citation(s) in RCA: 176] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2015] [Revised: 05/17/2015] [Accepted: 06/09/2015] [Indexed: 01/19/2023]
Abstract
The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.
Collapse
Affiliation(s)
- Andrey Karshikoff
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| | - Rudolf Ladenstein
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| |
Collapse
|
28
|
Aramini JM, Vorobiev SM, Tuberty LM, Janjua H, Campbell ET, Seetharaman J, Su M, Huang YJ, Acton TB, Xiao R, Tong L, Montelione GT. The RAS-Binding Domain of Human BRAF Protein Serine/Threonine Kinase Exhibits Allosteric Conformational Changes upon Binding HRAS. Structure 2015; 23:1382-1393. [PMID: 26165597 DOI: 10.1016/j.str.2015.06.003] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2015] [Revised: 05/27/2015] [Accepted: 06/04/2015] [Indexed: 02/07/2023]
Abstract
RAS binding is a critical step in the activation of BRAF protein serine/threonine kinase and stimulation of the mitogen-activated protein kinase signaling pathway. Mutations in both RAS and BRAF are associated with many human cancers. Here, we report the solution nuclear magnetic resonance (NMR) and X-ray crystal structures of the RAS-binding domain (RBD) from human BRAF. We further studied the complex between BRAF RBD and the GppNHp bound form of HRAS in solution. Backbone, side-chain, and (19)F NMR chemical shift perturbations reveal unexpected changes distal to the RAS-binding face that extend through the core of the RBD structure. Moreover, backbone amide hydrogen/deuterium exchange NMR data demonstrate conformational ensemble changes in the RBD core structure upon complex formation. These changes in BRAF RBD reveal a basis for allosteric regulation of BRAF structure and function, and suggest a mechanism by which RAS binding can signal the drastic domain rearrangements required for activation of BRAF kinase.
Collapse
Affiliation(s)
- James M Aramini
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA.
| | - Sergey M Vorobiev
- Department of Biological Sciences, Northeast Structural Genomics Consortium, Columbia University, New York, NY 10027, USA
| | - Lynda M Tuberty
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Haleema Janjua
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Elliot T Campbell
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Jayaraman Seetharaman
- Department of Biological Sciences, Northeast Structural Genomics Consortium, Columbia University, New York, NY 10027, USA
| | - Min Su
- Department of Biological Sciences, Northeast Structural Genomics Consortium, Columbia University, New York, NY 10027, USA
| | - Yuanpeng J Huang
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Thomas B Acton
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Rong Xiao
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Liang Tong
- Department of Biological Sciences, Northeast Structural Genomics Consortium, Columbia University, New York, NY 10027, USA
| | - Gaetano T Montelione
- Center for Advanced Biotechnology and Medicine, Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Northeast Structural Genomics Consortium, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Biochemistry and Molecular Biology, Robert Wood Johnson Medical School, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA.
| |
Collapse
|
29
|
Rathi PC, Jaeger KE, Gohlke H. Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis. PLoS One 2015; 10:e0130289. [PMID: 26147762 PMCID: PMC4493141 DOI: 10.1371/journal.pone.0130289] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Accepted: 05/18/2015] [Indexed: 11/24/2022] Open
Abstract
Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R2 = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.
Collapse
Affiliation(s)
- Prakash Chandra Rathi
- Institute of Pharmaceutical and Medical Chemistry, Heinrich-Heine-University, Düsseldorf, Germany
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich-Heine-University, Düsseldorf, Germany
- Institute of Bio- and Geosciences IBG-1: Biotechnology, Research Centre Jülich, Jülich, Germany
| | - Holger Gohlke
- Institute of Pharmaceutical and Medical Chemistry, Heinrich-Heine-University, Düsseldorf, Germany
- * E-mail:
| |
Collapse
|
30
|
Bhutani I, Loharch S, Gupta P, Madathil R, Parkesh R. Structure, dynamics, and interaction of Mycobacterium tuberculosis (Mtb) DprE1 and DprE2 examined by molecular modeling, simulation, and electrostatic studies. PLoS One 2015; 10:e0119771. [PMID: 25789990 PMCID: PMC4366402 DOI: 10.1371/journal.pone.0119771] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2014] [Accepted: 02/03/2015] [Indexed: 11/18/2022] Open
Abstract
The enzymes decaprenylphosphoryl-β-D-ribose oxidase (DprE1) and decaprenylphosphoryl-β-D-ribose-2-epimerase (DprE2) catalyze epimerization of decaprenylphosporyl ribose (DPR) todecaprenylphosporyl arabinose (DPA) and are critical for the survival of Mtb. Crystal structures of DprE1 so far reported display significant disordered regions and no structural information is known for DprE2. We used homology modeling, protein threading, molecular docking and dynamics studies to investigate the structural and dynamic features of Mtb DprE1 and DprE2 and DprE1-DprE2 complex. A three-dimensional model for DprE2 was generated using the threading approach coupled with ab initio modeling. A 50 ns simulation of DprE1 and DprE2 revealed the overall stability of the structures. Principal Component Analysis (PCA) demonstrated the convergence of sampling in both DprE1 and DprE2. In DprE1, residues in the 269–330 area showed considerable fluctuation in agreement with the regions of disorder observed in the reported crystal structures. In DprE2, large fluctuations were detected in residues 95–113, 146–157, and 197–226. The study combined docking and MD simulation studies to map and characterize the key residues involved in DprE1-DprE2 interaction. A 60 ns MD simulation for DprE1-DprE2 complex was also performed. Analysis of data revealed that the docked complex is stabilized by H-bonding, hydrophobic and ionic interactions. The key residues of DprE1 involved in DprE1-DprE2 interactions belong to the disordered region. We also examined the docked complex of DprE1-BTZ043 to investigate the binding pocket of DprE1 and its interactions with the inhibitor BTZ043. In summary, we hypothesize that DprE1-DprE2 interaction is crucial for the synthesis of DPA and DprE1-DprE2 complex may be a new therapeutic target amenable to pharmacological validation. The findings have important implications in tuberculosis (TB) drug discovery and will facilitate drug development efforts against TB.
Collapse
Affiliation(s)
- Isha Bhutani
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh 160036, India
| | - Saurabh Loharch
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh 160036, India
| | - Pawan Gupta
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh 160036, India
| | - Rethi Madathil
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh 160036, India
| | - Raman Parkesh
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh 160036, India
- * E-mail:
| |
Collapse
|
31
|
Effects of N-glycosylation on protein conformation and dynamics: Protein Data Bank analysis and molecular dynamics simulation study. Sci Rep 2015; 5:8926. [PMID: 25748215 PMCID: PMC4352867 DOI: 10.1038/srep08926] [Citation(s) in RCA: 165] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2014] [Accepted: 02/11/2015] [Indexed: 12/29/2022] Open
Abstract
N-linked glycosylation is one of the most important, chemically complex, and ubiquitous post-translational modifications in all eukaryotes. The N-glycans that are covalently linked to proteins are involved in numerous biological processes. There is considerable interest in developments of general approaches to predict the structural consequences of site-specific glycosylation and to understand how these effects can be exploited in protein design with advantageous properties. In this study, the impacts of N-glycans on protein structure and dynamics are systematically investigated using an integrated computational approach of the Protein Data Bank structure analysis and atomistic molecular dynamics simulations of glycosylated and deglycosylated proteins. Our study reveals that N-glycosylation does not induce significant changes in protein structure, but decreases protein dynamics, likely leading to an increase in protein stability. Overall, these results suggest not only a common role of glycosylation in proteins, but also a need for certain proteins to be properly glycosylated to gain their intrinsic dynamic properties.
Collapse
|
32
|
Fetics SK, Guterres H, Kearney BM, Buhrman G, Ma B, Nussinov R, Mattos C. Allosteric effects of the oncogenic RasQ61L mutant on Raf-RBD. Structure 2015; 23:505-516. [PMID: 25684575 PMCID: PMC7755167 DOI: 10.1016/j.str.2014.12.017] [Citation(s) in RCA: 175] [Impact Index Per Article: 19.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 12/07/2014] [Accepted: 12/17/2014] [Indexed: 10/24/2022]
Abstract
The Ras/Raf/MEK/ERK signal transduction pathway is a major regulator of cell proliferation activated by Ras-guanosine triphosphate (GTP). The oncogenic mutant RasQ61L is not able to hydrolyze GTP in the presence of Raf and thus is a constitutive activator of this mitogenic pathway. The Ras/Raf interaction is essential for the activation of the Raf kinase domain through a currently unknown mechanism. We present the crystal structures of the Ras-GppNHp/Raf-RBD and RasQ61L-GppNHp/Raf-RBD complexes, which, in combination with MD simulations, reveal differences in allosteric interactions leading from the Ras/Raf interface to the Ras calcium-binding site and to the remote Raf-RBD loop L4. In the presence of Raf, the RasQ61L mutant has a rigid switch II relative to the wild-type and increased flexibility at the interface with switch I, which propagates across Raf-RBD. We show that in addition to local perturbations on Ras, RasQ61L has substantial long-range effects on the Ras allosteric lobe and on Raf-RBD.
Collapse
Affiliation(s)
- Susan K Fetics
- Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA 02115, USA; Department of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, NC 27695, USA
| | - Hugo Guterres
- Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA 02115, USA
| | - Bradley M Kearney
- Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA 02115, USA; Department of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, NC 27695, USA
| | - Greg Buhrman
- Department of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, NC 27695, USA
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc., Cancer and Inflammation Program, National Cancer Institute, Frederick, MD 21702, USA
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc., Cancer and Inflammation Program, National Cancer Institute, Frederick, MD 21702, USA; Department of Human Genetics, Sackler Institute of Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Carla Mattos
- Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA 02115, USA; Department of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, NC 27695, USA.
| |
Collapse
|
33
|
Abstract
Riboswitches are noncoding regulatory elements that control gene expression in response to the presence of metabolites, which bind to the aptamer domain. Metabolite binding appears to occur through a combination of conformational selection and induced fit mechanism. This demands to characterize the structural dynamics of the apo state of aptamer domains. In principle, molecular dynamics (MD) simulations can give insights at the atomistic level into the dynamics of the aptamer domain. However, it is unclear to what extent contemporary force fields can bias such insights. Here, we show that the Amber force field ff99 yields the best agreement with detailed experimental observations on differences in the structural dynamics of wild type and mutant aptamer domains of the guanine-sensing riboswitch (Gsw), including a pronounced influence of Mg2+. In contrast, applying ff99 with parmbsc0 and parmχOL modifications (denoted ff10) results in strongly damped motions and overly stable tertiary loop-loop interactions. These results are based on 58 MD simulations with an aggregate simulation time>11 μs, careful modeling of Mg2+ ions, and thorough statistical testing. Our results suggest that the moderate stabilization of the χ-anti region in ff10 can have an unwanted damping effect on functionally relevant structural dynamics of marginally stable RNA systems. This suggestion is supported by crystal structure analyses of Gsw aptamer domains that reveal χ torsions with high-anti values in the most mobile regions. We expect that future RNA force field development will benefit from considering marginally stable RNA systems and optimization toward good representations of dynamics in addition to structural characteristics.
Collapse
Affiliation(s)
- Christian A Hanke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany.
| |
Collapse
|
34
|
Acharya C, Kufareva I, Ilatovskiy AV, Abagyan R. PeptiSite: a structural database of peptide binding sites in 4D. Biochem Biophys Res Commun 2014; 445:717-23. [PMID: 24406170 DOI: 10.1016/j.bbrc.2013.12.132] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Accepted: 12/26/2013] [Indexed: 12/11/2022]
Abstract
We developed PeptiSite, a comprehensive and reliable database of biologically and structurally characterized peptide-binding sites, in which each site is represented by an ensemble of its complexes with protein, peptide and small molecule partners. The unique features of the database include: (1) the ensemble site representation that provides a fourth dimension to the otherwise three dimensional data, (2) comprehensive characterization of the binding site architecture that may consist of a multimeric protein assembly with cofactors and metal ions and (3) analysis of consensus interaction motifs within the ensembles and identification of conserved determinants of these interactions. Currently the database contains 585 proteins with 650 peptide-binding sites. http://peptisite.ucsd.edu/ link allows searching for the sites of interest and interactive visualization of the ensembles using the ActiveICM web-browser plugin. This structural database for protein-peptide interactions enables understanding of structural principles of these interactions and may assist the development of an efficient peptide docking benchmark.
Collapse
Affiliation(s)
- Chayan Acharya
- UCSD, Skaggs School of Pharmacy and Pharmaceutical Sciences, La Jolla, CA 92093, USA
| | - Irina Kufareva
- UCSD, Skaggs School of Pharmacy and Pharmaceutical Sciences, La Jolla, CA 92093, USA
| | - Andrey V Ilatovskiy
- UCSD, Skaggs School of Pharmacy and Pharmaceutical Sciences, La Jolla, CA 92093, USA; Division of Molecular and Radiation Biophysics, Petersburg Nuclear Physics Institute, Gatchina 188300, Russia; Research and Education Center "Biophysics", PNPI and St. Petersburg State Polytechnical University, St. Petersburg 195251, Russia
| | - Ruben Abagyan
- UCSD, Skaggs School of Pharmacy and Pharmaceutical Sciences, La Jolla, CA 92093, USA.
| |
Collapse
|
35
|
Yu H, Huang H. Engineering proteins for thermostability through rigidifying flexible sites. Biotechnol Adv 2013; 32:308-15. [PMID: 24211474 DOI: 10.1016/j.biotechadv.2013.10.012] [Citation(s) in RCA: 163] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2013] [Revised: 09/04/2013] [Accepted: 10/29/2013] [Indexed: 01/06/2023]
Abstract
Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.
Collapse
Affiliation(s)
- Haoran Yu
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Key Laboratory of Systems Bioengineering, Ministry of Education, Tianjin University, Tianjin 300072, China
| | - He Huang
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Key Laboratory of Systems Bioengineering, Ministry of Education, Tianjin University, Tianjin 300072, China.
| |
Collapse
|
36
|
Fox N, Streinu I. Towards accurate modeling of noncovalent interactions for protein rigidity analysis. BMC Bioinformatics 2013; 14 Suppl 18:S3. [PMID: 24564209 PMCID: PMC3817810 DOI: 10.1186/1471-2105-14-s18-s3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Background Protein rigidity analysis is an efficient computational method for extracting flexibility information from static, X-ray crystallography protein data. Atoms and bonds are modeled as a mechanical structure and analyzed with a fast graph-based algorithm, producing a decomposition of the flexible molecule into interconnected rigid clusters. The result depends critically on noncovalent atomic interactions, primarily on how hydrogen bonds and hydrophobic interactions are computed and modeled. Ongoing research points to the stringent need for benchmarking rigidity analysis software systems, towards the goal of increasing their accuracy and validating their results, either against each other and against biologically relevant (functional) parameters. We propose two new methods for modeling hydrogen bonds and hydrophobic interactions that more accurately reflect a mechanical model, without being computationally more intensive. We evaluate them using a novel scoring method, based on the B-cubed score from the information retrieval literature, which measures how well two cluster decompositions match. Results To evaluate the modeling accuracy of KINARI, our pebble-game rigidity analysis system, we use a benchmark data set of 20 proteins, each with multiple distinct conformations deposited in the Protein Data Bank. Cluster decompositions for them were previously determined with the RigidFinder method from Gerstein's lab and validated against experimental data. When KINARI's default tuning parameters are used, an improvement of the B-cubed score over a crude baseline is observed in 30% of this data. With our new modeling options, improvements were observed in over 70% of the proteins in this data set. We investigate the sensitivity of the cluster decomposition score with case studies on pyruvate phosphate dikinase and calmodulin. Conclusion To substantially improve the accuracy of protein rigidity analysis systems, thorough benchmarking must be performed on all current systems and future extensions. We have measured the gain in performance by comparing different modeling methods for noncovalent interactions. We showed that new criteria for modeling hydrogen bonds and hydrophobic interactions can significantly improve the results. The two new methods proposed here have been implemented and made publicly available in the current version of KINARI (v1.3), together with the benchmarking tools, which can be downloaded from our software's website, http://kinari.cs.umass.edu.
Collapse
|
37
|
Abstract
Protein rigidity and flexibility can be analyzed accurately and efficiently using the program floppy inclusion and rigid substructure topography (FIRST). Previous studies using FIRST were designed to analyze the rigidity and flexibility of proteins using a single static (snapshot) structure. It is however well known that proteins can undergo spontaneous sub-molecular unfolding and refolding, or conformational dynamics, even under conditions that strongly favor a well-defined native structure. These (local) unfolding events result in a large number of conformers that differ from each other very slightly. In this context, proteins are better represented as a thermodynamic ensemble of 'native-like' structures, and not just as a single static low-energy structure. Working with this notion, we introduce a novel FIRST-based approach for predicting rigidity/flexibility of the protein ensemble by (i) averaging the hydrogen bonding strengths from the entire ensemble and (ii) by refining the mathematical model of hydrogen bonds. Furthermore, we combine our FIRST-ensemble rigidity predictions with the ensemble solvent accessibility data of the backbone amides and propose a novel computational method which uses both rigidity and solvent accessibility for predicting hydrogen-deuterium exchange (HDX). To validate our predictions, we report a novel site specific HDX experiment which characterizes the native structural ensemble of Acylphosphatase from hyperthermophile Sulfolobus solfataricus (Sso AcP). The sub-structural conformational dynamics that is observed by HDX data, is closely matched with the FIRST-ensemble rigidity predictions, which could not be attained using the traditional single 'snapshot' rigidity analysis. Moreover, the computational predictions of regions that are protected from HDX and those that undergo exchange are in very good agreement with the experimental HDX profile of Sso AcP.
Collapse
Affiliation(s)
- Adnan Sljoka
- Department of Mathematics and Statistics, York University, 4700 Keele Street, Toronto, M3J 1P3, Canada
| | | |
Collapse
|
38
|
Pfleger C, Gohlke H. Efficient and robust analysis of biomacromolecular flexibility using ensembles of network topologies based on fuzzy noncovalent constraints. Structure 2013; 21:1725-34. [PMID: 23994009 DOI: 10.1016/j.str.2013.07.012] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2013] [Revised: 07/04/2013] [Accepted: 07/17/2013] [Indexed: 11/19/2022]
Abstract
We describe an approach (ENT(FNC)) for performing rigidity analyses of biomacromolecules on ensembles of network topologies (ENT) generated from a single input structure. The ENT is based on fuzzy noncovalent constraints, which considers thermal fluctuations of biomacromolecules without actually sampling conformations. Definitions for fuzzy noncovalent constraints were derived from persistency data from molecular dynamics (MD) simulations. A very good agreement between local flexibility and rigidity characteristics from ENT(FNC) and MD simulations-generated ensembles is found. Regarding global characteristics, convincing results were obtained when relative thermostabilities of citrate synthase and lipase A structures were computed. The ENT(FNC) approach significantly improves the robustness of rigidity analyses, is highly efficient, and does not require a protein-specific parameterization. Its low computational demand makes it especially valuable for the analysis of large data sets, e.g., for data-driven protein engineering.
Collapse
Affiliation(s)
- Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität, 40225 Düsseldorf, Germany
| | | |
Collapse
|
39
|
Krüger DM, Rathi PC, Pfleger C, Gohlke H. CNA web server: rigidity theory-based thermal unfolding simulations of proteins for linking structure, (thermo-)stability, and function. Nucleic Acids Res 2013; 41:W340-8. [PMID: 23609541 PMCID: PMC3692064 DOI: 10.1093/nar/gkt292] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
The Constraint Network Analysis (CNA) web server provides a user-friendly interface to the CNA approach developed in our laboratory for linking results from rigidity analyses to biologically relevant characteristics of a biomolecular structure. The CNA web server provides a refined modeling of thermal unfolding simulations that considers the temperature dependence of hydrophobic tethers and computes a set of global and local indices for quantifying biomacromolecular stability. From the global indices, phase transition points are identified where the structure switches from a rigid to a floppy state; these phase transition points can be related to a protein’s (thermo-)stability. Structural weak spots (unfolding nuclei) are automatically identified, too; this knowledge can be exploited in data-driven protein engineering. The local indices are useful in linking flexibility and function and to understand the impact of ligand binding on protein flexibility. The CNA web server robustly handles small-molecule ligands in general. To overcome issues of sensitivity with respect to the input structure, the CNA web server allows performing two ensemble-based variants of thermal unfolding simulations. The web server output is provided as raw data, plots and/or Jmol representations. The CNA web server, accessible at http://cpclab.uni-duesseldorf.de/cna or http://www.cnanalysis.de, is free and open to all users with no login requirement.
Collapse
Affiliation(s)
- Dennis M Krüger
- Computational Pharmaceutical Chemistry Group, Department of Mathematics and Natural Sciences, Heinrich-Heine-University, 40225 Düsseldorf, Germany
| | | | | | | |
Collapse
|
40
|
Pfleger C, Rathi PC, Klein DL, Radestock S, Gohlke H. Constraint Network Analysis (CNA): a Python software package for efficiently linking biomacromolecular structure, flexibility, (thermo-)stability, and function. J Chem Inf Model 2013; 53:1007-15. [PMID: 23517329 DOI: 10.1021/ci400044m] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
For deriving maximal advantage from information on biomacromolecular flexibility and rigidity, results from rigidity analyses must be linked to biologically relevant characteristics of a structure. Here, we describe the Python-based software package Constraint Network Analysis (CNA) developed for this task. CNA functions as a front- and backend to the graph-based rigidity analysis software FIRST. CNA goes beyond the mere identification of flexible and rigid regions in a biomacromolecule in that it (I) provides a refined modeling of thermal unfolding simulations that also considers the temperature-dependence of hydrophobic tethers, (II) allows performing rigidity analyses on ensembles of network topologies, either generated from structural ensembles or by using the concept of fuzzy noncovalent constraints, and (III) computes a set of global and local indices for quantifying biomacromolecular stability. This leads to more robust results from rigidity analyses and extends the application domain of rigidity analyses in that phase transition points ("melting points") and unfolding nuclei ("structural weak spots") are determined automatically. Furthermore, CNA robustly handles small-molecule ligands in general. Such advancements are important for applying rigidity analysis to data-driven protein engineering and for estimating the influence of ligand molecules on biomacromolecular stability. CNA maintains the efficiency of FIRST such that the analysis of a single protein structure takes a few seconds for systems of several hundred residues on a single core. These features make CNA an interesting tool for linking biomacromolecular structure, flexibility, (thermo-)stability, and function. CNA is available from http://cpclab.uni-duesseldorf.de/software for nonprofit organizations.
Collapse
Affiliation(s)
- Christopher Pfleger
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University, Universitätsstr. 1, 40225, Düsseldorf, Germany
| | | | | | | | | |
Collapse
|
41
|
Li X, Ye L, Wang X, Wang X, Liu H, Qian X, Zhu Y, Yu H. Molecular docking, molecular dynamics simulation, and structure-based 3D-QSAR studies on estrogenic activity of hydroxylated polychlorinated biphenyls. THE SCIENCE OF THE TOTAL ENVIRONMENT 2012; 441:230-238. [PMID: 23137989 DOI: 10.1016/j.scitotenv.2012.08.072] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2012] [Revised: 08/21/2012] [Accepted: 08/24/2012] [Indexed: 06/01/2023]
Abstract
Hydroxylated polychlorinated biphenyls (HO-PCBs), major metabolites of PCBs, have been reported to present agonist or antagonist interactions with estrogen receptor α (ERα) and induce ER-mediated responses. In this work, a multistep framework combining molecular docking, molecular dynamics (MD) simulations, and structure-based three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were performed to explore the influence of structural features on the estrogenic activities of HO-PCBs, and to investigate the molecular mechanism of ERα-ligand interactions. The CoMSIA (comparative molecular similarity indices analysis) model was developed from the conformations obtained from molecular docking. The model exhibited statistically significant results as the cross-validated correlation coefficient q² was 0.648, the non-cross-validated correlation coefficient r² was 0.968, and the external predictive correlation coefficient r(pred)² was 0.625. The key amino acid residues were identified by molecular docking, and the detailed binding modes of the compounds with different activities were determined by MD simulations. The binding free energies correlated well with the experimental activity. An energetic analysis, MM-GBSA energy decomposition, revealed that the van der Waals interaction was the major driving force for the binding of compounds to ERα. The hydrogen bond interactions between the ligands and residue His524 help to stabilize the conformation of ligands at the binding pocket. These results are expected to be beneficial to predict estrogenic activities of other HO-PCB congeners and helpful for understanding the binding mechanism of HO-PCBs and ERα.
Collapse
Affiliation(s)
- Xiaolin Li
- State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment, Nanjing University, Nanjing 210046, PR China
| | | | | | | | | | | | | | | |
Collapse
|
42
|
Baussand J, Kleinjung J. Specific Conformational States of Ras GTPase upon Effector Binding. J Chem Theory Comput 2012; 9:738-749. [PMID: 23316125 PMCID: PMC3541755 DOI: 10.1021/ct3007265] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2012] [Indexed: 12/31/2022]
Abstract
![]()
To uncover the structural and dynamical determinants
involved in
the highly specific binding of Ras GTPase to its effectors, the conformational
states of Ras in uncomplexed form and complexed to the downstream
effectors Byr2, PI3Kγ, PLCε, and RalGDS were investigated
using molecular dynamics and cross-comparison of the trajectories.
The subtle changes in the dynamics and conformations of Ras upon effector
binding require an analysis that targets local changes independent
of global motions. Using a structural alphabet, a computational procedure
is proposed to quantify local conformational changes. Positions detected
by this approach were characterized as either specific for a particular
effector, specific for an effector domain type, or as effector unspecific.
A set of nine structurally connected residues (Ras residues 5–8,
32–35, 39–42, 55–59, 73–78, and 161–165),
which link the effector binding site to the distant C-terminus, changed
dynamics upon effector binding, indicating a potential effector-unspecific
signaling route within the Ras structure. Additional conformational
changes were detected along the N-terminus of the central β-sheet.
Besides the Ras residues at the effector interface (e.g., D33, E37,
D38, and Y40), which adopt effector-specific local conformations,
the binding signal propagates from the interface to distant hot-spot
residues, in particular to Y5 and D57. The results of this study reveal
possible conformational mechanisms for the stabilization of the active
state of Ras upon downstream effector binding and for the structural
determinants responsible for effector specificity.
Collapse
Affiliation(s)
- Julie Baussand
- Division of Mathematical Biology, MRC National Institute for Medical Research, The Ridgeway, Mill Hill, London NW7 1AA, United Kingdom
| | | |
Collapse
|
43
|
Pfleger C, Radestock S, Schmidt E, Gohlke H. Global and local indices for characterizing biomolecular flexibility and rigidity. J Comput Chem 2012; 34:220-33. [PMID: 23007873 DOI: 10.1002/jcc.23122] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2012] [Revised: 08/26/2012] [Accepted: 08/28/2012] [Indexed: 12/24/2022]
Abstract
Understanding flexibility and rigidity characteristics of biomolecules is a prerequisite for understanding biomolecular structural stability and function. Computational methods have been implemented that directly characterize biomolecular flexibility and rigidity by constraint network analysis. For deriving maximal advantage from these analyses, their results need to be linked to biologically relevant characteristics of a structure. Such links are provided by global and local measures ("indices") of biomolecular flexibility and rigidity. To date, more than 14 indices are available with sometimes overlapping or only vague definitions. We present concise definitions of these indices, analyze the relation between, and the scope and limitations of them, and compare their informative value. For this, we probe the structural stability of the calcium binding protein α-lactalbumin as a showcase, both in the "ground state" and after perturbing the system by changing the network topology. In addition, we introduce three indices for the first time that extend the application domain of flexibility and rigidity analyses. The results allow us to provide guidelines for future studies suggesting which of these indices could best be used for analyzing, understanding, and quantifying structural features that are important for biomolecular stability and function. Finally, we make suggestions for proper index notations in future studies to prevent the misinterpretation and to facilitate the comparison of results obtained from flexibility and rigidity analyses.
Collapse
Affiliation(s)
- Christopher Pfleger
- Department of Mathematics and Natural Sciences, Institute for Pharmaceutical and Medicinal Chemistry, Heinrich-Heine-University, Düsseldorf, Germany
| | | | | | | |
Collapse
|
44
|
Haglund E, Sułkowska JI, He Z, Feng GS, Jennings PA, Onuchic JN. The unique cysteine knot regulates the pleotropic hormone leptin. PLoS One 2012; 7:e45654. [PMID: 23029163 PMCID: PMC3454405 DOI: 10.1371/journal.pone.0045654] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2012] [Accepted: 08/20/2012] [Indexed: 11/30/2022] Open
Abstract
Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.
Collapse
Affiliation(s)
- Ellinor Haglund
- Department of Chemistry and Biochemistry and Center for theoretical Biological Physics (CTBP), University of California San Diego, La Jolla, California, United States of America
| | - Joanna I. Sułkowska
- Department of Chemistry and Biochemistry and Center for theoretical Biological Physics (CTBP), University of California San Diego, La Jolla, California, United States of America
| | - Zhao He
- Department of Pathology; School of Medicine and Molecular Biology Section, Division of Biological Sciences, University of California San Diego, La Jolla, California, United States of America
| | - Gen-Sheng Feng
- Department of Pathology; School of Medicine and Molecular Biology Section, Division of Biological Sciences, University of California San Diego, La Jolla, California, United States of America
| | - Patricia A. Jennings
- Department of Chemistry and Biochemistry and Center for theoretical Biological Physics (CTBP), University of California San Diego, La Jolla, California, United States of America
| | - José N. Onuchic
- Center for Theoretical Biological physics and Department of Physics and Astronomy, Chemistry, and Biochemistry and Cell Biology, Rice University, Houston, Texas, United States of America
| |
Collapse
|
45
|
Rader AJ, Yennamalli RM, Harter AK, Sen TZ. A rigid network of long-range contacts increases thermostability in a mutant endoglucanase. J Biomol Struct Dyn 2012; 30:628-37. [PMID: 22731517 DOI: 10.1080/07391102.2012.689696] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35 V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35 V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.
Collapse
Affiliation(s)
- A J Rader
- Department of Physics, Indiana University-Purdue University Indianapolis, Indianapolis, IN 46202, USA.
| | | | | | | |
Collapse
|
46
|
Heal JW, Jimenez-Roldan JE, Wells SA, Freedman RB, Römer RA. Inhibition of HIV-1 protease: the rigidity perspective. ACTA ACUST UNITED AC 2012; 28:350-7. [PMID: 22291339 DOI: 10.1093/bioinformatics/btr683] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
MOTIVATION HIV-1 protease is a key drug target due to its role in the life cycle of the HIV-1 virus. Rigidity analysis using the software First is a computationally inexpensive method for inferring functional information from protein crystal structures. We evaluate the rigidity of 206 high-resolution (2 Å or better) X-ray crystal structures of HIV-1 protease and compare the effects of different inhibitors binding to the enzyme. RESULTS Inhibitor binding has little effect on the overall rigidity of the protein homodimer, including the rigidity of the active site. The principal effect of inhibitor binding on rigidity is to constrain the flexibility of the β-hairpin flaps, which move to allow access to the active site of the enzyme. We show that commercially available antiviral drugs which target HIV-1 protease can be divided into two classes, those which significantly affect flap rigidity and those which do not. The non-peptidic inhibitor tipranavir is distinctive in its consistently strong effect on flap rigidity. CONTACT jack.heal@warwick.ac.uk; r.roemer@warwick.ac.uk SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- J W Heal
- MOAC Doctoral Training Centre, University of Warwick, Coventry CV4 7AL, UK.
| | | | | | | | | |
Collapse
|
47
|
Grimme D, González-ruiz D, Gohlke* H. Computational Strategies and Challenges for Targeting Protein–Protein Interactions with Small Molecules. PHYSICO-CHEMICAL AND COMPUTATIONAL APPROACHES TO DRUG DISCOVERY 2012. [DOI: 10.1039/9781849735377-00319] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
48
|
Maftei M, Tian X, Manea M, Exner TE, Schwanzar D, von Arnim CAF, Przybylski M. Interaction structure of the complex between neuroprotective factor humanin and Alzheimer's β-amyloid peptide revealed by affinity mass spectrometry and molecular modeling. J Pept Sci 2012; 18:373-82. [PMID: 22522311 DOI: 10.1002/psc.2404] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2011] [Revised: 01/19/2012] [Accepted: 01/20/2012] [Indexed: 02/02/2023]
Abstract
Humanin (HN) is a linear 24-aa peptide recently detected in human Alzheimer's disease (AD) brain. HN specifically inhibits neuronal cell death in vitro induced by ß-amyloid (Aß) peptides and by amyloid precursor protein and its gene mutations in familial AD, thereby representing a potential therapeutic lead structure for AD; however, its molecular mechanism of action is not well understood. We report here the identification of the binding epitopes between HN and Aß(1-40) and characterization of the interaction structure through a molecular modeling study. Wild-type HN and HN-sequence mutations were synthesized by SPPS and the HPLC-purified peptides characterized by MALDI-MS. The interaction epitopes between HN and Aß(1-40) were identified by affinity-MS using proteolytic epitope excision and extraction, followed by elution and mass spectrometric characterization of the affinity-bound peptides. The affinity-MS analyses revealed HN(5-15) as the epitope sequence of HN, whereas Aß(17-28) was identified as the Aß interaction epitope. The epitopes and binding sites were ascertained by ELISA of the complex of HN peptides with immobilized Aß(1-40) and by ELISA with Aß(1-40) and Aß-partial sequences as ligands to immobilized HN. The specificity and affinity of the HN-Aß interaction were characterized by direct ESI-MS of the HN-Aß(1-40) complex and by bioaffinity analysis using a surface acoustic wave biosensor, providing a K(D) of the complex of 610 nm. A molecular dynamics simulation of the HN-Aß(1-40) complex was consistent with the binding specificity and shielding effects of the HN and Aß interaction epitopes. These results indicate a specific strong association of HN and Aß(1-40) polypeptide and provide a molecular basis for understanding the neuroprotective function of HN.
Collapse
Affiliation(s)
- Madalina Maftei
- Laboratory of Analytical Chemistry and Biopolymer Structure Analysis, Department of Chemistry, University of Konstanz, 78457, Konstanz, Germany
| | | | | | | | | | | | | |
Collapse
|
49
|
Lemmon G, Kaufmann K, Meiler J. Prediction of HIV-1 protease/inhibitor affinity using RosettaLigand. Chem Biol Drug Des 2012; 79:888-96. [PMID: 22321894 DOI: 10.1111/j.1747-0285.2012.01356.x] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Predicting HIV-1 protease/inhibitor binding affinity as the difference between the free energy of the inhibitor bound and unbound state remains difficult as the unbound state exists as an ensemble of conformations with various degrees of flap opening. We improve computational prediction of protease/inhibitor affinity by invoking the hypothesis that the free energy of the unbound state while difficult to predict is less sensitive to mutation. Thereby the HIV-1 protease/inhibitor binding affinity can be approximated with the free energy of the bound state alone. Bound state free energy can be predicted from comparative models of HIV-1 protease mutant/inhibitor complexes. Absolute binding energies are predicted with R = 0.71 and SE = 5.91 kJ/mol. Changes in binding free energy upon mutation can be predicted with R = 0.85 and SE = 4.49 kJ/mol. Resistance mutations that lower inhibitor binding affinity can thereby be recognized early in HIV-1 protease inhibitor development.
Collapse
Affiliation(s)
- Gordon Lemmon
- Department of Chemistry, Center for Structural Biology, Institute of Chemical Biology, Vanderbilt University, Nashville, TN 37232, USA
| | | | | |
Collapse
|
50
|
Jiang W, Ghosh D. Motion and flexibility in human cytochrome p450 aromatase. PLoS One 2012; 7:e32565. [PMID: 22384274 PMCID: PMC3288111 DOI: 10.1371/journal.pone.0032565] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2011] [Accepted: 01/31/2012] [Indexed: 11/26/2022] Open
Abstract
The crystal structures of human placental aromatase in complex with the substrate androstenedione and exemestane have revealed an androgen-specific active site and the structural basis for higher order organization. However, X-ray structures do not provide accounts of movements due to short-range fluctuations, ligand binding and protein-protein association. In this work, we conduct normal mode analysis (NMA) revealing the intrinsic fluctuations of aromatase, deduce the internal modes in membrane-free and membrane-integrated monomers as well as the intermolecular modes in oligomers, and propose a quaternary organization for the endoplasmic reticulum (ER) membrane integration. Dynamics of the crystallographic oligomers from NMA is found to be in agreement with the isotropic thermal factors from the X-ray analysis. Calculations of the root mean square fluctuations of the C-alpha atoms from their equilibrium positions confirm that the rigid-core structure of aromatase is intrinsic regardless of the changes in steroid binding interactions, and that aromatase self-association does not deteriorate the rigidity of the catalytic cleft. Furthermore, NMA on membrane-integrated aromatase shows that the internal modes in all likelihood contribute to breathing of the active site access channel. The collective intermolecular hinge bending and twisting modes provide the flexibility in the quaternary association necessary for membrane integration of the aromatase oligomers. Taken together, fluctuations of the active site, the access channel, and the heme-proximal cavity, and a dynamic quaternary organization could all be essential components of the functional aromatase in its role as an ER membrane-embedded steroidogenic enzyme.
Collapse
Affiliation(s)
| | - Debashis Ghosh
- Department of Pharmacology, State University of New York Upstate Medical University, Syracuse, New York, United States of America
- * E-mail:
| |
Collapse
|