Identification of novel purine and pyrimidine cyclin-dependent kinase inhibitors with distinct molecular interactions and tumor cell growth inhibition profiles

A computational study of the protein-ligand interactions in CDK2 inhibitors: using quantum mechanics/molecular mechanics interaction energy as a predictor of the biological activity

We report a combined quantum mechanics/molecular mechanics (QM/MM) study to determine the protein-ligand interaction energy between CDK2 (cyclin-dependent kinase 2) and five inhibitors with the N(2)-substituted 6-cyclohexyl-methoxy-purine scaffold. The computational results in this work show that the QM/MM interaction energy is strongly correlated to the biological activity and can be used as a predictor, at least within a family of substrates. A detailed analysis of the protein-ligand structures obtained from molecular dynamics simulations shows specific interactions within the active site that, in some cases, have not been reported before to our knowledge. The computed interaction energy gauges the strength of protein-ligand interactions. Finally, energy decomposition and multiple regression analyses were performed to check the contribution of the electrostatic and van der Waals energies to the total interaction energy and to show the capabilities of the computational model to identify new potent inhibitors.

Design, synthesis, and antiproliferative and CDK2-cyclin a inhibitory activity of novel flavopiridol analogues

The design and synthesis of a small library of 8-amidoflavone, 8-sulfonamidoflavone, 8-amido-7-hydroxyflavone, and heterocyclic analogues of flavopiridol is reported. The potential activity of these compounds as kinase inhibitors was evaluated by cytotoxicity studies in MCF-7 and ID-8 cancer cell lines and inhibition of CDK2-Cyclin A enzyme activity in vitro. The antiproliferative and CDK2-Cyclin A inhibitory activity of these analogues was significantly lower than the activity of flavopiridol. Molecular docking simulations were carried out and these studies suggested a different binding orientation inside the CDK2 binding pocket for these analogues compared to flavopiridol.

A novel computational analysis of ligand-induced conformational changes in the ATP binding sites of cyclin dependent kinases

Protein kinases in general are known to be very flexible macromolecules. In this article, the conformational plasticity of the ATP binding site in cyclin dependent kinases is analyzed. Movement of the two lysine residues lining the ATP binding site are shown to play a major role in the conformational variability of the site. Linear models are developed to identify and quantify ligand properties that maximally influence the lysine side chain conformations. A few simple properties of the ligands are shown to account for more than 70% of the variation in the lysine conformations. The results are validated using test data and molecular simulation studies. Illustrative applications of the results of this analysis to finding the appropriate crystal structure for molecular docking and binding mode predictions of novel ligands are provided. This work provides a new approach to quantify ligand-induced conformational changes in the active sites of flexible proteins and to find the appropriate crystal structure for docking novel ligands.

Novel structural features of CDK inhibition revealed by an ab initio computational method combined with dynamic simulations

The rational development of specific inhibitors for the approximately 500 protein kinases encoded in the human genome is impeded by a poor understanding of the structural basis for the activity and selectivity of small molecules that compete for ATP binding. Combining classical dynamic simulations with a novel ab initio computational approach linear-scalable to molecular interactions involving thousands of atoms, we have investigated the binding of five distinct inhibitors to the cyclin-dependent kinase CDK2. We report here that polarization and dynamic hydrogen bonding effects, so far undetected by crystallography, affect both their activity and selectivity. The effects arise from the specific solvation patterns of water molecules in the ATP binding pocket or the intermittent formation of hydrogen bonds during the dynamics of CDK/inhibitor interactions and explain the unexpectedly high potency of certain inhibitors such as 3-(3H-imidazol-4-ylmethylene)-5-methoxy-1,3-dihydro-indol-2-one (SU9516). The Lys89 residue in the ATP-binding pocket of CDK2 is observed to form temporary hydrogen bonds with the three most potent inhibitors. This residue is replaced in CDK4 by Thr89, whose shorter side-chain cannot form similar bonds, explaining the relative selectivity of the inhibitors for CDK2. Our results provide a generally applicable computational method for the analysis of biomolecular structures and reveal hitherto unrecognized features of the interaction between protein kinases and their inhibitors.

Differential binding of inhibitors to active and inactive CDK2 provides insights for drug design

The cyclin-dependent kinases (CDKs) have been characterized in complex with a variety of inhibitors, but the majority of structures solved are in the inactive form. We have solved the structures of six inhibitors in both the monomeric CDK2 and binary CDK2/cyclinA complexes and demonstrate that significant differences in ligand binding occur depending on the activation state. The binding mode of two ligands in particular varies substantially in active and inactive CDK2. Furthermore, energetic analysis of CDK2/cyclin/inhibitors demonstrates that a good correlation exists between the in vitro potency and the calculated energies of interaction, but no such relationship exists for CDK2/inhibitor structures. These results confirm that monomeric CDK2 ligand complexes do not fully reflect active conformations, revealing significant implications for inhibitor design while also suggesting that the monomeric CDK2 conformation can be selectively inhibited.

Study of the inhibition of cyclin-dependent kinases with roscovitine and indirubin-3'-oxime from molecular dynamics simulations

Molecular dynamics simulations were performed to elucidate the interactions of CDK2 and CDK5 complexes with three inhibitors: R-roscovitine, S-roscovitine, and indirubin-3'-oxime. The preference of the two complexes for R-roscovitine over the S enantiomer, as reported by the experiment, was also found by the simulations. More importantly, the simulations showed that the cause of the stronger affinity for the R enantiomer is the presence of an important hydrogen bond between R-roscovitine and the kinases not found with S-roscovitine. The simulations also showed two amino acid mutations in the active site of CDK5/R-roscovitine that favor binding-enhanced electrostatic contributions, making the inhibitor more effective for CDK5 than for CDK2. This suggests that the effectiveness of roscovitine-like inhibitors can be improved by enhancing their electrostatic interaction with the kinases. Finally, molecular mechanics-Possion-Boltzmann/surface area calculations of the CDK5/indirubin-3'-oxime system in both water-excluded and water-included environments gave significantly different electrostatic contributions to the binding. The simulations detected the displacement of a water molecule in the active site of the water-included CDK/indirubin-3'-oxime system. This resulted in a more conserved binding pattern than the water-excluded structure. Hence, in the design of new indirubin-like inhibitors, it is important to include the water molecule in the analysis.

