Homology model of the CDK1/cyclin B complex

We describe a refined homology model of a CDK1/cyclin B complex that was previously used for the structure-based optimization of the Paullone class of inhibitors. The preliminary model was formed from the homologous regions of the deposited CDK2/cyclin A crystal structure. Further refinement of the CDK1/cyclin B complex was accomplished using molecular mechanics and hydropathic analysis with a protocol of constraints and local geometry searches. For the most part, our CKD1/cyclin B homology model is very similar to the high resolution CDK2/cyclin A crystal structure regarding secondary and tertiary features. However, minor discrepancies between the two kinase structures suggest the possibility that ligand design may be specifically tuned for either CDK1 or CDK2. Our examination of the CDK1/cyclin B model includes a comparison with the CDK2/cyclin A crystal structure in the PSTAIRE interface region, connecting portions to the ATP binding domain, as well as the ATP binding site itself.     

Homology Model of the CDK1/cyclin B Complex

Abstract

We describe a refined homology model of a CDK1/cyclin B complex that was previously used for the structure-based optimization of the Paullone class of inhibitors. The preliminary model was formed from the homologous regions of the deposited CDK2/cyclin A crystal structure. Further refinement of the CDK1/cyclin B complex was accomplished using molecular mechanics and hydropathic analysis with a protocol of constraints and local geometry searches. For the most part, our CKD1/cyclin B homology model is very similar to the high resolution CDK2/cyclin A crystal structure regarding secondary and tertiary features. However, minor discrepancies between the two kinase structures suggest the possibility that ligand design may be specifically tuned for either CDK1 or CDK2. Our examination of the CDK1/cyclin B model includes a comparison with the CDK2/cyclin A crystal structure in the PSTAIRE interface region, connecting portions to the ATP binding domain, as well as the ATP binding site itself.     

Structure-based design of 2-arylamino-4-cyclohexylmethyl-5-nitroso-6-aminopyrimidine inhibitors of cyclin-dependent kinases 1 and 2

A series of O(4)-cyclohexylmethyl-5-nitroso-6-aminopyrimidines bearing 2-arylamino substituents was synthesised and evaluated for CDK1 and CDK2 inhibitory activity. Consistent with analogous studies with O(6)-cyclohexylmethylpurines, 2-arylaminopyrimidines with a sulfonamide or carboxamide group at the 4'-position were potent inhibitors, with IC(50) values against CDK2 of 1.1+/-0.3 and 34+/-8 nM, respectively. The crystal structure of the 4'-carboxamide derivative, in complex with phospho-Thr160 CDK2/cyclin A, confirmed the expected binding mode of the inhibitor, and revealed an additional interaction between the carboxamide function and an aspartate residue.     

Cyclin A2 mutagenesis analysis: a new insight into CDK activation and cellular localization requirements

Cyclin A2 is essential at two critical points in the somatic cell cycle: during S phase, when it activates CDK2, and during the G2 to M transition when it activates CDK1. Based on the crystal structure of Cyclin A2 in association with CDKs, we generated a panel of mutants to characterize the specific amino acids required for partner binding, CDK activation and subcellular localization. We find that CDK1, CDK2, p21, p27 and p107 have overlapping but distinct requirements for association with this protein. Our data highlight the crucial importance of the N-terminal α helix, in conjunction with the α3 helix within the cyclin box, in activating CDK. Several Cyclin A2 mutants selectively bind to either CDK1 or CDK2. We demonstrate that association of Cyclin A2 to proteins such as CDK2 that was previously suggested as crucial is not a prerequisite for its nuclear localization, and we propose that the whole protein structure is involved.     

Structure-based design and protein X-ray analysis of a protein kinase inhibitor.

A 5-aryl-1H-pyrazole molecular scaffold was designed to ligate the ATP binding site of cyclin dependent kinase 2 (CDK2) on the basis of crystallographic information. A search of the compound collection of Novartis using this scaffold as substructure query led to the identification of PKF049-365 as a representative of a new class of inhibitors of the cell cycle kinases CDK1/2. The three-dimensional structure of CDK2 in complex with PKF049-365 was subsequently determined by protein crystallography and refined to 1.53 A resolution. The X-ray analysis confirmed the binding mode expected from the design hypothesis. In addition, it revealed an alternative binding orientation involving a second tautomeric form of the inhibitor that was not envisaged during the design stage.     

