Molecular Dynamics Simulations and Classical Multidimensional Scaling Unveil New Metastable States in the Conformational Landscape of CDK2

Protein kinases are key regulatory nodes in cellular networks and their function has been shown to be intimately coupled with their structural flexibility. However, understanding the key structural mechanisms of large conformational transitions remains a difficult task. CDK2 is a crucial regulator of cell cycle. Its activity is finely tuned by Cyclin E/A and the catalytic segment phosphorylation, whereas its deregulation occurs in many types of cancer. ATP competitive inhibitors have failed to be approved for clinical use due to toxicity issues raised by a lack of selectivity. However, in the last few years type III allosteric inhibitors have emerged as an alternative strategy to selectively modulate CDK2 activity. In this study we have investigated the conformational variability of CDK2. A low dimensional conformational landscape of CDK2 was modeled using classical multidimensional scaling on a set of 255 crystal structures. Microsecond-scale plain and accelerated MD simulations were used to populate this landscape by using an out-of-sample extension of multidimensional scaling. CDK2 was simulated in the apo-form and in complex with the allosteric inhibitor 8-anilino-1-napthalenesulfonic acid (ANS). The apo-CDK2 landscape analysis showed a conformational equilibrium between an Src-like inactive conformation and an active-like form. These two states are separated by different metastable states that share hybrid structural features with both forms of the kinase. In contrast, the CDK2/ANS complex landscape is compatible with a conformational selection picture where the binding of ANS in proximity of the αC helix causes a population shift toward the inactive conformation. Interestingly, the new metastable states could enlarge the pool of candidate structures for the development of selective allosteric CDK2 inhibitors. The method here presented should not be limited to the CDK2 case but could be used to systematically unmask similar mechanisms throughout the human kinome.     

Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop

Nanoseconds long molecular dynamics (MD) trajectories of differently active complexes of human cyclin-dependent kinase 2 (inactive CDK2/ATP, semiactive CDK2/Cyclin A/ATP, fully active pT160-CDK2/Cyclin A/ATP, inhibited pT14-; pY15-; and pT14,pY15,pT160-CDK2/Cyclin A/ATP) were compared. The MD simulations results of CDK2 inhibition by phosphorylation at T14 and/or Y15 sites provide insight into the structural aspects of CDK2 deactivation. The inhibitory sites are localized in the glycine-rich loop (G-loop) positioned opposite the activation T-loop. Phosphorylation of T14 and both inhibitory sites T14 and Y15 together causes ATP misalignment for phosphorylation and G-loop conformational change. This conformational change leads to the opening of the CDK2 substrate binding box. The phosphorylated Y15 residue negatively affects substrate binding or its correct alignment for ATP terminal phospho-group transfer to the CDK2 substrate. The MD simulations of the CDK2 activation process provide results in agreement with previous X-ray data.     

Aminoimidazo[1,2-a]pyridines as a new structural class of cyclin-dependent kinase inhibitors. Part 1: Design, synthesis, and biological evaluation

We have identified a novel structural class of protein serine/threonine kinase inhibitors comprised of an aminoimidazo[1,2-a]pyridine nucleus. Compounds from this family are shown to potently inhibit cyclin-dependent kinases by competing with ATP for binding to a catalytic subunit of the protein. Structure-based design approach was used to direct this chemical scaffold toward generating potent and selective CDK2 inhibitors. The discovery of this new class of ATP-site directed protein kinase inhibitors, aminoimidazo[1,2-a]pyridines, provides the basis of new medicinal chemistry tool in search for an effective treatment of cancer and other diseases that involve protein kinase signaling pathways.     

Fluorescent peptide biosensor for probing the relative abundance of cyclin-dependent kinases in living cells

