Binding free energy calculation with QM/MM hybrid methods for Abl-Kinase inhibitor

We report a Quantum mechanics/Molecular Mechanics–Poisson-Boltzmann/ Surface Area (QM/MM-PB/SA) method to calculate the binding free energy of c-Abl human tyrosine kinase by combining the QM and MM principles where the ligand is treated quantum mechanically and the rest of the receptor by classical molecular mechanics. To study the role of entropy and the flexibility of the protein ligand complex in a solvated environment, molecular dynamics calculations are performed using a hybrid QM/MM approach. This work shows that the results of the QM/MM approach are strongly correlated with the binding affinity. The QM/MM interaction energy in our reported study confirms the importance of electronic and polarization contributions, which are often neglected in classical MM-PB/SA calculations. Moreover, a comparison of semi-empirical methods like DFTB-SCC, PM3, MNDO, MNDO-PDDG, and PDDG-PM3 is also performed. The results of the study show that the implementation of a DFTB-SCC semi-empirical Hamiltonian that is derived from DFT gives better results than other methods. We have performed such studies using the AMBER molecular dynamic package for the first time. The calculated binding free energy is also in agreement with the experimentally determined binding affinity for c-Abl tyrosine kinase complex with Imatinib.     

A computational analysis of the binding affinities of FKBP12 inhibitors using the MM-PB/SA method

The FK506-binding proteins have been targets of pharmaceutical interests over years. We have studied the binding of a set of 12 nonimmunosuppressive small-molecule inhibitors to FKBP12 through molecular dynamics simulations. Each complex was subjected to 1-ns MD simulation conducted in an explicit solvent environment under constant temperature and pressure. The binding free energy of each complex was then computed by the MM-PB/SA method in the AMBER program. Our MM-PB/SA computation produced a good correlation between the experimentally determined and the computed binding free energies with a correlation coefficient (R(2)) of 0.93 and a standard deviation as low as 0.30 kcal/mol. The vibrational entropy term given by the normal mode analysis was found to be helpful for achieving this correlation. Moreover, an adjustment to one weight factor in the PB/SA model was essential to correct the absolute values of the final binding free energies to a reasonable range. A head-to-head comparison of our MM-PB/SA model with a previously reported Linear Response Approximation (LRA) model suggested that the MM-PB/SA method is more robust in binding affinity prediction for this class of compounds.     

In silico profiling of tyrosine kinases binding specificity and drug resistance using Monte Carlo simulations with the ensembles of protein kinase crystal structures

The molecular basis of the tyrosine kinases binding specificity and drug resistance against cancer drugs Imatinib and Dasatinib is elucidated using Monte Carlo simulations of the inhibitor-receptor binding with the ensembles of protein kinase crystal structures. In silico proteomics analysis unravels mechanisms by which structural plasticity of the tyrosine kinases is linked with the conformational preferences of Imatinib and Dasatinib in achieving effective drug binding with a distinct spectrum of the tyrosine kinome. The differences in the inhibitor sensitivities to the ABL kinase mutants are rationalized based on variations in the binding free energy profiles with the conformational states of the ABL kinase. While Imatinib binding is highly sensitive to the activation state of the enzyme, the computed binding profile of Dasatinib is remarkably tolerant to the conformational state of ABL. A comparative analysis of the inhibitor binding profiles with the clinically important ABL kinase mutants has revealed an excellent agreement with the biochemical and proteomics data. We have found that conformational adaptability of the kinase inhibitors to structurally different conformational states of the tyrosine kinases may have pharmacological relevance in acquiring a specific array of potent activities and regulating a scope of the inhibitor resistance mutations. This study outlines a useful approach for understanding and predicting the molecular basis of the inhibitor sensitivity against potential kinase targets and drug resistance.     

Role of pH on dimeric interactions for DENV envelope protein: an insight from molecular dynamics study

The entry of dengue viruses is mediated by pH triggering in the host cells. Here we have studied the DENV E protein stability and binding of its units at low and normal pH using MD and MM-PB/SA method for the first time. To investigate the role of pH in dissociation of dimeric protein, we have performed a concise study of hydrogen bonding and other interactions between units of dimer at low and normal pH. The Generalized Born calculation connotes that dimeric unit was relatively less stable and less proned for dimerisation at low pH. Our results provide a theoretical verification for previous assumptions of pH triggering mechanism of dengue envelope protein. During the pH alteration, we found a large decrement in salt bridges which were observed at normal pH. We also compared the flexibility of each unit and found that they exhibit different fluctuations during molecular dynamics simulations.     

