Unveiling Conformational States of CDK6 Caused by Binding of Vcyclin Protein and Inhibitor by Combining Gaussian Accelerated Molecular Dynamics and Deep Learning.

CDK6 plays a key role in the regulation of the cell cycle and is considered a crucial target for cancer therapy. In this work, conformational transitions of CDK6 were identified by using Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy landscapes (FELs). DL finds that the binding pocket as well as the T-loop binding to the Vcyclin protein are involved in obvious differences of conformation contacts. This result suggests that the binding pocket of inhibitors (LQQ and AP9) and the binding interface of CDK6 to the Vcyclin protein play a key role in the function of CDK6. The analyses of FELs reveal that the binding pocket and the T-loop of CDK6 have disordered states. The results from principal component analysis (PCA) indicate that the binding of the Vcyclin protein affects the fluctuation behavior of the T-loop in CDK6. Our QM/MM-GBSA calculations suggest that the binding ability of LQQ to CDK6 is stronger than AP9 with or without the binding of the Vcyclin protein. Interaction networks of inhibitors with CDK6 were analyzed and the results reveal that LQQ contributes more hydrogen binding interactions (HBIs) and hot interaction spots with CDK6. In addition, the binding pocket endures flexibility changes from opening to closing states and the Vcyclin protein plays an important role in the stabilizing conformation of the T-loop. We anticipate that this work could provide useful information for further understanding the function of CDK6 and developing new promising inhibitors targeting CDK6.

Molecular docking and MD simulation studies of 4-thiazol-N-(pyridin-2-yl)pyrimidin-2-amine derivatives as novel inhibitors targeted to CDK2/4/6.

PURPOSE: Nowadays, cyclin-dependent kinase 4/6 (CDK4/6) inhibitors have been approved for treating metastatic breast cancer and have achieved inspiring curative effects. But some discoveries have indicated that CDK 4/6 are not the requisite factors in some cell types because CDK2 partly compensates for the inhibition of CDK4/6. Thus, it is urgent to design CDK2/4/6 inhibitors for significantly enhancing their potency. This study aims to explore the mechanism of the binding of CDK2/4/6 kinases and their inhibitors to design novel CDK2/4/6 inhibitors for significantly enhancing their potency in different kinds of cancers. MATERIALS AND METHODS: A series of 72 disparately functionalized 4-substituted N-phenylpyrimidin-2-amine derivatives exhibiting potent inhibitor activities against CDK2, CDK4 and CDK6 were collected to apply to this research. The total set of these derivatives was divided into a training set (54 compounds) and a test set (18 compounds). The derivatives were constructed through the sketch molecule module in SYBYL 6.9 software. A Powell gradient algorithm and Tripos force field were used to calculate the minimal structural energy and the minimized structure was used as the initial conformation for molecular docking. By the means of 3D-QSAR models, partial least squares (PLS) analysis, molecular dynamics (MD) simulations and binding free energy calculations, we can find the relationship between structure and biological activity. RESULTS: In this study, we used molecular docking, 3D-QSAR and molecular dynamics simulation methods to comprehensively analyze the interaction and structure-activity relationships of 72 new CDK2/4/6 inhibitors. We used detailed statistical data to reasonably verify the constructed 3D-QSAR models for three receptors (q2 of CDK2 = 0.714, R2pred = 0.764, q2 = 0.815; R2pred of CDK4 = 0.681, q2 = 0.757; R2pred of CDK6 = 0.674). MD simulations and decomposition energy analysis validated the reasonability of the docking results and identified polar interactions as crucial factors that influence the different bioactivities of the studied inhibitors of CDK2/4/6 receptors, especially the electrostatic interactions of Lys33/35/43 and Asp145/158/163. The nonpolar interaction with Ile10/12/19 was also critical for the differing potencies of the CDK2/4/6 inhibitors. We concluded that the following probably enhanced the bioactivity against CDK2/4/6 kinases: (1) electronegative groups at the N1-position and electropositive and moderate-sized groups at ring E; (2) electrogroups featured at R2; (3) carbon atoms at the X-position or ring C replaced by a benzene ring; and (4) an electrogroup as R4. CONCLUSION: Previous studies, to our knowledge, only utilized a single approach of 3D-QSAR and did not integrate this method with other sophisticated techniques such as molecular dynamics simulations to discover new potential inhibitors of CDK2, CDK4, or CDK6. So we applied the intergenerational technology, such as 3D-QSAR technology, molecular docking simulation techniques, molecular dynamics simulations and MMPBSA19/MMGBSA20-binding free energy calculations to statistically explore the correlations between the structure with biological activities. The constructed 3D-QSAR models of the three receptors were reasonable and confirmed by the excellent statistical data. We hope the results obtained from this work will provide some useful references for the development of novel CDK2/4/6 inhibitors.