Cyclin-dependant kinases play a central role in coordinating cell growth and division, and in sustaining proliferation of cancer cells, thereby constituting attractive pharmacological targets. However, there are no direct means of assessing their relative abundance in living cells, current approaches being limited to antigenic and proteomic analysis of fixed cells. In order to probe the relative abundance of these kinases directly in living cells, we have developed a fluorescent peptide biosensor with biligand affinity for CDKs and cyclins in vitro, that retains endogenous CDK/cyclin complexes from cell extracts, and that bears an environmentally-sensitive probe, whose fluorescence increases in a sensitive fashion upon recognition of its targets. CDKSENS was introduced into living cells, through complexation with the cell-penetrating carrier CADY2 and applied to assess the relative abundance of CDK/Cyclins through fluorescence imaging and ratiometric quantification. This peptide biosensor technology affords direct and sensitive readout of CDK/cyclin complex levels, and reports on differences in complex formation when tampering with a single CDK or cyclin. CDKSENS further allows for detection of differences between different healthy and cancer cell lines, thereby enabling to distinguish cells that express high levels of these heterodimeric kinases, from cells that present decreased or defective assemblies. This fluorescent biosensor technology provides information on the overall status of CDK/Cyclin complexes which cannot be obtained through antigenic detection of individual subunits, in a non-invasive fashion which does not require cell fixation or extraction procedures. As such it provides promising perspectives for monitoring the response to therapeutics that affect CDK/Cyclin abundance, for cell-based drug discovery strategies and fluorescence-based cancer diagnostics.     

Structure-based discovery of the first allosteric inhibitors of cyclin-dependent kinase 2

Allosteric targeting of protein kinases via displacement of the structural αC helix with type III allosteric inhibitors is currently gaining a foothold in drug discovery. Recently, the first crystal structure of CDK2 with an open allosteric pocket adjacent to the αC helix has been described, prospecting new opportunities to design more selective inhibitors, but the structure has not yet been exploited for the structure-based design of type III allosteric inhibitors. In this work we report the results of a virtual screening campaign that resulted in the discovery of the first-in-class type III allosteric ligands of CDK2. Using a combination of docking and post-docking analyses made with our tool BEAR, 7 allosteric ligands (hit rate of 20%) with micromolar affinity for CDK2 were identified, some of them inhibiting the growth of breast cancer cell lines in the micromolar range. Competition experiments performed in the presence of the ATP-competitive inhibitor staurosporine confirmed that the 7 ligands are truly allosteric, in agreement with their design. Of these, compound 2 bound CDK2 with an EC₅₀ value of 3 μM and inhibited the proliferation of MDA-MB231 and ZR-75–1 breast cancer cells with IC₅₀ values of approximately 20 μM, while compound 4 had an EC₅₀ value of 71 μM and IC₅₀ values around 4 μM. Remarkably, the most potent compound 4 was able to selectively inhibit CDK2-mediated Retinoblastoma phosphorylation, confirming that its mechanism of action is fully compatible with a selective inhibition of CDK2 phosphorylation in cells. Finally, hit expansion through analog search of the most potent inhibitor 4 revealed an additional ligand 4g with similar in vitro potency on breast cancer cells.     

Design of a novel class of peptide inhibitors of cyclin-dependent kinase/cyclin activation

Cyclin dependent kinases (CDKs) are key regulators of the cell cycle progression and therefore constitute excellent targets for the design of anticancer agents. Most of the inhibitors identified to date inhibit kinase activity by interfering with the ATP-binding site of CDKs. We recently proposed that the protein/protein interface and conformational changes required in the molecular mechanism of CDK2-cyclin A activation were potential targets for the design of specific inhibitors of cell cycle progression. To this aim, we have designed and characterized a small peptide, termed C4, derived from amino acids 285-306 in the alpha5 helix of cyclin A. We demonstrate that this peptide does not interfere with complex formation but forms stable complexes with CDK2-cyclin A. The C4 peptide significantly inhibits kinase activity of complexes harboring CDK2 in a competitive fashion with respect to substrates but does not behave as an ATP antagonist. Moreover, when coupled with the protein transduction domain of Tat, the C4 peptide blocks the proliferation of tumor cell lines, thereby constituting a potent lead for the development of specific CDK-cyclin inhibitors.     

