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.     

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.     

Direct binding assay for the detection of type IV allosteric inhibitors of Abl

Abelson (Abl) tyrosine kinase is an important cellular enzyme that is rendered constitutively active in the breakpoint cluster region (BCR)-Abl fusion protein, contributing to several forms of leukemia. Although inhibiting BCR-Abl activity with imatinib shows great clinical success, many patients acquire secondary mutations that result in resistance to imatinib. Second-generation inhibitors such as dasatinib and nilotinib can overcome the majority of these mutations but fail to treat patients with an especially prevalent T315I mutation at the gatekeeper position of the kinase domain. However, a combination of nilotinib with an allosteric type IV inhibitor was recently shown to overcome this clinically relevant point mutation. In this study, we present the development of a direct binding assay that enables the straightforward detection of allosteric inhibitors which bind within the myristate pocket of Abl. The assay is amenable to high-throughput screening and exclusively detects the binding of ligands to this unique allosteric site.     

Computational systematic selectivity of the Fasalog inhibitors between ROCK-I and ROCK-II kinase isoforms in Alzheimer’s disease

Human Rho-associated coiled-coil forming kinase (ROCK) is a class of essential neurokinases that consists of two structurally conserved isoforms ROCK-I and ROCK-II; they have been revealed to play distinct roles in the pathogenesis of Alzheimer’s disease (AD) and other neurological disorders. Selective targeting of the two kinase isoforms with small-molecule inhibitors is a great challenge due to the surprisingly high homology in kinase domain (92 %) and the full identity in kinase active site (100 %). Here, we describe a computational protocol to systematically profile the selectivity of Fasudil and its 25 analogs (termed as Fasalogs) between the two kinase isoforms. It is suggested that the substitution of Fasudil’s 1,4-diazepane moiety with rigid ring such as Ripasudil and Dimehtylfasudil would render the resulting inhibitors of ROCK-II over ROCK-I (II-o-I) selectivity, while the substitution with long, flexible group such as H-89 and BDBM92607 tends to have I-o-II selectivity. Structural analysis reveals that the inhibitor affinity is not only determined by the identical active site, but also contributed from the non-identical first and second shells of the site as well as other non-conserved kinase regions, which can indirectly influence the active site and inhibitor binding through allosteric effect. A further kinase assay basically confirms the computational findings, which also exhibits a good consistence with theoretical selectivity over 10 tested samples (R_p = 0.89). In particular, the Fasalog compounds Dimehtylfasudil and H-89 are identified as II-o-I and I-o-II selective inhibitors. They can be considered as promising lead molecular entities to develop new specific ROCK isoform-selective Fasalog inhibitors.     

Sphingosine kinase 1 (SK1) allosteric inhibitors that target the dimerization site

The sphingosine kinase 1 (SK1)/sphingosine-1-phosphate (S1P) signaling pathway is a crucial target for numerous human diseases from cancer to cardiovascular diseases. However, available SK1 inhibitors that target the active site suffer from poor potency, selectivity and pharmacokinetic properties. The selectivity issue of the kinases, which share a highly-conserved ATP-pocket, can be overcome by targeting the less-conserved allosteric sites. SK1 is known to function minimally as a dimer; however, the crystal structure of the SK1 dimer has not been determined. In this study, a template-based algorithm implemented in PRISM was used to predict the SK1 dimer structure and then the possible allosteric sites at the dimer interface were determined via SiteMap. These sites were used in a virtual screening campaign that includes an integrated workflow of structure-based pharmacophore modeling, virtual screening, molecular docking, re-screening of common scaffolds to propose a series of compounds with different scaffolds as potential allosteric SK1 inhibitors. Finally, the stability of the SK1-ligand complexes was analyzed by molecular dynamics simulations. As a final outcome, ligand 7 having a 4,9-dihydro-1H-purine scaffold and ligand 12 having a 2,3,4,9-tetrahydro-1H-β-carboline scaffold were found to be potential selective inhibitors for SK1.     

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.     

An isoform-selective, small-molecule inhibitor targets the autoregulatory mechanism of p21-activated kinase

