c-Src binds to the cancer drug imatinib with an inactive Abl/c-Kit conformation and a distributed thermodynamic penalty

Sequence determinants of a specific inactive protein kinase conformation

Only a small percentage of protein kinases have been shown to adopt a distinct inactive ATP-binding site conformation, called the Asp-Phe-Gly-out (DFG-out) conformation. Given the high degree of homology within this enzyme family, we sought to understand the basis of this disparity on a sequence level. We identified two residue positions that sensitize mitogen-activated protein kinases (MAPKs) to inhibitors that stabilize the DFG-out inactive conformation. After characterizing the structure and dynamics of an inhibitor-sensitive MAPK mutant, we demonstrated the generality of this strategy by sensitizing a kinase (apoptosis signal-regulating kinase 1) not in the MAPK family to several DFG-out stabilizing ligands, using the same residue positions. The use of specific inactive conformations may aid the study of noncatalytic roles of protein kinases, such as binding partner interactions and scaffolding effects.

Conformation-selective inhibitors reveal differences in the activation and phosphate-binding loops of the tyrosine kinases Abl and Src

Over the past decade, an increasingly diverse array of potent and selective inhibitors that target the ATP-binding sites of protein kinases have been developed. Many of these inhibitors, like the clinically approved drug imatinib (Gleevec), stabilize a specific catalytically inactive ATP-binding site conformation of their kinases targets. Imatinib is notable in that it is highly selective for its kinase target, Abl, over other closely related tyrosine kinases, such as Src. In addition, imatinib is highly sensitive to the phosphorylation state of Abl's activation loop, which is believed to be a general characteristic of all inhibitors that stabilize a similar inactive ATP-binding site conformation. In this report, we perform a systematic analysis of a diverse series of ATP-competitive inhibitors that stabilize a similar inactive ATP-binding site conformation as imatinib with the tyrosine kinases Src and Abl. In contrast to imatinib, many of these inhibitors have very similar potencies against Src and Abl. Furthermore, only a subset of this class of inhibitors is sensitive to the phosphorylation state of the activation loop of these kinases. In attempting to explain this observation, we have uncovered an unexpected correlation between Abl's activation loop and another flexible active site feature, called the phosphate-binding loop (p-loop). These studies shed light on how imatinib is able to obtain its high target selectivity and reveal how the conformational preference of flexible active site regions can vary between closely related kinases.

iPfam: a database of protein family and domain interactions found in the Protein Data Bank

The database iPfam, available at http://ipfam.org, catalogues Pfam domain interactions based on known 3D structures that are found in the Protein Data Bank, providing interaction data at the molecular level. Previously, the iPfam domain–domain interaction data was integrated within the Pfam database and website, but it has now been migrated to a separate database. This allows for independent development, improving data access and giving clearer separation between the protein family and interactions datasets. In addition to domain–domain interactions, iPfam has been expanded to include interaction data for domain bound small molecule ligands. Functional annotations are provided from source databases, supplemented by the incorporation of Wikipedia articles where available. iPfam (version 1.0) contains >9500 domain–domain and 15 500 domain–ligand interactions. The new website provides access to this data in a variety of ways, including interactive visualizations of the interaction data.

Computational analysis of the binding specificity of Gleevec to Abl, c-Kit, Lck, and c-Src tyrosine kinases

Gleevec, a well-known cancer therapeutic agent, is an effective inhibitor of several tyrosine kinases, including Abl and c-Kit, but displays less potency to inhibit closely homologous tyrosine kinases, such as Lck and c-Src. Because many structural features of the binding site are highly conserved in these homologous kinases, the molecular determinants responsible for the binding specificity of Gleevec remain poorly understood. To address this issue, free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent was used to compute the binding affinity of Gleevec to Abl, c-Kit, Lck, and c-Src. The results of the FEP/MD calculations are in good agreement with experiments, enabling a detailed and quantitative dissection of the absolute binding free energy in terms of various thermodynamic contributions affecting the binding specificity of Gleevec to the kinases. Dominant binding free energy contributions arises from the van der Waals dispersive interaction, compensating about two-thirds of the unfavorable free energy penalty associated with the loss of translational, rotational, and conformational freedom of the ligand upon binding. In contrast, the contributions from electrostatic and repulsive interactions nearly cancel out due to solvent effects. Furthermore, the calculations show the importance of the conformation of the kinase activation loop. Among the kinases examined, Abl provides the most favorable binding environment for Gleevec via optimal protein-ligand interactions and a small free energy cost for loss of the translational, rotational, and conformational freedom upon ligand binding. The FEP/MD calculations additionally reveal that Lck and c-Src provide similar nonbinding interactions with the bound-Gleevec, but the former pays less entropic penalty for the ligand losing its translational, rotational, and conformational motions to bind, examining the empirically observed differential binding affinities of Gleevec between the two Src-family kinases.

Activation state-selective kinase inhibitor assay based on ion mobility-mass spectrometry