Towards a Molecular Understanding of the Link between Imatinib Resistance and Kinase Conformational Dynamics

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics.     

Accurate prediction of the binding free energy and analysis of the mechanism of the interaction of replication protein A (RPA) with ssDNA

The eukaryotic replication protein A (RPA) has several pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Moreover, RPA seems to have a crucial role in organizing the sequential assembly of DNA processing proteins along single stranded DNA (ssDNA). The strong RPA affinity for ssDNA, K(A) between 10(-9)-10(-10) M, is characterized by a low cooperativity with minor variation for changes on the nucleotide sequence. Recently, new data on RPA interactions was reported, including the binding free energy of the complex RPA70AB with dC(8) and dC(5), which has been estimated to be -10 ± 0.4 kcal mol(-1) and -7 ± 1 kcal mol(-1), respectively. In view of these results we performed a study based on molecular dynamics aimed to reproduce the absolute binding free energy of RPA70AB with the dC(5) and dC(8) oligonucleotides. We used several tools to analyze the binding free energy, rigidity, and time evolution of the complex. The results obtained by MM-PBSA method, with the use of ligand free geometry as a reference for the receptor in the separate trajectory approach, are in excellent agreement with the experimental data, with ±4 kcal mol(-1) error. This result shows that the MM-PB(GB)SA methods can provide accurate quantitative estimates of the binding free energy for interacting complexes when appropriate geometries are used for the receptor, ligand and complex. The decomposition of the MM-GBSA energy for each residue in the receptor allowed us to correlate the change of the affinity of the mutated protein with the ΔG(gas+sol) contribution of the residue considered in the mutation. The agreement with experiment is optimal and a strong change in the binding free energy can be considered as the dominant factor in the loss for the binding affinity resulting from mutation.     

Understanding the molecular basis of MK2-p38α signaling complex assembly: insights into protein-protein interaction by molecular dynamics and free energy studies

The formation of a p38 MAPK and MAPK-activated protein kinase 2 (MK2) signaling complex is physiologically relevant to cellular responses such as the proinflammatory cytokine production. The interaction between p38α isoform and MK2 is of great importance for this signaling. In this study, molecular dynamics simulation and binding free energy calculation were performed on the MK2-p38α signaling complex to investigate the protein-protein interaction between the two proteins. Dynamic domain motion analyses were performed to analyze the conformational changes between the unbound and bound states of proteins during the interaction. The activation loop, αF-I helices, and loops among α helices in the C-lobe of MK2 are found to be highly flexible and exhibit significant changes upon p38α binding. The results also show that after the binding of p38α, the N- and C-terminal domains of MK2 display an opening and twisting motion centered on the activation loop. The molecular mechanics Poisson-Boltzmann and generalized-Born surface area (MM-PB/GBSA) methods were used to calculate binding free energies between MK2 and p38α. The analysis of the components of binding free energy calculation indicates that the van der Waals interaction and the nonpolar solvation energy provide the driving force for the binding process, while the electrostatic interaction contributes critically to the specificity, rather than to MK2-p38α binding affinity. The contribution of each residue at the interaction interface to the binding affinity of MK2 with p38α was also analyzed by free energy decomposition. Several important residues responsible for the protein-protein interaction were also identified.     

Conformational flexibility and binding energy profile of c-Abl tyrosine kinase complexed with Imatinib: an insight from MD study

Structural biology of kinase and in particular of tyrosine kinase has given detailed insights into the intrinsic flexibility of the catalytic domain and has provided a rational basis for obtaining selective inhibitors. In this paper, we have studied the conformational flexibility of c-Abl tyrosine kinase complexed with Imatinib (STI), in the presence of TIP3P water in physiological conditions at neutral pH. The conformational studies suggest that the flexibility of activation loop is responsible to facilitate the nucleotide binding and release. Owing to the conformational adaptability, adenosine triphosphate (ATP) binds at a particular site in the loop region of the tyrosine kinase. The molecular mechanics Poisson–Boltzmann surface area methods are analysed, as is a free-energy pathways method, which shows the stable binding with free energy ? 6.04 kcal/mol for STI. The binding energy calculated by the Sietraj method is approximately the same as the experimental binding energy of STI with c-Abl kinase. It is suggested that the conserved glutamic acid and lysine residues are necessary for the stability and optimum activity of inhibitor. This study may be helpful in rational drug designing of new kinase inhibitors.     

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.     

Incorporating replacement free energy of binding‐site waters in molecular docking

