Molecular interactions of c-ABL mutants in complex with imatinib/nilotinib: a computational study using linear interaction energy (LIE) calculations

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.     

Interplay between kinase domain autophosphorylation and F-actin binding domain in regulating imatinib sensitivity and nuclear import of BCR-ABL

Background

The constitutively activated BCR-ABL tyrosine kinase of chronic myeloid leukemia (CML) is localized exclusively to the cytoplasm despite the three nuclear localization signals (NLS) in the ABL portion of this fusion protein. The NLS function of BCR-ABL is re-activated by a kinase inhibitor, imatinib, and in a kinase-defective BCR-ABL mutant. The mechanism of this kinase-dependent inhibition of the NLS function is not understood.

Methodology/Principal Findings

By examining the subcellular localization of mutant BCR-ABL proteins under conditions of imatinib and/or leptomycin B treatment to inhibit nuclear export, we have found that mutations of three specific tyrosines (Y232, Y253, Y257, according to ABL-1a numbering) in the kinase domain can inhibit the NLS function of kinase-proficient and kinase-defective BCR-ABL. Interestingly, binding of imatinib to the kinase-defective tyrosine-mutant restored the NLS function, suggesting that the kinase domain conformation induced by imatinib-binding is critical to the re-activation of the NLS function. The C-terminal region of ABL contains an F-actin binding domain (FABD). We examined the subcellular localization of several FABD-mutants and found that this domain is also required for the activated kinase to inhibit the NLS function; however, the binding to F-actin per se is not important. Furthermore, we found that some of the C-terminal deletions reduced the kinase sensitivity to imatinib.

Conclusions/Significance

Results from this study suggest that an autophosphorylation-dependent kinase conformation together with the C-terminal region including the FABD imposes a blockade of the BCR-ABL NLS function. Conversely, conformation of the C-terminal region including the FABD can influence the binding affinity of imatinib for the kinase domain. Elucidating the structural interactions among the kinase domain, the NLS region and the FABD may therefore provide insights on the design of next generation BCR-ABL inhibitors for the treatment of CML.     

Uncovering new aspects of protein interactions through analysis of specificity landscapes in peptide recognition domains

Protein interactions underlie all biological processes. An important class of protein interactions, often observed in signaling pathways, consists of peptide recognition domains binding short protein segments on the surface of their target proteins. Recent developments in experimental techniques have uncovered many such interactions and shed new lights on their specificity. To analyze these data, novel computational methods have been introduced that can accurately describe the specificity landscape of peptide recognition domains and predict new interactions. Combining large-scale analysis of binding specificity data with structure-based modeling can further reveal new biological insights into the molecular recognition events underlying signaling pathways.     

BCR-ABL tyrosine kinase inhibitor pharmacophore model derived from a series of phenylaminopyrimidine-based (PAP) derivatives

To reveal novel insights into the inhibition of BCR-ABL tyrosine kinase, pharmacophore mapping studies were performed for a series of phenylaminopyrimidine-based (PAP) derivatives, including imatinib (Gleevec). A seven-point pharmacophore model with one hydrophobic group (H), two hydrogen bond donors (D) and four aromatic rings (R) was developed using phase (pharmacophore alignment & scoring engine). The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of 0.886 and a survival score of 4.97 for training set molecules. The model showed excellent predictive power, with a correlation coefficient of Q² = 0.768 for an external test set of ten molecules. The results obtained from our studies provide a valuable tool for designing new lead molecules with potent activity.     

Ponatinib Is a Pan-BCR-ABL Kinase Inhibitor: MD Simulations and SIE Study