The discovery of activation state dependent kinase inhibitors, which bind specifically to the inactive conformation of the protein, is considered to be a promising pathway to improved cancer treatments. Identifying such inhibitors is challenging, however, because they can have Kd values similar to molecules known to inhibit kinase function by interacting with the active form. Further, while inhibitor induced changes within the kinase tertiary structure are significant, few technologies are able to correctly assign inhibitor binding modes in a high-throughput fashion based exclusively on protein-inhibitor complex formation and changes in local protein structure. We have developed a new assay, using ion mobility-mass spectrometry, capable of both rapidly detecting inhibitor binding and classifying the resultant kinase binding modes. Here, we demonstrate the ability of our approach to classify a broad set of kinase inhibitors, using micrograms of protein, without the need for protein modification or tagging.

Structural basis for activation of ZAP-70 by phosphorylation of the SH2-kinase linker

Serial activation of the tyrosine kinases Lck and ZAP-70 initiates signaling downstream of the T cell receptor. We previously reported the structure of an autoinhibited ZAP-70 variant in which two regulatory tyrosine residues (315 and 319) in the SH2-kinase linker were replaced by phenylalanine. We now present a crystal structure of ZAP-70 in which Tyr 315 and Tyr 319 are not mutated, leading to the recognition of a five-residue sequence register error in the SH2-kinase linker of the original crystallographic model. The revised model identifies distinct roles for these two tyrosines. As seen in a recently reported structure of the related tyrosine kinase Syk, Tyr 315 of ZAP-70 is part of a hydrophobic interface between the regulatory apparatus and the kinase domain, and the integrity of this interface would be lost upon engagement of doubly phosphorylated peptides by the SH2 domains. Tyr 319 is not necessarily dislodged by SH2 engagement, which activates ZAP-70 only ∼5-fold in vitro. In contrast, phosphorylation by Lck activates ZAP-70 ∼100-fold. This difference is due to the ability of Tyr 319 to suppress ZAP-70 activity even when the SH2 domains are dislodged from the kinase domain, providing stringent control of ZAP-70 activity downstream of Lck.

Targeting the unactivated conformations of protein kinases for small molecule drug discovery

BACKGROUND

The number of drugs in active clinical development or on the market that target the unactivated conformational states of protein kinases is growing and represents a significant portion of kinase research at biopharmaceutical companies. These non-classical kinase inhibitors have a mode of action which may overcome some of the liabilities of classical ATP-site inhibitors that substantially overlap the space that ATP occupies in the activated kinase.

OBJECTIVE

This review will discuss state-of-the-art methods of inhibiting protein kinases by targeting the unactivated conformations of the enzyme with small molecules directed to the ATP binding region.

METHODS

Biochemical and structural biology publications and public domain crystal structures were evaluated to identify key concepts in drug discovery for unactivated protein kinase inhibitors that target the ATP binding region.

CONCLUSION

The potential for enhanced selectivity, potency and duration of pharmacological action may allow non-classical kinase therapeutics to be used for chronic dosing in non-life-threatening indications. Moreover, by targeting additional conformational space on the kinase protein it is possible that new chemical matter will be discovered such that current intellectual property limitations on traditional ATP-site chemical scaffolds may be circumvented.

Chemogenomics in drug discovery: computational methods based on the comparison of binding sites

Novel computational methods for understanding relationships between ligands and all possible biological targets have emerged in recent years. Proteins are connected to each other based on the similarity of their ligands or based on the similarity of their binding sites. The assumption is that compounds sharing chemical similarity should share targets and that targets with a similar binding site should also share ligands. A large number of computational techniques have been developed to assess ligand and binding site similarity, which can be used to make quantitative predictions of the most probable biological target of a given compound. This review covers the recent advances in new computational methods for relating biological targets based on the similarity of their binding sites. Binding site comparisons are used for the prediction of their most likely ligands, their possible cross reactivity and selectivity. These comparisons can also be used to infer the function of novel uncharacterized proteins.

Fragment-based discovery of new highly substituted 1H-pyrrolo[2,3-b]- and 3H-imidazolo[4,5-b]-pyridines as focal adhesion kinase inhibitors

Focal adhesion kinase (FAK) is considered as an attractive target for oncology, and small-molecule inhibitors are reported to be in clinical testing. In a surface plasmon resonance (SPR)-mediated fragment screening campaign, we discovered bicyclic scaffolds like 1H-pyrazolo[3,4-d]pyrimidines binding to the hinge region of FAK. By an accelerated knowledge-based fragment growing approach, essential pharmacophores were added. The establishment of highly substituted unprecedented 1H-pyrrolo[2,3-b]pyridine derivatizations provided compounds with submicromolar cellular FAK inhibition potential. The combination of substituents on the bicyclic templates and the nature of the core structure itself have a significant impact on the compounds FAK selectivity. Structural analysis revealed that the appropriately substituted pyrrolo[2,3-b]pyridine induced a rare helical DFG-loop conformation. The discovered synthetic route to introduce three different substituents independently paves the way for versatile applications of the 7-azaindole core.

Explaining why Gleevec is a specific and potent inhibitor of Abl kinase

