1
|
Shadman H, Ziebarth JD, Gallops CE, Luo R, Li Z, Chen HF, Wang Y. Map conformational landscapes of intrinsically disordered proteins with polymer physics quantities. Biophys J 2024; 123:1253-1263. [PMID: 38615193 PMCID: PMC11140466 DOI: 10.1016/j.bpj.2024.04.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 02/20/2024] [Accepted: 04/10/2024] [Indexed: 04/15/2024] Open
Abstract
Disordered proteins are conformationally flexible proteins that are biologically important and have been implicated in devastating diseases such as Alzheimer's disease and cancer. Unlike stably folded structured proteins, disordered proteins sample a range of different conformations that needs to be accounted for. Here, we treat disordered proteins as polymer chains, and compute a dimensionless quantity called instantaneous shape ratio (Rs), as Rs = Ree2/Rg2, where Ree is end-to-end distance and Rg is radius of gyration. Extended protein conformations tend to have high Ree compared with Rg, and thus have high Rs values, whereas compact conformations have smaller Rs values. We use a scatter plot of Rs (representing shape) against Rg (representing size) as a simple map of conformational landscapes. We first examine the conformational landscape of simple polymer models such as Random Walk, Self-Avoiding Walk, and Gaussian Walk (GW), and we notice that all protein/polymer maps lie within the boundaries of the GW map. We thus use the GW map as a reference and, to assess conformational diversity, we compute the fraction of the GW conformations (fC) covered by each protein/polymer. Disordered proteins all have high fC scores, consistent with their disordered nature. Each disordered protein accesses a different region of the reference map, revealing differences in their conformational ensembles. We additionally examine the conformational maps of the nonviral gene delivery vector polyethyleneimine at various protonation states, and find that they resemble disordered proteins, with coverage of the reference map decreasing with increasing protonation state, indicating decreasing conformational diversity. We propose that our method of combining Rs and Rg in a scatter plot generates a simple, meaningful map of the conformational landscape of a disordered protein, which in turn can be used to assess conformational diversity of disordered proteins.
Collapse
Affiliation(s)
- Hossain Shadman
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Jesse D Ziebarth
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Caleb E Gallops
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Ray Luo
- Chemical and Materials Physics Graduate Program, Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California
| | - Zhengxin Li
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | - Yongmei Wang
- Department of Chemistry, The University of Memphis, Memphis, Tennessee.
| |
Collapse
|
2
|
LeMaster DM, Bashir Q, Hernández G. Propagation of conformational instability in FK506-binding protein FKBP12. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:140990. [PMID: 38142946 PMCID: PMC10939819 DOI: 10.1016/j.bbapap.2023.140990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 12/12/2023] [Accepted: 12/18/2023] [Indexed: 12/26/2023]
Abstract
FKBP12 is the archetype of the FK506 binding domains that define the family of FKBP proteins which participate in the regulation of various distinct physiological signaling processes. As the drugs FK506 and rapamycin inhibit many of these FKBP proteins, there is need to develop therapeutics which exhibit selectivity within this family. The long β4-β5 loop of the FKBP domain is known to regulate transcriptional activity for the steroid hormone receptors and appears to participate in regulating calcium channel activity for the cardiac and skeletal muscle ryanodine receptors. The β4-β5 loop of FKBP12 has been shown to undergo extensive conformational dynamics, and here we report hydrogen exchange measurements for a series of mutational variants in that loop which indicate deviations from a two-state kinetics for those dynamics. In addition to a previously characterized local transition near the tip of this loop, evidence is presented for a second site of conformational dynamics in the stem of this loop. These mutation-dependent hydrogen exchange effects extend beyond the β4-β5 loop, primarily by disrupting the hydrogen bond between the Gly 58 amide and the Tyr 80 carbonyl oxygen which links the two halves of the structural rim that surrounds the active site cleft. Mutationally-induced opening of the cleft between Gly 58 and Tyr 80 not only modulates the global stability of the protein, it promotes a conformational transition in the distant β2-β3a hairpin that modulates the binding affinity for a FKBP51-selective inhibitor previously designed to exploit a localized conformational transition at the homologous site.
Collapse
Affiliation(s)
- David M LeMaster
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America
| | - Qamar Bashir
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America
| | - Griselda Hernández
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America.
| |
Collapse
|
3
|
Anderson JS, LeMaster DM, Hernández G. Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states. J Biol Chem 2023; 299:105159. [PMID: 37579948 PMCID: PMC10514456 DOI: 10.1016/j.jbc.2023.105159] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 08/08/2023] [Accepted: 08/10/2023] [Indexed: 08/16/2023] Open
Abstract
Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, New York, USA
| | - David M LeMaster
- New York State Department of Health, Wadsworth Center, Albany, New York, USA
| | - Griselda Hernández
- New York State Department of Health, Wadsworth Center, Albany, New York, USA.
| |
Collapse
|
4
|
Schulze CJ, Seamon KJ, Zhao Y, Yang YC, Cregg J, Kim D, Tomlinson A, Choy TJ, Wang Z, Sang B, Pourfarjam Y, Lucas J, Cuevas-Navarro A, Santos CA, Vides A, Li C, Marquez A, Zhong M, Vemulapalli V, Weller C, Gould A, Whalen DM, Salvador A, Milin A, Saldajeno-Concar M, Dinglasan N, Chen A, Evans J, Knox JE, Koltun ES, Singh M, Nichols R, Wildes D, Gill AL, Smith JAM, Lito P. Chemical remodeling of a cellular chaperone to target the active state of mutant KRAS. Science 2023; 381:794-799. [PMID: 37590355 PMCID: PMC10474815 DOI: 10.1126/science.adg9652] [Citation(s) in RCA: 37] [Impact Index Per Article: 37.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2023] [Accepted: 06/28/2023] [Indexed: 08/19/2023]
Abstract
The discovery of small-molecule inhibitors requires suitable binding pockets on protein surfaces. Proteins that lack this feature are considered undruggable and require innovative strategies for therapeutic targeting. KRAS is the most frequently activated oncogene in cancer, and the active state of mutant KRAS is such a recalcitrant target. We designed a natural product-inspired small molecule that remodels the surface of cyclophilin A (CYPA) to create a neomorphic interface with high affinity and selectivity for the active state of KRASG12C (in which glycine-12 is mutated to cysteine). The resulting CYPA:drug:KRASG12C tricomplex inactivated oncogenic signaling and led to tumor regressions in multiple human cancer models. This inhibitory strategy can be used to target additional KRAS mutants and other undruggable cancer drivers. Tricomplex inhibitors that selectively target active KRASG12C or multiple RAS mutants are in clinical trials now (NCT05462717 and NCT05379985).
Collapse
Affiliation(s)
| | - Kyle J. Seamon
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Yulei Zhao
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Yu C. Yang
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Jim Cregg
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Dongsung Kim
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Aidan Tomlinson
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Tiffany J. Choy
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Zhican Wang
- Department of Non-clinical Development and Clinical Pharmacology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Ben Sang
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Yasin Pourfarjam
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Jessica Lucas
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Antonio Cuevas-Navarro
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Carlos Ayala Santos
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Alberto Vides
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Chuanchuan Li
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
| | - Abby Marquez
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Mengqi Zhong
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | | | - Caroline Weller
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Andrea Gould
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Daniel M. Whalen
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Anthony Salvador
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Anthony Milin
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Mae Saldajeno-Concar
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Nuntana Dinglasan
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Anqi Chen
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Jim Evans
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - John E. Knox
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Elena S. Koltun
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Mallika Singh
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Robert Nichols
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - David Wildes
- Department of Biology, Revolution Medicines, Inc., Redwood City, CA, 94063
| | - Adrian L. Gill
- Department of Discovery Chemistry, Revolution Medicines, Inc., Redwood City, CA, 94063
| | | | - Piro Lito
- Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer, New York, NY, 10065
- Department of Medicine, Memorial Sloan Kettering Cancer Center, New York, NY, 10065
- Department of Medicine, Weill Cornell Medical College, New York, NY, 10065
| |
Collapse
|
5
|
Jeong KB, Ryu M, Kim JS, Kim M, Yoo J, Chung M, Oh S, Jo G, Lee SG, Kim HM, Lee MK, Chi SW. Single-molecule fingerprinting of protein-drug interaction using a funneled biological nanopore. Nat Commun 2023; 14:1461. [PMID: 37015934 PMCID: PMC10073129 DOI: 10.1038/s41467-023-37098-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 03/01/2023] [Indexed: 04/06/2023] Open
Abstract
In drug discovery, efficient screening of protein-drug interactions (PDIs) is hampered by the limitations of current biophysical approaches. Here, we develop a biological nanopore sensor for single-molecule detection of proteins and PDIs using the pore-forming toxin YaxAB. Using this YaxAB nanopore, we demonstrate label-free, single-molecule detection of interactions between the anticancer Bcl-xL protein and small-molecule drugs as well as the Bak-BH3 peptide. The long funnel-shaped structure and nanofluidic characteristics of the YaxAB nanopore enable the electro-osmotic trapping of diverse folded proteins and high-resolution monitoring of PDIs. Distinctive nanopore event distributions observed in the two-dimensional (ΔI/Io-versus-IN) plot illustrate the ability of the YaxAB nanopore to discriminate individual small-molecule drugs bound to Bcl-xL from non-binders. Taken together, our results present the YaxAB nanopore as a robust platform for label-free, ultrasensitive, single-molecule detection of PDIs, opening up a possibility for low-cost, highly efficient drug discovery against diverse drug targets.
Collapse
Affiliation(s)
- Ki-Baek Jeong
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea
| | - Minju Ryu
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea
| | - Jin-Sik Kim
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea
| | - Minsoo Kim
- Department of Physics, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea
| | - Jejoong Yoo
- Department of Physics, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea
| | - Minji Chung
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Sohee Oh
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Gyunghee Jo
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
| | - Seong-Gyu Lee
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
| | - Ho Min Kim
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
- Graduate School of Medical Science and Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Mi-Kyung Lee
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea.
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea.
| | - Seung-Wook Chi
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea.
- School of Pharmacy, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea.
| |
Collapse
|
6
|
Sora V, Laspiur AO, Degn K, Arnaudi M, Utichi M, Beltrame L, De Menezes D, Orlandi M, Stoltze UK, Rigina O, Sackett PW, Wadt K, Schmiegelow K, Tiberti M, Papaleo E. RosettaDDGPrediction for high-throughput mutational scans: From stability to binding. Protein Sci 2023; 32:e4527. [PMID: 36461907 PMCID: PMC9795540 DOI: 10.1002/pro.4527] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 11/25/2022] [Accepted: 11/25/2022] [Indexed: 12/05/2022]
Abstract
Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein-protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction.
Collapse
Affiliation(s)
- Valentina Sora
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Adrian Otamendi Laspiur
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Kristine Degn
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Matteo Arnaudi
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Mattia Utichi
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Ludovica Beltrame
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Dayana De Menezes
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Matteo Orlandi
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Ulrik Kristoffer Stoltze
- Department of Clinical GeneticsCopenhagen University Hospital RigshospitaletCopenhagenDenmark
- Department of Pediatrics and Adolescent MedicineUniversity Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Olga Rigina
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Peter Wad Sackett
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Karin Wadt
- Department of Clinical GeneticsCopenhagen University Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Kjeld Schmiegelow
- Department of Pediatrics and Adolescent MedicineUniversity Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Matteo Tiberti
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
| | - Elena Papaleo
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| |
Collapse
|
7
|
Carius Y, Ries F, Gries K, Trentmann O, Lancaster CRD, Willmund F. Structural features of chloroplast trigger factor determined at 2.6 Å resolution. Acta Crystallogr D Struct Biol 2022; 78:1259-1272. [PMID: 36189745 PMCID: PMC9527764 DOI: 10.1107/s2059798322009068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Accepted: 09/11/2022] [Indexed: 11/21/2022] Open
Abstract
The folding of newly synthesized polypeptides requires the coordinated action of molecular chaperones. Prokaryotic cells and the chloroplasts of plant cells possess the ribosome-associated chaperone trigger factor, which binds nascent polypeptides at their exit stage from the ribosomal tunnel. The structure of bacterial trigger factor has been well characterized and it has a dragon-shaped conformation, with flexible domains responsible for ribosome binding, peptidyl-prolyl cis-trans isomerization (PPIase) activity and substrate protein binding. Chloroplast trigger-factor sequences have diversified from those of their bacterial orthologs and their molecular mechanism in plant organelles has been little investigated to date. Here, the crystal structure of the plastidic trigger factor from the green alga Chlamydomonas reinhardtii is presented at 2.6 Å resolution. Due to the high intramolecular flexibility of the protein, diffraction to this resolution was only achieved using a protein that lacked the N-terminal ribosome-binding domain. The eukaryotic trigger factor from C. reinhardtii exhibits a comparable dragon-shaped conformation to its bacterial counterpart. However, the C-terminal chaperone domain displays distinct charge distributions, with altered positioning of the helical arms and a specifically altered charge distribution along the surface responsible for substrate binding. While the PPIase domain shows a highly conserved structure compared with other PPIases, its rather weak activity and an unusual orientation towards the C-terminal domain points to specific adaptations of eukaryotic trigger factor for function in chloroplasts.
Collapse
Affiliation(s)
- Yvonne Carius
- Department of Structural Biology, Saarland University, Center of Human and Molecular Biology (ZHMB), Faculty of Medicine, Building 60, 66421 Homburg, Germany
| | - Fabian Ries
- Molecular Genetics of Eukaryotes, University of Kaiserslautern, Erwin-Schrödinger-Strasse 70, 67663 Kaiserslautern, Germany
| | - Karin Gries
- Molecular Genetics of Eukaryotes, University of Kaiserslautern, Erwin-Schrödinger-Strasse 70, 67663 Kaiserslautern, Germany
| | - Oliver Trentmann
- Molecular Botany, University of Kaiserslautern, Erwin-Schrödinger-Strasse 70, 67663 Kaiserslautern, Germany
| | - C. Roy D. Lancaster
- Department of Structural Biology, Saarland University, Center of Human and Molecular Biology (ZHMB), Faculty of Medicine, Building 60, 66421 Homburg, Germany
| | - Felix Willmund
- Molecular Genetics of Eukaryotes, University of Kaiserslautern, Erwin-Schrödinger-Strasse 70, 67663 Kaiserslautern, Germany
| |
Collapse
|
8
|
Leveraging Fungal and Human Calcineurin-Inhibitor Structures, Biophysical Data, and Dynamics To Design Selective and Nonimmunosuppressive FK506 Analogs. mBio 2021; 12:e0300021. [PMID: 34809463 PMCID: PMC8609367 DOI: 10.1128/mbio.03000-21] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Calcineurin is a critical enzyme in fungal pathogenesis and antifungal drug tolerance and, therefore, an attractive antifungal target. Current clinically accessible calcineurin inhibitors, such as FK506, are immunosuppressive to humans, so exploiting calcineurin inhibition as an antifungal strategy necessitates fungal specificity in order to avoid inhibiting the human pathway. Harnessing fungal calcineurin-inhibitor crystal structures, we recently developed a less immunosuppressive FK506 analog, APX879, with broad-spectrum antifungal activity and demonstrable efficacy in a murine model of invasive fungal infection. Our overarching goal is to better understand, at a molecular level, the interaction determinants of the human and fungal FK506-binding proteins (FKBP12) required for calcineurin inhibition in order to guide the design of fungus-selective, nonimmunosuppressive FK506 analogs. To this end, we characterized high-resolution structures of the Mucor circinelloides FKBP12 bound to FK506 and of the Aspergillus fumigatus, M. circinelloides, and human FKBP12 proteins bound to the FK506 analog APX879, which exhibits enhanced selectivity for fungal pathogens. Combining structural, genetic, and biophysical methodologies with molecular dynamics simulations, we identify critical variations in these structurally similar FKBP12-ligand complexes. The work presented here, aimed at the rational design of more effective calcineurin inhibitors, indeed suggests that modifications to the APX879 scaffold centered around the C15, C16, C18, C36, and C37 positions provide the potential to significantly enhance fungal selectivity. IMPORTANCE Invasive fungal infections are a leading cause of death in the immunocompromised patient population. The rise in drug resistance to current antifungals highlights the urgent need to develop more efficacious and highly selective agents. Numerous investigations of major fungal pathogens have confirmed the critical role of the calcineurin pathway for fungal virulence, making it an attractive target for antifungal development. Although FK506 inhibits calcineurin, it is immunosuppressive in humans and cannot be used as an antifungal. By combining structural, genetic, biophysical, and in silico methodologies, we pinpoint regions of the FK506 scaffold and a less immunosuppressive analog, APX879, centered around the C15 to C18 and C36 to C37 positions that could be altered with selective extensions and/or deletions to enhance fungal selectivity. This work represents a significant advancement toward realizing calcineurin as a viable target for antifungal drug discovery.
Collapse
|
9
|
Youssef N, Susko E, Bielawski JP. Consequences of Stability-Induced Epistasis for Substitution Rates. Mol Biol Evol 2020; 37:3131-3148. [DOI: 10.1093/molbev/msaa151] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
AbstractDo interactions between residues in a protein (i.e., epistasis) significantly alter evolutionary dynamics? If so, what consequences might they have on inference from traditional codon substitution models which assume site-independence for the sake of computational tractability? To investigate the effects of epistasis on substitution rates, we employed a mechanistic mutation-selection model in conjunction with a fitness framework derived from protein stability. We refer to this as the stability-informed site-dependent (S-SD) model and developed a new stability-informed site-independent (S-SI) model that captures the average effect of stability constraints on individual sites of a protein. Comparison of S-SI and S-SD offers a novel and direct method for investigating the consequences of stability-induced epistasis on protein evolution. We developed S-SI and S-SD models for three natural proteins and showed that they generate sequences consistent with real alignments. Our analyses revealed that epistasis tends to increase substitution rates compared with the rates under site-independent evolution. We then assessed the epistatic sensitivity of individual site and discovered a counterintuitive effect: Highly connected sites were less influenced by epistasis relative to exposed sites. Lastly, we show that, despite the unrealistic assumptions, traditional models perform comparably well in the presence and absence of epistasis and provide reasonable summaries of average selection intensities. We conclude that epistatic models are critical to understanding protein evolutionary dynamics, but epistasis might not be required for reasonable inference of selection pressure when averaging over time and sites.
Collapse
Affiliation(s)
- Noor Youssef
- Department of Biology, Dalhousie University, Halifax, Nova Scotia, Canada
- Centre for Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Edward Susko
- Centre for Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, Nova Scotia, Canada
- Department of Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Joseph P Bielawski
- Department of Biology, Dalhousie University, Halifax, Nova Scotia, Canada
- Centre for Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, Nova Scotia, Canada
- Department of Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia, Canada
| |
Collapse
|
10
|
Setny P. Conserved internal hydration motifs in protein kinases. Proteins 2020; 88:1578-1591. [DOI: 10.1002/prot.25977] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Revised: 05/18/2020] [Accepted: 07/05/2020] [Indexed: 01/14/2023]
Affiliation(s)
- Piotr Setny
- Centre of New Technologies University of Warsaw Warsaw Poland
| |
Collapse
|
11
|
Oshima H, Re S, Sugita Y. Prediction of Protein–Ligand Binding Pose and Affinity Using the gREST+FEP Method. J Chem Inf Model 2020; 60:5382-5394. [DOI: 10.1021/acs.jcim.0c00338] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Hiraku Oshima
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Suyong Re
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
12
|
Song D, Liu H, Luo R, Chen HF. Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins. J Chem Inf Model 2020; 60:2257-2267. [PMID: 32227937 PMCID: PMC10449432 DOI: 10.1021/acs.jcim.0c00059] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The need for accurate and efficient force fields for modeling 3D structures of macrobiomolecules and in particular intrinsically disordered proteins (IDPs) has increased with recent findings to associate IDPs and human diseases. However, most conventional protein force fields and recent IDP-specific force fields are limited in reproducing accurate structural features of IDPs. Here, we present an environmental specific precise force field (ESFF1) based on CMAP corrections of 71 different sequence environments to improve the accuracy and efficiency of MD simulation for both IDPs and folded proteins. MD simulations of 84 different short peptides, IDPs, and structured proteins show that ESFF1 can accurately reproduce spectroscopic properties for different peptides and proteins whether they are disordered or ordered. The successful ab initio folding of five fast-folding proteins further supports the reliability of ESFF1. The extensive analysis documented here shows that ESFF1 is able to achieve a reasonable balance between ordered and disordered states in protein simulations.
Collapse
Affiliation(s)
- Dong Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Molecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, California 92697-3900, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
- Shanghai Center for Bioinformation Technology, Shanghai 200235, China
| |
Collapse
|
13
|
Nerattini F, Figliuzzi M, Cardelli C, Tubiana L, Bianco V, Dellago C, Coluzza I. Identification of Protein Functional Regions. Chemphyschem 2020; 21:335-347. [PMID: 31944517 DOI: 10.1002/cphc.201900898] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2019] [Revised: 11/01/2019] [Indexed: 11/12/2022]
Abstract
Protein sequence stores the information relative to both functionality and stability, thus making it difficult to disentangle the two contributions. However, the identification of critical residues for function and stability has important implications for the mapping of the proteome interactions, as well as for many pharmaceutical applications, e. g. the identification of ligand binding regions for targeted pharmaceutical protein design. In this work, we propose a computational method to identify critical residues for protein functionality and stability and to further categorise them in strictly functional, structural and intermediate. We evaluate single site conservation and use Direct Coupling Analysis (DCA) to identify co-evolved residues both in natural and artificial evolution processes. We reproduce artificial evolution using protein design and base our approach on the hypothesis that artificial evolution in the absence of any functional constraint would exclusively lead to site conservation and co-evolution events of the structural type. Conversely, natural evolution intrinsically embeds both functional and structural information. By comparing the lists of conserved and co-evolved residues, outcomes of the analysis on natural and artificial evolution, we identify the functional residues without the need of any a priori knowledge of the biological role of the analysed protein.
Collapse
Affiliation(s)
- Francesca Nerattini
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Matteo Figliuzzi
- Sorbonne Universites, UPMC, Institut de Biologie Paris-Seine, CNRS, Laboratoire de Biologie Computationnelle et Quantitative UMR, 7238, Paris, France
| | - Chiara Cardelli
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Luca Tubiana
- Physics Department, Universitá degli studi di Trento, via Sommarive 14, 38123, Trento, IT
| | - Valentino Bianco
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria.,Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid, 28040, Spain
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria
| | - Ivan Coluzza
- CIC biomaGUNE, Paseo Miramon 182, 20014 San Sebastian, Spain, and IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain
| |
Collapse
|
14
|
Alavilli H, Lee H, Park M, Yun DJ, Lee BH. Enhanced multiple stress tolerance in Arabidopsis by overexpression of the polar moss peptidyl prolyl isomerase FKBP12 gene. PLANT CELL REPORTS 2018; 37:453-465. [PMID: 29247292 DOI: 10.1007/s00299-017-2242-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/29/2017] [Accepted: 12/05/2017] [Indexed: 06/07/2023]
Abstract
PaFKBP12 overexpression in Arabidopsis resulted in stress tolerance to heat, ABA, drought, and salt stress, in addition to growth promotion under normal conditions. Polytrichastrum alpinum (alpine haircap moss) is one of polar organisms that can withstand the severe conditions of the Antarctic. In this study, we report the isolation of a peptidyl prolyl isomerase FKBP12 gene (PaFKBP12) from P. alpinum collected in the Antarctic and its functional implications in development and stress responses in plants. In P. alpinum, PaFKBP12 expression was induced by heat and ABA. Overexpression of PaFKBP12 in Arabidopsis increased the plant size, which appeared to result from increased rates of cell cycle. Under heat stress conditions, PaFKBP12-overexpressing lines (PaFKBP12-OE) showed better growth and survival than the wild type. PaFKBP12-OE also showed higher root elongation rates, better shoot growth and enhanced survival at higher concentrations of ABA in comparison to the wild type. In addition, PaFKBP12-OE were more tolerant to drought and salt stress than the wild type. All these phenotypes were accompanied with higher induction of the stress responsive genes in PaFKBP12-OE than in the wild type. Taken together, our findings revealed important functions of PaFKBP12 in plant development and abiotic stress responses.
Collapse
Affiliation(s)
| | - Hyoungseok Lee
- Unit of Polar Genomics, Korea Polar Research Institute, Incheon, 21990, South Korea
| | - Mira Park
- Department of Life Science, Sogang University, Seoul, 04107, South Korea
- Unit of Polar Genomics, Korea Polar Research Institute, Incheon, 21990, South Korea
| | - Dae-Jin Yun
- Department of Biomedical Science and Engineering, Konkuk University, Seoul, 05029, South Korea
| | - Byeong-Ha Lee
- Department of Life Science, Sogang University, Seoul, 04107, South Korea.
| |
Collapse
|
15
|
Solomentsev G, Diehl C, Akke M. Conformational Entropy of FK506 Binding to FKBP12 Determined by Nuclear Magnetic Resonance Relaxation and Molecular Dynamics Simulations. Biochemistry 2018; 57:1451-1461. [DOI: 10.1021/acs.biochem.7b01256] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Gleb Solomentsev
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Carl Diehl
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mikael Akke
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
16
|
Spyrakis F, Ahmed MH, Bayden AS, Cozzini P, Mozzarelli A, Kellogg GE. The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery. J Med Chem 2017; 60:6781-6827. [PMID: 28475332 DOI: 10.1021/acs.jmedchem.7b00057] [Citation(s) in RCA: 98] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.
Collapse
Affiliation(s)
- Francesca Spyrakis
- Dipartimento di Scienza e Tecnologia del Farmaco, Università degli Studi di Torino , Via Pietro Giuria 9, 10125 Torino, Italy
| | - Mostafa H Ahmed
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| | - Alexander S Bayden
- CMD Bioscience , 5 Science Park, New Haven, Connecticut 06511, United States
| | - Pietro Cozzini
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Modellistica Molecolare, Università degli Studi di Parma , Parco Area delle Scienze 59/A, 43121 Parma, Italy
| | - Andrea Mozzarelli
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Biochimica, Università degli Studi di Parma , Parco Area delle Scienze 23/A, 43121 Parma, Italy.,Istituto di Biofisica, Consiglio Nazionale delle Ricerche , Via Moruzzi 1, 56124 Pisa, Italy
| | - Glen E Kellogg
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| |
Collapse
|
17
|
Weininger U, Modig K, Geitner AJ, Schmidpeter PAM, Koch JR, Akke M. Dynamics of Aromatic Side Chains in the Active Site of FKBP12. Biochemistry 2016; 56:334-343. [DOI: 10.1021/acs.biochem.6b01157] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Ulrich Weininger
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
- Institute
of Physics, Biophysics, Martin-Luther-University Halle-Wittenberg, D-06120 Halle (Saale), Germany
| | - Kristofer Modig
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
| | - Anne-Juliane Geitner
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Philipp A. M. Schmidpeter
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Johanna R. Koch
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Mikael Akke
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
18
|
Ubiquitylation Directly Induces Fold Destabilization of Proteins. Sci Rep 2016; 6:39453. [PMID: 27991582 PMCID: PMC5172356 DOI: 10.1038/srep39453] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2016] [Accepted: 11/23/2016] [Indexed: 12/24/2022] Open
Abstract
Ubiquitin is a common post-translational modifier and its conjugation is a key signal for proteolysis by the proteasome. Because the molecular mass of ubiquitin is larger than that of other modifiers such as phosphate, acetyl, or methyl groups, ubiquitylation not only influences biochemical signaling, but also may exert physical effects on its substrate proteins by increasing molecular volume and altering shape anisotropy. Here we show that ubiquitylation destabilizes the fold of two proteins, FKBP12 and FABP4, and that elongation of the conjugated ubiquitin chains further enhances this destabilization effect. Moreover, NMR relaxation analysis shows that ubiquitylation induces characteristic structural fluctuations in the backbone of both proteins. These results suggest that the ubiquitylation-driven structural fluctuations lead to fold destabilization of its substrate proteins. Thus, physical destabilization by ubiquitylation may facilitate protein degradation by the proteasome.
Collapse
|
19
|
LeMaster DM, Hernandez G. Conformational Dynamics in FKBP Domains: Relevance to Molecular Signaling and Drug Design. Curr Mol Pharmacol 2016; 9:5-26. [PMID: 25986571 DOI: 10.2174/1874467208666150519113146] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Revised: 02/25/2015] [Accepted: 05/17/2015] [Indexed: 01/05/2023]
Abstract
Among the 22 FKBP domains in the human genome, FKBP12.6 and the first FKBP domains (FK1) of FKBP51 and FKBP52 are evolutionarily and structurally most similar to the archetypical FKBP12. As such, the development of inhibitors with selectivity among these four FKBP domains poses a significant challenge for structure-based design. The pleiotropic effects of these FKBP domains in a range of signaling processes such as the regulation of ryanodine receptor calcium channels by FKBP12 and FKBP12.6 and steroid receptor regulation by the FK1 domains of FKBP51 and FKBP52 amply justify the efforts to develop selective therapies. In contrast to their close structural similarities, these four FKBP domains exhibit a substantial diversity in their conformational flexibility. A number of distinct conformational transitions have been characterized for FKBP12 spanning timeframes from 20 s to 10 ns and in each case these dynamics have been shown to markedly differ from the conformational behavior for one or more of the other three FKBP domains. Protein flexibilitybased inhibitor design could draw upon the transitions that are significantly populated in only one of the targeted proteins. Both the similarities and differences among these four proteins valuably inform the understanding of how dynamical effects propagate across the FKBP domains as well as potentially how such intramolecular transitions might couple to the larger scale transitions that are central to the signaling complexes in which these FKBP domains function.
Collapse
Affiliation(s)
| | - Griselda Hernandez
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, New York, 12201, USA; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, New York, 12201, USA.
| |
Collapse
|
20
|
Hydration of proteins and nucleic acids: Advances in experiment and theory. A review. Biochim Biophys Acta Gen Subj 2016; 1860:1821-35. [PMID: 27241846 DOI: 10.1016/j.bbagen.2016.05.036] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2015] [Revised: 05/20/2016] [Accepted: 05/26/2016] [Indexed: 11/21/2022]
Abstract
BACKGROUND Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics. SCOPE OF REVIEW This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties. MAJOR CONCLUSIONS In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration. GENERAL SIGNIFICANCE An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.
Collapse
|
21
|
Barman A, Smitherman C, Souffrant M, Gadda G, Hamelberg D. Conserved Hydration Sites in Pin1 Reveal a Distinctive Water Recognition Motif in Proteins. J Chem Inf Model 2015; 56:139-47. [DOI: 10.1021/acs.jcim.5b00560] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Arghya Barman
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Crystal Smitherman
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Michael Souffrant
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Giovanni Gadda
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Departments
of Chemistry and ‡Biology and the §Centers for Diagnostics and Therapeutics and ∥Biotechnology
and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
22
|
Setny P. Prediction of Water Binding to Protein Hydration Sites with a Discrete, Semiexplicit Solvent Model. J Chem Theory Comput 2015; 11:5961-72. [PMID: 26642995 DOI: 10.1021/acs.jctc.5b00839] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Buried water molecules are ubiquitous in protein structures and are found at the interface of most protein-ligand complexes. Determining their distribution and thermodynamic effect is a challenging yet important task, of great of practical value for the modeling of biomolecular structures and their interactions. In this study, we present a novel method aimed at the prediction of buried water molecules in protein structures and estimation of their binding free energies. It is based on a semiexplicit, discrete solvation model, which we previously introduced in the context of small molecule hydration. The method is applicable to all macromolecular structures described by a standard all-atom force field, and predicts complete solvent distribution within a single run with modest computational cost. We demonstrate that it indicates positions of buried hydration sites, including those filled by more than one water molecule, and accurately differentiates them from sterically accessible to water but void regions. The obtained estimates of water binding free energies are in fair agreement with reference results determined with the double decoupling method.
Collapse
Affiliation(s)
- Piotr Setny
- Centre of New Technologies, University of Warsaw , Banacha 2c, 02-097 Warsaw, Poland
| |
Collapse
|
23
|
Chakraborty D, Taly A, Sterpone F. Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins. J Phys Chem B 2015; 119:12760-70. [PMID: 26335353 DOI: 10.1021/acs.jpcb.5b05791] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.
Collapse
Affiliation(s)
- Debashree Chakraborty
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Antoine Taly
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
24
|
Statistical survey of the buried waters in the Protein Data Bank. Amino Acids 2015; 48:193-202. [PMID: 26315961 DOI: 10.1007/s00726-015-2064-4] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2015] [Accepted: 07/29/2015] [Indexed: 12/18/2022]
Abstract
The structures of buried water molecules were studied in an ensemble of high-quality and non-redundant protein crystal structures. Buried water molecules were clustered and classified in lake-like clusters, which are completely isolated from the bulk solvent, and bay-like clusters, which are in contact with the bulk solvent through a surface water molecule. Buried water molecules are extremely common: lake-like clusters are found in 89 % of the protein crystal structures and bay-like clusters in 93 %. Clusters with only one water molecule are much more common than larger clusters. Both cluster types incline to be surrounded by loop residues, and to a minor extent by residues in extended secondary structure. Helical residues on the contrary do not tend to surround clusters of buried water molecules. One buried water molecule is found every 30-50 amino acid residues, depending on the secondary structures that are more abundant in the protein. Both main- and side-chain atoms are in contact with buried waters; they form four hydrogen bonds with the first water and 1-1.5 additional hydrogen bond for each additional water in the cluster. Consequently, buried water molecules appear to be firmly packed and rigid like the protein atoms. In this regard, it is remarkable to observe that prolines often surround water molecules buried in the protein interior. Interestingly, clusters of buried water molecules tend to be just beneath the protein surface. Moreover, water molecules tend to form a one-dimensional wire rather than more compact arrangements. This agrees with recent evidence of the mechanisms of solvent exchange between internal cavities and bulk solvent.
Collapse
|
25
|
Hörner M, Kaufmann B, Cotugno G, Wiedtke E, Büning H, Grimm D, Weber W. A chemical switch for controlling viral infectivity. Chem Commun (Camb) 2015; 50:10319-22. [PMID: 25058661 DOI: 10.1039/c4cc03292f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Chemically triggered molecular switches for controlling the fate and function of biological systems are fundamental to the emergence of synthetic biology and the development of biomedical applications. We here present the first chemically triggered switch for controlling the infectivity of adeno-associated viral (AAV) vectors.
Collapse
Affiliation(s)
- Maximilian Hörner
- Faculty of Biology, University of Freiburg, Schänzlestrasse 1, 79104 Freiburg, Germany.
| | | | | | | | | | | | | |
Collapse
|
26
|
Chuntonov L, Pazos IM, Ma J, Gai F. Kinetics of exchange between zero-, one-, and two-hydrogen-bonded states of methyl and ethyl acetate in methanol. J Phys Chem B 2015; 119:4512-20. [PMID: 25738661 DOI: 10.1021/acs.jpcb.5b00745] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
It has recently been shown that the ester carbonyl stretching vibration can be used as a sensitive probe of local electrostatic field in molecular systems. To further characterize this vibrational probe and extend its potential applications, we studied the kinetics of chemical exchange between differently hydrogen-bonded (H-bonded) ester carbonyl groups of methyl acetate (MA) and ethyl acetate (EA) in methanol. We found that, while both MA and EA can form zero, one, or two H-bonds with the solvent, the population of the 2hb state in MA is significantly smaller than that in EA. Using a combination of linear and nonlinear infrared measurements and numerical simulations, we further determined the rate constants for the exchange between these differently H-bonded states. We found that for MA the chemical exchange reaction between the two dominant states (i.e., 0hb and 1hb states) has a relaxation rate constant of 0.14 ps(-1), whereas for EA the three-state chemical exchange reaction occurs in a predominantly sequential manner with the following relaxation rate constants: 0.11 ps(-1) for exchange between 0hb and 1hb states and 0.12 ps(-1) for exchange between 1hb and 2hb states.
Collapse
Affiliation(s)
- Lev Chuntonov
- †Ultrafast Optical Processes Laboratory and ‡Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Ileana M Pazos
- †Ultrafast Optical Processes Laboratory and ‡Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Jianqiang Ma
- †Ultrafast Optical Processes Laboratory and ‡Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Feng Gai
- †Ultrafast Optical Processes Laboratory and ‡Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
27
|
Ladani ST, Souffrant MG, Barman A, Hamelberg D. Computational perspective and evaluation of plausible catalytic mechanisms of peptidyl-prolyl cis-trans isomerases. Biochim Biophys Acta Gen Subj 2015; 1850:1994-2004. [PMID: 25585011 DOI: 10.1016/j.bbagen.2014.12.023] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2014] [Revised: 12/23/2014] [Accepted: 12/29/2014] [Indexed: 01/16/2023]
Abstract
BACKGROUND Peptidyl prolyl cis-trans isomerization of the protein backbone is involved in the regulation of many biological processes. Cis-trans isomerization is notoriously slow and is catalyzed by a family of cis-trans peptidyl prolyl isomerases (PPIases) that have been implicated in many diseases. A general consensus on how these enzymes speed up prolyl isomerization has not been reached after decades of both experimental and computational studies. SCOPE OF REVIEW Computational studies carried out to understand the catalytic mechanism of the prototypical FK506 binding protein 12, Cyclophilin A and peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) are reviewed. A summary and an evaluation of the implications of the proposed mechanisms from computational studies are presented. MAJOR CONCLUSIONS The analysis of computational studies and evaluation of the proposed mechanisms provide a general consensus and a better understanding of PPIase catalysis. The speedup of the rate of peptidyl-prolyl isomerization by PPIases can be best described by a catalytic mechanism in which the substrate in transition state configuration is stabilized. The enzymes preferentially bind the transition state configuration of the substrate relative to the cis conformation, which in most cases is bound better than the trans conformation of the substrate. Stabilization of the transition state configuration of the substrate leads to a lower free energy barrier and a faster rate of isomerization when compared to the uncatalyzed isomerization reaction. GENERAL SIGNIFICANCE Fully understanding the catalytic mechanism of PPIases has broad implications for drug design, elucidation of the molecular basis of many diseases, protein engineering, and enzyme catalysis in general. This article is part of a Special Issue entitled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.
Collapse
Affiliation(s)
- Safieh Tork Ladani
- Department of Chemistry and the Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30302-3965, USA
| | - Michael G Souffrant
- Department of Chemistry and the Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30302-3965, USA
| | - Arghya Barman
- Department of Chemistry and the Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30302-3965, USA
| | - Donald Hamelberg
- Department of Chemistry and the Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30302-3965, USA.
| |
Collapse
|
28
|
Structure of the rabbit ryanodine receptor RyR1 at near-atomic resolution. Nature 2014; 517:50-55. [PMID: 25517095 PMCID: PMC4338550 DOI: 10.1038/nature14063] [Citation(s) in RCA: 327] [Impact Index Per Article: 32.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Accepted: 11/12/2014] [Indexed: 12/12/2022]
Abstract
The ryanodine receptors (RyRs) are high-conductance intracellular Ca2+ channels that play a pivotal role in the excitation-contraction coupling of skeletal and cardiac muscles. RyRs are the largest known ion channels, with a homotetrameric organization and approximately 5000 residues in each protomer. Here we report the structure of the rabbit RyR1 in complex with its modulator FKBP12 at an overall resolution of 3.8 Å, determined by single-particle electron cryo-microscopy. Three previously uncharacterized domains, named Central, Handle, and Helical domains, display the armadillo repeat fold. These domains, together with the amino-terminal domain, constitute a network of superhelical scaffold for binding and propagation of conformational changes. The channel domain exhibits the voltage-gated ion channel superfamily fold with distinct features. A negative charge-enriched hairpin loop connecting S5 and the pore helix is positioned above the entrance to the selectivity filter vestibule. The four elongated S6 segments form a right-handed helical bundle that closes the pore at the cytoplasmic border of the membrane. Allosteric regulation of the pore by the cytoplasmic domains is mediated through extensive interactions between the Central domains and the channel domain. These structural features explain high ion conductance by RyRs and the long-range allosteric regulation of channel activities.
Collapse
|
29
|
Structure-affinity properties of a high-affinity ligand of FKBP12 studied by molecular simulations of a binding intermediate. PLoS One 2014; 9:e114610. [PMID: 25502559 PMCID: PMC4264844 DOI: 10.1371/journal.pone.0114610] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2014] [Accepted: 11/10/2014] [Indexed: 11/19/2022] Open
Abstract
With a view to explaining the structure-affinity properties of the ligands of the protein FKBP12, we characterized a binding intermediate state between this protein and a high-affinity ligand. Indeed, the nature and extent of the intermolecular contacts developed in such a species may play a role on its stability and, hence, on the overall association rate. To find the binding intermediate, a molecular simulation protocol was used to unbind the ligand by gradually decreasing the biasing forces introduced. The intermediate was subsequently refined with 17 independent stochastic boundary molecular dynamics simulations that provide a consistent picture of the intermediate state. In this state, the core region of the ligand remains stable, notably because of the two anchoring oxygen atoms that correspond to recurrent motifs found in all FKBP12 ligand core structures. Besides, the non-core regions participate in numerous transient intermolecular and intramolecular contacts. The dynamic aspect of most of the contacts seems important both for the ligand to retain at least a part of its configurational entropy and for avoiding a trapped state along the binding pathway. Since the transient and anchoring contacts contribute to increasing the stability of the intermediate, as a corollary, the dissociation rate constant of this intermediate should be decreased, resulting in an increase of the affinity constant . The present results support our previous conclusions and provide a coherent rationale for explaining the prevalence in high-affinity ligands of (i) the two oxygen atoms found in carbonyl or sulfonyl groups of dissimilar core structures and of (ii) symmetric or pseudo-symmetric mobile groups of atoms found as non-core moieties. Another interesting aspect of the intermediate is the distortion of the flexible 80 s loop of the protein, mainly in its tip region, that promotes the accessibility to the bound state.
Collapse
|
30
|
Rahaman O, Kalimeri M, Melchionna S, Hénin J, Sterpone F. Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains. J Phys Chem B 2014; 119:8939-49. [PMID: 25317828 DOI: 10.1021/jp507571u] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.
Collapse
Affiliation(s)
- Obaidur Rahaman
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Maria Kalimeri
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Simone Melchionna
- ‡CNR-IPCF, Consiglio Nazionale delle Ricerche, Physics Dept., Univ. La Sapienza, P.le A. Moro 2, 00185, Rome, Italy
| | - Jérôme Hénin
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Fabio Sterpone
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| |
Collapse
|
31
|
Anderson JS, Mustafi SM, Hernández G, LeMaster DM. Statistical allosteric coupling to the active site indole ring flip equilibria in the FK506-binding domain. Biophys Chem 2014; 192:41-8. [PMID: 25016286 DOI: 10.1016/j.bpc.2014.06.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2014] [Accepted: 06/16/2014] [Indexed: 10/25/2022]
Abstract
In solution, the Trp 59 indole ring at the base of the active site cleft in the FKBP domain protein FKBP12 is rotated by ~90° at a population level of 20%, relative to its canonical crystallographic orientation. NMR measurements on the homologous FK1 domains of human FKBP51 and FKBP52 indicate no observable indole ring flip conformation, while the V101I variant of FKBP12 decreases the population having a perpendicular indole orientation by 10-fold. A set of three parallel 400 ns CHARMM27 molecular simulations for both wild type FKBP12 and the V101I variant examined how this ring flip might be energetically coupled to a transition of the Glu 60 sidechain which interacts with the backbone of the 50's loop located ~12 Å from the indole nitrogen. Analysis of the transition matrix for the local dynamics of the Glu 60 sidechain, the Trp 59 sidechain, and of the structurally interposed α-helix hydrogen bonding pattern yielded a statistical allosteric coupling of 10 kJ/mol with negligible concerted dynamical coupling for the transitions of the two sidechains.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, NY 12308, United States
| | - Sourajit M Mustafi
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States
| | - Griselda Hernández
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, NY 12201, United States
| | - David M LeMaster
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, NY 12201, United States.
| |
Collapse
|
32
|
Structural basis of conformational transitions in the active site and 80's loop in the FK506-binding protein FKBP12. Biochem J 2014; 458:525-36. [PMID: 24405377 PMCID: PMC3940039 DOI: 10.1042/bj20131429] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The extensive set of NMR doublings exhibited by the immunophilin FKBP12 (FK506-binding protein 12) arose from a slow transition to the cis-peptide configuration at Gly89 near the tip of the 80′s loop, the site for numerous protein-recognition interactions for both FKBP12 and other FKBP domain proteins. The 80′s loop also exhibited linebroadening, indicative of microsecond to millisecond conformational dynamics, but only in the trans-peptide state. The G89A variant shifted the trans–cis peptide equilibrium from 88:12 to 33:67, whereas a proline residue substitution induced fully the cis-peptide configuration. The 80′s loop conformation in the G89P crystal structure at 1.50 Å resolution differed from wild-type FKBP12 primarily at residues 88, 89 and 90, and it closely resembled that reported for FKBP52. Structure-based chemical-shift predictions indicated that the microsecond to millisecond dynamics in the 80′s loop probably arose from a concerted main chain (ψ88 and ϕ89) torsion angle transition. The indole side chain of Trp59 at the base of the active-site cleft was reoriented ~90o and the adjacent backbone was shifted in the G89P crystal structure. NOE analysis of wild-type FKBP12 demonstrated that this indole populates the perpendicular orientation at 20%. The 15N relaxation analysis was consistent with the indole reorientation occurring in the nanosecond timeframe. Recollection of the G89P crystal data at 1.20 Å resolution revealed a weaker wild-type-like orientation for the indole ring. Differences in the residues that underlie the Trp59 indole ring and altered interactions linking the 50′s loop to the active site suggested that reorientation of this ring may be disfavoured in the other six members of the FKBP domain family that bear this active-site tryptophan residue. Extensive resonance doubling arises from a cis–trans peptide transition at Gly89, whereas linebroadening appears due to a concerted shift in the neighbouring torsion angles. The active site Trp59 ring adopts a perpendicular orientation at a population of 20%.
Collapse
|
33
|
Chen H, Mustafi SM, LeMaster DM, Li Z, Héroux A, Li H, Hernández G. Crystal structure and conformational flexibility of the unligated FK506-binding protein FKBP12.6. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2014; 70:636-46. [PMID: 24598733 PMCID: PMC3949516 DOI: 10.1107/s1399004713032112] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/09/2013] [Accepted: 11/25/2013] [Indexed: 12/15/2022]
Abstract
The primary known physiological function of FKBP12.6 involves its role in regulating the RyR2 isoform of ryanodine receptor Ca(2+) channels in cardiac muscle, pancreatic β islets and the central nervous system. With only a single previously reported X-ray structure of FKBP12.6, bound to the immunosuppressant rapamycin, structural inferences for this protein have been drawn from the more extensive studies of the homologous FKBP12. X-ray structures at 1.70 and 1.90 Å resolution from P2₁ and P3₁21 crystal forms are reported for an unligated cysteine-free variant of FKBP12.6 which exhibit a notable diversity of conformations. In one monomer from the P3₁21 crystal form, the aromatic ring of Phe59 at the base of the active site is rotated perpendicular to its typical orientation, generating a steric conflict for the immunosuppressant-binding mode. The peptide unit linking Gly89 and Val90 at the tip of the protein-recognition `80s loop' is flipped in the P2₁ crystal form. Unlike the >30 reported FKBP12 structures, the backbone conformation of this loop closely follows that of the first FKBP domain of FKBP51. The NMR resonances for 21 backbone amides of FKBP12.6 are doubled, corresponding to a slow conformational transition centered near the tip of the 80s loop, as recently reported for 31 amides of FKBP12. The comparative absence of doubling for residues along the opposite face of the active-site pocket in FKBP12.6 may in part reflect attenuated structural coupling owing to increased conformational plasticity around the Phe59 ring.
Collapse
Affiliation(s)
- Hui Chen
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - Sourajit M. Mustafi
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - David M. LeMaster
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| | - Zhong Li
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - Annie Héroux
- Department of Biology, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Hongmin Li
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| | - Griselda Hernández
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| |
Collapse
|
34
|
Xu Y, Leitner DM. Vibrational Energy Flow through the Green Fluorescent Protein–Water Interface: Communication Maps and Thermal Boundary Conductance. J Phys Chem B 2014; 118:7818-26. [DOI: 10.1021/jp412141z] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Yao Xu
- Department
of Chemistry and Chemical Physics Program, University of Nevada, Reno, Nevada 89557, United States
- Physikalische
Chemie II, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - David M. Leitner
- Department
of Chemistry and Chemical Physics Program, University of Nevada, Reno, Nevada 89557, United States
| |
Collapse
|
35
|
Abstract
The 1H-15N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80′s loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50′s loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~105-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80′s loop. Line-broadening was also decreased for the residues in the α-helix and 50′s loop, whereas line-broadening in the 40′s loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40′s loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.
Collapse
|
36
|
Shi Y, Xia Z, Zhang J, Best R, Wu C, Ponder JW, Ren P. The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins. J Chem Theory Comput 2013; 9:4046-4063. [PMID: 24163642 DOI: 10.1021/ct4003702] [Citation(s) in RCA: 437] [Impact Index Per Article: 39.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.
Collapse
Affiliation(s)
- Yue Shi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| | | | | | | | | | | | | |
Collapse
|
37
|
Mustafi SM, Chen H, Li H, Lemaster DM, Hernández G. Analysing the visible conformational substates of the FK506-binding protein FKBP12. Biochem J 2013; 453:371-380. [PMID: 23688288 DOI: 10.1042./bj20130276] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
The 1H-15N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80's loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50's loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~105-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80's loop. Line-broadening was also decreased for the residues in the α-helix and 50's loop, whereas line-broadening in the 40's loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40's loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.
Collapse
Affiliation(s)
- Sourajit M Mustafi
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | | | | | | | | |
Collapse
|
38
|
Demonte D, Drake EJ, Lim KH, Gulick AM, Park S. Structure-based engineering of streptavidin monomer with a reduced biotin dissociation rate. Proteins 2013; 81:1621-33. [PMID: 23670729 DOI: 10.1002/prot.24320] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2013] [Revised: 04/18/2013] [Accepted: 04/22/2013] [Indexed: 11/07/2022]
Abstract
We recently reported the engineering of monomeric streptavidin, mSA, corresponding to one subunit of wild type (wt) streptavidin tetramer. The monomer was designed by homology modeling, in which the streptavidin and rhizavidin sequences were combined to engineer a high affinity binding pocket containing residues from a single subunit only. Although mSA is stable and binds biotin with nanomolar affinity, its fast off rate (koff ) creates practical challenges during applications. We obtained a 1.9 Å crystal structure of mSA bound to biotin to understand their interaction in detail, and used the structure to introduce targeted mutations to improve its binding kinetics. To this end, we compared mSA to shwanavidin, which contains a hydrophobic lid containing F43 in the binding pocket and binds biotin tightly. However, the T48F mutation in mSA, which introduces a comparable hydrophobic lid, only resulted in a modest 20-40% improvement in the measured koff . On the other hand, introducing the S25H mutation near the bicyclic ring of bound biotin increased the dissociation half life (t½ ) from 11 to 83 min at 20°C. Molecular dynamics (MD) simulations suggest that H25 stabilizes the binding loop L3,4 by interacting with A47, and protects key intermolecular hydrogen bonds by limiting solvent entry into the binding pocket. Concurrent T48F or T48W mutation clashes with H25 and partially abrogates the beneficial effects of H25. Taken together, this study suggests that stabilization of the binding loop and solvation of the binding pocket are important determinants of the dissociation kinetics in mSA.
Collapse
Affiliation(s)
- Daniel Demonte
- Department of Chemical and Biological Engineering, University at Buffalo, Buffalo, New York, 14260
| | | | | | | | | |
Collapse
|
39
|
Guan JY, Keizers PHJ, Liu WM, Löhr F, Skinner SP, Heeneman EA, Schwalbe H, Ubbink M, Siegal G. Small-molecule binding sites on proteins established by paramagnetic NMR spectroscopy. J Am Chem Soc 2013; 135:5859-68. [PMID: 23509882 DOI: 10.1021/ja401323m] [Citation(s) in RCA: 80] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Determining the three-dimensional structure of a small molecule-protein complex with weak affinity can be a significant challenge. We present a paramagnetic NMR method to determine intermolecular structure restraints based on pseudocontact shifts (PCSs). Since the ligand must be in fast exchange between free and bound states and the fraction bound can be as low as a few percent, the method is ideal for ligands with high micromolar to millimolar dissociation constants. Paramagnetic tags are attached, one at a time, in a well-defined way via two arms at several sites on the protein surface. The ligand PCSs were measured from simple 1D (1)H spectra and used as docking restraints. An independent confirmation of the complex structure was carried out using intermolecular NOEs. The results show that structures derived from these two approaches are similar. The best results are obtained if the magnetic susceptibility tensors of the tags are known, but it is demonstrated that with two-armed probes, the magnetic susceptibility tensor can be predicted with sufficient accuracy to provide a low-resolution model of the ligand orientation and the location of the binding site in the absence of isotope-labeled protein. This approach can facilitate fragment-based drug discovery in obtaining structural information on the initial fragment hits.
Collapse
Affiliation(s)
- Jia-Ying Guan
- Gorlaeus Laboratories, Leiden Institute of Chemistry, Leiden University, Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | | | | | | | | | | | | | | | | |
Collapse
|
40
|
Anderson JS, Hernández G, LeMaster DM. Assessing the chemical accuracy of protein structures via peptide acidity. Biophys Chem 2012. [PMID: 23182463 DOI: 10.1016/j.bpc.2012.10.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Although the protein native state is a Boltzmann conformational ensemble, practical applications often require a representative model from the most populated region of that distribution. The acidity of the backbone amides, as reflected in hydrogen exchange rates, is exquisitely sensitive to the surrounding charge and dielectric volume distribution. For each of four proteins, three independently determined X-ray structures of differing crystallographic resolution were used to predict exchange for the static solvent-exposed amide hydrogens. The average correlation coefficients range from 0.74 for ubiquitin to 0.93 for Pyrococcus furiosus rubredoxin, reflecting the larger range of experimental exchange rates exhibited by the latter protein. The exchange prediction errors modestly correlate with the crystallographic resolution. MODELLER 9v6-derived homology models at ~60% sequence identity (36% identity for chymotrypsin inhibitor CI2) yielded correlation coefficients that are ~0.1 smaller than for the cognate X-ray structures. The most recently deposited NOE-based ubiquitin structure and the original NMR structure of CI2 fail to provide statistically significant predictions of hydrogen exchange. However, the more recent RECOORD refinement study of CI2 yielded predictions comparable to the X-ray and homology model-based analyses.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, New York 12308, USA
| | | | | |
Collapse
|
41
|
Gollan PJ, Bhave M, Aro EM. The FKBP families of higher plants: Exploring the structures and functions of protein interaction specialists. FEBS Lett 2012; 586:3539-47. [PMID: 22982859 DOI: 10.1016/j.febslet.2012.09.002] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2012] [Revised: 08/31/2012] [Accepted: 09/03/2012] [Indexed: 01/24/2023]
Abstract
The FK506-binding proteins (FKBPs) are known both as the receptors for immunosuppressant drugs and as prolyl isomerase (PPIase) enzymes that catalyse rotation of prolyl bonds. FKBPs are characterised by the inclusion of at least one FK506-binding domain (FKBd), the receptor site for proline and the active site for PPIase catalysis. The FKBPs form large and diverse families in most organisms, with the largest FKBP families occurring in higher plants. Plant FKBPs are molecular chaperones that interact with specific protein partners to regulate a diversity of cellular processes. Recent studies have found that plant FKBPs operate in intricate and coordinated mechanisms for regulating stress response and development processes, and discoveries of new interaction partners expand their cellular influences to gene expression and photosynthetic adaptations. This review presents an examination of the molecular and structural features and functional roles of the higher plant FKBP family within the context of these recent findings, and discusses the significance of domain conservation and variation for the development of a diverse, versatile and complex chaperone family.
Collapse
Affiliation(s)
- Peter J Gollan
- Environment and Biotechnology Centre, Faculty of Life and Social Sciences, Swinburne University of Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia.
| | | | | |
Collapse
|
42
|
Xiol J, Cora E, Koglgruber R, Chuma S, Subramanian S, Hosokawa M, Reuter M, Yang Z, Berninger P, Palencia A, Benes V, Penninger J, Sachidanandam R, Pillai RS. A role for Fkbp6 and the chaperone machinery in piRNA amplification and transposon silencing. Mol Cell 2012; 47:970-9. [PMID: 22902560 DOI: 10.1016/j.molcel.2012.07.019] [Citation(s) in RCA: 106] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2012] [Revised: 06/25/2012] [Accepted: 07/17/2012] [Indexed: 11/29/2022]
Abstract
Epigenetic silencing of transposons by Piwi-interacting RNAs (piRNAs) constitutes an RNA-based genome defense mechanism. Piwi endonuclease action amplifies the piRNA pool by generating new piRNAs from target transcripts by a poorly understood mechanism. Here, we identified mouse Fkbp6 as a factor in this biogenesis pathway delivering piRNAs to the Piwi protein Miwi2. Mice lacking Fkbp6 derepress LINE1 (L1) retrotransposon and display reduced DNA methylation due to deficient nuclear accumulation of Miwi2. Like other cochaperones, Fkbp6 associates with the molecular chaperone Hsp90 via its tetratricopeptide repeat (TPR) domain. Inhibition of the ATP-dependent Hsp90 activity in an insect cell culture model results in the accumulation of short antisense RNAs in Piwi complexes. We identify these to be byproducts of piRNA amplification that accumulate only in nuage-localized Piwi proteins. We propose that the chaperone machinery normally ejects these inhibitory RNAs, allowing turnover of Piwi complexes for their continued participation in piRNA amplification.
Collapse
Affiliation(s)
- Jordi Xiol
- European Molecular Biology Laboratory, 6 Rue Jules Horowitz, BP 181, 38042 Grenoble, France
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Application of Hydration Thermodynamics to the Evaluation of Protein Structures and Protein-Ligand Binding. ENTROPY 2012. [DOI: 10.3390/e14081443] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
44
|
Designing electrostatic interactions in biological systems via charge optimization or combinatorial approaches: insights and challenges with a continuum electrostatic framework. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1252-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
45
|
Molecular interactions of c-ABL mutants in complex with imatinib/nilotinib: a computational study using linear interaction energy (LIE) calculations. J Mol Model 2012; 18:4333-41. [PMID: 22570081 DOI: 10.1007/s00894-012-1436-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2011] [Accepted: 04/17/2012] [Indexed: 10/28/2022]
Abstract
In spite of the effectiveness of Imatinib for chronic myeloid leukemia (CML) treatment, resistance has repeatedly been reported and is associated with point mutations in the BCR-ABL chimeric gene. To overcome this resistance, several inhibitors of BCR-ABL tyrosine kinase activity were developed. In this context, computational simulations have become a powerful tool for understanding drug-protein interactions. Herein, we report a comparative molecular dynamics analysis of the interaction between two tyrosine kinase inhibitors (imatinib or nilotinib) against wild type c-ABL protein and 12 mutants, using the semi-empirical linear interaction energy (LIE) method, to assess the feasibility of this approach for studying resistance against the inhibitory activity of these drugs. In addition, to understand the structural changes that are associated with resistance, we describe the behavior of water molecules that interact simultaneously with specific residues (Glu286, Lys271 and Asp381) of c-ABL (wild type or mutant) and their relationship with drug resistance. Experimental IC50 values for the interaction between imatinib, wild type c-ABL, and 12 mutants were used to obtain the proper LIE coefficients (α, β and γ) to estimate the free energy of the binding of imatinib with wild-type and mutant proteins, and values were extrapolated for the analysis of the nilotinib/c-ABL interaction. Our results indicate that LIE was suitable to predict the superior inhibitory activity of nilotinib and the resistance to inhibition that was observed in c-ABL mutants. Additionally, for c-ABL mutants, the observed number of water molecules being turned over while interacting with amino acids Glu286, Lys271 and Asp381 was associated with resistance to imatinib, resulting in a less effective inhibition of the kinase activity.
Collapse
|
46
|
Gopalakrishnan R, Kozany C, Wang Y, Schneider S, Hoogeland B, Bracher A, Hausch F. Exploration of pipecolate sulfonamides as binders of the FK506-binding proteins 51 and 52. J Med Chem 2012; 55:4123-31. [PMID: 22455398 DOI: 10.1021/jm201747c] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
FK506-binding proteins (FKBP) 51 and 52 are cochaperones that modulate the signal transduction of steroid hormone receptors. Single nucleotide polymorphisms in the gene encoding FKBP51 have been associated with a variety of psychiatric disorders. Rapamycin and FK506 are two macrocyclic natural products, which tightly bind to most FKBP family members, including FKBP51 and FKBP52. A bioisosteric replacement of the α-ketoamide moiety of rapamycin and FK506 with a sulfonamide was envisaged with the retention of the conserved hydrogen bonds. A focused solid support-based synthesis protocol was developed, which led to ligands with submicromolar affinity for FKBP51 and FKBP52. The molecular binding mode for one sulfonamide analogue was confirmed by X-ray crystallography.
Collapse
|
47
|
Afonine PV, Grosse-Kunstleve RW, Echols N, Headd JJ, Moriarty NW, Mustyakimov M, Terwilliger TC, Urzhumtsev A, Zwart PH, Adams PD. Towards automated crystallographic structure refinement with phenix.refine. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2012; 68:352-67. [PMID: 22505256 PMCID: PMC3322595 DOI: 10.1107/s0907444912001308] [Citation(s) in RCA: 4052] [Impact Index Per Article: 337.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/27/2011] [Accepted: 01/11/2012] [Indexed: 11/10/2022]
Abstract
phenix.refine is a program within the PHENIX package that supports crystallographic structure refinement against experimental data with a wide range of upper resolution limits using a large repertoire of model parameterizations. It has several automation features and is also highly flexible. Several hundred parameters enable extensive customizations for complex use cases. Multiple user-defined refinement strategies can be applied to specific parts of the model in a single refinement run. An intuitive graphical user interface is available to guide novice users and to assist advanced users in managing refinement projects. X-ray or neutron diffraction data can be used separately or jointly in refinement. phenix.refine is tightly integrated into the PHENIX suite, where it serves as a critical component in automated model building, final structure refinement, structure validation and deposition to the wwPDB. This paper presents an overview of the major phenix.refine features, with extensive literature references for readers interested in more detailed discussions of the methods.
Collapse
Affiliation(s)
- Pavel V Afonine
- Lawrence Berkeley National Laboratory, One Cyclotron Road, MS64R0121, Berkeley, CA 94720, USA.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
48
|
Chiba S, Harano Y, Roth R, Kinoshita M, Sakurai M. Evaluation of protein-ligand binding free energy focused on its entropic components. J Comput Chem 2011; 33:550-60. [DOI: 10.1002/jcc.22891] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2011] [Revised: 10/18/2011] [Accepted: 10/24/2011] [Indexed: 11/12/2022]
|
49
|
The structure of a Burkholderia pseudomallei immunophilin-inhibitor complex reveals new approaches to antimicrobial development. Biochem J 2011; 437:413-22. [PMID: 21574961 DOI: 10.1042/bj20110345] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Mips (macrophage infectivity potentiators) are a subset of immunophilins associated with virulence in a range of micro-organisms. These proteins possess peptidylprolyl isomerase activity and are inhibited by drugs including rapamycin and tacrolimus. We determined the structure of the Mip homologue [BpML1 (Burkholderia pseudomallei Mip-like protein 1)] from the human pathogen and biowarfare threat B. pseudomallei by NMR and X-ray crystallography. The crystal structure suggests that key catalytic residues in the BpML1 active site have unexpected conformational flexibility consistent with a role in catalysis. The structure further revealed BpML1 binding to a helical peptide, in a manner resembling the physiological interaction of human TGFβRI (transforming growth factor β receptor I) with the human immunophilin FKBP12 (FK506-binding protein 12). Furthermore, the structure of BpML1 bound to the class inhibitor cycloheximide N-ethylethanoate showed that this inhibitor mimics such a helical peptide, in contrast with the extended prolyl-peptide mimicking shown by inhibitors such as tacrolimus. We suggest that Mips, and potentially other bacterial immunophilins, participate in protein-protein interactions in addition to their peptidylprolyl isomerase activity, and that some roles of Mip proteins in virulence are independent of their peptidylprolyl isomerase activity.
Collapse
|
50
|
Targeting FKBP isoforms with small-molecule ligands. Curr Opin Pharmacol 2011; 11:365-71. [PMID: 21803654 DOI: 10.1016/j.coph.2011.04.007] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2011] [Accepted: 04/13/2011] [Indexed: 11/21/2022]
Abstract
The FK506 binding protein (FKBP) family of proteins provide an interesting series of drug targets since different isoforms modulate diverse cellular pathways. There are therapeutic opportunities in the fields of cancer therapy, neurodegenerative conditions and psychiatric disorders. X-ray crystallographic or NMR data are available for eight of fourteen human FKBPs covering ten of the twenty-two different FKBP domains. We have made a detailed sequence and structural comparison of human FKBP domains. These data show that the chemical scaffolds common to the immunosuppressive inhibitors FK506 and rapamycin bind to the most conserved region of the binding site. This observation opens the way to the design of isoform specific inhibitors.
Collapse
|