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.     

Insights into the conformational switching mechanism of the human vascular endothelial growth factor receptor type 2 kinase domain

Human vascular endothelial growth factor receptor type 2 (h-VEFGR2) is a receptor tyrosine kinase involved in the angiogenesis process and regarded as an interesting target for the design of anticancer drugs. Its activation/inactivation mechanism is related to conformational changes in its cytoplasmatic kinase domain, involving first among all the αC-helix in N-lobe and the A-loop in C-lobe. Affinity of inhibitors for the active or inactive kinase form could dictate the open or closed conformation of the A-loop, thus making the different conformations of the kinase domain receptor (KDR) domain different drug targets in drug discovery. In this view, a detailed knowledge of the conformational landscape of KDR domain is of central relevance to rationalize the efficiency and selectivity of kinase inhibitors. Here, molecular dynamics simulations were used to gain insight into the conformational switching activity of the KDR domain and to identify intermediate conformations between the two limiting active and inactive conformations. Specific energy barriers have been selectively removed to induce, and hence highlight at the atomistic level, the regulation mechanism of the A-loop opening. The proposed strategy allowed to repeatedly observe the escape of the KDR domain from the DFG-out free energy basin and to identify rare intermediate conformations between the DFG-out and the DFG-in structures to be employed in a structure-based drug discovery process.     

Molecular dynamics simulation and free energy calculation studies of the binding mechanism of allosteric inhibitors with p38α MAP kinase

p38 MAP kinase is a promising target for anti-inflammatory treatment. The classical kinase inhibitors imatinib and sorafenib as well as BI-1 and BIRB-796 were reported to bind in the DFG-out form of human p38α, known as type II or allosteric kinase inhibitors. Although DFG-out conformation has attracted great interest in the design of type II kinase inhibitors, the structural requirements for binding and mechanism of stabilization of DFG-out conformation remain unclear. As allosteric inhibition is important to the selectivity of kinase inhibitor, herein the binding modes of imatinib, sorafenib, BI-1 and BIRB-796 to p38α were investigated by molecular dynamics simulation. Binding free energies were calculated by molecular mechanics/Poisson-Boltzmann surface area method. The predicted binding affinities can give a good explanation of the activity difference of the studied inhibitors. Furthermore, binding free energies decomposition analysis and further structural analysis indicate that the dominating effect of van der Waals interaction drives the binding process, and key residues, such as Lys53, Gly71, Leu75, Ile84, Thr106, Met109, Leu167, Asp168, and Phe169, play important roles by forming hydrogen bond, salt bridge, and hydrophobic interactions with the DFG-out conformation of p38α. Finally, we also conducted a detailed analysis of BI-1, imatinib, and sorafenib binding to p38α in comparison with BIRB-796 exploited for gaining potency as well as selectivity of p38 inhibitors. These results are expected to be useful for future rational design of novel type II p38 inhibitors.     

Dynamic Equilibrium of the Aurora A Kinase Activation Loop Revealed by Single‐Molecule Spectroscopy

The conformation of the activation loop (T‐loop) of protein kinases underlies enzymatic activity and influences the binding of small‐molecule inhibitors. By using single‐molecule fluorescence spectroscopy, we have determined that phosphorylated Aurora A kinase is in dynamic equilibrium between a DFG‐in‐like active T‐loop conformation and a DFG‐out‐like inactive conformation, and have measured the rate constants of interconversion. Addition of the Aurora A activating protein TPX2 shifts the equilibrium towards an active T‐loop conformation whereas addition of the inhibitors MLN8054 and CD532 favors an inactive T‐loop. We show that Aurora A binds TPX2 and MLN8054 simultaneously and provide a new model for kinase conformational behavior. Our approach will enable conformation‐specific effects to be integrated into inhibitor discovery across the kinome, and we outline some immediate consequences for structure‐based drug discovery.     

Cell-based proteome profiling of potential dasatinib targets by use of affinity-based probes