Binding‐site water molecules play a crucial role in protein‐ligand recognition, either being displaced upon ligand binding or forming water bridges to stabilize the complex. However, rigorously treating explicit binding‐site waters is challenging in molecular docking, which requires to fully sample ensembles of waters and to consider the free energy cost of replacing waters. Here, we describe a method to incorporate structural and energetic properties of binding‐site waters into molecular docking. We first developed a solvent property analysis (SPA) program to compute the replacement free energies of binding‐site water molecules by post‐processing molecular dynamics trajectories obtained from ligand‐free protein structure simulation in explicit water. Next, we implemented a distance‐dependent scoring term into DOCK scoring function to take account of the water replacement free energy cost upon ligand binding. We assessed this approach in protein targets containing important binding‐site waters, and we demonstrated that our approach is reliable in reproducing the crystal binding geometries of protein‐ligand‐water complexes, as well as moderately improving the ligand docking enrichment performance. In addition, SPA program (free available to academic users upon request) may be applied in identifying hot‐spot binding‐site residues and structure‐based lead optimization. Proteins 2014; 82:1765–1776. © 2014 Wiley Periodicals, Inc.     

Imatinib can act as an Allosteric Activator of Abl Kinase

Imatinib is an ATP-competitive inhibitor of Bcr-Abl kinase and the first drug approved for chronic myelogenous leukemia (CML) treatment. Here we show that imatinib binds to a secondary, allosteric site located in the myristoyl pocket of Abl to function as an activator of the kinase activity. Abl transitions between an assembled, inhibited state and an extended, activated state. The equilibrium is regulated by the conformation of the αΙ helix, which is located nearby the allosteric pocket. Imatinib binding to the allosteric pocket elicits an αΙ helix conformation that is not compatible with the assembled state, thereby promoting the extended state and stimulating the kinase activity. Although in wild-type Abl the catalytic pocket has a much higher affinity for imatinib than the allosteric pocket does, the two binding affinities are comparable in Abl variants carrying imatinib-resistant mutations in the catalytic site. A previously isolated imatinib-resistant mutation in patients appears to be mediating its function by increasing the affinity of imatinib for the allosteric pocket, providing a hitherto unknown mechanism of drug resistance. Our results highlight the benefit of combining imatinib with allosteric inhibitors to maximize their inhibitory effect on Bcr-Abl.     

Computational screening of fatty acid synthase inhibitors against thioesterase domain

Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein–ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.     

Salt bridges: geometrically specific, designable interactions

Salt bridges occur frequently in proteins, providing conformational specificity and contributing to molecular recognition and catalysis. We present a comprehensive analysis of these interactions in protein structures by surveying a large database of protein structures. Salt bridges between Asp or Glu and His, Arg, or Lys display extremely well-defined geometric preferences. Several previously observed preferences are confirmed, and others that were previously unrecognized are discovered. Salt bridges are explored for their preferences for different separations in sequence and in space, geometric preferences within proteins and at protein-protein interfaces, co-operativity in networked salt bridges, inclusion within metal-binding sites, preference for acidic electrons, apparent conformational side chain entropy reduction on formation, and degree of burial. Salt bridges occur far more frequently between residues at close than distant sequence separations, but, at close distances, there remain strong preferences for salt bridges at specific separations. Specific types of complex salt bridges, involving three or more members, are also discovered. As we observe a strong relationship between the propensity to form a salt bridge and the placement of salt-bridging residues in protein sequences, we discuss the role that salt bridges might play in kinetically influencing protein folding and thermodynamically stabilizing the native conformation. We also develop a quantitative method to select appropriate crystal structure resolution and B-factor cutoffs. Detailed knowledge of these geometric and sequence dependences should aid de novo design and prediction algorithms.     

Binding of the proline-rich segment of myelin basic protein to SH3 domains: spectroscopic, microarray, and modeling studies of ligand conformation and effects of posttranslational modifications