BCR-ABL kinase domain inhibition can be used to treat chronic myeloid leukemia. The inhibitors such as imatinib, dasatinib and nilotinib are effective drugs but are resistant to some BCR-ABL mutations. The pan-BCR-ABL kinase inhibitor ponatinib exhibits potent activity against native, T315I, and all other clinically relevant mutants, and showed better inhibition than the previously known inhibitors. We have studied the molecular dynamics simulations and calculated solvated interaction energies of native and fourteen mutant BCR-ABL kinases (M244V, G250E, Q252H, Y253F, Y253H, E255K, E255V, T315A, T315I, F317L, F317V, M351T, F359V and H396P) complexed with ponatinib. These studies revealed that the interactions between ponatinib and individual residues in BCR-ABL kinase are also affected due to the remote residue mutations. We report that some residues, Met244, Lys245, Gln252, Gly254, Leu370 and Leu298 do not undergo any conformational changes, while the fluctuations in residues from P-loop, β3-, β5- strands and αC- helix are mainly responsible for ponatinib binding to native and all mutant BCR-ABL kinases. Our work provides the molecular mechanisms of native and mutant BCR-ABL kinases inhibition by ponatinib at atomic level that has not been studied before.     

How does the novel T315L mutation of breakpoint cluster region-abelson (BCR-ABL) kinase confer resistance to ponatinib: a comparative molecular dynamics simulation study

Abstract

Acute lymphocytic leukemia (ALL) is one of the most dangerous types of leukemia, and about 40% of them is Philadelphia chromosome-positive acute lymphocytic leukemia (Ph?+?ALL). Ph?+?ALL is caused by the fusion of the breakpoint cluster region (BCR) and the Ableson (ABL) genes, named the BCR-ABL fused gene that codes for an autonomously active tyrosine kinase. Tyrosine kinase inhibitors (TKIs) are among the first-line therapeutic agents for the treatment of Ph?+?ALL. Drug resistance are the major obstacle, limiting their clinical utility. The latest third-generation TKIs, ponatinib, can tackle most abnormal BCR-ABL kinases, including the T315I mutant that is resistant to first- and second-generations TKIs such as imatinib. However, drug resistance still emerges with the novel T315L mutation and the underlying mechanisms remain elusive. Here, using molecular dynamics (MD) simulations, we explored into the detailed interactions between ponatinib and BCR-ABL in the wild-type (WT), T315I, and T315L systems. The simulations revealed the significant conformational changes of ponatinib in its binding site due to the T315L mutation and the underlying structural mechanisms. Binding free energy analysis unveiled that the affinity of ponatinib to BCR-ABL decreased upon T315L mutation, which resulted in its unfavorable binding and drug resistance. Key residues responsible for the unfavored unbinding were also identified. This study elucidates the detailed mechanisms for the resistance of ponatinib in Ph?+?ALL triggered by the T315L mutation and will provide insights for future drug development and optimization.     

Structure-activity relationship of C6-C3 phenylpropanoids on xanthine oxidase-inhibiting and free radical-scavenging activities

We employed the techniques of DNA relaxation, DPPH (1,1-diphenyl-2-picrylhydrazyl hydrate), and DMPO (5,5-dimethyl-1-pyrroline-N-oxide)-electron spin resonance (ESR), to study the effects of reactive oxygen species (ROS) suppression by 11 selected C6-C3 phenylpropanoid derivatives under oxidative conditions. We also investigated the effects of the derivatives on the inhibition of xanthine oxidase (XO) activity, and the structure-activity relationships (SARs) of these derivatives against XO activity were further examined using computer-aided molecular modeling. Caffeic acid was the most potent radical scavenger among the 11 test compounds. Our results suggest that the chemical structure and number of hydroxyl groups on the benzene ring of phenylpropanoids are correlated with the effects of ROS suppression. All test derivatives were competitive inhibitors of XO. The results of the structure-based molecular modeling exhibited interactions between phenylpropanoid derivatives and the molybdopterin region of XO. The para-hydroxyl of phenylpropanoid derivatives was pointed toward the guanidinium group of Arg 880. The phenylpropanoid derivatives containing the meta-or ortho-hydroxyl formed hydrogen bonds with Thr 1010. In addition, meta-hydroxyl formed hydrogen bonds with the peptide bond between the residues of Thr1010 and Phe1009. CAPE, the phenylenethyl ester of phenylpropanoids, had the highest affinity toward the binding site of XO, and we speculated that this was due to hydrophobic interactions of the phenylethyl ester with several hydrophobic residues surrounding the active site. The hypoxanthine/XO reaction in the DMPO-ESR technique was used to correlate the effects of these phenylpropanoid derivatives on enzyme inhibition and ROS suppression, and the results showed that caffeic acid and CAPE were the two most potent agents among the tested compounds. We further assessed the effects of the test compounds on living cells, and CAPE was the most potent agent for protecting cells against ROS-mediated damage among the tested phenylpropanoids.     