Protein kinases (PKs) play an important role in the development and progression of cancer by regulating cell growth, survival, invasion, metastasis, and angiogenesis. Dasatinib (BMS-354825), a dual Src/Abl inhibitor, is a promising therapeutic agent with oral bioavailability. It has been used for the treatment of imatinib-resistant chronic myelogenous leukemia (CML). Most kinase inhibitors, including Dasatinib, inhibit multiple cellular targets and do not possess exquisite cellular specificity. Recent efforts in kinase research thus focus on the development of large-scale, proteome-wide chemical profiling methods capable of rapid identification of potential cellular (on- and off-) targets of kinase inhibitors. Most existing approaches, however, are still problematic and in many cases not compatible with live-cell studies. In this work, we have successfully developed a cell-permeable kinase probe (DA-2) capable of proteome-wide profiling of potential cellular targets of Dasatinib. In this way, highly regulated, compartmentalized kinase-drug interactions were maintained. By comparing results obtained from different proteomic setups (live cells, cell lysates, and immobilized affinity matrix), we found DA-2 was able to identify significantly more putative kinase targets. In addition to Abl and Src family tyrosine kinases, a number of previously unknown Dasatinib targets have been identified, including several serine/threonine kinases (PCTK3, STK25, eIF-2A, PIM-3, PKA C-α, and PKN2). They were further validated by pull-down/immunoblotting experiments as well as kinase inhibition assays. Further studies are needed to better understand the exact relevance of Dasatinib and its pharmacological effects in relation to these newly identified cellular targets. The approach developed herein should be amenable to the study of many of the existing reversible drugs/drug candidates.     

Inverse screening of Simvastatin kinase targets from glioblastoma druggable kinome

The statin drug Simvastatin is a HMG-CoA reductase inhibitor that has been widely used to lower blood lipid. However, the drug is clinically observed to reposition a significant suppressing potency on glioblastoma (GBM) by unexpectedly targeting diverse kinase pathways involved in GBM tumorigensis. Here, an inverse screening strategy is described to discover potential kinase targets of Simvastatin. Various human protein kinases implicated in GBM are enriched to define a druggable kinome; the binding behavior of Simvastatin to the kinome is profiled systematically via an integrative computational approach, from which most kinases have only low or moderate binding potency to Simvastatin, while only few are identified as promising kinase hits. It is revealed that Simvastatin can potentially interact with certain known targets or key regulators of GBM such as ErbB, c-Src and FGFR signaling pathways, but exhibit low affinity to the well-established GBM target of PI3K/Akt/mTOR pathway. Further assays determine that Simvastatin can inhibit kinase hits EGFR, MET, SRC and HER2 at nanomolar level, which are comparable with those of cognate kinase inhibitors. Structural analyses reveal that the sophisticated T790 M gatekeeper mutation can considerably reduce Simvastatin sensitivity to EGFR by inducing the ligand change between different binding modes.     

Molecular insight into the T798M gatekeeper mutation-caused acquired resistance to tyrosine kinase inhibitors in ErbB2-positive breast cancer

Human epidermal growth factor receptor 2 (ErbB2) is an attractive therapeutic target for metastatic breast cancer. The kinase has been clinically observed to harbor a gatekeeper mutation T798M in its active site, which causes acquired resistance to the first-line targeted breast cancer therapy with small-molecule tyrosine kinase inhibitors. Previously, several theories have been proposed to explain the molecular mechanism of gatekeeper mutation-caused drug resistance, such as blocking of inhibitor binding and increasing of ATP affinity. In the current study, the direct binding of three wild type-selective inhibitors (Lapatinib, AEE788 and TAK-285) and two wild type-sparing inhibitors (Staurosporine and Bosutinib) to the wild-type ErbB2 and its T798M mutant are investigated in detail by using rigorous computational analysis and binding affinity assay. Substitution of the polar threonine with a bulky methionine at residue 798 can impair and improve the direct binding affinity of wild type-selective and wild type-sparing inhibitors, respectively. Hindrance effect is responsible for the affinity decrease of wild type-selective inhibitors, while additional nonbonded interactions contribute to the affinity increase of wild type-sparing inhibitors, thus conferring selectivity to the inhibitors for mutant over wild type. The binding affinity of Staurosporine and Bosutinib to ErbB2 kinase domain is improved by 11.9-fold and 2.1-fold upon T798M mutation, respectively. Structural analysis reveals that a nonbonded network of S–π contact interactions (for Staurosporine) or an S-involving halogen bond (for Bosutinib) forms with the sulfide group of mutant Met798 residue.     

Sequence, structure, and active site analyses of p38 MAP kinase: exploiting DFG-out conformation as a strategy to design new type II leads