Tyrosine kinases present attractive drug targets for specific types of cancers. Gleevec, a well-known therapeutic agent against chronic myelogenous leukemia, is an effective inhibitor of Abl tyrosine kinase. However, Gleevec fails to inhibit closely homologous tyrosine kinases, such as c-Src. Because many structural features of the binding site are conserved, the molecular determinants responsible for binding specificity are not immediately apparent. Some have attributed the difference in binding specificity of Gleevec to subtle variations in ligand-protein interactions (binding affinity control), whereas others have proposed that it is the conformation of the DFG motif, in which ligand binding is only accessible to Abl and not to c-Src (conformational selection control). To address this issue, the absolute binding free energy was computed using all-atom molecular dynamics simulations with explicit solvent. The results of the free energy simulations are in good agreement with experiments, thereby enabling a meaningful decomposition of the binding free energy to elucidate the factors controlling Gleevec's binding specificity. The latter is shown to be controlled by a conformational selection mechanism and also by differences in key van der Waals interactions responsible for the stabilization of Gleevec in the binding pocket of Abl.

Structure and dynamic regulation of Abl kinases

The c-abl proto-oncogene encodes a unique protein-tyrosine kinase (Abl) distinct from c-Src, c-Fes, and other cytoplasmic tyrosine kinases. In normal cells, Abl plays prominent roles in cellular responses to genotoxic stress as well as in the regulation of the actin cytoskeleton. Abl is also well known in the context of Bcr-Abl, the oncogenic fusion protein characteristic of chronic myelogenous leukemia. Selective inhibitors of Bcr-Abl, of which imatinib is the prototype, have had a tremendous impact on clinical outcomes in chronic myelogenous leukemia and revolutionized the field of targeted cancer therapy. In this minireview, we focus on the structural organization and dynamics of Abl kinases and how these features influence inhibitor sensitivity.

The ins and outs of bcr-abl inhibition

The development of inhibitors against Abl has changed the landscape for the treatment of chronic myelogenous leukemia (CML) and cancer in general. Beginning with the monumental discovery and approval of imatinib for CML, a second generation of inhibitors, nilotinib and dasatinib, has now gained approval for the treatment of CML. Notably, these second-generation inhibitors are active against many of the mutations in the Abl kinase that confer resistance to imatinib. However, resistance remains a major problem, and new inhibitors such as ponatinib and GNF2/GNF5 have been developed, with activity towards the common gatekeeper T315I mutation. We review here the mechanisms of Abl inhibition with an emphasis on structural elements that are important for the selectivity and design of new molecules. In particular, we focus on how changes in the conformation of the P-loop, the activation loop, the DFG motif, and other structural elements of Abl have been instrumental in developing an understanding of inhibitor binding.

Identification of pyrimidine derivatives as hSMG-1 inhibitors

hSMG-1 kinase plays a dual role in a highly conserved RNA surveillance pathway termed nonsense-mediated RNA decay (NMD) and in cellular genotoxic stress response. Since deregulation of cellular responses to stress contributes to tumor growth and resistance to chemotherapy, hSMG-1 is a potential target for cancer treatment. From our screening efforts, we have identified pyrimidine derivatives as hSMG-1 kinase inhibitors. We report structure-based optimization of this pan-kinase scaffold to improve its biochemical profile and overall kinome selectivity, including mTOR and CDK, to generate the first reported selective hSMG-1 tool compound.

Virtual target screening: validation using kinase inhibitors

Computational methods involving virtual screening could potentially be employed to discover new biomolecular targets for an individual molecule of interest (MOI). However, existing scoring functions may not accurately differentiate proteins to which the MOI binds from a larger set of macromolecules in a protein structural database. An MOI will most likely have varying degrees of predicted binding affinities to many protein targets. However, correctly interpreting a docking score as a hit for the MOI docked to any individual protein can be problematic. In our method, which we term "Virtual Target Screening (VTS)", a set of small drug-like molecules are docked against each structure in the protein library to produce benchmark statistics. This calibration provides a reference for each protein so that hits can be identified for an MOI. VTS can then be used as tool for: drug repositioning (repurposing), specificity and toxicity testing, identifying potential metabolites, probing protein structures for allosteric sites, and testing focused libraries (collection of MOIs with similar chemotypes) for selectivity. To validate our VTS method, twenty kinase inhibitors were docked to a collection of calibrated protein structures. Here, we report our results where VTS predicted protein kinases as hits in preference to other proteins in our database. Concurrently, a graphical interface for VTS was developed.

Development of a highly selective c-Src kinase inhibitor

Generating highly selective probes to interrogate protein kinase function in biological studies remains a challenge, and new strategies are required. Herein, we describe the development of the first highly selective and cell-permeable inhibitor of c-Src, a key signaling kinase in cancer. Our strategy involves extension of traditional inhibitor design by appending functionality proposed to interact with the phosphate-binding loop of c-Src. Using our selective inhibitor, we demonstrate that selective inhibition is significantly more efficacious than pan-kinase inhibition in slowing the growth of cancer cells. We also show that inhibition of c-Abl kinase, an off-target of most c-Src inhibitors, promotes oncogenic cell growth.

An amino-indazole scaffold with spectrum selective kinase inhibition of FLT3, PDGFRα and kit