Autoregulatory domains found within kinases may provide more unique targets for chemical inhibitors than the conserved ATP-binding pocket targeted by most inhibitors. The kinase Pak1 contains an autoinhibitory domain that suppresses the catalytic activity of its kinase domain. Pak1 activators relieve this autoinhibition and initiate conformational rearrangements and autophosphorylation events leading to kinase activation. We developed a screen for allosteric inhibitors targeting Pak1 activation and identified the inhibitor IPA-3. Remarkably, preactivated Pak1 is resistant to IPA-3. IPA-3 also inhibits activation of related Pak isoforms regulated by autoinhibition, but not more distantly related Paks, nor >200 other kinases tested. Pak1 inhibition by IPA-3 in live cells supports a critical role for Pak in PDGF-stimulated Erk activation. These studies illustrate an alternative strategy for kinase inhibition and introduce a highly selective, cell-permeable chemical inhibitor of Pak.     

Computational simulation of ligand docking to L-type pyruvate kinase subunit

Computational blind docking approach was used for mapping of possible binding sites in L-type pyruvate kinase subunit for peptides, RRASVA and the phosphorylated derivative RRAS(P_i)VA, which model the phosphorylatable N-terminal regulatory domain of the enzyme. In parallel, the same docking analysis was done for both substrates of this enzyme, phosphoenolpyruvate (PEP) and adenosine diphosphate (ADP), and for docking of fructose 1,6-bisphosphate (FBP), which is the allosteric activator of the enzyme. The binding properties of the entire surface of the protein were scanned and several possible binding sites were identified in domains A and C of the protein, while domain B revealed no docking sites for peptides or for substrates or the allosteric regulator. It was found that the docking sites of different ligands were partially overlapping, pointing to the possibility that some regulatory effects, observed in the case of L-type pyruvate kinase, may be caused by the competition of different ligands for the same binding sites.     

QSAR‐based targeting anaplastic lymphoma kinase (ALK) variants with noncognate inhibitors in pediatric acute lymphoblastic leukemia

Human anaplastic lymphoma kinase (ALK) is a potential target for the treatment of pediatric acute lymphoblastic leukemia. However, a number of residue mutations in ALK kinase domain have been observed to cause drug resistance in pediatric acute lymphoblastic leukemia chemotherapy. Here, a chemometrics quantitative structure‐activity relationship predictor was developed using a structure‐based panel of kinase‐inhibitor activity data. The predictor was validated rigorously through internal cross‐validation and external blind test to ensure its statistical reliability, which was then used to computationally construct a systematic activity profile of 13 noncognate kinase inhibitors against both wild‐type ALK^wt and cancer‐related variants ALK^vt. It is revealed that most noncognate inhibitors exhibit weak potency on ALK^wt, but some of them are able to selectively target ALK^vt over ALK^wt. The chemometrics findings were then evaluated by using a kinase inhibition protocol; results showed that few noncognate inhibitors are 2‐ to 5‐fold higher potent against ALK variant than wild‐type kinase.     

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.     

Markov State Models and Molecular Dynamics Simulations Provide Understanding of the Nucleotide-Dependent Dimerization-Based Activation of LRRK2 ROC Domain

Mutations in leucine-rich repeat kinase 2 (LRRK2) are recognized as the most frequent cause of Parkinson’s disease (PD). As a multidomain ROCO protein, LRRK2 is characterized by the presence of both a Ras-of-complex (ROC) GTPase domain and a kinase domain connected through the C-terminal of an ROC domain (COR). The bienzymatic ROC–COR–kinase catalytic triad indicated the potential role of GTPase domain in regulating kinase activity. However, as a functional GTPase, the detailed intrinsic regulation of the ROC activation cycle remains poorly understood. Here, combining extensive molecular dynamics simulations and Markov state models, we disclosed the dynamic structural rearrangement of ROC’s homodimer during nucleotide turnover. Our study revealed the coupling between dimerization extent and nucleotide-binding state, indicating a nucleotide-dependent dimerization-based activation scheme adopted by ROC GTPase. Furthermore, inspired by the well-known R1441C/G/H PD-relevant mutations within the ROC domain, we illuminated the potential allosteric molecular mechanism for its pathogenetic effects through enabling faster interconversion between inactive and active states, thus trapping ROC in a prolonged activated state, while the implicated allostery could provide further guidance for identification of regulatory allosteric pockets on the ROC complex. Our investigations illuminated the thermodynamics and kinetics of ROC homodimer during nucleotide-dependent activation for the first time and provided guidance for further exploiting ROC as therapeutic targets for controlling LRRK2 functionality in PD treatment.     

Identification of Flavonoids as Putative ROS-1 Kinase Inhibitors Using Pharmacophore Modeling for NSCLC Therapeutics