A new knowledge, structure, and sequence based strategy involving the effective exploitation of the DFG-out conformation is delineated. A comprehensive analysis of the structure, sequence, cocrystals, and active sites of p38 MAP kinase crystal structures present in Protein Data Bank (PDB) and the FDA approved MAP kinase drugs has been done, and the information is used for the design of type II leads. The 98 crystal structures, 138 cocrystals, and 31 FDA drugs comprise of 7 different sequences of 2 organisms viz., Homo sapiens and Mus musculus differing in sequence length, constituting both homo- and heterochains. Multiple sequence alignment with ClustalW showed >95% sequence similarity with highly conserved domains and a high propensity for mutations in the activation loop. The bound ligands were extracted, and their interactions with DFG in and out conformations were studied. These cocrystals and FDA drugs were fragmented on the basis of their binding interactions and their affinity to ATP and allosteric sites. The fragment library thus generated contains 106 fragments with overlapping drug fragments. A blue print constituting three main parts viz., head (ATP region), linker (DFG region), and tail (allosteric region) has thus been formulated and used to design 64 type II p38 MAP kinase inhibitors. The above strategy has been employed to design potent type II p38 MAP kinase inhibitors, which are shown to be very promising.     

Protein Alpha Shape Similarity Analysis (PASSA): a new method for mapping protein binding sites. Application in the design of a selective inhibitor of tyrosine kinase 2

We have developed a method that we have called Protein Alpha Shape Similarity Analysis (PASSA), that identifies interaction sites that can be utilised to achieve selectivity towards a protein. We have shown that this method is able to identify residues of tyrosine kinases that interact with known selective inhibitors using the following test cases: Abelson (Abl) kinase in complex with STI-571 and Janus kinase 2 (Jak2) in complex with AG-490. The 3D structures of the tyrosine kinase domains of Tyrosine kinase 2 (Tyk2) and Jak2 have been predicted by homology modelling. Computational docking of AG-490 and a set of tyrphostins known not to inhibit Jak2 indicated that our homology models are able to separate inhibitors from non-inhibitors. PASSA has also been used to identify unique properties of Tyk2. According to our results, interactions with hydrogen acceptors and donors on the following residues can be utilised to achieve selectivity towards Tyk2: Y955, E1053, D1062 and S1063. These residues are placed close to non-conserved hydrophobic pockets. The PASSA results, together with results from Multiple Copy Simultaneous Search (MCSS) were used to suggest functional groups of a selective Tyk2 inhibitor.     

Structural basis for selective inhibition of Src family kinases by PP1.

BACKGROUND: Small-molecule inhibitors that can target individual kinases are powerful tools for use in signal transduction research. It is difficult to find such compounds because of the enormous number of protein kinases and the highly conserved nature of their catalytic domains. Recently, a novel, potent, Src family selective tyrosine kinase inhibitor was reported (PP1). Here, we study the structural basis for this inhibitor's specificity for Src family kinases. RESULTS: A single residue corresponding to Ile338 (v-Src numbering; Thr338 in c-Src) in Src family tyrosine kinases largely controls PP1's ability to inhibit protein kinases. Mutation of Ile338 to a larger residue such as methionine or phenylalanine in v-Src makes this inhibitor less potent. Conversely, mutation of Ile338 to alanine or glycine increases PP1's potency. PP1 can inhibit Ser/Thr kinases if the residue corresponding to Ile338 in v-Src is mutated to glycine. We have accurately predicted several non-Src family kinases that are moderately (IC(50) approximately 1 microM) inhibited by PP1, including c-Abl and the MAP kinase p38. CONCLUSIONS: Our mutagenesis studies of the ATP-binding site in both tyrosine kinases and Ser/Thr kinases explain why PP1 is a specific inhibitor of Src family tyrosine kinases. Determination of the structural basis of inhibitor specificity will aid in the design of more potent and more selective protein kinase inhibitors. The ability to desensitize a particular kinase to PP1 inhibition of residue 338 or conversely to sensitize a kinase to PP1 inhibition by mutation should provide a useful basis for chemical genetic studies of kinase signal transduction.     

Design and Synthesis of a Series of Novel Macrocycle Janus Kinase 2 Inhibitors

Macrocycle has attracted the attention of many researchers in the field of medicinal chemistry due to its unique advantages and good prospects, but the difficulties in drug design and synthesis of macrocycle limit its applications. In this study, a series of macrocyclic derivatives designed from anaplastic lymphoma kinase (ALK) inhibitor lorlatinib were synthesized as Janus kinase 2 (JAK2) selective inhibitors. Among them, 17f had the best inhibitory activity (IC₅₀ = 0.177 μmol·L^–1) and selectivity for JAK2 over JAK1 and JAK3, which indicated that design of the macrocyclic derivatives might be a feasible strategy for the discovery of novel selective JAK2 inhibitors.     

