Distinguishing the optimal binding mechanism through reversible and irreversible inhibition analysis of HSP72 protein in cancer therapy

Over the past two decades, covalent inhibitors have gained much interest and are living up to their reputation as a powerful tool in drug discovery. Covalent inhibitors possess several significant advantages, including increased biochemical efficiency, prolonged duration and the ability to target shallow, solvent-exposed substrate-binding domains. One of the enzymes that have been both covalently and non-covalently targeted is the heat shock protein 72 (HSP72). This elevated enzyme expression in cancer cells may be responsible for tumorigenesis and tumor progression by providing chemotherapy resistance. A critical gap remains in the molecular understanding of the structural mechanism's covalent and non-covalent binding to HSP72. In this study, we explore the most optimal binding mechanism in the inhibition of the HSP72. Based on the molecular dynamic analyses, it was evident that the non-covalent complex showed more stability than the covalent complex. The covalent ligand, however, was more able to induce and stabilize the sealed conformation of the HSP72-NBD ATP binding domain throughout the. Also, the non-covalent ligand does not induce any significant conformational change as it remained close to the shape of the unbound complex; and the affinity is only dependent on the multiple hydrogen bonds in contrast to the covalent ligand. This is supported by the secondary structure elements and principal component analysis that was more dominant in the covalently inhibited complex. Covalent bond induced the α-helices sealed conformation of the HSP72-NBD; based on our findings, this will prevent other small molecules from interacting at the ATP binding site domain. Moreover, inhibition of the ATP binding domain can directly affect the ATPs protein folding mechanism of the HSP72 enzyme. The essential dynamic analysis presented in this report compliments the binding mechanism of HSP72, establishing covalent inhibition as the preferred method of inhibiting the HSP72 protein. The findings from this study may assist in the design of more target-specific HSP72 covalent inhibitors exploring the surface-exposed lysine residues.

[1]  R. V. van Montfort,et al.  An Irreversible Inhibitor of HSP72 that Unexpectedly Targets Lysine‐56 , 2017, Angewandte Chemie.

[2]  Yuping Wei,et al.  Noncovalent interaction-assisted drug delivery system with highly efficient uptake and release of paclitaxel for anticancer therapy , 2017, International journal of nanomedicine.

[3]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[4]  R. V. van Montfort,et al.  Exploiting Protein Conformational Change to Optimize Adenosine-Derived Inhibitors of HSP70 , 2016, Journal of medicinal chemistry.

[5]  J. C. Patterson,et al.  Molecular Dynamics Study of Zn(Aβ) and Zn(Aβ)2 , 2013, PloS one.

[6]  H. Hotelling.  Analysis of a complex of statistical variables into principal components. , 1933 .

[7]  R. A. Bauer.  Covalent inhibitors in drug discovery: from accidental discoveries to avoided liabilities and designed therapies. , 2015, Drug discovery today.

[8]  Y. Ghasemi,et al.  MARK4 protein can explore the active-like conformations in its non-phosphorylated state , 2019, Scientific Reports.

[9]  Matthew P. Repasky,et al.  Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. , 2006, Journal of medicinal chemistry.

[10]  M. Mayer,et al.  Functional Analysis of Hsp70 Inhibitors , 2013, PloS one.

[11]  Ivan Merelli,et al.  Molecular Mechanism of Allosteric Communication in Hsp70 Revealed by Molecular Dynamics Simulations , 2012, PLoS Comput. Biol..

[12]  Jianzhong Chen,et al.  Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations , 2018, RSC advances.

[13]  Caterina Barillari,et al.  Targeting HSP70: The second potentially druggable heat shock protein and molecular chaperone? , 2010, Cell cycle.

[14]  M. Soliman,et al.  Dynamics of the thumb-finger regions in a GH11 xylanase Bacillus circulans : comparison between the Michaelis and covalent intermediate , 2015 .

[15]  E. Seifert,et al.  OriginPro 9.1: Scientific Data Analysis and Graphing Software - Software Review , 2014, J. Chem. Inf. Model..

[16]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[17]  Haruki Nakamura,et al.  A Novel Approach of Dynamic Cross Correlation Analysis on Molecular Dynamics Simulations and Its Application to Ets1 Dimer–DNA Complex , 2014, PloS one.

[18]  T. Richmond,et al.  Solvent accessible surface area and excluded volume in proteins. Analytical equations for overlapping spheres and implications for the hydrophobic effect. , 1984, Journal of molecular biology.

[19]  Monique Laberge,et al.  Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics. , 2008, Biophysical journal.

[20]  Lila M. Gierasch,et al.  An Interdomain Energetic Tug-of-War Creates the Allosterically Active State in Hsp70 Molecular Chaperones , 2012, Cell.