Mechanisms of resistance to BCR-ABL and other kinase inhibitors

In this article, we are reviewing the molecular mechanisms that lead to kinase inhibitor resistance. As the oncogenic BCR-ABL kinase is the target of the first approved small-molecule kinase inhibitor imatinib, we will first focus on the structural and mechanistic basis for imatinib resistance. We will then show ways how next generations of BCR-ABL inhibitors and alternative targeting strategies have helped to offer effective treatment options for imatinib-resistant patients. Based on these insights, we discuss commonalities and further mechanisms that lead to resistance to other kinase inhibitors in solid tumors. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

Intracellular retention of ABL kinase inhibitors determines commitment to apoptosis in CML cells

Clinical development of imatinib in CML established continuous target inhibition as a paradigm for successful tyrosine kinase inhibitor (TKI) therapy. However, recent reports suggested that transient potent target inhibition of BCR-ABL by high-dose TKI (HD-TKI) pulse-exposure is sufficient to irreversibly commit cells to apoptosis. Here, we report a novel mechanism of prolonged intracellular TKI activity upon HD-TKI pulse-exposure (imatinib, dasatinib) in BCR-ABL-positive cells. Comprehensive mechanistic exploration revealed dramatic intracellular accumulation of TKIs which closely correlated with induction of apoptosis. Cells were rescued from apoptosis upon HD-TKI pulse either by repetitive drug wash-out or by overexpression of ABC-family drug transporters. Inhibition of ABCB1 restored sensitivity to HD-TKI pulse-exposure. Thus, our data provide evidence that intracellular drug retention crucially determines biological activity of imatinib and dasatinib. These studies may refine our current thinking on critical requirements of TKI dose and duration of target inhibition for biological activity of TKIs.     

Molecular Probing of the HPV-16 E6 Protein Alpha Helix Binding Groove with Small Molecule Inhibitors

The human papillomavirus (HPV) HPV E6 protein has emerged as a central oncoprotein in HPV-associated cancers in which sustained expression is required for tumor progression. A majority of the E6 protein interactions within the human proteome use an alpha-helix groove interface for binding. The UBE3A/E6AP HECT domain ubiquitin ligase binds E6 at this helix-groove interface. This enables formation of a trimeric complex with p53, resulting in destruction of this tumor suppressor. While recent x-ray crystal structures are useful, examples of small molecule probes that can modulate protein interactions at this interface are limited. To develop insights useful for potential structure-based design of ligands for HPV E6, a series of 2,6-disubstituted benzopyranones were prepared and tested as competitive antagonists of E6-E6AP helix-groove interactions. These small molecule probes were used in both binding and functional assays to evaluate recognition features of the E6 protein. Evidence for an ionic functional group interaction within the helix groove was implicated by the structure-activity among the highest affinity ligands. The molecular topographies of these protein-ligand interactions were evaluated by comparing the binding and activities of single amino acid E6 mutants with the results of molecular dynamic simulations. A group of arginine residues that form a rim-cap over the E6 helix groove offer compensatory roles in binding and recognition of the small molecule probes. The flexibility and impact on the overall helix-groove shape dictated by these residues offer new insights for structure-based targeting of HPV E6.     

Real-Time Analysis of Imatinib- and Dasatinib-Induced Effects on Chronic Myelogenous Leukemia Cell Interaction with Fibronectin