Here we describe the synthesis and characterization of a number of 3-amino-1H-indazol-6-yl-benzamides that were designed to target the 'DFG-out' conformation of the kinase activation loop. Several compounds such as 4 and 11 exhibit single-digit nanomolar EC(50)s against FLT3, c-Kit and the gatekeeper T674M mutant of PDGFRα.

Targeted kinase selectivity from kinase profiling data

Kinase selectivity plays a major role in the design strategy of lead series and in the ultimate success of kinase drug discovery programs. Although profiling compounds against a large panel of protein kinases has become a standard part of modern drug discovery, data accumulated from these kinase panels may be underutilized for new kinase projects. We present a method that can be used to optimize the selectivity profile of a compound using historical kinase profiling data. This method proposes chemical transformations based on pairs of very similar compounds, which are both active against a desired target kinase and differ in activity against another kinase. We show that these transformations are transferable across scaffolds, thus making this tool valuable to exploit kinase profiling data for unrelated series of compounds.

Structure, regulation, signaling, and targeting of abl kinases in cancer

Abl kinases are prototypic cytoplasmic tyrosine kinases and are involved in a variety of chromosomal aberrations in different cancers. This causes the expression of Abl fusion proteins, such as Bcr-Abl, that are constitutively activated and drivers of tumorigenesis. Over the past decades, biochemical and functional studies on the molecular mechanisms of Abl regulation have gone hand in hand with progression of our structural understanding of autoinhibited and active Abl conformations. In parallel, Abl oncoproteins have become prime molecular targets for cancer therapy, using adenosine triphosphate (ATP)-competitive kinase inhibitors, such as imatinib. Abl-targeting drugs serve as a paradigm for our understanding of kinase inhibitor action, specificity, and resistance development. In this review article, I will review the molecular mechanisms that are responsible for the regulation of Abl kinase activity and how oncogenic Abl fusions signal. Furthermore, past and ongoing efforts to target Abl oncoproteins using ATP-competitive and allosteric inhibitors, as well as future possibilities using combination therapy, will be discussed.

Structural and spectroscopic analysis of the kinase inhibitor bosutinib and an isomer of bosutinib binding to the Abl tyrosine kinase domain

Chronic myeloid leukemia (CML) is caused by the kinase activity of the BCR-Abl fusion protein. The Abl inhibitors imatinib, nilotinib and dasatinib are currently used to treat CML, but resistance to these inhibitors is a significant clinical problem. The kinase inhibitor bosutinib has shown efficacy in clinical trials for imatinib-resistant CML, but its binding mode is unknown. We present the 2.4 Å structure of bosutinib bound to the kinase domain of Abl, which explains the inhibitor's activity against several imatinib-resistant mutants, and reveals that similar inhibitors that lack a nitrile moiety could be effective against the common T315I mutant. We also report that two distinct chemical compounds are currently being sold under the name "bosutinib", and report spectroscopic and structural characterizations of both. We show that the fluorescence properties of these compounds allow inhibitor binding to be measured quantitatively, and that the infrared absorption of the nitrile group reveals a different electrostatic environment in the conserved ATP-binding sites of Abl and Src kinases. Exploiting such differences could lead to inhibitors with improved selectivity.

Highly specific, bisubstrate-competitive Src inhibitors from DNA-templated macrocycles

Protein kinases are attractive therapeutic targets, but their high sequence and structural conservation complicates the development of specific inhibitors. We recently discovered from a DNA-templated macrocycle library inhibitors with unusually high selectivity among Src-family kinases. Starting from these compounds, we developed and characterized in molecular detail potent macrocyclic inhibitors of Src kinase and its cancer-associated gatekeeper mutant. We solved two co-crystal structures of macrocycles bound to Src kinase. These structures reveal the molecular basis of the combined ATP- and substrate peptide-competitive inhibitory mechanism and the remarkable kinase specificity of the compounds. The most potent compounds inhibit Src activity in cultured mammalian cells. Our work establishes that macrocycles can inhibit protein kinases through a bi-substrate competitive mechanism with high potency and exceptional specificity, reveals the precise molecular basis for their desirable properties, and provides new insights into the development of Src-specific inhibitors with potential therapeutic relevance.

The different flexibility of c-Src and c-Abl kinases regulates the accessibility of a druggable inactive conformation

c-Src and c-Abl are two closely related protein kinases that constitute important anticancer targets. Despite their high sequence identity, they show different sensitivities to the anticancer drug imatinib, which binds specifically to a particular inactive conformation in which the Asp of the conserved DFG motif points outward (DFG-out). We have analyzed the DFG conformational transition of the two kinases using massive molecular dynamics simulations, free energy calculations, and isothermal titration calorimetry. On the basis of the reconstruction of the free energy surfaces for the DFG-in to DFG-out conformational changes of c-Src and c-Abl, we propose that the different flexibility of the two kinases results in a different stability of the DFG-out conformation and might be the main determinant of imatinib selectivity.

Targeting cancer with small-molecular-weight kinase inhibitors