Palbociclib and Michael-acceptor hybrid compounds as CDK4/6 covalent inhibitors: Improved potency, broad anticancer spectrum and overcoming drug resistance.

To search for potent CDK4/6 covalent inhibitors, total 14 compounds have been designed and synthesized by connecting different Michael-acceptor to the piperazine moiety of palbociclib. All the compounds displayed good antiproliferative activity against human hepatoma cell (HepG2), non-small cell lung cancer (A549), and breast cancer (MDA-MB-231 and MCF-7) cell lines. In particular, compound A4 showed the highest inhibitory activity to MDA-MB-231 and MCF-7 cells with IC50 values of 0.51 µM and 0.48 µM, respectively. More importantly, A4 also showed strong inhibition against MDA-MB-231/palbociclib cells, indicating that A4 could effectively avoid the resistance of palbociclib. In the enzyme test, A4 showed selective inhibitory activity against CDK4/6, with the IC50 value of 18 nM and 13 nM, respectively. It was also found that A4 could efficiently induce apoptosis and arrest the cell cycle at G0/G1 phase. Moreover, A4 could significantly decrease the phosphorylation level of CDK4 and CDK6. HPLC and molecular modeling studies suggested that A4 could form a covalent bond with the target protein.

Naringenin as a potential inhibitor of human cyclin-dependent kinase 6: Molecular and structural insights into anti-cancer therapeutics.

Cancer is one of the major causes of global deaths and needs immediate therapeutic development. So far, several strategies have been undertaken to prevent cancer, including kinase targeting by small-molecule inhibitors. Cyclin dependent kinase 6 (CDK6) plays an essential role in cancer progression and development as its overexpression is associated with tumor development and progression. The present study demonstrated that Naringenin (NAG) binds strongly to CDK6 with a binding affinity of -7.51 kcal/mol. ATPase assay of CDK6 in the presence of NAG shows that it inhibits CDK6 with an IC50 = 3.13 µM. Fluorescence and isothermal titration calorimetry studies demonstrated that NAG binds to CDK6 with the binding constant (K) values of 3.55 × 106 M-1 and 7.06 ± 2.70 × 106 M-1, respectively. The cell-based functional studies showed that NAG decreases the cell viability of human cancer cell lines, induces apoptosis, and reduces their colonization ability. Outcomes of the present in silico and in vitro studies highlighted the significance of NAG for the development of anti-cancer leads in terms of CDK6 inhibitors and provided future implications for combinatorial anti-cancer therapies.

Effect of pH on the structure and function of cyclin-dependent kinase 6.