Attachment of stem leukemic cells to the bone marrow extracellular matrix increases their resistance to chemotherapy and contributes to the disease persistence. In chronic myelogenous leukemia (CML), the activity of the fusion BCR-ABL kinase affects adhesion signaling. Using real-time monitoring of microimpedance, we studied in detail the kinetics of interaction of human CML cells (JURL-MK1, MOLM-7) and of control BCR-ABL-negative leukemia cells (HEL, JURKAT) with fibronectin-coated surface. The effect of two clinically used kinase inhibitors, imatinib (a relatively specific c-ABL inhibitor) and dasatinib (dual ABL/SRC family kinase inhibitor), on cell binding to fibronectin is described. Both imatinib and low-dose (several nM) dasatinib reinforced CML cell interaction with fibronectin while no significant change was induced in BCR-ABL-negative cells. On the other hand, clinically relevant doses of dasatinib (100 nM) had almost no effect in CML cells. The efficiency of the inhibitors in blocking the activity of BCR-ABL and SRC-family kinases was assessed from the extent of phosphorylation at autophosphorylation sites. In both CML cell lines, SRC kinases were found to be transactivated by BCR-ABL. In the intracellular context, EC50 for BCR-ABL inhibition was in subnanomolar range for dasatinib and in submicromolar one for imatinib. EC50 for direct inhibition of LYN kinase was found to be about 20 nM for dasatinib and more than 10 µM for imatinib. Cells pretreated with 100 nM dasatinib were still able to bind to fibronectin and SRC kinases are thus not necessary for the formation of cell-matrix contacts. However, a minimal activity of SRC kinases might be required to mediate the increase in cell adhesivity induced by BCR-ABL inhibition. Indeed, active (autophosphorylated) LYN was found to localize in cell adhesive structures which were visualized using interference reflection microscopy.     

Comparison between molecularly defined and conventional therapeutics in a conditional BCR-ABL cell culture model

Accumulating knowledge about the molecular mechanisms causing human diseases can support the development of targeted therapies such as imatinib, a BCR-ABL-specific tyrosine kinase inhibitor to treat chronic myeloid leukemia (CML). Here, we use lentivirus-mediated RNA interference (RNAi) targeting BCR-ABL and the downstream signaling molecules SHP2, STAT5, and Gab2 to compare the efficacy and specificity of molecularly defined therapeutics with that of conventional cytotoxic drugs (cytarabine, doxorubicin, etoposide) in a conditional BCR-ABL cell culture model. IC(50) values were determined for each drug in TonB cells cultured either with interleukin-3 (IL-3) or BCR-ABL, and molecularly defined therapies were studied using lentivirally expressed shRNAs. We demonstrate that conventional anti-leukemic drugs have small or no differential effects under different cell culture conditions, whereas both imatinib and specific RNAi significantly inhibit proliferation of TonB cells in the presence of BCR-ABL but not IL-3. To study molecularly defined combination therapy, we evaluated either imatinib in TonB cells with target-specific RNAi or we used lentiviral vectors to induce combinatorial RNAi through simultaneous expression of two shRNAs. These combination therapies result in increased efficacy without loss in specificity. Interestingly, combinatorial RNAi can specifically deplete TonB cell cultures in the presence of BCR-ABL, even without targeting the oncogene itself. This model provides a tool to evaluate potential therapeutic targets and to quantify efficacy and specificity preclinically of new combination therapies in BCR-ABL-positive cells.     

BCR-ABL Gene Expression Is Required for Its Mutations in a Novel KCL-22 Cell Culture Model for Acquired Resistance of Chronic Myelogenous Leukemia