Discovery of novel furanylbenzamide inhibitors that target oncogenic tyrosine phosphatase SHP2 in leukemia cells

Disturbance of the dynamic balance between tyrosine phosphorylation and dephosphorylation of signaling molecules, controlled by protein tyrosine kinases and protein tyrosine phosphatases (PTPs), is known to lead to the development of cancer. While most approved targeted cancer therapies are tyrosine kinase inhibitors, PTPs have long been stigmatized as undruggable and have only recently gained renewed attention in drug discovery. One PTP target is the Src-homology 2 domain–containing phosphatase 2 (SHP2). SHP2 is implicated in tumor initiation, progression, metastasis, and treatment resistance, primarily because of its role as a signaling nexus of the extracellular signal–regulated kinase pathway, acting upstream of the small GTPase Ras. Efforts to develop small molecules that target SHP2 are ongoing, and several SHP2 allosteric inhibitors are currently in clinical trials for the treatment of solid tumors. However, while the reported allosteric inhibitors are highly effective against cells expressing WT SHP2, none have significant activity against the most frequent oncogenic SHP2 variants that drive leukemogenesis in several juvenile and acute leukemias. Here, we report the discovery of novel furanylbenzamide molecules as inhibitors of both WT and oncogenic SHP2. Importantly, these inhibitors readily cross cell membranes, bind and inhibit SHP2 under physiological conditions, and effectively decrease the growth of cancer cells, including triple-negative breast cancer cells, acute myeloid leukemia cells expressing either WT or oncogenic SHP2, and patient-derived acute myeloid leukemia cells. These novel compounds are effective chemical probes of active SHP2 and may serve as starting points for therapeutics targeting WT or mutant SHP2 in cancer.     

Optimization of norbornyl‐based carbocyclic nucleoside analogs as cyclin‐dependent kinase 2 inhibitors

We report on the discovery of norbornyl moiety as a novel structural motif for cyclin‐dependent kinase 2 (CDK2) inhibitors which was identified by screening a carbocyclic nucleoside analogue library. Three micromolar hits were expanded by the use of medicinal chemistry methods into a series of 16 novel compounds. They had prevailingly micromolar activities against CDK2 and the best compound of the series attained IC₅₀ of 190 nM. The binding modes were explored in molecular details by modeling and docking. Quantum mechanics‐based scoring was used to rationalize the affinities. In conclusion, the discovered 9‐hydroxymethylnorbornyl moiety was shown by joint experimental‐theoretical efforts to be able to serve as a novel substituent for CDK2 inhibitors. This finding opens door to the exploration of chemical space towards more effective derivatives targeting this important class of protein kinases.     

Evolutionary variation and adaptation in a conserved protein kinase allosteric network: Implications for inhibitor design

The activation of protein kinases involves conformational changes in key functional regions of the kinase domain, a detailed understanding of which is essential for the design of selective protein kinase inhibitors. Through statistical analysis of protein kinase sequences and crystal structures from diverse organisms, we recently proposed that the activation of protein kinases involves a hidden strain switch in the catalytic loop. Specifically, we demonstrated that the backbone torsion-angles of residues in the catalytic loop switch from a “relaxed” to “strained” conformation upon kinase activation and the strained geometry results in a network of hydrogen bonds involving conserved non-catalytic residues in the ATP and substrate binding lobes. Here, we further explore this activation mechanism by analyzing families that lack the canonical hydrogen bonding interactions with the strained backbone. We find that alternative mechanisms have evolved to maintain catalytic loop strain. In PIM kinase, for example, two water molecules account for the lack of a conserved aspartate in the substrate binding by hydrogen bonds to the strained backbone. We discuss the relevance of these findings in the design of family-specific allosteric inhibitors, and in predicting the structural and functional impact of cancer mutations that alter the strain associated hydrogen bonding network. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

