Three-dimensional model of the cyclin-dependent kinase 1 (CDK1): Ab initio active site parameters for molecular dynamics studies of CDKS

Structural insights and influence of V599 mutations on the overall dynamics of BRAF protein against its kinase domains

Mutations in the BRAF gene are well known for their oncogenic effects. Point mutations in V599 are particularly oncogenic and are considered important for therapeutic purposes. Along with wild type, other V599 mutated BRAF variants viz. V599E, V599D and V599R are reported and crystals of the former two with inhibitor (BAY43-9006) are further detailed. Both wild-type and mutated BRAF forms show similar interaction patterns with BAY43-9006, but the 599th residue did not show any involvement in the interactions. Upon BAY43-9006 binding, kinase domains of both forms were found adopting essentially identical conformations. However, BAY43-9006 shows a varied activity profile in the case of the wild and V599E variant of the BRAF protein. Furthermore, MMGBSA binding energy results for all four BRAF variants, further revealed the importance of the 599th residue. In-depth analysis viz. molecular dynamics, residue correlation studies and residue interaction network (RIN) analyses were conducted, providing a deep insight into the 599th residue and its impact on the overall dynamics of BRAF protein. Our findings reveal that the mutated residue at the 599th position not only changed the BAY43-9006-BRAF binding behaviour but also produced a massive impact on the overall dynamic behaviour of the protein. The insights obtained herein could be of great relevance for designing new BRAF inhibitors aimed at getting ideal activity against all BRAF forms.

Canonical and Alternative Pathways in Cyclin-Dependent Kinase 1/Cyclin B Inactivation upon M-Phase Exit in Xenopus laevis Cell-Free Extracts

Cyclin-Dependent Kinase 1 (CDK1) is the major M-phase kinase known also as the M-phase Promoting Factor or MPF. Studies performed during the last decade have shown many details of how CDK1 is regulated and also how it regulates the cell cycle progression. Xenopus laevis cell-free extracts were widely used to elucidate the details and to obtain a global view of the role of CDK1 in M-phase control. CDK1 inactivation upon M-phase exit is a primordial process leading to the M-phase/interphase transition during the cell cycle. Here we discuss two closely related aspects of CDK1 regulation in Xenopus laevis cell-free extracts: firstly, how CDK1 becomes inactivated and secondly, how other actors, like kinases and phosphatases network and/or specific inhibitors, cooperate with CDK1 inactivation to assure timely exit from the M-phase.

Pharmacological targeting of CDK9 in cardiac hypertrophy

Cardiac hypertrophy allows the heart to adapt to workload, but persistent or unphysiological stimulus can result in pump failure. Cardiac hypertrophy is characterized by an increase in the size of differentiated cardiac myocytes. At the molecular level, growth of cells is linked to intensive transcription and translation. Several cyclin-dependent kinases (CDKs) have been identified as principal regulators of transcription, and among these CDK9 is directly associated with cardiac hypertrophy. CDK9 phosphorylates the C-terminal domain of RNA polymerase II and thus stimulates the elongation phase of transcription. Chronic activation of CDK9 causes not only cardiac myocyte enlargement but also confers predisposition to heart failure. Due to the long interest of molecular oncologists and medicinal chemists in CDKs as potential targets of anticancer drugs, a portfolio of small-molecule inhibitors of CDK9 is available. Recent determination of CDK9's crystal structure now allows the development of selective inhibitors and their further optimization in terms of biochemical potency and selectivity. CDK9 may therefore constitute a novel target for drugs against cardiac hypertrophy.

A computational study of a recreated G protein-GEF reaction intermediate competent for nucleotide exchange: fate of the Mg ion

Small G-proteins of the superfamily Ras function as molecular switches, interacting with different cellular partners according to their activation state. G-protein activation involves the dissociation of bound GDP and its replacement by GTP, in an exchange reaction that is accelerated and regulated in the cell by guanine-nucleotide exchange factors (GEFs). Large conformational changes accompany the exchange reaction, and our understanding of the mechanism is correspondingly incomplete. However, much knowledge has been derived from structural studies of blocked or inactive mutant GEFs, which presumably closely represent intermediates in the exchange reaction and yet which are by design incompetent for carrying out the nucleotide exchange reaction. In this study we have used comparative modelling to recreate an exchange-competent form of a late, pre-GDP-ejection intermediate species in Arf1, a well-characterized small G-protein. We extensively characterized three distinct models of this intermediate using molecular dynamics simulations, allowing us to address ambiguities related to the mutant structural studies. We observed in particular the unfavorable nature of Mg associated forms of the complex and the establishment of closer Arf1-GEF contacts in its absence. The results of this study shed light on GEF-mediated activation of this small G protein and on predicting the fate of the Mg ion at a critical point in the exchange reaction. The structural models themselves furnish additional targets for interfacial inhibitor design, a promising direction for exploring potentially druggable targets with high biological specificity.