Structure and regulation of the CDK5-p25(nck5a) complex

CDK5 plays an indispensable role in the central nervous system, and its deregulation is involved in neurodegeneration. We report the crystal structure of a complex between CDK5 and p25, a fragment of the p35 activator. Despite its partial structural similarity with the cyclins, p25 displays an unprecedented mechanism for the regulation of a cyclin-dependent kinase. p25 tethers the unphosphorylated T loop of CDK5 in the active conformation. Residue Ser159, equivalent to Thr160 on CDK2, contributes to the specificity of the CDK5-p35 interaction. Its substitution with threonine prevents p35 binding, while the presence of alanine affects neither binding nor kinase activity. Finally, we provide evidence that the CDK5-p25 complex employs a distinct mechanism from the phospho-CDK2-cyclin A complex to establish substrate specificity.     

A flavonoid gossypin binds to cyclin-dependent kinase 2

Flavonoids have low toxicity and mild activity. In order to find flavonoids showing cyclin-dependent kinase 2 (CDK2) binding effects, 347 flavonoid derivatives were docked into the crystal structure of the CDK2. The docking study showed that gossypin has a good conformational match with CDK2, which was confirmed by the binding affinity assay using NMR experiments.     

Crystallographic approach to identification of cyclin-dependent kinase 4 (CDK4)-specific inhibitors by using CDK4 mimic CDK2 protein

Genetic alteration of one or more components of the p16(INK4A)-CDK4,6/cyclin D-retinoblastoma pathway is found in more than half of all human cancers. Therefore, CDK4 is an attractive target for the development of a novel anticancer agent. However, it is difficult to make CDK4-specific inhibitors that do not possess activity for other kinases, especially CDK2, because the CDK family has high structural homology. The three-dimensional structure of CDK2, particularly that bound with the inhibitor, has provided useful information for the synthesis of CDK2-specific inhibitors. The same approach used to make CDK4-specific inhibitors was hindered by the failure to obtain a crystal structure of CDK4. To overcome this problem, we synthesized a CDK4 mimic CDK2 protein in which the ATP binding pocket of CDK2 was replaced with that of CDK4. This CDK4 mimic CDK2 was crystallized both in the free and inhibitor-bound form. The structural information thus obtained was found to be useful for synthesis of a CDK4-specific inhibitor that does not have substantial CDK2 activity. Namely, the data suggest that CDK4 has additional space that will accommodate a large substituent such as the CDK4 selective inhibitor. Inhibitors designed to bind into this large cavity should be selective for CDK4 without having substantial CDK2 activity. This design principle was confirmed in the x-ray crystal structure of the CDK4 mimic CDK2 with a new CDK4 selective inhibitor bound.     

Discovery of cyclin-dependent kinase inhibitor, CR229, using structurebased drug screening

To generate new scaffold candidates as highly selective and potent cyclin-dependent kinase (CDK) inhibitors, structure-based drug screening was performed utilizing 3D pharmacophore conformations of known potent inhibitors. As a result, CR229 (6-bromo-2,3,4,9-tetrahydro-carbolin-1-one) was generated as the hit-compound. A computational docking study using the X-ray crystallographic structure of CDK2 in complex with CR229 was evaluated. This predicted binding mode study of CR229 with CDK2 demonstrated that CR229 interacted effectively with the Leu83 and Glu81 residues in the ATP-binding pocket of CDK2 for the possible hydrogen bond formation. Furthermore, biochemical studies on inhibitory effects of CR229 on various kinases in the human cervical cancer HeLa cells demonstrated that CR229 was a potent inhibitor of CDK2 (IC50: 3 microM), CDK1 (IC50: 4.9 microM), and CDK4 (IC50: 3 microM), yet had much less inhibitory effect (IC50: >20 microM) on other kinases, such as casein kinase 2-1 (CK2- alpha1), protein kinase A (PKA), and protein kinase C (PKC). Accordingly, these data demonstrate that CR229 is a potent CDK inhibitor with anticancer efficacy.     