The role of residue stability in transient protein-protein interactions involved in enzymatic phosphate hydrolysis. A computational study

Finding why protein-protein interactions (PPIs) are so specific can provide a valuable tool in a variety of fields. Statistical surveys of so-called transient complexes (like those relevant for signal transduction mechanisms) have shown a tendency of polar residues to participate in the interaction region. Following this scheme, residues in the unbound partners have to compete between interacting with water or interacting with other residues of the protein. On the other hand, several works have shown that the notion of active site electrostatic preorganization can be used to interpret the high efficiency in enzyme reactions. This preorganization can be related to the instability of the residues important for catalysis. In some enzymes, in addition, conformational changes upon binding to other proteins lead to an increase in the activity of the enzymatic partner. In this article the linear response approximation version of the semimacroscopic protein dipoles Langevin dipoles (PDLD/S-LRA) model is used to evaluate the stability of several residues in two phosphate hydrolysis enzymes upon complexation with their activating partners. In particular, the residues relevant for PPI and for phosphate hydrolysis in the CDK2/Cyclin A and Ras/GAP complexes are analyzed. We find that the evaluation of the stability of residues in these systems can be used to identify not only active site regions but it can also be used as a guide to locate "hot spots" for PPIs. We also show that conformational changes play a major role in positioning interfacing residues in a proper "energetic" orientation, ready to interact with the residues in the partner protein surface. Thus, we extend the preorganization theory to PPIs, extrapolating the results we obtained from the above-mentioned complexes to a more general case. We conclude that the correlation between stability of a residue in the surface and the likelihood that it participates in the interaction can be a general fact for transient PPIs.

Inactivation of O6-alkylguanine DNA alkyltransferase as a means to enhance chemotherapy

DNA adducts at the O6-position of guanine are a result of the carcinogenic, mutagenic and cytotoxic actions of methylating and chloroethylating agents. The presence of the DNA repair protein O6-alkylguanine-DNA alkyltransferase (AGT) renders cells resistant to the biological effects induced by agents that attack at this position. O6-Benzylguanine (O6-BG) is a low molecular weight substrate of AGT and therefore, results in sensitizing cells and tumors to alkylating agent-induced cytotoxicity and antitumor activity. Presently, chemotherapy regimens of O6-BG in combination with BCNU, temozolomide and Gliadel are in clinical development. Other ongoing clinical trials include expression of mutant AGT proteins that confer resistance to O6-BG in bone marrow stem cells, in an effort to reduce the potential enhanced toxicity and mutagenicity of alkylating agents in the bone marrow. O6-BG has also been found to enhance the cytotoxicity of agents that do not form adducts at the O6-position of DNA, including platinating agents. O6-BG's mechanism of action with these agents is not fully understood; however, it is independent of AGT activity or AGT inactivation. A better understanding of the effects of this agent will contribute to its clinical usefulness and the design of better analogs to further improve cancer chemotherapy.

2D autocorrelation modelling of the inhibitory activity of cytokinin-derived cyclin-dependent kinase inhibitors

The inhibitory activity towards p34(cdc2)/cyclin b kinase (CBK) enzyme of 30 cytokinin-derived compounds has been successfully modelled using 2D spatial autocorrelation vectors. Predictive linear and non-linear models were obtained by forward stepwise multi-linear regression analysis (MRA) and artificial neural network (ANN) approaches respectively. A variable selection routine that selected relevant non-linear information from the data set was employed prior to networks training. The best ANN with three input variables was able to explain about 87% data variance in comparison with 80% by the linear equation using the same number of descriptors. Similarly, the neural network had higher predictive power. The MRA model showed a linear dependence between the inhibitory activities and the spatial distributions of masses, electronegativities and van der Waals volumes on the inhibitors molecules. Meanwhile, ANN model evidenced the occurrence of non-linear relationships between the inhibitory activity and the mass distribution at different topological distance on the cytokinin-derived compounds. Furthermore, inhibitors were well distributed regarding its activity levels in a Kohonen self-organizing map (SOM) built using the input variables of the best neural network.

Selectivity and potency of cyclin-dependent kinase inhibitors

Members of the cyclin-dependent kinase (CDK) family play key roles in various cellular processes. There are 11 members of the CDK family known till now. CDKs are activated by forming noncovalent complexes with cyclins such as A-, B-, C-, D- (D1, D2, and D3), and E-type cyclins. Each isozyme of this family is responsible for particular aspects (cell signaling, transcription, etc) of the cell cycle, and some of the CDK isozymes are specific to certain kinds of tissues. Aberrant expression and overexpression of these kinases are evidenced in many disease conditions. Inhibition of isozymes of CDKs specifically can yield beneficiary treatment modalities with minimum side effects. More than 80 3-dimensional structures of CDK2, CDK5, and CDK6 complexed with inhibitors have been published. This review provides an understanding of the structural aspects of CDK isozymes and binding modes of various known CDK inhibitors so that these kinases can be better targeted for drug discovery and design. The amino acid residues that constitute the cyclin binding region, the substrate binding region, and the area around the adenosine triphosphate (ATP) binding site have been compared for CDK isozymes. Those amino acids at the ATP binding site that could be used to improve the potency and subtype specificity have been described.

Protein structures in virtual screening: a case study with CDK2

The influence of protein structure on the successful reproduction of known ligand poses by high-throughput docking programs is rarely discussed. Two commonly used programs, Glide and GOLD, were used to dock a set of CDK2 inhibitors of known bound pose into 20 different CDK2 structures. The numbers of docked poses that reproduced the known pose are reported. Depending on the program and protein structure, 0.3%-96.2% of the ligands docked with the correct pose. Although it is not possible to say that any one structure is "the best" for virtual screening, there are some structures that are clearly better than others. The main determinants of this are the volume of the binding site into which the ligands are docked and the exact orientation of the residues forming the binding site.