Protein conformational transitions: the closure mechanism of a kinase explored by atomistic simulations

Kinase large-scale conformational rearrangement is an issue of enormous biological and pharmacological relevance. Atomistic simulations able to capture the dynamics and the energetics of kinase large-scale motions are still in their infancy. Here, we present a computational study in which the atomistic dynamics of the "open-to-closed" movement of the cyclin-dependent kinase 5 (CDK5) have been simulated. Simulations were carried out using a new sampling method that is able to find the lowest free-energy channel between an initial state and a final state. This large-scale movement has a two-step mechanism: first, the alphaC-helix rotates by approximately 45 degrees , allowing the interaction between Glu51 and Arg149; then the CDK5 activation loop refolds to assume the closed conformation. We have also estimated the free-energy profile associated with the global motion and identified a CDK5 intermediate, which could be exploited for drug-design purposes. Our new sampling method turned out to be well-suited for investigating at an atomistic level the energetics and dynamics of kinase large-scale conformational motions.

Computational molecular biology approaches to ligand-target interactions

Binding of small molecules to their targets triggers complex pathways. Computational approaches are keys for predictions of the molecular events involved in such cascades. Here we review current efforts at characterizing the molecular determinants in the largest membrane-bound receptor family, the G-protein-coupled receptors (GPCRs). We focus on odorant receptors, which constitute more than half GPCRs. The work presented in this review uncovers structural and energetic aspects of components of the cellular cascade. Finally, a computational approach in the context of radioactive boron-based antitumoral therapies is briefly described.

Development of anaplastic lymphoma kinase (ALK) small-molecule inhibitors for cancer therapy

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase (RTK) involved in the genesis of several human cancers; indeed, ALK was initially identified in constitutively activated and oncogenic fusion forms--the most common being nucleophosmin (NPM)-ALK--in a non-Hodgkin's lymphoma (NHL) known as anaplastic large-cell lymphoma (ALCL) and subsequent studies identified ALK fusions in the human sarcomas called inflammatory myofibroblastic tumors (IMTs). In addition, two recent reports have suggested that the ALK fusion, TPM4-ALK, may be involved in the genesis of a subset of esophageal squamous cell carcinomas. While the cause-effect relationship between ALK fusions and malignancies such as ALCL and IMT is very well established, more circumstantial links implicate the involvement of the full-length, normal ALK receptor in the genesis of additional malignancies including glioblastoma, neuroblastoma, breast cancer, and others; in these instances, ALK is believed to foster tumorigenesis following activation by autocrine and/or paracrine growth loops involving the reported ALK ligands, pleiotrophin (PTN) and midkine (MK). There are no currently available ALK small-molecule inhibitors approved for clinical cancer therapy; however, recognition of the variety of malignancies in which ALK may play a causative role has recently begun to prompt developmental efforts in this area. This review provides a succinct summary of normal ALK biology, the confirmed and putative roles of ALK fusions and the full-length ALK receptor in the development of human cancers, and efforts to target ALK using small-molecule kinase inhibitors.

New alternative phosphorylation sites on the cyclin dependent kinase 1/cyclin a complex in p53-deficient human cells treated with etoposide: possible association with etoposide-induced apoptosis

Cell cycle arrest is a major cellular response to DNA damage preceding the decision to repair or die. Many malignant cells have non-functional p53 rendering them more "aggressive" in nature. Arrest in p53-negative cells occurs at the G2M cell cycle checkpoint. Failure of DNA damaged cells to arrest at G2 results in entry into mitosis and potential death through aberrant mitosis and/or apoptosis. The pivotal kinase regulating the G2M checkpoint is Cdk1/cyclin B whose activity is controlled by phosphorylation. The p53-negative myeloid leukemia cell lines K562 and HL-60 were used to determine Cdk1 phosphorylation status during etoposide treatment. Cdk1 tyrosine 15 phosphorylation was associated with G2M arrest, but not with cell death. Cdk1 tyrosine 15 phosphorylation also led to suppression of nuclear cyclin B-associated Cdk1 kinase activity. However cell death, associated with broader tyrosine phosphorylation of Cdk1 was not attributed to tyrosine 15 alone. This broader phosphoryl isoform of Cdk1 was associated with cyclin A and not cyclin B. Alternative phosphorylations sites were predicted as tyrosines 4, 99 and 237 by computer analysis. No similar pattern was found on Cdk2. These findings suggest novel Cdk1 phosphorylation sites, which appear to be associated with p53-independent cell death following etoposide treatment.