Acquired resistance through genetic mutations is a common phenomenon in several cancer therapies using molecularly targeted drugs, best exemplified by the BCR-ABL inhibitor imatinib in treating chronic myelogenous leukemia (CML). Overcoming acquired resistance is a daunting therapeutic challenge, and little is known about how these mutations evolve. To facilitate understanding the resistance mechanisms, we developed a novel culture model for CML acquired resistance in which the CML cell line KCL-22, following initial response to imatinib, develops resistant T315I BCR-ABL mutation. We demonstrate that the emergence of BCR-ABL mutations do not require pre-existing BCR-ABL mutations derived from the original patient as the subclones of KCL-22 cells can form various BCR-ABL mutations upon imatinib treatment. BCR-ABL mutation rates vary from cell clone to clone and passages, in contrast to the relatively stable mutation rate of the hypoxanthine-guanine phosphoribosyltransferase gene. Strikingly, development of BCR-ABL mutations depends on its gene expression because BCR-ABL knockdown completely blocks KCL-22 cell relapse on imatinib and acquisition of mutations. We further show that the endogenous BCR-ABL locus has significantly higher mutagenesis potential than the transduced randomly integrated BCR-ABL cDNA. Our study suggests important roles of BCR-ABL gene expression and its native chromosomal locus for acquisition of BCR-ABL mutations and provides a new tool for further studying resistance mechanisms.     

A computational docking study for prediction of binding mode of diospyrin and derivatives: Inhibitors of human and leishmanial DNA topoisomerase-I

A computational approach was utilized to study the relative binding modes of diospyrin (bisnaphthoquinonoid) with the crystal structure of human DNA-TopoI and the recently reported Leishmania donavani DNA-TopoI. Additionally, the binding site interactions of amino derivatives of diospyrin with human TopoI were studied extensively. Based on the docking results, binding modes of diospyrin with the human and leishmanial TopoI catalytic core were predicted. The parallel use of two efficient and predictive docking programs, GOLD and Ligandfit, allowed mutual validation of the predicted binding poses. A reasonably good correlation coefficient between the calculated docking scores and the experimentally determined cytotoxicity helped in validating the docking method. Furthermore, a structure-based pharmacophore model was developed for L. donavani DNA-TopoI inhibition which helped in elucidating the topological and spatial requirements of the ligand-receptor interactions. This study provides an understanding of the structural basis of ligand binding to the topoisomerase receptor, which may be used for the structure-based design of potent and novel ligands for anticancer and antileishmanial therapy. To our knowledge, this is the first report of a binding mode exploration study for diospyrin and its derivatives as inhibitors of the leishmanial and human TopoI enzymes.     

Different roles of p53 in the regulation of DNA damage caused by 1,2-heteroannelated anthraquinones and doxorubicin

The anthracyclin antibiotic agent doxorubicin (DXR) has been widely used as a chemotherapeutic drug for more than 40 years, but its clinical use has been limited by its cardiotoxicity. The mechanism of action of DXR remains uncertain and controversial. A series of 1,2-heteroannelated anthraquinones and anthra[1,2-d]imidazole-6,11-dione compounds were synthesized and their cytotoxicity profiles were analyzed using the National Cancer Institute 60 (NCI 60) platform and human telomerase inhibition assays. In the current study, three of the 1,2-heteroannelated anthraquinones, NSC745795, NSC745885 and NSC745887, were found to differ from each other with respect to their effects on cell cycle regulation, apoptosis, autophagy, senescence and their abilities to induce DNA damage. The differences depended on the presence or absence of a heterocyclic moiety, which suggested that the differences were due, at least in part, to differential effects on specific cellular targets, such as p53. In contrast to DXR, which induced p53 expression, treatment with NSC745885 resulted in the degradation of several proteins, including p53, via proteasome-dependent and -independent pathways in HeLa cells. These results provide insights into the molecular mechanisms governing cell inhibition by 1,2-heteroannelated anthraquinone derivatives and suggest that these mechanisms could serve as the basis for new structure-based drug designs.     

Expression and Activity of Fyn Mediate Proliferation and Blastic Features of Chronic Myelogenous Leukemia