Alternative binding modes of an inhibitor to two different kinases.

Protein kinases are targets for therapeutic agents designed to intervene in signaling processes in the diseased state. Most kinase inhibitors are directed towards the conserved ATP binding site. Because the essential features of this site are conserved in all eukaryotic protein kinases, it is generally assumed that the same compound will bind in a similar manner to different protein kinases. The inhibitor 4,5,6,7-tetrabromobenzotriazole (TBB) is a selective inhibitor for the protein kinase CK2 (IC50 1.6 micro m) (Sarno et al. (2001) FEBS Letts.496, 44-48). Three other kinases [cyclin-dependent protein kinase 2 (CDK2), phosphorylase kinase and glycogen synthase kinase 3beta] exhibit approximately 10-fold weaker affinity for TBB than CK2. We report the crystal structure of TBB in complex with phospho-CDK2-cyclin A at 2.2 A resolution and compare the interactions with those observed for TBB bound to CK2. TBB binds at the ATP binding site of both kinases. In CDK2, each of the four bromine atoms makes polar contacts either to main chain oxygens in the hinge region of the kinase or to water molecules, in addition to several van der Waals contacts. The mode of binding of TBB to CDK2 is different from that to CK2. TBB in CDK2 is displaced more towards the hinge region between the N- and C-terminal lobes and rotated relative to TBB in CK2. The ATP binding pocket is wider in CDK2 than in CK2 resulting in fewer van der Waals contacts but TBB in CK2 does not contact the hinge. The structures show that, despite the conservation of the ATP binding pocket, the inhibitor is able to exploit different recognition features so that the same compound can bind in different ways to the two different kinases.     

Homology modeling,molecular dynamic simulation and docking studies of cyclin dependent kinase 1

In order to develop promising cyclin dependent kinase 1 inhibitors, homology modeling, docking and molecular dynamic simulation techniques were applied to get insight into the functional and structural properties of cyclin dependent kinase 1 (CDK1). Since there is no reported CDK1 crystal structural data, the three dimensional structure of CDK1 was constructed based on homology modeling. An extensive dynamic simulation was also performed on a Flavopiridol-CDK1 complex for probing the binding pattern of Flavopiridol in the active site of CDK1. The binding modes of other inhibitors to CDK1 were also proposed by molecular docking. The structural requirement for developing more potent CDK1 inhibitors was obtained by the above-mentioned molecular simulations and pharmacophore modeling.     

Discovery of a novel family of CDK inhibitors with the program LIDAEUS: structural basis for ligand-induced disordering of the activation loop

A family of 4-heteroaryl-2-amino-pyrimidine CDK2 inhibitor lead compounds was discovered with the new database-mining program LIDAEUS through in silico screening. Four compounds with IC(50) values ranging from 17 to 0.9 microM were selected for X-ray crystal analysis. Two distinct binding modes are observed, one of which resembles the hydrogen bonding pattern of bound ATP. In the second binding mode, the ligands trigger a conformational change in the activation T loop by inducing movement of Lys(33) and Asp(145) side chains. The family of molecules discovered provides an excellent starting point for the design and synthesis of tight binding inhibitors, which may lead to a new class of antiproliferative drugs.     

Synthesis, biological evaluation, and molecular docking studies of N-((1,3-diphenyl-1H-pyrazol-4-yl)methyl)aniline derivatives as novel anticancer agents