Non-small cell lung cancer (NSCLC) is a lethal non-immunogenic malignancy and proto-oncogene ROS-1 tyrosine kinase is one of its clinically relevant oncogenic markers. The ROS-1 inhibitor, crizotinib, demonstrated resistance due to the Gly2032Arg mutation. To curtail this resistance, researchers developed lorlatinib against the mutated kinase. In the present study, a receptor-ligand pharmacophore model exploiting the key features of lorlatinib binding with ROS-1 was exploited to identify inhibitors against the wild-type (WT) and the mutant (MT) kinase domain. The developed model was utilized to virtually screen the TimTec flavonoids database and the retrieved drug-like hits were subjected for docking with the WT and MT ROS-1 kinase. A total of 10 flavonoids displayed higher docking scores than lorlatinib. Subsequent molecular dynamics simulations of the acquired flavonoids with WT and MT ROS-1 revealed no steric clashes with the Arg2032 (MT ROS-1). The binding free energy calculations computed via molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) demonstrated one flavonoid (Hit) with better energy than lorlatinib in binding with WT and MT ROS-1. The Hit compound was observed to bind in the ROS-1 selectivity pocket comprised of residues from the β-3 sheet and DFG-motif. The identified Hit from this investigation could act as a potent WT and MT ROS-1 inhibitor.     

Structural Insight and Development of EGFR Tyrosine Kinase Inhibitors

Lung cancer has a high prevalence, with a growing number of new cases and mortality every year. Furthermore, the survival rate of patients with non-small-cell lung carcinoma (NSCLC) is still quite low in the majority of cases. Despite the use of conventional therapy such as tyrosine kinase inhibitor for Epidermal Growth Factor Receptor (EGFR), which is highly expressed in most NSCLC cases, there was still no substantial improvement in patient survival. This is due to the drug’s ineffectiveness and high rate of resistance among individuals with mutant EGFR. Therefore, the development of new inhibitors is urgently needed. Understanding the EGFR structure, including its kinase domain and other parts of the protein, and its activation mechanism can accelerate the discovery of novel compounds targeting this protein. This study described the structure of the extracellular, transmembrane, and intracellular domains of EGFR. This was carried out along with identifying the binding pose of commercially available inhibitors in the ATP-binding and allosteric sites, thereby clarifying the research gaps that can be filled. The binding mechanism of inhibitors that have been used clinically was also explained, thereby aiding the structure-based development of new drugs.     

Creating chemical diversity to target protein kinases

Protein kinases play crucial roles in regulating virtually every cellular process and are currently attracting tremendous interest as drug targets from the pharmaceutical industry. The major challenges facing the development of the potential kinase inhibitor drugs are: selectivity, physical properties (solubility, molecular weight), and pharmacological properties (bioavailability, half life, toxicity, etc.) This review focuses on how selective protein kinase inhibitors that target the ATP and allosteric binding sites are currently being identified and optimized.     

Identification of Potent,Selective Protein Kinase C Inhibitors Based on a Phorbol Skeleton

The elucidation of specific functions of protein kinase C (PKC) subtypes in physiological processes is an important challenge for the future development of new drug targets. Subtype‐selective PKC agonists and antagonists are useful biological tools for this purpose. Most of the currently used PKC modulators elicit their activities through binding to the ATP binding site of PKC, which shares many features with other kinases. PKC modulators that target the PKC regulatory domain are considered to be advantageous in terms of selectivity, because the structure of the regulatory domain is intrinsic to each PKC subtype. In this paper, we describe the identification of new potent and conventional PKC‐selective inhibitors that target the regulatory domain. The inhibitors contain a phorbol skeleton, a naturally occurring potent and selective PKC regulatory domain binder, with a perfluorinated alkyl group and a polyether hydrophilic chain on a terephthaloyl aromatic ring at the C12 position. Both of these substituents are essential for the potent inhibitory activity. Specifically, the binding affinity between PKC and the phorbol ester analogues was improved by an electron‐deficient aromatic ring at C12. This finding cannot be explained by the previously proposed binding model and suggests a new binding mode between phorbol esters and PKC.     

When Kinases Meet PROTACs