Structural model of the Plasmodium CDK, Pfmrk, a novel target for malaria therapeutics

Malaria, with 300-500 million clinical cases resulting in 1-3 million fatalities a year, is one of the most deadly tropical diseases. As current antimalarial therapeutics become increasingly ineffective due to parasitic resistance, there exists an urgent need to develop and pursue new therapeutic strategies. Recent genome sequencing and molecular cloning projects have identified several enzymes from Plasmodium (P.) falciparum that may represent novel drug targets, including a family of proteins that are homologous to the mammalian cyclin-dependent kinases (CDKs). CDKs are essential for the control of the mammalian cell cycle and, based on the conservation of the CDKs across species, the plasmodial CDKs are expected to play a crucial role in parasitic growth. Here we present a 3D structural model of Pfmrk, a putative human CDK activating kinase (CAK) homolog in P. falciparum. Notable features of the present structural model include: (1) parameterization of the Mg2+ hexacoordination system using ab initio quantum chemical calculations to accurately represent the ATP-kinase interaction; and (2) comparison between the docking scores and measured binding affinities for a series of oxindole-based Pfmrk inhibitors of known activity. Detailed analysis of inhibitor-Pfmrk binding interactions enabled us to identify specific residues (viz. Met66, Met75, Met91, Met94 and Phe143) within the Pfmrk binding pocket that may play an important role in inhibitor binding affinity and selectivity. The availability of this Pfmrk structural model, together with insights gained from analysis of ligand-receptor interactions, should promote the rational design of potent and selective Pfmrk inhibitors as antimalarial therapeutics.

Significance of water molecules in the inhibition of cylin-dependent kinase 2 and 5 complexes

Interest in CDK2 and CDK5 has stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. In the present paper, eight Molecular Dynamics (MD) simulations are carried out to examine the importance of structure and dynamics of water in the active site of both CDK2 and CDK5 complexes with roscovitine and indirubin analogues. Together with previous results, the current work shows a highly conserved water-involved hydrogen bonding (HB) network in both CDK2- and CDK5-indirubin combinations to complete information from the X-ray crystallography. The simulations suggest the importance of such a network for combining the inhibitor to the host protein as well as the significance of using an activated CDK as a template when designing new inhibitors. Different binding patterns of roscovitine in CDK2 and CDK5 are detected during the simulations because of the different binding conformations of the group on the C2 side chain, which might offer a clue toward finding highly selective inhibitors with regards to CDK2 and CDK5.

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.

Structural insights into the ATP binding pocket of the anaplastic lymphoma kinase by site-directed mutagenesis, inhibitor binding analysis, and homology modeling

Anaplastic lymphoma kinase (ALK) is a valid target for anticancer therapy; however, potent ALK inhibitors suitable for clinical use are lacking. Because the majority of described kinase inhibitors bind in the ATP pocket of the kinase domain, we have characterized this pocket in ALK using site-directed mutagenesis, inhibition studies, and molecular modeling. Mutation of the gatekeeper residue, a key structural determinant influencing inhibitor binding, rendered the fusion protein, NPM/ALK, sensitive to inhibition by SKI-606 in the nanomolar range, while PD173955 inhibited the NPM/ALK mutant at micromolar concentrations. In contrast, both wild type and mutant NPM/ALK were insensitive to imatinib. Computer modeling indicated that docking solutions obtained with a homology model representing the intermediate conformation of the ALK kinase domain reflected closely experimental data. The good agreement between experimental and virtual results indicate that the ALK molecular models described here are useful tools for the rational design of ALK selective inhibitors. In addition, 4-phenylamino-quinoline compounds may have potential as templates for ALK inhibitors.

Molecular dynamics simulations on the inhibition of cyclin-dependent kinases 2 and 5 in the presence of activators