Molecular basis for nucleotide-binding specificity: role of the exocyclic amino group "N2" in recognition by a guanylyl-ribonuclease

Proteins interact with nucleotides to perform a multitude of functions within cells. These interactions are highly specific; however, the molecular basis for this specificity is not well understood. To identify factors critical for protein-guanine nucleotide recognition the binding of two closely related ligands, guanosine 3'-monophosphate (3'GMP) and inosine 3'-monophosphate (3'IMP), to Ribonuclease Sa (RNase Sa), a small, guanylyl-endoribonuclease from Streptomyces aureofaciens, was compared using isothermal titration calorimetry, NMR, X-ray crystallography and molecular dynamics simulations. This comparison has allowed for the determination of the contribution of the exocyclic amino group "N2" of the guanine base to nucleotide binding specificity. Calorimetric measurements indicate that RNase Sa has a higher affinity for 3'GMP (K=(1.5+/-0.2)x10(5)) over 3'IMP (K=(3.1+/-0.2)x10(4)) emphasizing the importance of N2 as a key determinant of RNase Sa guanine binding specificity. This result was unexpected as the published structural data for RNase Sa in complex with 3'GMP showed only a potential long-range interaction (>3.3A) between N2 and the side-chain of Glu41 of RNase Sa. The observed difference in affinity is largely due to a reduction in the favorable enthalpy change by 10 kJ/mol for 3'IMP binding as compared to 3'GMP that is accompanied by a significant difference in the heat capacity changes observed for binding the two ligands. To aid interpretation of the calorimetric data, the first crystal structure of a small, guanylyl ribonuclease bound to 3'IMP was determined to 2.0 A resolution. This structure has revealed small yet unexpected changes in the ligand conformation and differences in the conformations of the side-chains contacting the sugar and phosphate moieties as compared to the 3'GMP complex. The structural data suggest the less favorable enthalpy change is due to an overall lengthening of the contacts between RNase Sa and 3'IMP as compared to 3'GMP. The long-range interaction between N2 and Glu41 is critical for positioning of the nucleotide in the binding cleft for optimal contact formation. Thus, combined, the data demonstrate how a long-range interaction can have a significant impact on nucleotide binding affinity and energetics.

Modeling of Cyclin-Dependent Kinase Inhibition by 1H-Pyrazolo[3,4-d]Pyrimidine Derivatives Using Artificial Neural Network Ensembles

Artificial neural network ensembles were used for modeling the cyclin-dependent kinase inhibition of 1H-pyrazolo[3,4-d]pyrimidine derivatives. The structural characteristics of these inhibitors were encoded in relevant 3D-spatial descriptors extracted by genetic algorithm feature selection. Bayesian-regularized multilayer neural networks, trained by the back-propagation algorithm, were developed using these variables as inputs. The predictive power of the model was tested by leave-one-out cross validation. In addition, for a more rigorous measure of the predictive capacity, multiple validation sets were randomly generated as members of neural network ensembles, which makes doing averaged predictions feasible. In this way, the predictive power was analyzed accounting for the averaged test set R values and test set mean-square errors. Otherwise, Kohonen self-organizing maps were used as an additional tool for the same modeling. The location of the inhibitors in a map facilitates the analysis of the connection between compounds and serves as a useful tool for qualitative predictions.

Potentiation of paclitaxel-induced apoptosis by the novel cyclin-dependent kinase inhibitor NU6140: a possible role for survivin down-regulation

Cyclin-dependent kinases (CDK) play a crucial role in the control of the cell cycle. Aberrations in the control of cell cycle progression occur in the majority of human malignancies; hence, CDKs are promising targets for anticancer therapy. Here, we define the cellular effects of the novel CDK inhibitor NU6140, alone or in association with paclitaxel, with respect to inhibition of cell proliferation and cell cycle progression and induction of apoptosis in HeLa cervical carcinoma cells and in comparison with purvalanol A. Both CDK inhibitors induced a concentration-dependent cell cycle arrest at the G2-M phase and an increase in the apoptotic rate, with a concomitant down-regulation of the antiapoptotic protein survivin, a member of the inhibitors of apoptosis protein family. Notably, the addition of NU6140 to paclitaxel-treated cells resulted in markedly increased cytotoxic effect and apoptotic response in comparison with the paclitaxel-purvalanol A combination (86 ± 11% and 37 ± 8%, respectively). Similarly, the extent of caspase-9 and caspase-3 activation in paclitaxel-NU6140–treated cells was ∼4-fold higher than after the paclitaxel-purvalanol A combination. Moreover, an almost complete abrogation of the expression of the active, Thr34-phosphorylated form of survivin was observed in cells exposed to the paclitaxel-NU6140 combination. A synergistic effect of the paclitaxel-NU6140 combination, as a consequence of survivin inhibition and increased activation of caspase-9 and caspase-3, was also observed in OAW42/e ovarian cancer line but not in the derived OAW42/Surv subline ectopically expressing survivin. Results from this study indicate that NU6140 significantly potentiates the apoptotic effect of paclitaxel, with inhibition of survivin expression/phosphorylation as the potential mechanism.

8-Azapurines as new inhibitors of cyclin-dependent kinases

Purine inhibitors of cyclin-dependent kinases (CDK) seem to be a potential anticancer drug candidate as one of the first representatives, roscovitine, is passing Phase II clinical trials for cancer and glomerulonephritis. In this article, we describe a novel modification of the purine scaffold influencing CDK2 inhibitory activities as well as anticancer properties in cell lines of different histopathological origin. The introduced N at position 8 of the purine ring generally lowered CDK2 inhibitory activity of new 8-azapurines (1,2,3-triazolo[4,5-d]pyrimidines) in comparison to the model trisubstituted purines, whereas the antiproliferative potential of some derivatives remained very high, reflecting their ability to activate p53 tumor suppressor.

Synthetic, spectral, magnetic and in vitro cytotoxic activity studies of cobalt(II) complexes with cytokinin derivatives: X-ray structure of 6-(3-methoxybenzylamino)purinium chloride monohydrate