Enhanced selectivity for inhibition of analog-sensitive protein kinases through scaffold optimization

The ability to inhibit any protein kinase of interest with a small molecule is enabled by a combination of genetics and chemistry. Genetics is used to modify the active site of a single kinase to render it distinct from all naturally occurring kinases. Next, organic synthesis is used to develop a small molecule, which does not bind to wild-type kinases but is a potent inhibitor of the engineered kinase. This approach, termed chemical genetics, has been used to generate highly potent mutant kinase-specific inhibitors based on a pyrazolopyrimidine scaffold. Here, we asked if the selectivity of the resulting pyrazolopyrimidines could be improved, as they inhibit several wild-type kinases with low-micromolar IC₅₀ values. Our approach to improve the selectivity of allele-specific inhibitors was to explore a second kinase inhibitor scaffold. A series of 6,9-disubstituted purines was designed, synthesized, and evaluated for inhibitory activity against several kinases in vitro and in vivo. Several purines proved to be potent inhibitors against the analog-sensitive kinases and exhibited greater selectivity than the existing pyrazolopyrimidines.     

Anti-cancer activity of novel dibenzo[b,f]azepine tethered isoxazoline derivatives

Background  

Nef is an HIV-1 accessory protein essential for viral replication and AIDS progression. Nef interacts with a multitude of host cell signaling partners, including members of the Src kinase family. Nef preferentially activates Hck, a Src-family kinase (SFK) strongly expressed in macrophages and other HIV target cells, by binding to its regulatory SH3 domain. Recently, we identified a series of kinase inhibitors that preferentially inhibit Hck in the presence of Nef. These compounds also block Nef-dependent HIV replication, validating the Nef-SFK signaling pathway as an antiretroviral drug target. Our findings also suggested that by binding to the Hck SH3 domain, Nef indirectly affects the conformation of the kinase active site to favor inhibitor association.     

The selectivity and promiscuity of brain-neuroregenerative inhibitors between ROCK1 and ROCK2 isoforms: An integration of SB-QSSR modelling,QM/MM analysis and in vitro kinase assay

The Rho-associated kinases (ROCKs) have long been recognized as an attractive therapeutic target for various neurological diseases; selective inhibition of ROCK1 and ROCK2 isoforms would result in distinct biological effects on neurogenesis, neuroplasticity and neuroregeneration after brain surgery and traumatic brain injury. However, the discovery and design of isoform-selective inhibitors remain a great challenge due to the high conservation and similarity between the kinase domains of ROCK1 and ROCK2. Here, a structure-based quantitative structure–selectivity relationship (SB-QSSR) approach was used to correlate experimentally measured selectivity with the difference in inhibitor binding to the two kinase isoforms. The resulting regression models were examined rigorously through both internal cross-validation and external blind validation; a nonlinear predictor was found to have high fitting stability and strong generalization ability, which was then employed to perform virtual screening against a structurally diverse, drug-like compound library. Consequently, five and seven hits were identified as promising candidates of 1-o-2 and 2-o-1 selective inhibitors, respectively, from which seven purchasable compounds were tested in vitro using a standard kinase assay protocol to determine their inhibitory activity against and selectivity between ROCK1 and ROCK2. The structural basis, energetic property and biological implication underlying inhibitor selectivity and promiscuity were also investigated systematically using a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme.     

Expanding the scope of native chemical ligation – templated small molecule drug synthesis via benzanilide formation

We describe a reaction system that enables the synthesis of Bcr–Abl tyrosine kinase inhibitors (TKI) via benzanilide formation in water. The reaction is based on native chemical ligation (NCL). In contrast to previous applications, we used the NCL chemistry to establish aromatic rather than aliphatic amide bonds in coupling reactions between benzoyl and o-mercaptoaniline fragments. The method was applied for the synthesis of thiolated ponatinib and GZD824 derivatives. Acid treatment provided benzothiazole structures, which opens opportunities for diversification. Thiolation affected the affinity for Abl1 kinase only moderately. Of note, a ponatinib-derived benzothiazole also showed nanomolar affinity. NCL-enabled benzanilide formation may prove useful for fragment-based drug discovery. To show that benzanilide synthesis can be put under the control of a template, we connected the benzoyl and o-mercaptoaniline fragments to DNA and peptide nucleic acid (PNA) oligomers. Complementary RNA templates enabled adjacent binding of reactive conjugates triggering a rapid benzoyl transfer from a thioester-linked DNA conjugate to an o-mercaptoaniline-DNA or -PNA conjugate. We evaluated the influence of linker length and unpaired spacer nucleotides within the RNA template on the product yield. The data suggest that nucleic acid-templated benzanilide formation could find application in the establishment of DNA-encoded combinatorial libraries (DEL).