Protein and lipid kinases fulfill essential roles in many signaling pathways that regulate normal cell functions. Deregulation of these kinase activities lead to a variety of pathologies ranging from cancer to inflammatory diseases, diabetes, infectious diseases, cardiovascular disorders, cell growth and survival. 518 protein kinases and about 20 lipid-modifying kinases are encoded by the human genome, and a much larger proportion of additional kinases are present in parasite, bacterial, fungal, and viral genomes that are susceptible to exploitation as drug targets. Since many human diseases result from overactivation of protein and lipid kinases due to mutations and/or overexpression, this enzyme class represents an important target for the pharmaceutical industry. Approximately one third of all protein targets under investigation in the pharmaceutical industry are protein or lipid kinases.The kinase inhibitors that have been launched, thus far, are mainly in oncology indications and are directed against a handful of protein and lipid kinases. With one exception, all of these registered kinase inhibitors are directed toward the ATP-site and display different selectivities, potencies, and pharmacokinetic properties. At present, about 150 kinase-targeted drugs are in clinical development and many more in various stages of preclinical development. Kinase inhibitor drugs that are in clinical trials target all stages of signal transduction from the receptor protein tyrosine kinases that initiate intracellular signaling, through second-messenger-dependent lipid and protein kinases, and protein kinases that regulate the cell cycle.This review provides an insight into protein and lipid kinase drug discovery with respect to achievements, binding modes of inhibitors, and novel avenues for the generation of second-generation kinase inhibitors to treat cancers.

The evolution of protein kinase inhibitors from antagonists to agonists of cellular signaling

Kinases are highly regulated enzymes with diverse mechanisms controlling their catalytic output. Over time, chemical discovery efforts for kinases have produced ATP-competitive compounds, allosteric regulators, irreversible binders, and highly specific inhibitors. These distinct classes of small molecules have revealed many novel aspects about kinase-mediated signaling, and some have progressed from simple tool compounds into clinically validated therapeutics. This review explores several small-molecule inhibitors for kinases highlighting elaborate mechanisms by which kinase function is modulated. A complete surprise of targeted kinase drug discovery has been the finding of ATP-competitive inhibitors that behave as agonists, rather than antagonists, of their direct kinase target. These studies hint at a connection between ATP-binding site occupancy and networks of communication that are independent of kinase catalysis. Indeed, kinase inhibitors that induce changes in protein localization, protein-protein interactions, and even enhancement of catalytic activity of the target kinase have been found. The relevance of these findings to the therapeutic efficacy of kinase inhibitors and to the future identification of new classes of drug targets is discussed.

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.

Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.

Backbone assignment of the tyrosine kinase Src catalytic domain in complex with imatinib

The Src tyrosine kinase is the paradigm of an oncogenic kinase, and of regulation by intramolecular inhibitory interactions, as well as an important anticancer target due to its roles in cell proliferation and metastasis. The assignment of backbone (1)H(N), (13)C(α), (13)CO, and (15)N, and sidechain (13)C(β) resonances of the catalytic domain of Src (283 residues) in complex with the anticancer drug Imatinib is reported here. Consistent with previous X-ray studies of the same complex, most signals from the activation loop are not detected, indicating that, even in the presence of the drug, it probably adopts highly heterogeneous conformations in intermediate exchange. For the rest of the polypeptide backbone, assignments have been completed for ~88% of residues, with only a few solvent-exposed amides remaining unassigned.

Protein kinase biochemistry and drug discovery

Protein kinases are fascinating biological catalysts with a rapidly expanding knowledge base, a growing appreciation in cell regulatory control, and an ascendant role in successful therapeutic intervention. To better understand protein kinases, the molecular underpinnings of phosphoryl group transfer, protein phosphorylation, and inhibitor interactions are examined. This analysis begins with a survey of phosphate group and phosphoprotein properties which provide context to the evolutionary selection of phosphorylation as a central mechanism for biological regulation of most cellular processes. Next, the kinetic and catalytic mechanisms of protein kinases are examined with respect to model aqueous systems to define the elements of catalysis. A brief structural biology overview further delves into the molecular basis of catalysis and regulation of catalytic activity. Concomitant with a prominent role in normal physiology, protein kinases have important roles in the disease state. To facilitate effective kinase drug discovery, classic and emerging approaches for characterizing kinase inhibitors are evaluated including biochemical assay design, inhibitor mechanism of action analysis, and proper kinetic treatment of irreversible inhibitors. As the resulting protein kinase inhibitors can modulate intended and unintended targets, profiling methods are discussed which can illuminate a more complete range of an inhibitor's biological activities to enable more meaningful cellular studies and more effective clinical studies. Taken as a whole, a wealth of protein kinase biochemistry knowledge is available, yet it is clear that a substantial extent of our understanding in this field remains to be discovered which should yield many new opportunities for therapeutic intervention.

Predicting inactive conformations of protein kinases using active structures: conformational selection of type-II inhibitors