Myelin basic protein (MBP) is a multifunctional protein involved in maintaining the stability and integrity of the myelin sheath by a variety of interactions with membranes and with cytoskeletal and other proteins. A central segment of MBP is highly conserved in mammals and consists of a membrane surface-associated amphipathic alpha-helix, immediately followed by a proline-rich segment that we hypothesize is an SH3 ligand. We show by circular dichroic spectroscopy that this proline-rich segment forms a polyproline type II helix in vitro under physiological conditions and that phosphorylation at a constituent threonyl residue has a stabilizing effect on its conformation. Using SH3 domain microarrays, we observe that the unmodified recombinant murine 18.5 kDa MBP isoform (rmC1 component) binds the following SH3 domains: Yes1 > PSD95 > cortactin = PexD = Abl = Fyn = c-Src = Itk in order of decreasing affinity. A quasi-deiminated form of the protein (rmC8) binds the SH3 domains Yes1 > Fyn > cortactin = c-Src > PexD = Abl. Phosphorylation of rmC1 at 1-2 threonines within the proline-rich segment by mitogen-activated protein kinase in vitro has no effect on the binding specificity to the SH3 domains on the array. An SH3 domain of chicken Fyn is also demonstrated to bind to lipid membrane-associated C1, phosphorylated C1, and rmC8. Molecular docking simulations of the interaction of the putative SH3 ligand of classic MBP with the human Fyn SH3 domain indicate that the strength of the interaction is of the same order of magnitude as with calmodulin and that the molecular recognition and association is mediated by some weak CH...pi interactions between the ligand prolyl residues and the aromatic ones of the SH3 binding site. One such interaction is well-conserved and involves the stacking of an MBP-peptide prolyl and an SH3 domain tryptophanyl residue, as in most other SH3-ligand complexes. Lysyl and arginyl residues in the peptide canonically interact via salt bridges and cation-pi interactions with negatively charged and aromatic residues in the SH3 domain binding site. Posttranslational modifications (phosphorylation or methylation) of the ligand cause noticeable shifts in the conformation of the flexible peptide and its side chains but do not predict any major inhibition of the binding beyond somewhat less favorable interactions for peptides with phosphorylated seryl or threonyl residues.     

Mutation D816V alters the internal structure and dynamics of c-KIT receptor cytoplasmic region: implications for dimerization and activation mechanisms

The type III receptor tyrosine kinase (RTK) KIT plays a crucial role in the transmission of cellular signals through phosphorylation events that are associated with a switching of the protein conformation between inactive and active states. D816V KIT mutation is associated with various pathologies including mastocytosis and cancers. D816V-mutated KIT is constitutively active, and resistant to treatment with the anti-cancer drug Imatinib. To elucidate the activating molecular mechanism of this mutation, we applied a multi-approach procedure combining molecular dynamics (MD) simulations, normal modes analysis (NMA) and binding site prediction. Multiple 50-ns MD simulations of wild-type KIT and its mutant D816V were recorded using the inactive auto-inhibited structure of the protein, characteristic of type III RTKs. Computed free energy differences enabled us to quantify the impact of D816V on protein stability in the inactive state. We evidenced a local structural alteration of the activation loop (A-loop) upon mutation, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. A thorough normal mode analysis of several MD conformations led to a plausible molecular rationale to propose that JMR is able to depart its auto-inhibitory position more easily in the mutant than in wild-type KIT and is thus able to promote kinase mutant dimerization without the need for extra-cellular ligand binding. Pocket detection at the surface of NMA-displaced conformations finally revealed that detachment of JMR from the kinase domain in the mutant was sufficient to open an access to the catalytic and substrate binding sites.     

Multiscale computational study of ligand binding pathways: Case of p38 MAP kinase and its inhibitors

Protein kinases are one of the most important drug targets in the past 10 years. Understanding the inhibitor association processes will profoundly impact new binder designs with preferred binding kinetics. However, after more than a decade of effort, a complete atomistic-level study of kinase inhibitor binding pathways is still lacking. As all kinases share a similar scaffold, we used p38 kinase as a model system to investigate the conformational dynamics and free energy transition of inhibitor binding toward kinases. Two major kinase conformations, Asp-Phe-Gly (DFG)-in and DFG-out, and three types of inhibitors, type I, II, and III, were thoroughly investigated in this work. We performed Brownian dynamics simulations and up to 340 μs Gaussian-accelerated molecular dynamics simulations to capture the inhibitor binding paths and a series of conformational transitions of the p38 kinase from its apo to inhibitor-bound form. Eighteen successful binding trajectories, including all types of inhibitors, are reported herein. Our simulations suggest a mechanism of inhibitor recruitment, a faster ligand association step to a pre-existing DFG-in/DFG-out p38 protein, followed by a slower molecular rearrangement step to adjust the protein-ligand conformation followed by a shift in the energy landscape to reach the final bound state. The ligand association processes also reflect the energetic favor of type I and type II/III inhibitor binding through ATP and allosteric channels, respectively. These different binding routes are directly responsible for the fast (type I binders) and slow (type II/III binders) kinetics of different types of p38 inhibitors. Our findings also echo the recent study of p38 inhibitor dissociation, implying that ligand unbinding could undergo a reverse path of binding, and both processes share similar metastates. This study deepens the understanding of molecular and energetic features of kinase inhibitor-binding processes and will inspire future drug development from a kinetic point of view.     