Cobalt(II) complexes with 6-(2-hydroxybenzylamino)purine (HL1), 6-(2-methoxybenzylamino)purine (HL2), 6-(3-methoxybenzylamino)purine (HL3) and 6-(4-methoxybenzylamino)purine (HL4) of the composition [Co(L1)Cl(H2O)2].H2O (1), [Co(L2)Cl(H2O)2] (2), [Co(L3)2(H2O)2].2H2O (3), [Co(L4)2(H2O)2].2H2O (4) have been synthesized. The compounds have been characterized by elemental analysis, FT-IR, ES+ MS (electrospray mass spectra in the positive ion mode) and electronic spectroscopies, magnetic and conductivity data as tetrahedral high-spin cobalt(II) complexes. The thermal stability of the complexes has also been studied. The cytotoxicity of the complexes (1-4) was determined by a Calcein acetoxymethyl (AM) assay. Human malignant melanoma (G361), human chronic myelogenous erythroleukemia (K562), human osteogenic sarcoma (HOS) and human breast adenocarcinoma (MCF7) cell lines were used for the testing. The molecular structure of 6-(3-methoxybenzylamino)purinium chloride monohydrate, H2L3+.Cl.H2O, i.e. a protonated form of the free HL(3) ligand, has been determined by a single crystal X-ray analysis. The geometry optimisation and infrared frequencies calculations of HL1, HL2, and H2L3+ and H2L4+ were performed using density-functional theory (DFT) calculations at the B3LYP/6-31G* level of the theory. The geometry of complex (1) was optimised at the same level of the theory.

Identifying the binding mode of a molecular scaffold

We describe a method for docking of a scaffold-based series and present its advantages over docking of individual ligands, for determining the binding mode of a molecular scaffold in a binding site. The method has been applied to eight different scaffolds of protein kinase inhibitors (PKI). A single analog of each of these eight scaffolds was previously crystallized with different protein kinases. We have used FlexX to dock a set of molecules that share the same scaffold, rather than docking a single molecule. The main mode of binding is determined by the mode of binding of the largest cluster among the docked molecules that share a scaffold. Clustering is based on our 'nearest single neighbor' method [J. Chem. Inf. Comput. Sci., 43 (2003) 208-217]. Additional criteria are applied in those cases in which more than one significant binding mode is found. Using the proposed method, most of the crystallographic binding modes of these scaffolds were reconstructed. Alternative modes, that have not been detected yet by experiments, could also be identified. The method was applied to predict the binding mode of an additional molecular scaffold that was not yet reported and the predicted binding mode has been found to be very similar to experimental results for a closely related scaffold. We suggest that this approach be used as a virtual screening tool for scaffold-based design processes.

Molecular models of cyclin-dependent kinase 1 complexed with inhibitors

Roscovitine and flavopiridol have been shown to potently inhibit cyclin-dependent kinase 1 and 2 (CDK1 and 2). The structures of CDK2 complexed with roscovitine and deschoroflavopiridol have been reported, however no crystallographic structure is available for complexes of CDK1 with inhibitors. The present work describes two molecular models for the binary complexes CDK1:roscovitine and CDK1:flavopiridol. These structural models indicate that both inhibitors strongly bind to the ATP-binding pocket of CDK1 and structural comparison of the CDK complexes correlates the structures with differences in inhibition of these CDKs by flavopiridol and roscovitine. This article explains the structural basis for the observed differences in activity of these inhibitors.

N2-substituted O6-cyclohexylmethylguanine derivatives: potent inhibitors of cyclin-dependent kinases 1 and 2

The adenosine 5'-triphosphate (ATP) competitive cyclin-dependent kinase inhibitor O(6)-cyclohexylmethylguanine (NU2058, 1) has been employed as the lead in a structure-based drug discovery program resulting in the discovery of the potent CDK1 and -2 inhibitor NU6102 (3, IC(50) = 9.5 nM and 5.4 nM vs CDK1/cyclinB and CDK2/cyclinA3, respectively). The SAR for this series have been explored further by the synthesis and evaluation of 45 N(2)-substituted analogues of NU2058. These studies have confirmed the requirement for the hydrogen bonding N(2)-NH group and the requirement for an aromatic N(2)-substituent to confer potency in the series. Additional potency is conferred by the presence of a group capable of donating a hydrogen bond at the 4'-position, for example, the 4'-hydroxy derivative (25, IC(50) = 94 nM and 69 nM vs CDK1/cyclinB and CDK2/cyclinA3, respectively), 4'-monomethylsulfonamide derivative (28, IC(50) = 9 nM and 7.0 nM vs CDK1/cyclinB and CDK2/cyclinA3, respectively), and 4'-carboxamide derivative (34, IC(50) = 67 nM and 64 nM vs CDK1/cyclinB and CDK2/cyclinA3, respectively). X-ray crystal structures have been obtained for key compounds and have been used to explain the observed trends in activity.

Effect of cell cycle inhibition on Cisplatin-induced cytotoxicity

Pharmacological inhibitors of cyclin-dependent kinase (CDK)2 are currently in preclinical and clinical development. The purpose of our work was to evaluate a series of guanine derivatives for their ability to inhibit CDK2, affect cell cycle progression, and enhance the cytotoxic and apoptotic effects of cisplatin. A panel of guanine derivatives, including O(6)-benzylguanine (O(6)-BG), S(6)-benzyl-6-thioguanine (S(6)-BG), S(6)-[(cyclohexyl)methyl]-6-thioguanine (S(6)-CMG), O(6)-[(cyclohexyl)methyl]guanine (O(6)-CMG), O(6)-benzyl-9-methylguanine (9-CH(3)-BG), O(6)-[(cyclohexyl)methyl]-9-methyl-guanine (9-CH(3)-CMG), and 7-benzylguanine (N7-BG), exhibited varying degrees of CDK2 inhibition with O(6)-CMG being the most potent and 9-CH(3)-BG, 9-CH(3)-CMG, and N7-BG the least potent compounds. Treatment with S(6)-CMG and O(6)-CMG significantly decreased the percentage of cells in S phase. In SQ20b and SCC61 head and neck cancer cell lines, the most potent CDK2 inhibitor, O(6)-CMG, was also the most effective at enhancing cisplatin-induced cytotoxicity and apoptosis. Cisplatin-induced DNA platination increased in SQ20b cells pretreated with S(6)-BG, S(6)-CMG, and O(6)-CMG. Treatment with both O(6)-BG and trichostatin A, an indirect cell cycle inhibitor, demonstrated additive effects on cisplatin-induced cytotoxicity. In summary, we have identified a group of guanine derivatives that were effective modulators of cisplatin-induced cytotoxicity and apoptosis.

