ChaperISM: improved chaperone binding prediction using position-independent scoring matrices

MOTIVATION Understanding the mechanisms of client protein interaction with Hsp70 chaperones is essential to analyze the complex dynamics in the context of normal or dysregulated metabolism. Because Hsp70 can bind millions of proteins, including key molecules involved in processes of stemness, tumorigenesis, and survival, in silico prediction of Hsp70 interactions has great value in validating possible new clients. Currently, two algorithms are available to predict binding to DnaK - the bacterial Hsp70 - but both are based on amino acid sequence and energy calculations of qualitative information - binders and non-binders. RESULTS We introduce a new algorithm to identify Hsp70 binding sequences in proteins - ChaperISM - a position-independent scoring matrix trained on either qualitative or quantitative chemiluminescence data previously published, which were obtained from the interaction between DnaK and different ligands. Both versions of ChaperISM, qualitative or quantitative, resulted in an improved performance in comparison to other state-of-the-art chaperone binding predictors. AVAILABILITY ChaperISM is implemented in Python version 3. The source code of ChaperISM is freely available for download at https://github.com/BioinfLab/ChaperISM. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Craig M. Ogata,et al.  Structural Analysis of Substrate Binding by the Molecular Chaperone DnaK , 1996, Science.

[2]  Joost J. J. van Durme,et al.  Accurate Prediction of DnaK-Peptide Binding via Homology Modelling and Experimental Data , 2009, PLoS Comput. Biol..

[3]  Bjoern Peters,et al.  Identifying MHC Class I Epitopes by Predicting the TAP Transport Efficiency of Epitope Precursors , 2003, The Journal of Immunology.

[4]  Thomas Lengauer,et al.  DynaPred: A structure and sequence based method for the prediction of MHC class I binding peptide sequences and conformations , 2006, ISMB.

[5]  P. Goloubinoff.  Editorial: The HSP70 Molecular Chaperone Machines , 2017, Front. Mol. Biosci..

[6]  Alessandro Sette,et al.  Generating quantitative models describing the sequence specificity of biological processes with the stabilized matrix method , 2005, BMC Bioinformatics.

[7]  M. Gething,et al.  BiP Binding Sequences in Antibodies (*) , 1995, The Journal of Biological Chemistry.

[8]  J. F. De Mesquita,et al.  Structural and Functional Analysis of Human SOD1 in Amyotrophic Lateral Sclerosis , 2013, PloS one.

[9]  E. Zuiderweg,et al.  Non-canonical Interactions between Heat Shock Cognate Protein 70 (Hsc70) and Bcl2-associated Anthanogene (BAG) Co-Chaperones Are Important for Client Release* , 2016, The Journal of Biological Chemistry.

[10]  P. K. Suresh,et al.  Computational identification and analysis of functional polymorphisms involved in the activation and detoxification genes implicated in endometriosis. , 2014, Gene.

[11]  M. Chawla-Sarkar,et al.  Cell death regulation during influenza A virus infection by matrix (M1) protein: a model of viral control over the cellular survival pathway , 2011, Cell Death and Disease.

[12]  M. Mayer,et al.  Hsp70 chaperones: Cellular functions and molecular mechanism , 2005, Cellular and Molecular Life Sciences.

[13]  L. Hendershot,et al.  The endoplasmic reticulum (ER) chaperone BiP is a master regulator of ER functions: Getting by with a little help from ERdj friends , 2018, The Journal of Biological Chemistry.

[14]  H. Sitte,et al.  A Cytosolic Relay of Heat Shock Proteins HSP70-1A and HSP90β Monitors the Folding Trajectory of the Serotonin Transporter* , 2014, The Journal of Biological Chemistry.

[15]  Roman Kityk,et al.  Insights into the molecular mechanism of allostery in Hsp70s , 2015, Front. Mol. Biosci..

[16]  A. M. Goswami.  Structural modeling and in silico analysis of non-synonymous single nucleotide polymorphisms of human 3β-hydroxysteroid dehydrogenase type 2 , 2015, Meta gene.

[17]  Bernd Bukau,et al.  Substrate specificity of the DnaK chaperone determined by screening cellulose‐bound peptide libraries , 1997, The EMBO journal.

[18]  P. Kloetzel,et al.  Modeling the MHC class I pathway by combining predictions of proteasomal cleavage,TAP transport and MHC class I binding , 2005, Cellular and Molecular Life Sciences CMLS.

[19]  L. Gierasch,et al.  How hsp70 molecular machines interact with their substrates to mediate diverse physiological functions. , 2015, Journal of molecular biology.

[20]  Roman Kityk,et al.  Nucleotide exchange factors Fes1 and HspBP1 mimic substrate to release misfolded proteins from Hsp70 , 2018, Nature Structural & Molecular Biology.

[21]  J. Reinstein,et al.  Bap (Sil1) regulates the molecular chaperone BiP by coupling release of nucleotide and substrate , 2018, Nature Structural & Molecular Biology.

[22]  Md. Solayman,et al.  In silico analysis of nonsynonymous single nucleotide polymorphisms of the human adiponectin receptor 2 (ADIPOR2) gene , 2017, Comput. Biol. Chem..

[23]  C. Hughes,et al.  Heterogeneous binding of the SH3 client protein to the DnaK molecular chaperone , 2015, Proceedings of the National Academy of Sciences.

[24]  François Stricher,et al.  The FoldX web server: an online force field , 2005, Nucleic Acids Res..

[25]  Takaya Saito,et al.  The Precision-Recall Plot Is More Informative than the ROC Plot When Evaluating Binary Classifiers on Imbalanced Datasets , 2015, PloS one.

[26]  J. Buchner,et al.  BiPPred: Combined sequence‐ and structure‐based prediction of peptide binding to the Hsp70 chaperone BiP , 2016, Proteins: Structure, Function, and Bioinformatics.