Protein kinases have been found to possess two characteristic conformations in their activation-loops: the active DFG-in conformation and the inactive DFG-out conformation. Recently, it has been very interesting to develop type-II inhibitors which target the DFG-out conformation and are more specific than the type-I inhibitors binding to the active DFG-in conformation. However, solving crystal structures of kinases with the DFG-out conformation remains a challenge, and this seriously hampers the application of the structure-based approaches in development of novel type-II inhibitors. To overcome this limitation, here we present a computational approach for predicting the DFG-out inactive conformation using the DFG-in active structures, and develop related conformational selection protocols for the uses of the predicted DFG-out models in the binding pose prediction and virtual screening of type-II ligands. With the DFG-out models, we predicted the binding poses for known type-II inhibitors, and the results were found in good agreement with the X-ray crystal structures. We also tested the abilities of the DFG-out models to recognize their specific type-II inhibitors by screening a database of small molecules. The AUC (area under curve) results indicated that the predicted DFG-out models were selective toward their specific type-II inhibitors. Therefore, the computational approach and protocols presented in this study are very promising for the structure-based design and screening of novel type-II kinase inhibitors.

Molecular mechanisms of drug resistance in tyrosine kinases cAbl and cKit

The inhibition of protein kinases has gained general acceptance as an effective approach to treat a wide range of cancers. However, in many cases, prolonged administration of kinase inhibitors often leads to acquired resistance, and the therapeutic effect is subsequently diminished. The wealth of recent studies using biochemical, kinetic, and structural approaches have revealed the molecular basis for the clinically observed resistance. In this review, we highlight several of the most common molecular mechanisms that lead to acquired resistance to kinase inhibitors observed with the cAbl (cellular form of the Abelson leukemia virus tyrosine kinase) and the type III receptor tyrosine kinase cKit, including a newly identified mechanism resulting from accelerated kinase activation caused by mutations in the activation loop. Strategies to overcome the loss of drug sensitivity that represents a challenge currently facing the field and the emerging approaches to circumvent resistance are discussed.

Activation state-dependent binding of small molecule kinase inhibitors: structural insights from biochemistry

Interactions between kinases and small molecule inhibitors can be activation state dependent. A detailed understanding of inhibitor binding therefore requires characterizing interactions across multiple activation states. We have systematically explored the effects of ABL1 activation loop phosphorylation and PDGFR family autoinhibitory juxtamembrane domain docking on inhibitor binding affinity. For a diverse compound set, the affinity patterns correctly classify inhibitors as having type I or type II binding modes, and we show that juxtamembrane domain docking can have dramatic negative effects on inhibitor affinity. The results have allowed us to associate ligand-induced conformational changes observed in cocrystal structures with specific energetic costs. The approach we describe enables investigation of the complex relationship between kinase activation state and compound binding affinity and should facilitate strategic inhibitor design.

Allostery and population shift in drug discovery

Proteins can exist in a large number of conformations around their native states that can be characterized by an energy landscape. The landscape illustrates individual valleys, which are the conformational substates. From the functional standpoint, there are two key points: first, all functionally relevant substates pre-exist; and second, the landscape is dynamic and the relative populations of the substates will change following allosteric events. Allosteric events perturb the structure, and the energetic strain propagates and shifts the population. This can lead to changes in the shapes and properties of target binding sites. Here we present an overview of dynamic conformational ensembles focusing on allosteric events in signaling. We propose that combining equilibrium fluctuation concepts with genomic screens could help drug discovery.

Effect of the SH3-SH2 domain linker sequence on the structure of Hck kinase

The coordination of activity in biological systems requires the existence of different signal transduction pathways that interact with one another and must be precisely regulated. The Src-family tyrosine kinases, which are found in many signaling pathways, differ in their physiological function despite their high overall structural similarity. In this context, the differences in the SH3-SH2 domain linkers might play a role for differential regulation, but the structural consequences of linker sequence remain poorly understood. We have therefore performed comparative molecular dynamics simulations of wildtype Hck and of a mutant Hck in which the SH3-SH2 domain linker is replaced by the corresponding sequence from the homologous kinase Lck. These simulations reveal that linker replacement not only affects the orientation of the SH3 domain itself, but also leads to an alternative conformation of the activation segment in the Hck kinase domain. The sequence of the SH3-SH2 domain linker thus exerts a remote effect on the active site geometry and might therefore play a role in modulating the structure of the inactive kinase or in fine-tuning the activation process itself.

Coupling of receptor conformation and ligand orientation determine graded activity

Small molecules stabilize specific protein conformations from a larger ensemble, enabling molecular switches that control diverse cellular functions. We show here that the converse also holds true: the conformational state of the estrogen receptor can direct distinct orientations of the bound ligand. 'Gain-of-allostery' mutations that mimic the effects of ligand in driving protein conformation allowed crystallization of the partial agonist ligand WAY-169916 with both the canonical active and inactive conformations of the estrogen receptor. The intermediate transcriptional activity induced by WAY-169916 is associated with the ligand binding differently to the active and inactive conformations of the receptor. Analyses of a series of chemical derivatives demonstrated that altering the ensemble of ligand binding orientations changes signaling output. The coupling of different ligand binding orientations to distinct active and inactive protein conformations defines a new mechanism for titrating allosteric signaling activity.

A molecular mechanics model for imatinib and imatinib:kinase binding