A minimalist approach to fragment-based ligand design using common rings and linkers: Application to kinase inhibitors

We present a novel method for stepwise scaffold assembly that integrates fragment-by-fragment ligand design approaches with high-throughput virtual library screening (COREGEN). As an extension of our earlier studies of common features present in drug molecules, we investigate the hypothesis that most pharmaceutically interesting ligands can be expressed in terms of the ring-linker frameworks that comprise them. Analysis of 119 published kinase inhibitors from at least 18 different targets illustrates that a basis set of 4 rings and 8 linkers is sufficient to describe approximately 90% of ring and linker occurrences, respectively. A similar result was derived from a larger set of approximately 40,000 kinase inhibitors from curated patents. A method for ring-linker-based assembly of scaffold libraries that uses experimental information to guide the placement of anchor fragments is validated using a set of reported kinase inhibitors of Bcr-Abl, Cdk2, and Src. In every case, the predominant structural motif of reported ligand cores is reproduced and variations are suggested. To underscore generality of this approach, a novel scaffold for a cyclooxygenase-2 (COX-2) selective ligand is proposed.

General Method for the Synthesis of 8-Arylsulfanyl Adenine Derivatives

We report a general method for the synthesis of 8-arylsulfanyl adenine derivatives using a mild protocol of coupling 8-mercaptoadenine with a variety of aryl iodides.

Pediatric neuroblastomas: genetic and epigenetic 'danse macabre'

Neuroblastomas are the most frequently occurring solid tumors in children under 5 years. Spontaneous regression is more common in neuroblastomas than in any other tumor type, especially in young patients under 12 months. Unfortunately, the full clinical spectrum of neuroblastomas also includes very aggressive tumors, unresponsive to multi-modality treatment and accounting for most of the pediatric cancer mortalities under 5 years of age. It is generally emphasized that more than one biological entity of neuroblastoma exists. Structural genetic defects such as amplification of MYCN, gain of chromosome 17q and LOH of 1p and several other chromosomal regions have proven to be valuable as prognostic factors and will be discussed in relation to their clinical relevance. Recent research is starting to uncover important molecular pathways involved in the pathogenesis of neuroblastomas. The aim of this review is to discuss several important aspects of the biology of the neuroblast, such as the role of overexpressed oncogenes like MYCN and cyclin D1, the mechanisms leading to decreased apoptosis, like overexpression of BCL-2, survivin, NM23, epigenetic silencing of caspase 8 and the role of tumor suppressor genes, like p53, p73 and RASSF1A. In addition, the role of specific proteins overexpressed in neuroblastomas, such as the neurotrophin receptors TrkA, B and C in relation to spontaneous regression and anti-angiogenesis will be discussed. Finally, we will try to relate these pathways to the embryonal origin of neuroblastomas and discuss possible new avenues in the therapeutic approach of future neuroblastoma patients.

Inhibitors of JAKs/STATs and the kinases: a possible new cluster of drugs

JAK(s)/STAT(s) relay cytokine signals through tyrosine site-specific phosphorylation of the proteins involved in cellular responses for the activation and proliferation of bone marrow-derived cells. In recent years, the constitutive or elevated expression of JAK/STAT has been found in cancer cells and oncogene transfected cells, and has been shown to be involved in the immune rejection of allografts and the inflammatory processes of autoimmune diseases. This review discusses the strategies for screening and rational design of selective, potent JAK/STAT and kinase inhibitors that are either ATP-competitive or non-ATP competitive, naturally derived or synthetic, as well as other unique inhibitors and analogues for different therapeutic indications.

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.

Oxindole-based compounds are selective inhibitors of Plasmodium falciparum cyclin dependent protein kinases

Cyclin dependent protein kinases (CDKs) have become attractive drug targets in an effort to identify effective inhibitors of the parasite Plasmodium falciparum, the causative agent of the most severe form of human malaria. We tested known CDK inhibitors for their ability to inhibit two malarial CDKs: Pfmrk and PfPK5. Many broad spectrum CDK inhibitors failed to inhibit Pfmrk suggesting that the active site differs from other CDKs in important ways. By screening compounds in the Walter Reed chemical database, we identified oxindole-based compounds as effective inhibitors of Pfmrk (IC(50) = 1.5 microM). These compounds have low cross-reactivity against PfPK5 and human CDK1 demonstrating selectivity for Pfmrk. Amino acid comparison of the active sites of Pfmrk and PfPK5 identified unique amino acid differences that may explain this selectivity and be exploited for further drug development efforts.

The cell cycle: a review of regulation, deregulation and therapeutic targets in cancer

The cell cycle is controlled by numerous mechanisms ensuring correct cell division. This review will focus on these mechanisms, i.e. regulation of cyclin-dependent kinases (CDK) by cyclins, CDK inhibitors and phosphorylating events. The quality checkpoints activated after DNA damage are also discussed. The complexity of the regulation of the cell cycle is also reflected in the different alterations leading to aberrant cell proliferation and development of cancer. Consequently, targeting the cell cycle in general and CDK in particular presents unique opportunities for drug discovery. This review provides an overview of deregulation of the cell cycle in cancer. Different families of known CDK inhibitors acting by ATP competition are also discussed. Currently, at least three compounds with CDK inhibitory activity (flavopiridol, UCN-01, roscovitine) have entered clinical trials.

Solution- and solid-phase parallel synthesis of 4-alkoxy-substituted pyrimidines with high molecular diversity

