Molecular Dynamic Behavior and Binding Affinity of Flavonoid Analogues to the Cyclin Dependent Kinase 6/cyclin D Complex

The cyclin dependent kinases (CDKs), each with their respective regulatory partner cyclin that are involved in the regulation of the cell cycle, apoptosis, and transcription, are potentially interesting targets for cancer therapy. The CDK6 complex with cyclin D (CDK6/cycD) drives cellular proliferation by phosphorylation of specific key target proteins. To understand the flavonoids that inhibit the CDK6/cycD functions, molecular dynamics simulations (MDSs) were performed on three inhibitors, fisetin (FST), apigenin (AGN), and chrysin (CHS), complexed with CDK6/cycD, including the two different binding orientations of CHS: FST-like (CHS_A) and deschloro-flavopiridol-like (CHS_B). For all three inhibitors, including both CHS orientations, the conserved interaction between the 4-keto group of the flavonoid and the backbone V101 nitrogen of CDK6 was strongly detected. The 3'- and 4'-OH groups on the flavonoid phenyl ring and the 3-OH group on the benzopyranone ring of inhibitor were found to significantly increase the binding and inhibitory efficiency. Besides the electrostatic interactions, especially through hydrogen bond formation, the van der Waals (vdW) interactions with the I19, V27, F98, H100, and L152 residues of CDK6 are also important factors in the binding efficiency of flavonoids against the CDK6/cycD complex. On the basis of the docking calculation and MM-PBSA method, the order of the predicted inhibitory affinities of these three inhibitors toward the CDK6/cycD was FST > AGN > CHS, which is in good agreement with the experimental data. In addition, CHS preferentially binds to the active CDK6 in a different orientation to FST and AGN but similar to its related analog, deschloro-flavopiridol. The obtained results are useful as the basic information for the further design of potent anticancer drugs specifically targeting the CDK6 enzyme.

[1]  Jian Shen,et al.  Calculation of binding energy differences for receptor–ligand systems using the Poisson‐Boltzmann method , 1995, J. Comput. Chem..

[2]  K. Vermeulen,et al.  The cell cycle: a review of regulation, deregulation and therapeutic targets in cancer , 2003, Cell proliferation.

[3]  Thanyada Rungrotmongkol,et al.  Source of oseltamivir resistance in avian influenza H5N1 virus with the H274Y mutation , 2009, Amino Acids.

[4]  J. Harper,et al.  Cyclin-dependent kinases. , 2001, Chemical reviews.

[5]  D O Morgan,et al.  Cyclin-dependent kinases: engines, clocks, and microprocessors. , 1997, Annual review of cell and developmental biology.

[6]  Nopporn Kaiyawet,et al.  Molecular insight into the specific binding of ADP-ribose to the nsP3 macro domains of chikungunya and Venezuelan equine encephalitis viruses: molecular dynamics simulations and free energy calculations. , 2010, Journal of molecular graphics & modelling.

[7]  Prachi Mishra,et al.  DNA hypermethylation as a chemotherapy target. , 2011, Cellular signalling.

[8]  Thanyada Rungrotmongkol,et al.  Understanding of known drug‐target interactions in the catalytic pocket of neuraminidase subtype N1 , 2008, Proteins.

[9]  S H Kim,et al.  Structural basis for specificity and potency of a flavonoid inhibitor of human CDK2, a cell cycle kinase. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Supot Hannongbua,et al.  On the lower susceptibility of oseltamivir to influenza neuraminidase subtype N1 than those in N2 and N9. , 2007, Biophysical journal.

[11]  Z Topcu,et al.  DNA topoisomerases as targets for anticancer drugs , 2001, Journal of clinical pharmacy and therapeutics.

[12]  Shunichi Takeda,et al.  Cyclin-dependent kinases and cell-cycle transitions: does one fit all? , 2008, Nature Reviews Molecular Cell Biology.

[13]  S. Kornbluth,et al.  Regulatory roles of cyclin dependent kinase phosphorylation in cell cycle control. , 1996, Current opinion in cell biology.

[14]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[15]  T. Rungrotmongkol,et al.  Source of high pathogenicity of an avian influenza virus H5N1: why H5 is better cleaved by furin. , 2008, Biophysical journal.

[16]  Laurent Meijer,et al.  Crystal structure of a human cyclin-dependent kinase 6 complex with a flavonol inhibitor, fisetin. , 2005, Journal of medicinal chemistry.

[17]  W. Dunphy,et al.  The decision to enter mitosis. , 1994, Trends in cell biology.

[18]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[19]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[20]  E. Nigg,et al.  Targets of cyclin-dependent protein kinases. , 1993, Current opinion in cell biology.

[21]  M. Malumbres,et al.  Targeting cell cycle kinases for cancer therapy. , 2007, Current medicinal chemistry.

[22]  A. Jemal,et al.  Global Cancer Statistics , 2011 .

[23]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[24]  M. Dorée,et al.  The cyclin‐dependent protein kinases and the control of cell division , 1994, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[25]  Thanyada Rungrotmongkol,et al.  How do carbon nanotubes serve as carriers for gemcitabine transport in a drug delivery system? , 2011, Journal of molecular graphics & modelling.

[26]  E. Karsenti,et al.  The role of cdc25 in checkpoints and feedback controls in the eukaryotic cell cycle , 1994, Journal of Cell Science.

[27]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[28]  Kornelia Polyak,et al.  Mechanism of CDK activation revealed by the structure of a cyclinA-CDK2 complex , 1995, Nature.

[29]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[30]  Supot Hannongbua,et al.  Accurate prediction of protonation state as a prerequisite for reliable MM‐PB(GB)SA binding free energy calculations of HIV‐1 protease inhibitors , 2008, J. Comput. Chem..

[31]  Thanyada Rungrotmongkol,et al.  Dynamic Behavior of Avian Influenza A Virus Neuraminidase Subtype H5N1 in Complex with Oseltamivir, Zanamivir, Peramivir, and Their Phosphonate Analogues , 2009, J. Chem. Inf. Model..