Imatinib is an important anticancer drug, which binds specifically to the Abl kinase and blocks its signalling activity. To model imatinib:protein interactions, we have developed a molecular mechanics force field for imatinib and four close analogues, which is consistent with the CHARMM force field for proteins and nucleic acids. Atomic charges and Lennard-Jones parameters were derived from a supermolecule ab initio approach. We considered the ab initio energies and geometries of a probe water molecule interacting with imatinib fragments at 32 different positions. We considered both a neutral and a protonated imatinib. The final RMS deviation between the ab initio and force field energies, averaged over both forms, was 0.2 kcal/mol. The model also reproduces the ab initio geometry and flexibility of imatinib. To apply the force field to imatinib:Abl simulations, it is also necessary to determine the most likely imatinib protonation state when it binds to Abl. This was done using molecular dynamics free energy simulations, where imatinib is reversibly protonated during a series of MD simulations, both in solution and in complex with Abl. The simulations indicate that imatinib binds to Abl in its protonated, positively-charged form. To help test the force field and the protonation prediction, we did MD free energy simulations that compare the Abl binding affinities of two imatinib analogs, obtaining good agreement with experiment. Finally, two new imatinib variants were considered, one of which is predicted to have improved Abl binding. This variant could be of interest as a potential drug.

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.Electronic supplementary material The online version of this article (doi:10.1007/s10867-010-9199-z) contains supplementary material, which is available to authorized users.

Fluorophore labeling of the glycine-rich loop as a method of identifying inhibitors that bind to active and inactive kinase conformations

Targeting protein kinases with small organic molecules is a promising strategy to regulate unwanted kinase activity in both chemical biology and medicinal chemistry research. Traditionally, kinase inhibitors are identified in activity-based screening assays using enzymatically active kinase preparations to measure the perturbation of substrate phosphorylation, often resulting in the enrichment of classical ATP competitive (Type I) inhibitors. However, addressing enzymatically incompetent kinase conformations offers new opportunities for targeted therapies and is moving to the forefront of kinase inhibitor research. Here we report the development of a new FLiK (Fluorescent Labels in Kinases) binding assay to detect small molecules that induce changes in the conformation of the glycine-rich loop. Due to cross-talk between the glycine-rich loop and the activation loop in kinases, this alternative labeling approach can also detect ligands that stabilize inactive kinase conformations, including slow-binding Type II and Type III kinase inhibitors. Protein X-ray crystallography validated the assay results and identified a novel DFG-out binding mode for a quinazoline-based inhibitor in p38alpha kinase. We also detected the high-affinity binding of a clinically relevant and specific VEGFR2 inhibitor, and we provide structural details of its binding mode in p38alpha, in which it stabilizes the DFG-out conformation. Last, we demonstrate the power of this new FLiK labeling strategy to detect the binding of Type I ligands that induce conformational changes in the glycine-rich loop as a means of gaining affinity for the target kinase. This approach may be a useful alternative to develop direct binding assays for kinases that do not adopt the DFG-out conformation while also avoiding the use of expensive kits, detection reagents, or radioactivity frequently employed with activity-based assays.

Discovery of a small-molecule type II inhibitor of wild-type and gatekeeper mutants of BCR-ABL, PDGFRalpha, Kit, and Src kinases: novel type II inhibitor of gatekeeper mutants

Many clinically validated kinases, such as BCR-ABL, c-Kit, PDGFR, and EGFR, become resistant to adenosine triphosphate-competitive inhibitors through mutation of the so-called gatekeeper amino acid from a threonine to a large hydrophobic amino acid, such as an isoleucine or methionine. We have developed a new class of adenosine triphosphate competitive inhibitors, exemplified by HG-7-85-01, which is capable of inhibiting T315I- BCR-ABL (clinically observed in chronic myeloid leukemia), T670I-c-Kit (clinically observed in gastrointestinal stromal tumors), and T674I/M-PDGFRalpha (clinically observed in hypereosinophilic syndrome). HG-7-85-01 is unique among all currently reported kinase inhibitors in having the ability to accommodate either a gatekeeper threonine, present in the wild-type forms of these kinases, or a large hydrophobic amino acid without becoming a promiscuous kinase inhibitor. The distinctive ability of HG-7-85-01 to simultaneously inhibit both wild-type and mutant forms of several kinases of clinical relevance is an important step in the development of the next generation of tyrosine kinase inhibitors.

Analysis of imatinib and sorafenib binding to p38alpha compared with c-Abl and b-Raf provides structural insights for understanding the selectivity of inhibitors targeting the DFG-out form of protein kinases

Protein kinases c-Abl, b-Raf, and p38alpha are recognized as important targets for therapeutic intervention. c-Abl and b-Raf are major targets of marketed oncology drugs Imatinib (Gleevec) and Sorafenib (Nexavar), respectively, and BIRB-796 is a p38alpha inhibitor that reached Phase II clinical trials. A shared feature of these drugs is the fact that they bind to the DFG-out forms of their kinase targets. Although the discovery of this class of kinase inhibitors has increased the level of emphasis on the design of DFG-out inhibitors, the structural determinants for their binding and stabilization of the DFG-out conformation remain unclear. To improve our understanding of these determinants, we determined cocrystal structures of Imatinib and Sorafenib with p38alpha. We also conducted a detailed analysis of Imatinib and Sorafenib binding to p38alpha in comparison with BIRB-796, including binding kinetics, binding interactions, the solvent accessible surface area (SASA) of the ligands, and stabilization of key structural elements of the protein upon ligand binding. Our results yield an improved understanding of the structural requirements for stabilizing the DFG-out form and a rationale for understanding the genesis of ligand selectivity among DFG-out inhibitors of protein kinases.

