Structure-based pharmacophore modeling, virtual screening and simulation studies for the identification of potent anticancerous phytochemical lead targeting cyclin-dependent kinase 2

J Biomol Struct Dyn. 2022;40(20):9815-9832. doi: 10.1080/07391102.2021.1936178. Epub 2021 Jun 21.

Abstract

Cyclin-dependent kinases are of critical importance in directing various cell cycle phases making them as potential tumor targets. Cyclin-dependent kinase 2 (CDK2) in particular plays a significant part during cell cycle events and its imbalance roots out tumorogenic environment. Herein, we built a structure-based pharmacophore model complementing the ATP pocket site of CDK2 with four pharmacophoric features, using a series of structures obtained from cluster analysis during MD simulation assessment. This was followed by its validation and further database screening against Taiwan indigenous plants database (5284 compounds). The screened compounds were subjected toward Lipinski's rule (RO5) and ADMET filter followed by docking analysis and simulation study. In filtering hits (10 compounds) via molecular docking against CDK2, Schinilenol with -8.1 kcal/mol fetched out as a best lead phytoinhibitor in the presence of standard drug (Dinaciclib). Additionally, pharmacophore mapping analysis also indicated relative fit values of dinaciclib and schinilenol as 2.37 and 2.31, respectively. Optimization, flexibility prediction and the stability of CDK2 in complex with the ligands were also ascertained by means of molecular dynamics for 50 ns, which further proposed schinilenol having better binding stability than dinaciclib with RMSD values ranging from 0.31 to 0.34 nm. Reactivity site, biological activity detection and cardiotoxicity assessment also proposed schinilenol as a better phytolead inhibitor than the existing dinaciclib. Abbreviations: CDK2: Cyclin dependent kinase2; ATP: Adenosine triphosphate; MD: Molecular dynamics, RO5: Rule of five; ADMET: Absorption, distribution, metabolism, and excretion; RMSD: Root mean square deviation; DS: Discovery Studio; SOM: Site of metabolism; RBPM: receptor based pharmacophore model; TIP: Schinilenol; hERG: human Ether-à-go-go - Related GeneCommunicated by Ramaswamy H. Sarma.

Keywords: CDK2; MD simulation; anticancer phytoinhibitor; pharmacophore modeling; schinilenol; screening.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Adenosine Triphosphate
  • Cyclin-Dependent Kinase 2* / antagonists & inhibitors
  • Humans
  • Hydrogen Bonding
  • Ligands
  • Molecular Docking Simulation
  • Molecular Dynamics Simulation
  • Pharmacophore*
  • Phytochemicals / pharmacology
  • Protein Kinase Inhibitors* / chemistry
  • Quantitative Structure-Activity Relationship

Substances

  • Adenosine Triphosphate
  • Cyclin-Dependent Kinase 2
  • Ligands
  • Phytochemicals
  • Protein Kinase Inhibitors