[21]  David L. Penkler,et al.  Perturbation-Response Scanning Reveals Key Residues for Allosteric Control in Hsp70 , 2017, J. Chem. Inf. Model..

[22]  Robert J. Ono,et al.  Role of non-covalent and covalent interactions in cargo loading capacity and stability of polymeric micelles. , 2014, Journal of controlled release : official journal of the Controlled Release Society.

[23]  P. Zhou,et al.  Specific noncovalent interactions at protein-ligand interface: implications for rational drug design. , 2012, Current medicinal chemistry.

[24]  B. Cravatt,et al.  Strategies for discovering and derisking covalent, irreversible enzyme inhibitors. , 2010, Future medicinal chemistry.

[25]  M. M. Hoque,et al.  Non-covalent interactions involving halogenated derivatives of capecitabine and thymidylate synthase: a computational approach , 2016, SpringerPlus.

[26]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[27]  D. Masison,et al.  Hsp70 structure, function, regulation and influence on yeast prions. , 2009, Protein and peptide letters.

[28]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

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

[30]  Yun-Xin Fu,et al.  Insights into the molecular mechanism underlying CD4-dependency and neutralization sensitivity of HIV-1: a comparative molecular dynamics study on gp120s from isolates with different phenotypes , 2018, RSC advances.

[31]  P E Bourne,et al.  The Protein Data Bank. , 2002, Nucleic acids research.

[32]  Bernd Bukau,et al.  Multistep mechanism of substrate binding determines chaperone activity of Hsp70 , 2000, Nature Structural Biology.

[33]  P. Christen,et al.  Sequence-specific rates of interaction of target peptides with the molecular chaperones DnaK and DnaJ. , 1998, Biochemistry.

[34]  Woody Sherman,et al.  Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments , 2013, Journal of Computer-Aided Molecular Design.

[35]  M. Soliman,et al.  Delving into Zika virus structural dynamics – a closer look at NS3 helicase loop flexibility and its role in drug discovery , 2017 .

[36]  M. Drysdale,et al.  Novel adenosine-derived inhibitors of 70 kDa heat shock protein, discovered through structure-based design. , 2009, Journal of medicinal chemistry.

[37]  Carlos Simmerling,et al.  Grid-based backbone correction to the ff12SB protein force field for implicit-solvent simulations. , 2015, Journal of chemical theory and computation.

[38]  D. McConkey,et al.  Bortezomib inhibits PKR-like endoplasmic reticulum (ER) kinase and induces apoptosis via ER stress in human pancreatic cancer cells. , 2005, Cancer research.

[39]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[40]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[41]  U. Manne,et al.  Targeting Hsp70: A possible therapy for cancer. , 2016, Cancer letters.

[42]  Suranjana Das,et al.  Optimized Hydrophobic Interactions and Hydrogen Bonding at the Target-Ligand Interface Leads the Pathways of Drug-Designing , 2010, PloS one.

[43]  J Osipiuk,et al.  Human Hsp70 molecular chaperone binds two calcium ions within the ATPase domain. , 1997, Structure.

[44]  E. Budianto,et al.  The molecular docking of 1,4-naphthoquinone derivatives as inhibitors of Polo-like kinase 1 using Molegro Virtual Docker - , 2014 .

[45]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[46]  Pedro Gonnet,et al.  P-SHAKE: A quadratically convergent SHAKE in O(n2) , 2007, J. Comput. Phys..

[47]  A. Frontera,et al.  Covalent and Non-covalent Noble Gas Bonding Interactions in XeFn Derivatives (n = 2–6): A Combined Theoretical and ICSD Analysis , 2020, Frontiers in Chemistry.

[48]  P. Christen,et al.  Kinetics of molecular chaperone action. , 1994, Science.

[49]  S. Sugano,et al.  Direct inter-subdomain interactions switch between the closed and open forms of the Hsp70 nucleotide-binding domain in the nucleotide-free state. , 2010, Acta crystallographica. Section D, Biological crystallography.

[50]  David van der Spoel,et al.  The Origin of Layer Structure Artifacts in Simulations of Liquid Water. , 2006 .

[51]  Jianzhong Chen,et al.  Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations , 2019, Nucleic acids research.

[52]  H. Abdi,et al.  Principal component analysis , 2010 .

[53]  M J Harvey,et al.  An Implementation of the Smooth Particle Mesh Ewald Method on GPU Hardware. , 2009, Journal of chemical theory and computation.

[54]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[55]  A. R. Srinivasan,et al.  Quasi‐harmonic method for studying very low frequency modes in proteins , 1984, Biopolymers.

[56]  Reversible versus irreversible inhibition modes of ERK2: a comparative analysis for ERK2 protein kinase in cancer therapy. , 2018, Future medicinal chemistry.