Cyclin-dependent kinase 6 (CDK6) is an important protein kinase that regulates cell growth, development, cell metabolism, inflammation, and apoptosis. Its overexpression is associated with reprogramming glucose metabolism through alternative pathways and apoptosis, which ultimately plays a significant role in cancer development. In the present study, we have investigated the structural and conformational changes in CDK6 at varying pH employing a multi-spectroscopic approach. Circular dichroism (CD) spectroscopy revealed at extremely acidic conditions (pH 2.0-4.0), the secondary structure of CDK6 got significantly disrupted, leading to aggregates formation. These aggregates were further characterized by employing Thioflavin T (ThT) fluorescence. No significant secondary structural changes were observed over the alkaline pH range (pH 7.0-11.0). Further, fluorescence and UV spectroscopy revealed that the tertiary structure of CDK6 was disrupted under extremely acidic conditions, with slight alteration occurring in mild acidic conditions. The tertiary structure remains intact over the entire alkaline range. Additionally, enzyme assay provided an insight into the functional aspect of CDK at varying pH; CDK6 activity was optimal in the pH range of 7.0-8.0. This study will provide a platform that provides newer insights into the pH-dependent dynamics and functional behavior of CDK6 in different CDK6 directed diseased conditions, viz. different types of cancers where changes in pH contribute to cancer development.

A Search for Cyclin-Dependent Kinase 4/6 Inhibitors by Pharmacophore-Based Virtual Screening, Molecular Docking, and Molecular Dynamic Simulations.

The G1 phase of cell cycle progression is regulated by Cyclin-Dependent Kinase 4 (CDK4) as well as Cyclin-Dependent Kinase 6 (CDK6), and the acivities of these enzymes are regulated by the catalytic subunit, cyclin D. Cell cycle control through selective pharmacological inhibition of CDK4/6 has proven to be beneficial in the treatment of estrogen receptor-positive (ER-positive) breast cancer, particularly improving the progression-free survival of patients. Thus, targeting specific inhibition on CDK4/6 is bound to increase therapeutic efficiency. This study aimed to obtain CDK4/6 inhibitors through a pharmacophore-based virtual screening of the ZINC15 purchasable compound database using the in silico method. The pharmacophore model was designed based on the FDA-approved cdk4/6 inhibitor structures, and molecular docking was performed to further screen the hit compounds obtained. A total of eight compounds were selected based on docking results and interactions with CDK4 and CDK6, using palbociclib as the reference drug. According to the results, the compounds of ZINC585292724 and ZINC585291674 were the best compounds based on free binding energy, as well as hydrogen bond stability, and, therefore, exhibit potential as starting points in the development of CDK4/6 inhibitors.

Targeting cyclin-dependent kinase 6 by vanillin inhibits proliferation of breast and lung cancer cells: Combined computational and biochemical studies.

Cyclin-dependent kinase 6 (CDK6) is a member of serine/threonine kinase family, and its overexpression is associated with cancer development. Thus, it is considered as a potential drug target for anticancer therapies. This study showed the CDK6 inhibitory potential of vanillin using combined experimental and computational methods. Structure-based docking and 200 ns molecular dynamics simulation studies revealed that the binding of vanillin stabilizes the CDK6 structure and provides mechanistic insights into the binding mechanism. Enzyme inhibition and fluorescence-binding studies showed that vanillin inhibits CDK6 with an half maximal inhibitory concentration = 4.99 µM and a binding constant (K) 4.1 × 107 M-1 . Isothermal titration calorimetry measurements further complemented our observations. Studies on human cancer cell lines (MCF-7 and A549) showed that vanillin decreases cell viability and colonization properties. The protein expression studies have further revealed that vanillin reduces the CDK6 expression and induces apoptosis in the cancer cells. In conclusion, our study presents the CDK6-mediated therapeutic implications of vanillin for anticancer therapies.

Discovery of new small-molecule cyclin-dependent kinase 6 inhibitors through computational approaches.