Small molecule drugs targeting kinases have revolutionized treatment options for millions of patients worldwide, especially in oncology. These targeted treatments have less side effects because they inhibit a specific dysfunctional kinase usually with relatively narrow selectivity. However, kinase inhibitors do have well‐established liabilities, most prominently the emergence of drug resistance. Moreover, the majority of kinases are multidomain and multifunctional proteins that in addition to their enzymatic activity have scaffolding and other roles, and inhibitors seldom address these alternative functions. Recently, small molecule mediated targeted protein degradation emerged as a new pharmacological strategy. The majority of small molecule degraders are bispecific molecules called proteolysis targeting chimeras (PROTACs), and their mechanism of action is based on simultaneous recruitment of the target of interest and an E3 ligase, resulting in target polyubiquitination and eventual destruction by the proteasome. Over the last couple of years, PROTAC strategy has been developed and validated for a range of targets, including kinases. Here, we introduce the targeted protein degradation strategy, provide an overview of representative kinase PROTACs, and describe design rationales, efficacy and specificity. We also discuss their potential advantages, as well as comment on some of the limitations of this emerging pharmacological modality.     

Selective Inhibitors of FKBP51 Employ Conformational Selection of Dynamic Invisible States

The recently discovered SAFit class of inhibitors against the Hsp90 co‐chaperone FKBP51 show greater than 10 000‐fold selectivity over its closely related paralogue FKBP52. However, the mechanism underlying this selectivity remained unknown. By combining NMR spectroscopy, biophysical and computational methods with mutational analysis, we show that the SAFit molecules bind to a transient pocket in FKBP51. This represents a weakly populated conformation resembling the inhibitor‐bound state of FKBP51, suggesting conformational selection rather than induced fit as the major binding mechanism. The inhibitor‐bound conformation of FKBP51 is stabilized by an allosteric network of residues located away from the inhibitor‐binding site. These residues stabilize the Phe67 side chain in a dynamic outward conformation and are distinct in FKBP52, thus rationalizing the basis for the selectivity of SAFit inhibitors. Our results represent a paradigm for the selective inhibition of transient binding pockets.     

Two-track virtual screening approach to identify both competitive and allosteric inhibitors of human small C-terminal domain phosphatase 1

Despite a wealth of persuasive evidence for the involvement of human small C-terminal domain phosphatase 1 (Scp1) in the impairment of neuronal differentiation and in Huntington’s disease, small-molecule inhibitors of Scp1 have been rarely reported so far. This study aims to the discovery of both competitive and allosteric Scp1 inhibitors through the two-track virtual screening procedure. By virtue of the improvement of the scoring function by implementing a new molecular solvation energy term and by reoptimizing the atomic charges for the active-site Mg²⁺ ion cluster, we have been able to identify three allosteric and five competitive Scp1 inhibitors with low-micromolar inhibitory activity. Consistent with the results of kinetic studies on the inhibitory mechanisms, the allosteric inhibitors appear to be accommodated in the peripheral binding pocket through the hydrophobic interactions with the nonpolar residues whereas the competitive ones bind tightly in the active site with a direct coordination to the central Mg²⁺ ion. Some structural modifications to improve the biochemical potency of the newly identified inhibitors are proposed based on the binding modes estimated with docking simulations.     

Molecular simulation studies to reveal the binding mechanisms of shikonin derivatives inhibiting VEGFR-2 kinase

Traditional vascular endothelial growth factor receptor 2 (VEGFR-2) inhibitors can manage angiogenesis; however, severe toxicity and resistance limit their long-term applications in clinical therapy. Shikonin (SHK) and its derivatives could be promising to inhibit the VEGFR-2 mediated angiogenesis, as they are reported to bind in the catalytic kinase domain with low affinity. However, the detailed molecular insights and binding dynamics of these natural inhibitors are unknown, which is crucial for potential SHK based lead design. Therefore, the present study employed molecular modeling and simulations techniques to get insight into the binding behaviors of SHK and its two derivates, β-hydroxyisovalerylshikonin (β-HIVS) and acetylshikonin (ACS). Here the intermolecular interactions between protein and ligands were studied by induced fit docking approach, which were further evaluated by treating QM/MM (quantum mechanics/molecular mechanics) and molecular dynamics (MD) simulation. The result showed that the naphthazarin ring of the SHK derivates is vital for strong binding to the catalytic domain; however, the binding stability can be modulated by the side chain modification. Because of having electrostatic potential, this ring makes essential interactions with the DFG (Asp1046 and Phe1047) motif and also allows interacting with the allosteric binding site. Taken together, the studies will advance our knowledge and scope for the development of new selective VEGFR-2 inhibitors based on SHK and its analogs.     

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.