A simple and straightforward methodology toward the synthesis of novel 2,6-disubstituted-4-alkoxypyrimidine derivatives of type 16 and 19 has been developed. This methodology, initially developed in solution, can be perfectly adapted to the solid support under analogous conditions, taking full advantage of automated parallel synthesis systems. This successful methodology benefits from the key role played by the thioether linkage placed at the 2-position in 3, 9, or 13 in a double manner: on one side, the steric effect exerted by the thioether linkage is likely to be responsible for the very high observed selectivity toward the formation of the O-alkylation products. On the other side, this sulfur linkage can serve not only as a robust point of attachment for the heterocycle, stable to a number of reaction conditions, but also as a means of introducing a new element of diversity through activation to the corresponding sulfone (safety-catch linker concept) and subsequent ipso-substitution reaction with a variety of different N-nucleophiles.

Synthesis of quinolinyl/isoquinolinyl[a]pyrrolo [3,4-c] carbazoles as cyclin D1/CDK4 inhibitors

A novel series of pyrrolo[3,4-c] carbazoles fused with a quinolinyl/isoquinolinyl moiety were synthesized and their D1/CDK4 inhibitory and antiproliferative activity were evaluated. Compound 8H, 14H-isoquinolinyl[6,5-a]-pyrrolo[3,4-c]carbazole-7,9-dione (1d) was found to be a highly potent D1/CDK4 inhibitor with an IC(50) of 69 nM. Compound 1d also inhibited tumor cell growth, arrested tumor cells in G1 phase and inhibited pRb phosphorylation.

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.

Rearrangements and increased expression of cyclin D1 (CCND1) in neuroblastoma

Cyclin D1 regulates G1 cell cycle progression by controlling the phosphorylation of the retinoblastoma protein. This pathway is frequently deregulated in many malignancies. In neuroblastoma, however, no consistent G1 cell cycle checkpoint aberrations have been found. We examined the possible deregulation of cyclin D1 (CCND1) in this tumor. mRNA expression profiles of neuroblastoma generated by SAGE (Serial Analysis of Gene Expression) revealed a high expression of CCND1 in a subset of neuroblastoma cell lines and tumors. The CCND1 expression level can be 0.3% of the total cellular mRNA. Northern blot analysis of CCND1 expression showed a relative overexpression in 16 of 23 neuroblastoma cell lines and 10 of 15 tumor samples. In the majority of cases, the high CCND1 mRNA levels also led to high CCND1 protein levels. In the search for mechanisms causing this relative overexpression, we screened for amplifications and rearrangements of CCND1. Five amplifications were found in 202 neuroblastoma tumors and cell lines. Analysis of the 3'-UTR of CCND1 showed a rearrangement in 1 of 96 tumors. These clonal aberrations of CCND1 together with the high expression suggest a role for deregulated CCND1 activity in neuroblastoma tumorigenesis.

4-Alkoxy-2,6-diaminopyrimidine derivatives: inhibitors of cyclin dependent kinases 1 and 2

The cyclin dependent kinase (cdk) inhibitor NU6027, 4-cyclohexylmethoxy-5-nitroso-pyrimidine-2,6-diamine (IC(50) vs cdk1/cyclinB1=2.9+/-0.1 microM and IC(50) vs cdk2/cyclinA3=2.2+/-0.6 microM), was used as the basis for the design of a series of 4-alkoxy-2,6-diamino-5-nitrosopyrimidine derivatives. The synthesis and evaluation of 21 compounds as potential inhibitors of cyclin-dependent kinases 1 and 2 is described and the structure-activity relationships relating to NU6027 have been probed. Simple alkoxy- or cycloalkoxy-groups at the O(4)-position were tolerated, with the 4-(2-methylbutoxy)-derivative (IC(50) vs cdk1/cyclinB1=12+/-2 microM and cdk2/cyclinA3=13+/-4 microM) retaining significant activity. Substitutions at the N(6) position were not tolerated. Replacement of the 5-nitroso substituent with ketone, oxime and semicarbazone groups essentially abolished activity. However, the derivative bearing an isosteric 5-formyl group, 2,6-diamino-4-cyclohexylmethoxy-pyrimidine-5-carbaldehyde, showed modest activity (IC(50) vs cdk1/cyclinB1=35+/-3 microM and cdk2/cyclinA3=43+/-3 microM). The X-ray crystal structure of the 5-formyl compound bound to cdk2 has been determined to 2.3A resolution. The intramolecular H-bond deduced from the structure with NU6027 bound to cdk2 is not evident in the structure with the corresponding formyl compound. Thus the parent compound, 4-cyclohexylmethoxy-5-nitrosopyrimidine-2,6-diamine (NU6027), remains the optimal basis for future structure-activity studies for cyclin-dependent kinase inhibitors in this series.

Aloisines, a new family of CDK/GSK-3 inhibitors. SAR study, crystal structure in complex with CDK2, enzyme selectivity, and cellular effects

Cyclin-dependent kinases (CDKs) regulate the cell cycle, apoptosis, neuronal functions, transcription, and exocytosis. The observation of CDK deregulations in various pathological situations suggests that CDK inhibitors may have a therapeutic value. In this article, we report on the identification of 6-phenyl[5H]pyrrolo[2,3-b]pyrazines (aloisines) as a novel potent CDK inhibitory scaffold. A selectivity study performed on 26 kinases shows that aloisine A is highly selective for CDK1/cyclin B, CDK2/cyclin A-E, CDK5/p25, and GSK-3 alpha/beta; the two latter enzymes have been implicated in Alzheimer's disease. Kinetic studies, as well as the resolution of a CDK2-aloisine cocrystal structure, demonstrate that aloisines act by competitive inhibition of ATP binding to the catalytic subunit of the kinase. As observed with all inhibitors reported so far, aloisine interacts with the ATP-binding pocket through two hydrogen bonds with backbone nitrogen and oxygen atoms of Leu 83. Aloisine inhibits cell proliferation by arresting cells in both G1 and G2.

Facilitation of addition–elimination reactions in pyrimidines and purines using trifluoroacetic acid in trifluoroethanol