A bead-based activity screen for small-molecule inhibitors of signal transduction in chronic myelogenous leukemia cells

Chronic myelogenous leukemia is characterized by the presence of the chimeric BCR-ABL gene, which is expressed as the constitutively active Bcr-Abl kinase. Although kinase activity is directly responsible for the clinical phenotype, current diagnostic and prognostic methods focus on a genetic classification system in which molecularly distinct subcategories are used to predict patient responses to small-molecule inhibitors of the Bcr-Abl kinase. Point mutations in the kinase domain are a central factor regulating inhibitor resistance; however, compensatory signaling caused by the activation of unrelated kinases can influence inhibitor efficacy. Kinase activity profiling can be used as a complementary approach to genetic screening and allows direct screening of small-molecule inhibitors. We developed a quantitative assay to monitor tyrosine kinase activities and inhibitor sensitivities in a model of chronic myelogenous leukemia using peptide reporters covalently immobilized on Luminex beads. Kinase activity is quantified by nonlinear regression from well-specific internal standard curves. Using optimized synthetic substrates and peptides derived from native substrates as probes, we measured kinase inhibition in cell lysates by the signal transduction inhibitors imatinib and dasatinib. Taking advantage of a convenient 96-well plate format, this assay also allows a straightforward and quantitative analysis of the differential effects of ATP and inhibitors on kinase activity. This method for analyzing a focused signaling network benefits from rigorous statistical analysis and short processing times, thereby offering a powerful tool for drug discovery and clinical testing.

Crystal structure of an Aurora-A mutant that mimics Aurora-B bound to MLN8054: insights into selectivity and drug design

The production of selective protein kinase inhibitors is often frustrated by the similarity of the enzyme active sites. For this reason, it is challenging to design inhibitors that discriminate between the three Aurora kinases, which are important targets in cancer drug discovery. We have used a triple-point mutant of Aurora-A (AurAx3) which mimics the active site of Aurora-B to investigate the structural basis of MLN8054 selectivity. The bias toward Aurora-A inhibition by MLN8054 is fully recapitulated by AurAx3 in vitro. X-ray crystal structures of the complex suggest that the basis for the discrimination is electrostatic repulsion due to the T217E substitution, which we have confirmed using a single-point mutant. The activation loop of Aurora-A in the AurAx3-MLN8054 complex exhibits an unusual conformation in which Asp274 and Phe275 side chains point into the interior of the protein. There is to our knowledge no documented precedent for this conformation, which we have termed DFG-up. The sequence requirements of the DFG-up conformation suggest that it might be accessible to only a fraction of kinases. MLN8054 thus circumvents the problem of highly homologous active sites. Binding of MLN8054 to Aurora-A switches the character of a pocket within the active site from polar to a hydrophobic pocket, similar to what is observed in the structure of Aurora-A bound to a compound that induces DFG-out. We propose that targeting this pocket may be a productive route in the design of selective kinase inhibitors and describe the structural basis for the rational design of these compounds.

Molecular dynamics simulations show that conformational selection governs the binding preferences of imatinib for several tyrosine kinases

Tyrosine kinases transmit cellular signals through a complex mechanism, involving their phosphorylation and switching between inactive and active conformations. The cancer drug imatinib binds tightly to several homologous kinases, including Abl, but weakly to others, including Src. Imatinib specifically targets the inactive, so-called "DFG-out" conformation of Abl, which differs from the preferred, "DFG-in" conformation of Src in the orientation of a conserved Asp-Phe-Gly (DFG) activation loop. However, recent x-ray structures showed that Src can also adopt the DFG-out conformation and uses it to bind imatinib. The Src/Abl-binding free energy difference can thus be decomposed into two contributions. Contribution i measures the different protein-imatinib interactions when either kinase is in its DFG-out conformation. Contribution ii depends on the ability of imatinib to select or induce this conformation, i.e. on the relative stabilities of the DFG-out and DFG-in conformations of each kinase. Neither contribution has been measured experimentally. We use molecular dynamics simulations to show that contribution i is very small, 0.2 +/- 0.6 kcal/mol; imatinib interactions are very similar in the two kinases, including long range electrostatic interactions with the imatinib positive charge. Contribution ii, deduced using the experimental binding free energy difference, is much larger, 4.4 +/- 0.9 kcal/mol. Thus, conformational selection, easy in Abl, difficult in Src, underpins imatinib specificity. Contribution ii has a simple interpretation; it closely approximates the stability difference between the DFG-out and DFG-in conformations of apo-Src. Additional calculations show that conformational selection also governs the relative binding of imatinib to the kinases c-Kit and Lck. These results should help clarify the current framework for engineering kinase signaling.