ABSTRACT
Recent experimental evidence suggests that mebendazole, a popular antiparasitic drug, binds to heat shock protein 90 (Hsp90) and inhibits acute myeloid leukemia cell growth. In this study we use quantum mechanics (QM), molecular similarity, and molecular dynamics (MD) calculations to predict possible binding poses of mebendazole to the adenosine triphosphate (ATP) binding site of Hsp90. Extensive conformational searches and minimization of the five mebendazole tautomers using the MP2/aug-cc-pVTZ theory level resulted in 152 minima. Mebendazole-Hsp90 complex models were subsequently created using the QM optimized conformations and protein coordinates obtained from experimental crystal structures that were chosen through similarity calculations. Nine different poses were identified from a total of 600 ns of explicit solvent, all-atom MD simulations using two different force fields. All simulations support the hypothesis that mebendazole is able to bind to the ATP binding site of Hsp90.
Subject(s)
Mebendazole , Molecular Dynamics Simulation , Adenosine Triphosphate , HSP90 Heat-Shock Proteins/chemistry , Humans , Mebendazole/pharmacology , Molecular Conformation , Protein Binding , Protein ConformationABSTRACT
The prognosis of elderly AML patients is still poor due to chemotherapy resistance. The Hedgehog (HH) pathway is important for leukemic transformation because of aberrant activation of GLI transcription factors. MBZ is a well-tolerated anthelmintic that exhibits strong antitumor effects. Herein, we show that MBZ induced strong, dose-dependent anti-leukemic effects on AML cells, including the sensitization of AML cells to chemotherapy with cytarabine. MBZ strongly reduced intracellular protein levels of GLI1/GLI2 transcription factors. Consequently, MBZ reduced the GLI promoter activity as observed in luciferase-based reporter assays in AML cell lines. Further analysis revealed that MBZ mediates its anti-leukemic effects by promoting the proteasomal degradation of GLI transcription factors via inhibition of HSP70/90 chaperone activity. Extensive molecular dynamics simulations were performed on the MBZ-HSP90 complex, showing a stable binding interaction at the ATP binding site. Importantly, two patients with refractory AML were treated with MBZ in an off-label setting and MBZ effectively reduced the GLI signaling activity in a modified plasma inhibitory assay, resulting in a decrease in peripheral blood blast counts in one patient. Our data prove that MBZ is an effective GLI inhibitor that should be evaluated in combination to conventional chemotherapy in the clinical setting.
Subject(s)
Leukemia, Myeloid, Acute/metabolism , Mebendazole/pharmacology , Proteasome Endopeptidase Complex/metabolism , Tubulin Modulators/pharmacology , Zinc Finger Protein GLI1/metabolism , Case-Control Studies , Heat-Shock Proteins/chemistry , Heat-Shock Proteins/metabolism , Humans , Leukemia, Myeloid, Acute/etiology , Leukemia, Myeloid, Acute/pathology , Leukemia, Myeloid, Acute/therapy , Proteolysis , Signal Transduction/drug effects , Structure-Activity Relationship , Transcription Factors/chemistry , Transcription Factors/metabolism , Zinc Finger Protein GLI1/antagonists & inhibitors , Zinc Finger Protein GLI1/chemistryABSTRACT
Aberrant activation of the hedgehog (HH) pathway is observed in many neoplasms, including acute myeloid leukemia (AML). The glioma-associated oncogene homolog (GLI) transcription factors are the main downstream effectors of the HH signaling cascade and are responsible for the proliferation and maintenance of leukemic stem cells, which support chemotherapy resistance and leukemia relapse. Cytarabine (Ara-C)-resistant variants of AML cell lines were established through long-term cultivation with successively increasing Ara-C concentrations. Subsequently, differences in GLI expression were analyzed by RT-qPCR. GLI3 mRNA levels were detectable in parental Kasumi-1, OCI-AML3, and OCI-AML5 cells, whereas GLI3 expression was completely silenced in all resistant counterparts. Therefore, we generated GLI3-knockdown cell lines using small hairpin RNAs (shRNA) and evaluated their sensitivity to Ara-C in vitro. The knockdown of GLI3 partly abolished the effect of Ara-C on colony formation and induction of apoptosis, indicating that GLI3 downregulation results in Ara-C resistance. Moreover, we analyzed the expression of several genes involved in Ara-C metabolism and transport. Knockdown of GLI3 resulted in the upregulation of SAM and HD domain-containing protein 1 (SAMHD1), cytidine deaminase (CDA), and ATP-binding cassette C11 (ABCC11)/multidrug resistance-associated protein 8 (MRP8), each of which has been identified as a predictive marker for Ara-C response in acute myeloid leukemia. Our results demonstrate that GLI3 downregulation is a potential mechanism to induce chemotherapy resistance in AML.