Excessive cell proliferation due to cell cycle disorders is one of the hallmarks of breast cancer. Cyclin-dependent kinases (CDKs), which are involved in the transition of the cell cycle from G1 phase to S phase by combining CDKs with cyclin, are considered promising targets with broad therapeutic potential based on their critical role in cell cycle regulation. Pharmacological evidence has shown that abnormal cell cycle due to the overexpression of CDK6 is responsible for the hyperproliferation of cancer cells. Blocking CDK6 expression inhibits tumour survival and growth. Therefore, CDK6 can be regarded as a potential target for anticancer therapeutics. Thus, small molecules that can be considered CDK inhibitors have been developed into promising anticancer drugs. In this study, combined structure-based and ligand-based in silicon models were created to identify new chemical entities against CDK6 with the appropriate pharmacokinetic properties. The database used to screen drug-like compounds in this thesis was based on the best E-pharmacophore hypothesis and the best ligand-based drug hypothesis. As a result, 147 common compounds were identified by further molecular docking. Surprisingly, the in vitro evaluation results of 20 of those compounds showed that the two had good CDK6 inhibitory effects. The best compound was subjected to kinase panel screening, followed by molecular dynamic simulations. The 50-ns MD studies revealed the pivotal role of VAL101 in the binding of inhibitors to CDK6. Overall, the identification of two new chemical entities with CDK6 inhibitory activity demonstrated the feasibility and potential of the new method.

A bioinformatics investigation into molecular mechanism of Yinzhihuang granules for treating hepatitis B by network pharmacology and molecular docking verification.

Yinzhihuang granules (YZHG) is a patented Chinese medicine for the treatment of hepatitis B. This study aimed to investigate the intrinsic mechanisms of YZHG in the treatment of hepatitis B and to provide new evidence and insights for its clinical application. The chemical compounds of YZHG were searched in the CNKI and PUBMED databases, and their putative targets were then predicted through a search of the SuperPred and Swiss Target Prediction databases. In addition, the targets of hepatitis B were obtained from TTD, PharmGKB and DisGeNET. The abovementioned data were visualized using Cytoscape 3.7.1, and network construction identified a total of 13 potential targets of YZHG in the treatment of hepatitis B. Molecular docking verification showed that CDK6, CDK2, TP53 and BRCA1 might be strongly correlated with hepatitis B treatment. Furthermore, GO and KEGG analyses indicated that the treatment of hepatitis B by YZHG might be related to positive regulation of transcription, positive regulation of gene expression, the hepatitis B pathway and the viral carcinogenesis pathway. Network pharmacology intuitively shows the multicomponent, multitarget and multichannel pharmacological effects of YZHG in the treatment of hepatitis B and provides a scientific basis for its mechanism of action.

Ellagic Acid Controls Cell Proliferation and Induces Apoptosis in Breast Cancer Cells via Inhibition of Cyclin-Dependent Kinase 6.

Cyclin-Dependent Kinase 6 (CDK6) plays an important role in cancer progression, and thus, it is considered as an attractive drug target in anticancer therapeutics. This study presents an evaluation of dietary phytochemicals, capsaicin, tocopherol, rosmarinic acid, ursolic acid, ellagic acid (EA), limonene, caffeic acid, and ferulic acid for their potential to inhibit the activity of CDK6. Molecular docking and fluorescence binding studies revealed appreciable binding affinities of these compounds to the CDK6. Among them, EA shows the highest binding affinity for CDK6, and thus a molecular dynamics simulation study of 200 ns was performed to get deeper insights into the binding mechanism and stability of the CDK6-EA complex. Fluorescence binding studies revealed that EA binds to the CDK6 with a binding constant of K = 107 M-1 and subsequently inhibits its enzyme activity with an IC50 value of 3.053 µM. Analysis of thermodynamic parameters of CDK6-EA complex formation suggested a hydrophobic interaction driven process. The treatment of EA decreases the colonization of cancer cells and induces apoptosis. Moreover, the expression of CDK6 has been downregulated in EA-treated human breast cancer cell lines. In conclusion, this study establishes EA as a potent CDK6 inhibitor that can be further evaluated in CDK6 directed anticancer therapies.

Discovery and SARs of 5-Chloro-N4-phenyl-N2-(pyridin-2-yl)pyrimidine-2,4-diamine Derivatives as Oral Available and Dual CDK 6 and 9 Inhibitors with Potent Antitumor Activity.