SNAr displacement reactions of 6-cyclohexylmethoxy-2-fluoropurine, 6-amino-2-butylsulfonyl-4-cyclohexylmethoxypyrimidine and 2-amino-6-chloropurine with substituted anilines (e.g. the weakly nucleophilic 4-aminobenzenesulfonamide) are dramaticallyaccelerated in the presence of trifluoroacetic acid and occur especially efficiently in 2,2,2-trifluoroethanol solvent.

Structural bioinformatics study of cyclin-dependent kinases complexed with inhibitors

The present work describes molecular models for the binary complexes CDK9, CDK5 and CDK1 complexed with Flavopiridol and Roscovitine. These structural models indicate that the inhibitors strongly bind to the ATP-binding pocket of CDKs and the structural comparison with the complexes CDK2:Flavopiridol and CDK2:Roscovitine correlates the structural differences with differences in inhibition of these CDKs by the inhibitors. These structures open the possibility of testing new inhibitor families, in addition to new substituents for the already known lead structures such as flavones and adenine derivatives.

Synthesis and evaluation of indenopyrazoles as cyclin-dependent kinase inhibitors. 2. Probing the indeno ring substituent pattern

We disclose a novel series of indenopyrazole-based cyclin-dependent kinase (CDK) inhibitors. Kinetic experiments confirmed our initial molecular modeling studies that the compounds are competitive with respect to adenosine 5'-triphosphate (ATP) and bind in the kinase ATP pocket. A unique combination of active pharmacophores led us to a series of semicarbazide-based inhibitors that are highly potent against CDK2 and CDK4 while maintaining selectivity against other relevant serine/threonine kinases. These compounds were active against a transformed human colon cancer cell line (HCT116) while maintaining an acceptable margin of activity against a normal fibroblast cell line. The compounds were found to be highly protein bound in our cell-based assay with the exception of 11k, which maintained a reasonable level of activity in the presence of human plasma proteins.

Dynamics and binding modes of free cdk2 and its two complexes with inhibitors studied by computer simulations

This article presents a molecular dynamics (MD) study of the cdk2 enzyme and its two complexes with the inhibitors isopentenyladenine and roscovitine using the Cornell et al. force field from the AMBER software package. The results show that inserting an inhibitor into the enzyme active site does not considerably change enzyme structure but it seemingly changes the distribution of internal motions. The inhibitor causes differences in the domain motions in free cdk2 and in its complexes. It was found out that repulsion of roscovitine N9 substituent causes conformational change on Lys 33 side chain. Isopentenyladenine forms with Lys 33 side chain terminal amino group a hydrogen bond. It implies that the cavity, where N9 substituent of roscovitine is buried, can adopt larger substituent due to Lys 33 side chain flexibility. The composition of electrostatic and van der Waals interactions between the inhibitor and the enzyme were also calculated along both cdk2/inhibitor MD trajectories together with MM-PB/GBSA analysis. These results show that isopentenyladenine-like inhibitors could be more effective after modifications leading to an increase in their van der Waals contact with the enzyme. We suggest that a way leading to better inhibitors occupying isopentenyladenine binding mode could be: to keep N9 and N7 purine positions free, to keep 3,3-dimethylallylamino group at C6 position, and to add, e.g., benzylamino group at C2 position. The results support the idea that the isopentenyladenine binding mode can be used for cdk2 inhibitors design and that all possibilities to improve this binding mode were not uncovered yet.

SU9516, a cyclin-dependent kinase 2 inhibitor, promotes accumulation of high molecular weight E2F complexes in human colon carcinoma cells

The E2F family plays a critical role in the expression of genes required for entry into and progression through S phase. E2F-mediated transcription is repressed by the tumor suppressor retinoblastoma protein (pRb), which results in sequestration of E2F in a multiprotein complex that includes pRb. Derepression of E2F results from a series of complex phosphorylation events mediated by cyclin D/cdk4 and cyclin E/cdk2. We have employed a novel 3-substituted indolinone compound, 3-[1-(3H-imidazol-4-yl)-meth-(Z)-ylidene]-5-methoxy-1,3-dihydro-indol-2-one (SU9516), which selectively inhibits cdk2 activity (Lane et al., Cancer Res 2001;61:6170-7) to investigate these events. Electrophoretic mobility gel shift assays were performed on SU9516-treated and -untreated HT-29, SW480, and RKO human colon cancer cell extracts. Treatment with 5 microM SU9516 prevented dissociation of pRb from E2F1 in all cell lines (HT-29>RKO>SW480). Treatment effects were time-dependent, demonstrating greater inhibition at 48 hr versus 24hr in HT-29 cells. Furthermore, E2F species were sequestered in complexes with p107, p130, DP-1, and cyclins A and E. After a 24-hr treatment with 5 microM SU9516, cyclin D1 and cdk2 levels decreased by 10-60%. These findings delineate a previously undescribed mechanism for SU9516-mediated cell growth arrest through down-regulation of cyclin D1, inhibition of cdk2 levels and activity, and pan-sequestration of E2F.

p42/p44 MAPKs are intracellular targets of the CDK inhibitor purvalanol

Chemical inhibitors of cyclin-dependent kinases (CDKs) have a great therapeutic potential against various proliferative and neurodegenerative disorders. Intensive screening of a combinatorial chemistry library of 2,6,9-trisubstituted purines has led to the identification of purvalanol, one of the most potent and selective CDK inhibitors to date. In preliminary studies, this compound demonstrates definite anti-mitotic properties, consistent with its nanomolar range efficiency towards purified CDK1 and CDK2. However, the actual intracellular targets of purvalanol remain to be identified, and a method for the determination of its in vivo selectivity was developed. In this technique, cell extracts were screened for purvalanol-interacting proteins by affinity chromatography on immobilized inhibitor. In addition to CDK1, p42/p44 MAPK were found to be two major purvalanol-interacting proteins in five different mammalian cell lines (CCL39, PC12, HBL100, MCF-7 and Jurkat cells), suggesting the generality of the purvalanol/p42/p44 MAPK interaction. The Chinese hamster lung fibroblast cell line CCL39 was used as a model to investigate the anti-proliferative properties of purvalanol. The compound inhibited cell growth with a GI(50) value of 2.5 microM and induced a G2/M block when added to exponentially growing cells. It did not appear to trigger massive activation of caspase. We next tested whether CDKs and p42/p44 MAPK were actually targeted by the compound in vivo. p42/p44 MAPK activity was visualized using an Elk-Gal4 luciferase reporter system and CDK1 activity was detected by the phosphonucleolin level. When cells were treated with purvalanol, p42/p44 MAPK and CDK1 activities were inhibited in a dose-dependent manner. Furthermore, purvalanol inhibited the nuclear accumulation of p42/p44 MAPK, an event dependent on the catalytic activity of these kinases. We conclude that the anti-proliferative properties of purvalanol are mediated by inhibition of both p42/p44 MAPK and CDKs. These observations highlight the potency of moderate selectivity compounds and encourage the search for new therapeutics which simultaneously target distinct but relevant pathways of cell proliferation.