Interests in CDK2 and CDK5 have stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. Molecular dynamics (MD) simulations have not only become a viable approach to drug design because of advances in computer technology but are increasingly an integral part of drug discovery processes. It is common in MD simulations of inhibitor/CDK complexes to exclude the activator of the CDKs in the structural models to keep computational time tractable. In this paper, we present simulation results of CDK2 and CDK5 with roscovitine using models with and without their activators (cyclinA and p25). While p25 was found to induce slight changes in CDK5, the calculations support that cyclinA leads to significant conformational changes near the active site of CDK2. This suggests that detailed and structure-based inhibitor design targeted at these CDKs should employ activator-included models of the kinases. Comparisons between P/CDK2/cyclinA/roscovitine and CDK5/p25/roscovitine complexes reveal differences in the conformations of the glutamine around the active sites, which may be exploited to find highly selective inhibitors with respect to CDK2 and CDK5.

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.

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.

Three-dimensional model of the human aromatase enzyme and density functional parameterization of the iron-containing protoporphyrin IX for a molecular dynamics study of heme-cysteinato cytochromes

Mammalian cytochromes P450 (CYP) are enzymes of great biological and pharmaco-toxicological relevance. Due to their membrane-bound nature, the structural characterization of these proteins is extremely difficult, and therefore computational techniques, such as comparative modeling, may help obtaining reliable structures of members of this family. An important feature of CYP is the presence of an iron-containing porphyrin group at the enzyme active site. This calls for quantum chemical calculations to derive charges and parameters suitable for classical force field-based investigations of this proteins family. In this report, we first carried out density functional theory (DFT) computations to derive suitable charges for the Fe2+-containing heme group of P450 enzymes. Then, by means of the homology modeling technique, and taking advantage of the recently published crystal structure of the human CYP2C9, we built a new model of the human aromatase (CYP19) enzyme. Furthermore, to study the thermal stability of the new model as well as to test the suitability of the new DFT-based heme parameters, molecular dynamics (MD) simulations were carried out on both CYP2C9 and CYP19. Finally, the last few ns of aromatase MD trajectories were investigated following the essential dynamics protocol that allowed the detection of some correlated motions among some protein domains.

Role of phosphorylated Thr160 for the activation of the CDK2/Cyclin A complex

The enzymatic activity of the CDK2/Cyclin A complex increases upon the specific phosphorylation of Thr160@CDK2. In the present study, we have performed a comparative molecular dynamics (MD) study of models of the complex CDK2/Cyclin A/Substrate, which differ for the presence or absence of the phosphate group bound to Thr160. The models are based on two X-ray structures available for CDK2/CyclinA and pCDK2/CyclinA/Substrate complexes. In this way, we analyze the influence of the phosphorylated Thr160 (pThr160) on both the flexibility of CDK2 activation loop (AL) and substrate binding in CDK2. Our calculations point to a decreased flexibility of the AL in the phosphorylated model, in fairly good agreement with experimental data, and to a key role of pThr160 for substrate recognition and stability. Multiple alignments of the CDKs sequences point to the very high conservation of the AL sequence among the CDKs, thus extending our results to all CDKs.

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.

Loop Flexibility and Solvent Dynamics as Determinants for the Selective Inhibition of Cyclin-Dependent Kinase 4: Comparative Molecular Dynamics Simulation Studies of CDK2 and CDK4

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.

Homology modeling and molecular dynamics study of GSK3/SHAGGY-like kinase

Although the GSK3/SHAGGY-like kinase is a highly conserved serine/threonine kinase implicated in many signaling pathways in eukaryotes, the lack of knowledge of its three-dimensional (3D) structure has hindered efforts to understand the binding specificities of substrate and catalytic mechanism. To understand the structure-activity relationships, the protein 3D structure was built by using homology modeling based on the known X-ray diffraction structure of Glycogen synthase kinase-3beta (Gsk3beta) and the model structure was further refined using unrestrained molecular dynamics simulations. The research indicates that the general 3D organization of the GSK3/SHAGGY-like kinase is a typical kinase family and comprises an N-terminal domain of beta-sheet and a larger C-terminal domain mainly constituted by alpha-helix. In order to understand the molecular interactions between the natural substrate-ATP and GSK3/SHAGGY-like kinase, a 3D model of the complex ATP-GSK3/SHAGGY-like kinase is developed by molecular docking program, which is helpful to guide the experimental realization and the new mutant designs as well. One important finding is that the identification of the key binding-site residue of Lys69 which plays an important role in the catalysis of GSK3/SHAGGY-like kinase and this is in consistent with experimental observation.

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.