A novel mode of protein kinase inhibition exploiting hydrophobic motifs of autoinhibited kinases: discovery of ATP-independent inhibitors of fibroblast growth factor receptor

Protein kinase inhibitors with enhanced selectivity can be designed by optimizing binding interactions with less conserved inactive conformations because such inhibitors will be less likely to compete with ATP for binding and therefore may be less impacted by high intracellular concentrations of ATP. Analysis of the ATP-binding cleft in a number of inactive protein kinases, particularly in the autoinhibited conformation, led to the identification of a previously undisclosed non-polar region in this cleft. This ATP-incompatible hydrophobic region is distinct from the previously characterized hydrophobic allosteric back pocket, as well as the main pocket. Generalized hypothetical models of inactive kinases were constructed and, for the work described here, we selected the fibroblast growth factor receptor (FGFR) tyrosine kinase family as a case study. Initial optimization of a FGFR2 inhibitor identified from a library of commercial compounds was guided using structural information from the model. We describe the inhibitory characteristics of this compound in biophysical, biochemical, and cell-based assays, and have characterized the binding mode using x-ray crystallographic studies. The results demonstrate, as expected, that these inhibitors prevent activation of the autoinhibited conformation, retain full inhibitory potency in the presence of physiological concentrations of ATP, and have favorable inhibitory activity in cancer cells. Given the widespread regulation of kinases by autoinhibitory mechanisms, the approach described herein provides a new paradigm for the discovery of inhibitors by targeting inactive conformations of protein kinases.     

Identification of novel CDK2 inhibitors by a multistage virtual screening method based on SVM,pharmacophore and docking model

Abstract

Cyclin-dependent kinase 2 (CDK2) is the family of Ser/Thr protein kinases that has emerged as a highly selective with low toxic cancer therapy target. A multistage virtual screening method combined by SVM, protein-ligand interaction fingerprints (PLIF) pharmacophore and docking was utilised for screening the CDK2 inhibitors. The evaluation of the validation set indicated that this method can be used to screen large chemical databases because it has a high hit-rate and enrichment factor (80.1% and 332.83 respectively). Six compounds were screened out from NCI, Enamine and Pubchem database. After molecular dynamics and binding free energy calculation, two compounds had great potential as novel CDK2 inhibitors and they also showed selective inhibition against CDK2 in the kinase activity assay.     

Synthesis of 6-(het) ary Xylocydine analogues and evaluating their inhibitory activities of CDK1 and CDK2 in vitro

A series of purine nucleoside analogues bearing an aryl and hetaryl group in position 6 were prepared and their biological activities were assessed by in vitro CDK1/Cyclin B1 and CDK2/Cyclin A2 kinase assay. From the synthesized chemicals, three Xylocydine derivatives 3h, 3i, and 3j exhibited specific inhibitory activities on CDK2/Cyclin A2 with IC(50) values of 4.6, 4.8, and 55 μM, respectively. Those three compounds all induced G1/S phase arrest in Human epithelial carcinoma cell line (HeLa), and the results suggested they may inhibit CDK2 activity in vitro. Furthermore, molecular modeling study, their docking into Cyclin Dependant Kinase 2 (CDK2) active site showed high docking scores. Taken together, these data suggest that, those three compounds are good inhibitors of CDK2 for studying this kinase signal transduction pathway in cell system.     

An in silico approach for the discovery of CDK5/p25 interaction inhibitors