The templated native chemical ligation between benzoyl thioesters and o-mercaptoaniline fragments proceeds in water and provides benzanilides that have nanomolar affinity for Abl1 kinase.     

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. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Discovery of EGFR selective 4,6-disubstituted pyrimidines from a combinatorial kinase-directed heterocycle library

The epidermal growth factor receptor (EGFR) tyrosine kinase was one of the first receptor tyrosine kinases to be targeted for drug development by the pharmaceutical industry due to its ubiquitous overexpression in a variety of tumors. Despite the validation of several quinazoline-based scaffolds in the clinic, there is a dearth of alternative chemical structure classes that are capable of inhibiting EGFR kinase activity selectively. Here we describe the discovery of potent and highly selective 4,6-disubstituted pyrimidine inhibitors of enzymatic and cellular EGFR activity and provide an explanation for their exceptional degree of kinase selectivity.     

Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Developing chemicals that inhibit checkpoint kinase 1 (Chk1) is a promising adjuvant therapeutic to improve the efficacy and selectivity of DNA-targeting agents. Reliable prediction of binding-free energy and binding affinity of Chk1 inhibitors can provide a guide for rational drug design. In this study, multiple docking strategies and Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) calculation were applied to predict the binding mode and free energy for a series of benzoisoquinolinones as Chk1 inhibitors. Reliable docking results were obtained using induced-fit docking and quantum mechanics/molecular mechanics (QM/MM) docking, which showed superior performance on both ligand binding pose and docking score accuracy to the rigid-receptor docking. Then, the Prime/MM-GBSA method based on the docking complex was used to predict the binding-free energy. The combined use of QM/MM docking and Prime/MM-GBSA method could give a high correlation between the predicted binding-free energy and experimentally determined pIC(50) . The molecular docking combined with Prime/MM-GBSA simulation can not only be used to rapidly and accurately predict the binding-free energy of novel Chk1 inhibitors but also provide a novel strategy for lead discovery and optimization targeting Chk1.     

Approaches to solving the rigid receptor problem by identifying a minimal set of flexible residues during ligand docking

BACKGROUND: Using fixed receptor sites derived from high-resolution crystal structures in structure-based drug design does not properly account for ligand-induced enzyme conformational change and imparts a bias into the discovery and design of novel ligands. We sought to facilitate the design of improved drug leads by defining residues most likely to change conformation, and then defining a minimal manifold of possible conformations of a target site for drug design based on a small number of identified flexible residues. RESULTS: The crystal structure of thymidylate synthase from an important pathogenic target Pneumocystis carinii (PcTS) bound to its substrate and the inhibitor, BW1843U89, is reported here and reveals a new conformation with respect to the structure of PcTS bound to substrate and the more conventional antifolate inhibitor, CB3717. We developed an algorithm for determining which residues provide 'soft spots' in the protein, regions where conformational adaptation suggests possible modifications for a drug lead that may yield higher affinity. Remodeling the active site of thymidylate synthase with new conformations for only three residues that were identified with this algorithm yields scores for ligands that are compatible with experimental kinetic data. CONCLUSIONS: Based on the examination of many protein/ligand complexes, we develop an algorithm (SOFTSPOTS) for identifying regions of a protein target that are more likely to accommodate plastically to regions of a drug molecule. Using these indicators we develop a second algorithm (PLASTIC) that provides a minimal manifold of possible conformations of a protein target for drug design, reducing the bias in structure-based drug design imparted by structures of enzymes co-crystallized with inhibitors.     

Covalent‐Allosteric Kinase Inhibitors

Targeting and stabilizing distinct kinase conformations is an instrumental strategy for dissecting conformation‐dependent signaling of protein kinases. Herein the structure‐based design, synthesis, and evaluation of pleckstrin homology (PH) domain‐dependent covalent‐allosteric inhibitors (CAIs) of the kinase Akt is reported. These inhibitors bind covalently to a distinct cysteine of the kinase and thereby stabilize the inactive kinase conformation. These modulators exhibit high potency and selectivity, and represent an innovative approach for chemical biology and medicinal chemistry research.