The BCR-ABL1 oncogene is a tyrosine kinase that activates many signaling pathways, resulting in the induction of chronic myeloid leukemia (CML). Kinase inhibitors, such as imatinib, have been developed for the treatment of CML; however, the terminal, blast crisis phase of the disease remains a clinical challenge. Blast crisis CML is difficult to treat due to resistance to tyrosine kinase inhibitors, increased genomic instability and acquired secondary mutations. Our recent studies uncovered a role for Fyn in promoting BCR-ABL1 mediated cell growth and sensitivity to imatinib. Here we demonstrate that Fyn contributes to BCR-ABL1 induced genomic instability, a feature of blast crisis CML. Bone marrow cells and mouse embryonic fibroblasts derived from Fyn knockout mice transduced with BCR-ABL1 display slowed growth and clonogenic potential as compared to Fyn wild-type BCR-ABL1 expressing counterparts. K562 cells overexpressing constitutively active Fyn kinase were larger in size and displayed an accumulation of genomic abnormalities such as chromosomal aberrations and polyploidy. Importantly, loss of Fyn protected mouse embryonic fibroblast cells from increased number of chromosomal aberrations and fragments induced by BCR-ABL1. Together, these results reveal a novel role for Fyn in regulating events required for genomic maintenance and suggest that Fyn kinase activity plays a role in the progression of CML to blast crisis.     

Computational studies of the binding mechanisms of fullerenes to human serum albumin

Fullerene and its derivatives show promising prospects for applications in a vast array of biological systems. A key aspect concerning their biomedical applications is how they interact with proteins from molecular levels, which is still poorly understood. In the current study, we investigated the structural and thermodynamic basis of the interactions between two pharmacologically relevant fullerene derivatives and human serum albumin (HSA) using molecular docking, molecular dynamics simulations, and binding free energy calculations. Our results demonstrate that fullerenes steadily bind with HSA at the interfacial cavity formed by subdomains IIA and IIIA. In agreement with available experimental data, our simulations show that the global structure of HSA becomes more compact in the presence of fullerene, while local structural dynamics of the binding cavity behaves diversely depending on the chemical properties of bound fullerenes. Binding free energy calculations confirmed that the interactions between fullerenes and HSA are dominantly stabilized by van der Waals forces and they further allowed the identification of key residues involved in fullerene binding. The structural and energetic insights obtained from this work may help for the development of fullerene-based drug delivery devices and therapeutic agents with improved biological profile.     

Nilotinib based pharmacophore models for BCRABL

Tyrosine kinase inhibitors have revolutionized the treatment of several malignancies, converting lethal diseases in a manageable aspect. Imitanib, a small molecule ABL kinase inhibitor is a highly effective therapy for early phase chronic myeloid leukemia (CML), which has constitutively active ABL kinase activity owing to the over expression of the BCR-ABL fusion protein. But some patients develop imatinib resistance, particularly in the advanced phases of CML.The discovery of resistance mechanisms of imitanib; urge forward the development of second generation drugs. Nilotinib, a second generation drug is more potent inhibitor of BCR-ABL than imatinib. But nilotinib also develops dermatologic events and headache in patients. Large information about BCR-ABL structure and its inhibitors are now available. Based on the pharmacophore modeling approaches, it is possible to decipher the molecular determinants to inhibit BCR-ABL. We conducted a structure based and ligand based study to identify potent natural compounds as BCR-ABL inhibitor. First kinase inhibitors were docked with the receptor (BCR-ABL) and nilotinib was selected as a pharmacophore due its high binding efficiency. Eleven compounds were selected out of 1457 substances which have mutual pharmacopohre features with nilotinib. These eleven compounds were validated and used for docking study to find the drug like molecules. The best molecules from the final set of screening candidates can be evaluated in cell lines and may represent a novel class of BCR-ABL inhibitors.

Abbreviations

CML - Chronic myeloid leukemia, PDGFR - Platelet derived growth factor receptor, TKI - Tyrosine kinase inhibitors.