Allosteric interactions between the myristate- and ATP-site of the Abl kinase

Abl kinase inhibitors targeting the ATP binding pocket are currently employed as potent anti-leukemogenic agents but drug resistance has become a significant clinical limitation. Recently, a compound that binds to the myristate pocket of Abl (GNF-5) was shown to act cooperatively with nilotinib, an ATP-competitive inhibitor to target the recalcitrant “T315I” gatekeeper mutant of Bcr-Abl. To uncover an explanation for how drug binding at a distance from the kinase active site could lead to inhibition and how inhibitors could combine their effects, hydrogen exchange mass spectrometry (HX MS) was employed to monitor conformational effects in the presence of both dasatinib, a clinically approved ATP-site inhibitor, and GNF-5. While dasatinib binding to wild type Abl clearly influenced Abl conformation, no binding was detected between dasatinib and T315I. GNF-5, however, elicited the same conformational changes in both wild type and T315I, including changes to dynamics within the ATP site located approximately 25 Å from the site of GNF-5 interaction. Simultaneous binding of dasatinib and GNF-5 to T315I caused conformational and/or dynamics changes in Abl such that effects of dasatinib on T315I were the same as when it bound to wild type Abl. These results provide strong biophysical evidence that allosteric interactions play a role in Abl kinase downregulation and that targeting sites outside the ATP binding site can provide an important pharmacological tool to overcome mutations that cause resistance to ATP-competitive inhibitors.     

Abi enhances Abl-mediated Cdc2 phosphorylation and inactivation

Abelson tyrosine kinase (Abl) is a non-receptor tyrosine kinase which is frequently coupled with adaptor proteins to interact with its substrates for the regulation of cytoskeleton rearrangement, cell growth and apoptosis in response to a variety of biological stimuli. The Abl interactor (Abi) family members were first identified as adaptor proteins of Abl for regulating Abl transforming and kinase activity. In the present study, we used a yeast two-hybrid screen to identify Cdc2 as a novel Abi-binding protein. This finding led us to investigate the role of Abi in linking Abl and Cdc2. These three proteins formed a trimeric complex inDrosophila and mammalian cells. The expression of Abi in cells greatly enhanced the formation of the Abl-Cdc2 complex, suggesting that Abi functions as an adaptor protein facilitating the binding between Abl and Cdc2. We show that Abi promotes Abl-mediated phosphorylation of Cdc2 at tyrosine 15 and inactivation of Cdc2 kinase activity. Furthermore, coexpression of Abl and Abi inDrosophila S2 cells led to suppression of cell growth. These data suggest that Abl signaling may be involved in the downregulation of Cdc2 kinase in cell cycle control.     

A novel mode of Gleevec binding is revealed by the structure of spleen tyrosine kinase

Spleen tyrosine kinase (Syk) is a non-receptor tyrosine kinase required for signaling from immunoreceptors in various hematopoietic cells. Phosphorylation of two tyrosine residues in the activation loop of the Syk kinase catalytic domain is necessary for signaling, a phenomenon typical of tyrosine kinase family members. Syk in vitro enzyme activity, however, does not depend on phosphorylation (activation loop tyrosine --> phenylalanine mutants retain catalytic activity). We have determined the x-ray structure of the unphosphorylated form of the kinase catalytic domain of Syk. The enzyme adopts a conformation of the activation loop typically seen only in activated, phosphorylated tyrosine kinases, explaining why Syk does not require phosphorylation for activation. We also demonstrate that Gleevec (STI-571, Imatinib) inhibits the isolated kinase domains of both unphosphorylated Syk and phosphorylated Abl with comparable potency. Gleevec binds Syk in a novel, compact cis-conformation that differs dramatically from the binding mode observed with unphosphorylated Abl, the more Gleevec-sensitive form of Abl. This finding suggests the existence of two distinct Gleevec binding modes: an extended, trans-conformation characteristic of tight binding to the inactive conformation of a protein kinase and a second compact, cis-conformation characteristic of weaker binding to the active conformation. Finally, the Syk-bound cis-conformation of Gleevec bears a striking resemblance to the rigid structure of the nonspecific, natural product kinase inhibitor staurosporine.