Cyclin-dependent kinases (CDKs) are promising therapeutic targets for cancer therapy. Herein, we describe our efforts toward the discovery of a series of 5-chloro-N4-phenyl-N2-(pyridin-2-yl)pyrimidine-2,4-diamine derivatives as dual CDK6 and 9 inhibitors. Intensive structural modifications lead to the identification of compound 66 as the most active dual CDK6/9 inhibitor with balancing potency against these two targets and good selectivity over CDK2. Further biological studies revealed that compound 66 was directly bound to CDK6/9, resulting in suppression of their downstream signaling pathway and inhibition of cell proliferation by blocking cell cycle progression and inducing cellular apoptosis. More importantly, compound 66 significantly inhibited tumor growth in a xenograft mouse model with no obvious toxicity, indicating the promising therapeutic potential of CDK6/9 dual inhibitors for cancer treatment. Therefore, the above results are of great importance in the development of dual CDK6/9 inhibitors for cancer therapy.

Potent and Preferential Degradation of CDK6 via Proteolysis Targeting Chimera Degraders.

A focused PROTAC library hijacking cancer therapeutic target CDK6 was developed. A design principle as "match/mismatch" was proposed for understanding the degradation profile differences in these PROTACs. Notably, potent PROTACs with specific and remarkable CDK6 degradation potential were generated by linking CDK6 inhibitor palbociclib and E3 ligase CRBN recruiter pomalidomide. The PROTAC strongly inhibited proliferation of hematopoietic cancer cells including multiple myeloma and robustly degraded copy-amplified/mutated forms of CDK6, indicating future potential clinical applications.

Artificial Intelligence Approach to Find Lead Compounds for Treating Tumors.

It has been demonstrated that MMP13 enzyme is related to most cancer cell tumors. The world's largest traditional Chinese medicine database was applied to screen for structure-based drug design and ligand-based drug design. To predict drug activity, machine learning models (Random Forest (RF), AdaBoost Regressor (ABR), Gradient Boosting Regressor (GBR)), and Deep Learning models were utilized to validate the Docking results, and we obtained an R2 of 0.922 on the training set and 0.804 on the test set in the RF algorithm. For the Deep Learning algorithm, R2 of the training set is 0.90, and R2 of the test set is 0.810. However, these TCM compounds fly away during the molecular dynamics (MD) simulation. We seek another method: peptide design. All peptide database were screened by the Docking process. Modification peptides were optimized the interaction modes, and the affinities were assessed with ZDOCK protocol and Refine Docked protein protocol. The 300 ns MD simulation evaluated the stability of receptor-peptide complexes. The double-site effect appeared on S2, a designed peptide based on a known inhibitor, when complexed with BCL2. S3, a designed peptide referred from endogenous inhibitor P16, competed against cyclin when binding with CDK6. The MDM2 inhibitors S5 and S6 were derived from the P53 structure and stable binding with MDM2. A flexible region of peptides S5 and S6 may enhance the binding ability by changing its own conformation, which was unforeseen. These peptides (S2, S3, S5, and S6) are potentially interesting to treat cancer; however, these findings need to be affirmed by biological testing, which will be conducted in the near future.

Cyclin D-Cdk4,6 Drives Cell-Cycle Progression via the Retinoblastoma Protein's C-Terminal Helix.

The cyclin-dependent kinases Cdk4 and Cdk6 form complexes with D-type cyclins to drive cell proliferation. A well-known target of cyclin D-Cdk4,6 is the retinoblastoma protein Rb, which inhibits cell-cycle progression until its inactivation by phosphorylation. However, the role of Rb phosphorylation by cyclin D-Cdk4,6 in cell-cycle progression is unclear because Rb can be phosphorylated by other cyclin-Cdks, and cyclin D-Cdk4,6 has other targets involved in cell division. Here, we show that cyclin D-Cdk4,6 docks one side of an alpha-helix in the Rb C terminus, which is not recognized by cyclins E, A, and B. This helix-based docking mechanism is shared by the p107 and p130 Rb-family members across metazoans. Mutation of the Rb C-terminal helix prevents its phosphorylation, promotes G1 arrest, and enhances Rb's tumor suppressive function. Our work conclusively demonstrates that the cyclin D-Rb interaction drives cell division and expands the diversity of known cyclin-based protein docking mechanisms.