The lack of selectivity of all existing ATP competitive inhibitors for a single cyclin-dependent kinase (CDK) has led us to redirect the structure-based molecule design from targeting the classic ATP-binding pocket in CDK5 toward the CDK5/p25 interface. The aim was to seek novel inhibition mechanisms to interrupt protein-protein interactions. A combined strategy of alanine-scanning calculations for locating binding sites, virtual screening for small molecules, molecular dynamics simulations for examining the binding stability of virtual screening hits and bio-assays for testing the level of inhibition was set up and used to explore novel inhibitors capable of interrupting the interactions between the proteins, and consequently of inhibiting the kinase activity. Two compounds were shown to inhibit the complex formation between CDK5 and p25 through p25 binding. They could open avenues for the discovery of new types of structures that prevent interactions between CDK5 and p25 or other CDK and activator proteins, and, more importantly, provide leads in the development of selective inhibitors among CDKs.     

Structural features of the protein kinase domain and targeted binding by small-molecule inhibitors

Protein kinases are key components in cellular signaling pathways as they carry out the phosphorylation of proteins, primarily on Ser, Thr, and Tyr residues. The catalytic activity of protein kinases is regulated, and they can be thought of as molecular switches that are controlled through protein–protein interactions and post-translational modifications. Protein kinases exhibit diverse structural mechanisms of regulation and have been fascinating subjects for structural biologists from the first crystal structure of a protein kinase over 30 years ago, to recent insights into kinase assemblies enabled by the breakthroughs in cryo-EM. Protein kinases are high-priority targets for drug discovery in oncology and other disease settings, and kinase inhibitors have transformed the outcomes of specific groups of patients. Most kinase inhibitors are ATP competitive, deriving potency by occupying the deep hydrophobic pocket at the heart of the kinase domain. Selectivity of inhibitors depends on exploiting differences between the amino acids that line the ATP site and exploring the surrounding pockets that are present in inactive states of the kinase. More recently, allosteric pockets outside the ATP site are being targeted to achieve high selectivity and to overcome resistance to current therapeutics. Here, we review the key regulatory features of the protein kinase family, describe the different types of kinase inhibitors, and highlight examples where the understanding of kinase regulatory mechanisms has gone hand in hand with the development of inhibitors.     

Detection of allosteric kinase inhibitors by displacement of active site probes

Non-adenosine triphosphate (ATP) competitive, allosteric inhibitors provide a promising avenue to develop highly selective small-molecule kinase inhibitors. Although this class of compounds is growing, detection of such inhibitors can be challenging as standard kinase activity assays preferentially detect compounds that bind to active kinases in an ATP competitive manner. We have previously described a time-resolved fluorescence resonance energy transfer (TR-FRET)-based kinase binding assay using the competitive displacement of ATP competitive active site fluorescent probes ("tracers"). Although this format has gained acceptance, published data with this and related formats are almost entirely without examples of non-ATP competitive compounds. Thus, this study addresses whether this format is useful for non-ATP competitive inhibitors. To this end, 15 commercially available non-ATP competitive inhibitors were tested for their ability to displace ATP competitive probes. Despite the diversity of both compound structures and their respective targets, 14 of the 15 compounds displaced the tracers with IC(50) values comparable to literature values. We conclude that such binding assays are well suited for the study of non-ATP competitive inhibitors. In addition, we demonstrate that allosteric inhibitors of BCR-Abl and MEK bind preferentially to the nonphosphorylated (i.e., inactive) form of the kinase, indicating that binding assays may be a preferred format in some cases.     

Structural basis for the ORC1‐Cyclin A association