A series of N-((1,3-diphenyl-1H-pyrazol-4-yl)methyl)aniline derivatives (5a-8d) have been designed and synthesized, and their biological activities were also evaluated as potential antitumor and cyclin dependent kinase 2 (CDK2) inhibitors. Among all the compounds, compound 5a displayed the most potent CDK2/cyclin E inhibitory activity in vitro, with an IC(50) of 0.98±0.06μM. Antitumor assays indicated that compound 5a owned high antiproliferative activity against MCF-7 and B16-F10 cancer cell lines with IC(50) values of 1.88±0.11 and 2.12±0.15μM, respectively. Docking simulation was performed to insert compound 5a into the crystal structure of CDK2 at active site to determine the probable binding model. Based on the preliminary results, compound 5a with potent inhibitory activity in tumor growth may be a potential anticancer agent.     

Quinazolines as cyclin dependent kinase inhibitors

Quinazolines have been identified as inhibitors of CDK4/D1 and CDK2/E. Aspects of the SAR were investigated using solution-phase, parallel synthesis. An X-ray crystal structure was obtained of quinazoline 51 bound in CDK2 and key interactions within the ATP binding pocket are defined.     

Imidazo[1,2-b]pyridazines: a potent and selective class of cyclin-dependent kinase inhibitors

Modification of imidazo[1,2-a]pyridine CDK inhibitors lead to identification of less lipophilic imidazo[1,2-b]pyridazine series of CDK inhibitors. Although several equivalent compounds from these two series have similar structure and show similar CDK activity, the SAR of the two series differs significantly. Protein inhibitor structure determination has confirmed differences in binding mode and given some understanding of these differences in SAR. Potent and selective imidazo[1,2-b]pyridazine inhibitors of CDK2 have been identified, which show >1 microM plasma levels following a 2mg/kg oral dose to mice.     

Binding of the potential antitumour agent indirubin-5-sulphonate at the inhibitor site of rabbit muscle glycogen phosphorylase b. Comparison with ligand binding to pCDK2-cyclin A complex.

The binding of indirubin-5-sulphonate (E226), a potential anti-tumour agent and a potent inhibitor (IC(50) = 35 nm) of cyclin-dependent kinase 2 (CDK2) and glycogen phosphorylase (GP) has been studied by kinetic and crystallographic methods. Kinetic analysis revealed that E226 is a moderate inhibitor of GPb (K(i) = 13.8 +/- 0.2 micro m) and GPa (K(i) = 57.8 +/- 7.1 micro m) and acts synergistically with glucose. To explore the molecular basis of E226 binding we have determined the crystal structure of the GPb/E226 complex at 2.3 A resolution. Structure analysis shows clearly that E226 binds at the purine inhibitor site, where caffeine and flavopiridol also bind [Oikonomakos, N.G., Schnier, J.B., Zographos, S.E., Skamnaki, V.T., Tsitsanou, K.E. & Johnson, L.N. (2000) J. Biol. Chem.275, 34566-34573], by intercalating between the two aromatic rings of Phe285 and Tyr613. The mode of binding of E226 to GPb is similar, but not identical, to that of caffeine and flavopiridol. Comparative structural analyses of the GPb-E226, GPb-caffeine and GPb-flavopiridol complex structures reveal the structural basis of the differences in the potencies of the three inhibitors and indicate binding residues in the inhibitor site that can be exploited to obtain more potent inhibitors. Structural comparison of the GPb-E226 complex structure with the active pCDK2-cyclin A-E226 complex structure clearly shows the different binding modes of the ligand to GPb and CDK2; the more extensive interactions of E226 with the active site of CDK2 may explain its higher affinity towards the latter enzyme.     

EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.     

Synthesis and biological activity of N-aryl-2-aminothiazoles: potent pan inhibitors of cyclin-dependent kinases

N-Aryl aminothiazoles 6-9 were prepared from 2-bromothiazole 5 and found to be CDK inhibitors. In cells they act as potent cytotoxic agents. Selectivity for CDK1, CDK2, and CDK4 was dependent of the nature of the N-aryl group and distinct from the CDK2 selective N-acyl analogues. The N-2-pyridyl analogues 7 and 19 showed pan CDK inhibitory activity. Elaborated analogues 19 and 23 exhibited anticancer activity in mice against P388 murine leukemia. The solid-state structure of 7 bound to CDK2 shows a similar binding mode to the N-acyl analogues.