Direct CDKN2 Modulation of CDK4 Alters Target Engagement of CDK4 Inhibitor Drugs.

The interaction of a drug with its target is critical to achieve drug efficacy. In cases where cellular environment influences target engagement, differences between individuals and cell types present a challenge for a priori prediction of drug efficacy. As such, characterization of environments conducive to achieving the desired pharmacologic outcome is warranted. We recently reported that the clinical CDK4/6 inhibitor palbociclib displays cell type-specific target engagement: Palbociclib engaged CDK4 in cells biologically sensitive to the drug, but not in biologically insensitive cells. Here, we report a molecular explanation for this phenomenon. Palbociclib target engagement is determined by the interaction of CDK4 with CDKN2A, a physiologically relevant protein inhibitor of CDK4. Because both the drug and CDKN2A prevent CDK4 kinase activity, discrimination between these modes of inhibition is not possible by traditional kinase assays. Here, we describe a chemo-proteomics approach that demonstrates high CDK4 target engagement by palbociclib in cells without functional CDKN2A and attenuated target engagement when CDKN2A (or related CDKN2/INK4 family proteins) is abundant. Analysis of biological sensitivity in engineered isogenic cells with low or absent CDKN2A and of a panel of previously characterized cell lines indicates that high levels of CDKN2A predict insensitivity to palbociclib, whereas low levels do not correlate with sensitivity. Therefore, high CDKN2A may provide a useful biomarker to exclude patients from CDK4/6 inhibitor therapy. This work exemplifies modulation of kinase target engagement by endogenous proteinaceous regulators and highlights the importance of cellular context in predicting inhibitor efficacy.

Homolog-Selective Degradation as a Strategy to Probe the Function of CDK6 in AML.

The design of selective small molecules is often stymied by similar ligand binding pockets. Here, we report BSJ-03-123, a phthalimide-based degrader that exploits protein-interface determinants to achieve proteome-wide selectivity for the degradation of cyclin-dependent kinase 6 (CDK6). Pharmacologic CDK6 degradation targets a selective dependency of acute myeloid leukemia cells, and transcriptomics and phosphoproteomics profiling of acute degradation of CDK6 enabled dynamic mapping of its immediate role in coordinating signaling and transcription.

2-Methylpyridine-1-ium-1-sulfonate as an Inducer of Apoptosis and Cell Cycle Arrest: A comparative in vitro and Computational Study.

"Let food be thy medicine and thy medicine be thy food" was expressed by Hippocrates and the health benefits of medicinal plants and natural products have been considered by humans since historic times. The current study aims to investigate the anti-cancer activity of 2-Methylpyridine-1-ium-1-sulfonate (MPS) isolated from bulbs of Allium hirtifolium. The MPS compound (in a dose-dependent manner) induced arrest the AGS cells in G1 and G2/M phases, and Caco-2 cells in G1 and S phases. These findings were associated with the down-regulation of cyclin D1, CDK4, and up-regulation of p21, p27 and p53. According to the morphological observations and DNA fragmentation assay, the MPS compound induced apoptosis in both cell lines, and also cause a significant increase in the expression of Bax/Bcl-2. In this context, our molecular docking results unveiled that the MPS compound has considerable affinity to interact with the minor groove of ctDNA and also with cell cycle kinases. To approve and find the accurate MPS mode of action against cancer cell lines (especially in gastrointestinal cancer) further studies is highly recommended.

Synthesis, cytotoxic evaluation, and molecular docking studies of the semi-synthetic "triterpenoid-steroid" hybrids.