Pharmacological inhibitors of cyclin-dependent kinases

Cyclin-dependent kinases (CDKs) regulate the cell division cycle, apoptosis, transcription and differentiation in addition to functions in the nervous system. Deregulation of CDKs in various diseases has stimulated an intensive search for selective pharmacological inhibitors of these kinases. More than 50 inhibitors have been identified, among which >20 have been co-crystallized with CDK2. These inhibitors all target the ATP-binding pocket of the catalytic site of the kinase. The actual selectivity of most known CDK inhibitors, and thus the underlying mechanism of their cellular effects, is poorly known. Pharmacological inhibitors of CDKs are currently being evaluated for therapeutic use against cancer, alopecia, neurodegenerative disorders (e.g. Alzheimer's disease, amyotrophic lateral sclerosis and stroke), cardiovascular disorders (e.g. atherosclerosis and restenosis), glomerulonephritis, viral infections (e.g. HCMV, HIV and HSV) and parasitic protozoa (Plasmodium sp. and Leishmania sp.).

Complexities in the development of cyclin-dependent kinase inhibitor drugs

Abnormalities in the normal regulation of the cell cycle are a hallmark of neoplasia. Drugs directed against the cyclin-dependent kinases (CDKs), which govern the normal orderly progression through the cell cycle, have been proposed to address the pathogenic defect in tumors. Recently, CDK family members that do not regulate the cell cycle directly but instead influence transcription (CDK7, CDK8, and CDK9) and neuronal and secretory cell function (CDK5) have been described. Continued synthetic chemistry efforts have defined important new selective inhibitors of CDKs, and strategies directed at newly described CDK-related targets, such as transcription control, can now be envisaged. CDKs remain important and novel targets whose potential needs to be more fully explored, albeit in light of the newly emerging complexities of their cellular physiology.

Probing the ATP ribose-binding domain of cyclin-dependent kinases 1 and 2 with O(6)-substituted guanine derivatives

O(6)-substituted guanines are adenosine 5'-triphosphate (ATP) competitive inhibitors of CDK1/cyclin B1 and CDK2/cyclin A, the O(6) substituent occupying the kinase ribose binding site. Fifty-eight O(6)-substituted guanines were prepared to probe the ribose pocket, and the structures of four representative compounds bound to monomeric CDK2 were determined by X-ray crystallography. Optimum binding occurs with a moderately sized aliphatic O(6) substituent that packs tightly against the hydrophobic patch presented by the glycine loop, centered on Val18, an interaction promoted by the conformational restraints imposed in a cyclohexylmethyl or cyclohexenylmethyl ring. Structure-based design generated (R)-(2-amino-9H-purin-6-yloxymethyl)pyrrolidin-2-one (56), which reproduces the reported hydrogen bonds formed between ATP and Asp86 and Gln131 but failed to improve inhibitory potency. Thus, the parent compound O(6)-cyclohexylmethylguanine (NU2058, 25) is the preferred starting point for exploring other areas of the kinase active site.

A comparison of the pharmacophore identification programs: Catalyst, DISCO and GASP

Three commercially available pharmacophore generation programs, Catalyst/HipHop, DISCO and GASP, were compared on their ability to generate known pharmacophores deduced from protein-ligand complexes extracted from the Protein Data Bank. Five different protein families were included Thrombin, Cyclin Dependent Kinase 2, Dihydrofolate Reductase, HIV Reverse Transcriptase and Thermolysin. Target pharmacophores were defined through visual analysis of the data sets. The pharmacophore models produced were evaluated qualitatively through visual inspection and according to their ability to generate the target pharmacophores. Our results show that GASP and Catalyst outperformed DISCO at reproducing the five target pharmacophores.

Structural classification of protein kinases using 3D molecular interaction field analysis of their ligand binding sites: target family landscapes

Protein kinases are critical components of signaling pathways and trigger various biological events. Several members of this superfamily are interesting targets for novel therapeutic approaches. All known eukaryotic protein kinases exhibit a conserved catalytic core domain with an adenosine 5'-triphosphate (ATP) binding site, which often is targeted in drug discovery programs. However, as ATP is common to kinases and other proteins, specific protein-ligand interactions are crucial prerequisites for valuable ATP site-directed ligands. In the present study, a set of 26 X-ray structures of eukaryotic protein kinases were classified into subfamilies with similar protein-ligand interactions in the ATP binding site using a chemometrical approach based on principal component analysis (PCA) and consensus PCA. This classification does not rely on protein sequence similarities, as descriptors are derived from three-dimensional (3D) binding site information only computed using GRID molecular interaction fields. The resulting classification, which we refer to as "target family landscape", lead to the identification of common binding pattern and specific interaction sites for particular kinase subfamilies. Moreover, those findings are in good agreement with experimental selectivity profiles for a series of 2,6,9-substituted purines as CDK inhibitors. Their interpretation in structural terms unveiled favorable substitutions toward selective CDK inhibitors and thus allowed for a rational design of specific ligands with minimized side effects. Additional 3D-quantitative structure-activity relationship (QSAR) analyses of a larger set of CDK-directed purines lead to the identification of essential structural requirements for affinity in this CDK ATP binding site. The combined interpretation of 3D-QSAR and the kinase target family landscape provides a consistent view to protein-ligand interactions, which are both favorable for affinity and selectivity in this important subfamily.