Progression of cell cycle is regulated by sequential expression of cyclins, which associate with distinct cyclin kinases to drive the transition between different cell cycle phases. The complex of Cyclin A with cyclin‐dependent kinase 2 (CDK2) controls the DNA replication activity through phosphorylation of a set of chromatin factors, which critically influences the S phase transition. It has been shown that the direct interaction between the Cyclin A‐CDK2 complex and origin recognition complex subunit 1 (ORC1) mediates the localization of ORC1 to centrosomes, where ORC1 inhibits cyclin E‐mediated centrosome reduplication. However, the molecular basis underlying the specific recognition between ORC1 and cyclins remains elusive. Here we report the crystal structure of Cyclin A‐CDK2 complex bound to a peptide derived from ORC1 at 2.54 å resolution. The structure revealed that the ORC1 peptide interacts with a hydrophobic groove, termed cyclin binding groove (CBG), of Cyclin A via a KXL motif. Distinct from other identified CBG‐binding sequences, an arginine residue flanking the KXL motif of ORC1 inserts into a neighboring acidic pocket, contributing to the strong ORC1‐Cyclin A association. Furthermore, structural and sequence analysis of cyclins reveals divergence on the ORC1‐binding sites, which may underpin their differential ORC1‐binding activities. This study provides a structural basis of the specific ORC1‐cyclins recognition, with implication in development of novel inhibitors against the cyclin/CDK complexes.     

Proteome-wide identification of cellular targets affected by bisindolylmaleimide-type protein kinase C inhibitors

Bisindolylmaleimide compounds such as GF109203X are potent inhibitors of protein kinase C (PKC) activity. Although bisindolylmaleimides are not entirely selective for PKC and are known to inhibit a few other protein kinases, these reagents have been extensively used to study the functional roles of PKC family enzymes in cellular signal transduction for more than a decade. Here, we establish a proteomics approach to gain further insights into the cellular effects of this compound class. Functional immobilization of suitable bisindolylmaleimide analogues in combination with the specific purification of cellular binding proteins by affinity chromatography led to the identification of several known and previously unknown enzyme targets. Subsequent in vitro binding and activity assays confirmed the protein kinases Ste20-related kinase and cyclin-dependent kinase 2 (CDK2) and the non-protein kinases adenosine kinase and quinone reductase type 2 as novel targets of bisindolylmaleimide inhibitors. As observed specifically for CDK2, minor chemical variation of the ligand by immobilizing the closely related bisindolylmaleimides III, VIII, and X dramatically affected target binding. These observed changes in affinity correlated with both the measured IC(50) values for in vitro CDK2 inhibition and results from molecular docking into the CDK2 crystal structure. Moreover, the conditions for affinity purification could be adapted in a way that immobilized bisindolylmaleimide III selectively interacted with either PKC alpha or ribosomal S6 protein kinase 1 only after activation of these kinases. Thus, we have established an efficient technique for the rapid identification of cellular bisindolylmaleimide targets and further demonstrate the comparative selectivity profiling of closely related kinase inhibitors within a cellular proteome.     

Allosteric modulation of kinases and GPCRs: design principles and structural diversity

Allosteric binding sites, as opposed to traditional orthosteric binding sites, offer unparalleled opportunities for drug discovery by providing high levels of selectivity, mimicking physiological conditions, affording fewer side effects because of desensitization/downregulation, and engendering ligands with chemotypes divergent from orthosteric ligands. For kinases, allosteric mechanisms described to date include alteration of protein kinase conformation blocking productive ATP binding which appear 'ATP competitive' or blocking kinase activation by conformational changes that are 'ATP non-competitive'. For GPCRs, allosteric mechanisms impart multiple modes of target modulation (positive allosteric modulation (PAM), negative allosteric modulation (NAM), neutral cooperativity, partial antagonism (PA), allosteric agonism and allosteric antagonism). Here, we review recent developments in the design principles and structural diversity of allosteric ligands for kinases and GPCRs.     

Pyrazolobenzodiazepines: Part I. Synthesis and SAR of a potent class of kinase inhibitors

A novel series of pyrazolobenzodiazepines 3 has been identified as potent inhibitors of cyclin-dependent kinase 2 (CDK2). Their synthesis and structure–activity relationships (SAR) are described. Representative compounds from this class reversibly inhibit CDK2 activity in vitro, and block cell cycle progression in human tumor cell lines. Further exploration has revealed that this class of compounds inhibits several kinases that play critical roles in cancer cell growth and division as well as tumor angiogenesis. Together, these properties suggest a compelling basis for their use as antitumor agents.