Synthetic transformations of steroids for drug discovery and improvement of drug effectiveness have been an important part of modern medicinal chemistry and pharmaceutical sciences. Pentacyclic triterpenoids, being represented in the nature by various structures and biogenetically related to steroids, can largely expand the spectrum of biologically active steroidal agents via synthesis of the so-called "triterpenoid-steroid" hybrids. In the presented work, the nitrile anion cyclizations of 3,4-secolupane and 3,4-seco-oleanane nitriles and follow-up synthetic transformations of the cyclized products with formation of the gemm-dimethyl-free A ring "triterpenoid-steroid" hybrids were studied. Furthermore, the resulting cyclic compounds were modified at C3, C4, and/or C5 positions of ring A, as well as at C20, C28, and C30 positions of the isopropylidene moiety in the case of lupane triterpenoids. The cytotoxic effect of the synthesized compounds against seven cancer cell lines HEp-2, HCT 116, MS, RD TE32, A549, MCF7, and PC3 was evaluated. The in silico identification of potential anticancer protein targets with regard to the compounds, which were active at micromolar concentrations against tested cell lines, was carried out. The molecular docking studies showed that compound 19, which demonstrated most pronounced cytotoxicity (IC50 0.64-3.17 µM) against all tested cell lines, fits well the active sites of CDK6 and HER2/neu.

Reduced state transition barrier of CDK6 from open to closed state induced by Thr177 phosphorylation and its implication in binding modes of inhibitors.

BACKGROUND: CDK6 is considered as a highly validated anticancer drug target due to its essential role in regulating cell cycle progression at G1 restriction point. Activation of CDK6 requires the phosphorylation of Thr177 on A-loop, but the structural insights of the activation mechanism remain unclear. METHODS: Herein, all-atoms molecular dynamics (MD) simulations were used to study the effects of Thr177 phosphorylation on the dynamic structure of CDK6-Vcyclin complex. RESULTS: MD results indicated that the free energy barrier of the transition from open to closed state decreased ~47.2% after Thr177 phosphorylation. Key steps along the state transition process were obtained from a cluster analysis. Binding preference of ten different inhibitors to open or closed state were also investigated through molecular docking along with MD simulations methods. CONCLUSIONS: Our results indicated that Thr177 phosphorylation increased the flexibility around the ATP-binding pocket. The transition of the ATP-binding pocket between open and closed states should be considered for understanding the binding of CDK6 inhibitors. GENERAL SIGNIFICANCE: This work could deepen the understanding of CDKs activation mechanism, and provide useful information for the discovery of new CDKs inhibitors with high affinity and specificity.

Palbociclib can overcome mutations in cyclin dependent kinase 6 that break hydrogen bonds between the drug and the protein.

Inhibition of cyclin dependent kinases (CDKs) 4 and 6 prevent cells from entering the synthesis phase of the cell cycle. CDK4 and 6 are therefore important drug targets in various cancers. The selective CDK4/6 inhibitor palbociclib is approved for the treatment of breast cancer and has shown activity in a cellular model of mixed lineage leukaemia (MLL)-rearranged acute myeloid leukaemia (AML). We studied the interactions of palbociclib and CDK6 using molecular dynamics simulations. Analysis of the simulations suggested several interactions that stabilized the drug in its binding site and that were not observed in the crystal structure of the protein-drug complex. These included a hydrogen bond to His 100 that was hitherto not reported and several hydrophobic contacts. Evolutionary-based bioinformatic analysis was used to suggest two mutants, D163G and H100L that would potentially yield drug resistance, as they lead to loss of important protein-drug interactions without hindering the viability of the protein. One of the mutants involved a change in the glycine of the well-conserved DFG motif of the kinase. Interestingly, CDK6-dependent human AML cells stably expressing either mutant retained sensitivity to palbociclib, indicating that the protein-drug interactions are not affected by these. Furthermore, the cells were proliferative in the absence of palbociclib, indicating that the Asp to Gly mutation in the DFG motif did not interfere with the catalytic activity of the protein.