Next Article in Journal
Identification of the Biosynthetic Gene Cluster of New Piperazic Acid-Containing Lipopeptides with Cytotoxic Activity in the Genome of Marine Streptomyces PHM034
Next Article in Special Issue
Curcumin Mitigates Malathion-Induced Renal Injury: Suppression of Apoptosis and Modulation of NF-κβ/TNF-α and Nrf2, and HO-1 Signaling
Previous Article in Journal
Brain Regulation of Cardiac Function during Hypoglycemia
Previous Article in Special Issue
Ultra-Performance Liquid Chromatography Coupled with Mass Metabolic Profiling of Ammi majus Roots as Waste Product with Isolation and Assessment of Oral Mucosal Toxicity of Its Psoralen Component Xanthotoxin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Marine-Derived Compounds for CDK5 Inhibition in Cancer: Integrating Multi-Stage Virtual Screening, MM/GBSA Analysis and Molecular Dynamics Investigations

by
Tagyedeen H. Shoaib
1,†,
Mohammed A. Almogaddam
1,†,
Yusra Saleh Andijani
2,
Samaher Ahmad Saib
3,
Najwa Mahmoud Almaghrabi
4,
Abdulaziz Fahad Elyas
5,
Rahmah Yasin Azzouni
6,
Ehda Ahmad Awad
7,
Shaimaa G. A. Mohamed
8,
Gamal A. Mohamed
9,
Sabrin R. M. Ibrahim
10,11,*,
Hazem G. A. Hussein
12,
Wadah Osman
13,14,
Ahmed Ashour
13,15,
Asmaa E. Sherif
13,15 and
Abdulrahim A. Alzain
1,*,†
1
Department of Pharmaceutical Chemistry, Faculty of Pharmacy, University of Gezira, Wad Madani 21111, Sudan
2
Department of Pharmacology and Toxicology, College of Pharmacy, Taibah University, Al-Madinah Al-Munawwarah 30078, Saudi Arabia
3
College of Pharmacy, Taibah University, Medina 42353, Saudi Arabia
4
Pharmaceutical Care, King Abdullah Medical Complex, MOH, Jeddah 23816, Saudi Arabia
5
Emergency Medical Services Department, Madinah National Hospital, Madinah 11461, Saudi Arabia
6
King Faisal Specialist Hospital & Research Center, Al-Madinah Al-Munawwarah 42523, Saudi Arabia
7
Prince Mohammed Bin Abdulaziz Hospital-Al Madinah Al Munawarah-NGHA, Ministry of National Guard Health Affairs, Kingdom of Saudi Arabia, Riyadh 41511, Saudi Arabia
8
Faculty of Dentistry, British University, El Sherouk City, Suez Desert Road, Cairo 11837, Egypt
9
Department of Natural Products and Alternative Medicine, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia
10
Preparatory Year Program, Department of Chemistry, Batterjee Medical College, Jeddah 21442, Saudi Arabia
11
Department of Pharmacognosy, Faculty of Pharmacy, Assiut University, Assiut 71526, Egypt
12
Preparatory Year Program, Batterjee Medical College, Jeddah 21442, Saudi Arabia
13
Department of Pharmacognosy, Faculty of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-kharj 11942, Saudi Arabia
14
Department of Pharmacognosy, Faculty of Pharmacy, University of Khartoum, Al-Qasr Ave., Khartoum 11111, Sudan
15
Department of Pharmacognosy, Faculty of Pharmacy, Mansoura University, Mansoura 35516, Egypt
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Metabolites 2023, 13(10), 1090; https://doi.org/10.3390/metabo13101090
Submission received: 4 September 2023 / Revised: 7 October 2023 / Accepted: 17 October 2023 / Published: 18 October 2023

Abstract

:
Cyclin-dependent kinase 5 (CDK5) plays a crucial role in various biological processes, including immune response, insulin secretion regulation, apoptosis, DNA (deoxyribonucleic acid) damage response, epithelial−mesenchymal transition (EMT), cell migration and invasion, angiogenesis, and myogenesis. Overactivation of CDK5 is associated with the initiation and progression of cancer. Inhibiting CDK5 has shown potential in suppressing cancer development. Despite advancements in CDK5-targeted inhibitor research, the range of compounds available for clinical and preclinical trials remains limited. The marine environment has emerged as a prolific source of diverse natural products with noteworthy biological activities, including anti-cancer properties. In this study, we screened a library of 47,450 marine natural compounds from the comprehensive marine natural product database (CMNPD) to assess their binding affinity with CDK5. Marine compounds demonstrating superior binding affinity compared to a reference compound were identified through high-throughput virtual screening, standard precision and extra-precision Glide docking modes. Refinement of the selected molecules involved evaluating molecular mechanics–generalized born surface area (MM/GBSA) free binding energy. The three most promising compounds, (excoecariphenol B, excoecariphenol A, and zyzzyanone B), along with the reference, exhibiting favorable binding characteristics were chosen for molecular dynamics (MD) simulations for 200 nanoseconds. These compounds demonstrated interaction stability with the target during MD simulations. The marine compounds identified in this study hold potential as effective CDK5 inhibitors and warrant subsequent experimental validation.

1. Introduction

Cancer is the world’s second most prevalent cause of death after coronary heart disease (CHD), and it is expected that it will become the first in the upcoming four decades [1]. Worldwide, the number of cancer-related deaths was 10 million out of 19.3 million according to updated WHO GLOBOCAN estimates. In descending order of fatality, lung, colorectal, liver, stomach, and breast cancers ranked as the most deadly forms of cancer [2].
CDK5 is a 292-amino-acids sequence protein that descends from the cyclin protein kinase (CDK) family alongside another 20 members [3,4]. For decades, CDKs have been projected as cell cycle regulators; they control G1 to S phase transition depending on their activation by cyclin [5]. CDK5 structure is 60% similar to other CDKs (i.e., CDK1 and CDK2 [3]); however, CDK5 is non-cyclin dependent since it lacks the amino acid motif that interacts with cyclin, and is not considered a direct cell cycle regulator [6,7]. CDK5 consists of an ATP (adenosine triphosphate) binding pocket intercalated between two lobes: the N-lobe which contains a 40 s subunit that is responsible for interaction with CDK5 activators, and the C-lobe that moderates ATP and substrate binding to the N-terminus [8,9].
Up-to-date research established an evident correlation between CDK5 aberrant activation and/or overexpression and tumorigenesis, mutagenesis dissemination, and cancer disease progression in a variety of cancers [10,11]. CDK5 has unignorable role in tumor progression and metastasis in numerous cancers, including head and neck squamous cell carcinomas, hepatocellular carcinoma, prostate cancer, medullary thyroid carcinoma, ovarian cancer, and chemotherapy-resistant cervical tumors of the uterine cervix and the bladder [12,13,14,15,16,17,18,19,20,21].
Despite the extensive studies that correlating CDK5 expression and cancer occurrence, drug therapies that selectively inhibit CDK5 were not available until very recently [22]. Selective targeting of CDK5 is extremely difficult due to the high degree of similarity with other CDKs [23]. Pan CDK inhibitors remain the chief therapeutic option for the treatment of CDK5-related illnesses [22]. First-generation pan CDK inhibitors, Roscovitine and flavopiridol, have entered Phase II clinical trials. However, they exhibited limited efficacy, and their side effects were deleterious. The second-generation pan CDKI, Dinaciclib, in phase II/III is far more potent than the first generation, but it displays low inhibitory activity against CDK5 and disastrous off-target toxicity [24,25]. At present, GFB-12811, modified thiazoles, and pyrido (3,2, d) pyridines have proven their selective CDK5 inhibitory properties. Nonetheless, they are still in the experimental phase [23,26,27].
Marine-based drug discovery has emerged as an exciting field in modern drug research and development since the discovery of the first marine derivative, MNP spongothymidine [28]. Currently, a lot of marine-derived medications have gained FDA approval (Ziconotide, Trabectedin, Eribulin, etc.), and many of them have entered clinical trials. Furthermore, a comprehensive marine natural product database (CMNPD) with more than 30,000 derivatives has become available for research purposes [29].
Computer-aided drug design is a multidisciplinary approach that intends to apply the fundamentals of computer sciences and software algorithms to assist the drug discovery process in terms of time and cost reduction [30,31,32].
The primary aim of this research endeavor is to exploit the various computer-aided drug design methodologies to systematically identify novel lead compounds demonstrating exclusive inhibitory activity against CDK5 within the context of comprehensive marine natural product database. Through comprehensive computational analyses and simulations, we seek to unravel promising candidate compounds with the potential to modulate CDK5, thus advancing our understanding of their therapeutic utility in CMNPD.

2. Materials and Methods

In this study, the Maestro version 12.8 software suite of Schrödinger and Desmond v6.5 module were simultaneously utilized to implement computational work to CDK5 protein structure analyses and the collected marine natural products.

2.1. Protein Preparation

The 3D structure of the crystalized CDK5 enzyme (PDB ID: 7VDQ) bound to a known potent selective inhibitor was retrieved from the RCSB Protein Data Bank [33]. The structure was submitted to the protein preparation wizard tool in Maestro 12.8 software suite. Correspondingly, the protein structure was preprocessed. Throughout this step, hydrogen atoms were added, bond orders for the constitutional amino acid residues were assigned, heteroatoms and solvent molecules beyond 5 Ǻ were deleted, hydrogen bonds were optimized, and zero-order bonds with metals were generated. In addition, the function of prime module Maestro was recruited to add the omitted side chains and loops. Eventually, the energy of the preprocessed and optimized structure was minimized based on OPLS force field energy parameters [34,35].

2.2. Ligand Preparation

A set of 47,450 2D structures was obtained from a comprehensive marine natural product database. Upon that, the default settings of LigPrep tool in Maestro 12.8 software suite were employed to prepare the collected ligands and decrease probable computational errors [36]. Subsequently, precise, energy-minimized, and refined 3D chemical assemblies were obtained, probable mistakes in the structures of the ligands library were removed, and the possible ionization states were generated at the physiological pH (7+/−2).

2.3. Molecular Docking

Following preparation of the collected dataset, the default settings of the Glide receptor grid generator tool in Maestro 12.8 were employed to create a cubic grid of 20 Ǻ, taking the centroid of the co-crystalized inhibitor of CDK5 as a reference place for ligand–receptor fitting [37].
The Glide module in Maestro 12.8 was implemented. The Glide module is a well-established docking platform that applies the Monte Carlo technique for ligand–receptor sampling and glide docking score (GScore), which is a modified form of ChemScore, as an energy scoring function retrieving the energy functions from the OPLS force field. The Glide module provides three modes of docking with different degrees of thoroughness and time labor [38]. High-throughput virtual screening (HTVS) represents the fastest mode; however, it uses a more forgiving scoring function and less extensive sampling techniques. It has been exploited to exclude compounds that possess poor receptor complementarity and reduce the number of tested molecules in the subsequent steps of the docking process. Thereafter, the top scoring molecules were subjected to Glide standard-precision docking (SP) followed by extra-precision (XP) docking that incorporates more extensive sampling and reliable scoring functions, and the top ten compounds were selected for further analysis based on their Glide docking score values [39,40,41].

2.4. MM/GBSA Binding Free Energy Calculation

The prime molecular mechanics–generalized born area (MM/GBSA) tool of Maestro is designed to calculate the binding free energy of the system based on the energy difference between the unbound minimized ligand–protein complex and the minimized complex after ligand binding, as depicted by the equation using OPLS_2005 as the force field and the VSGB2.0 solvent model [42].
∆G(bind) = Ecomplex(minimized) − Eligand(minimized) + Ereceptor(minimized)
The ten poses that have been preserved from the Glide XP docking output file underwent MM/GBSA calculations using the Prime module in Maestro version 12.8 software package and binding free energy (expressed as MM/GBSA ∆Gbind), and the energetic contribution in the ligand–receptor complex was estimated.

2.5. ADME and Toxicity Prediction

Early ADMET (absorption, distribution, metabolism, excretion, and toxicity) pro-filing is of a paramount essence within the drug discovery journey. Among the clinical trial data between 2010 and 2017, poor drug-like properties aborted the discovery of approximately 10% to 15%, whereas undesirable toxicity prohibited the development of 30% of newly discovered lead molecules [43]. In this study, we exploited the QikProp tool in the Maestro 12.8 software suite of Schrödinger to forecast the pharmacokinetic profile of the refined molecules in MM/GBSA calculations. Afterward, Pro Tox II webserver was recruited to get a deeper look into the toxicological tendencies of the selected compounds. The selection of these tools was guided with their light computational requirements and tremendous simulation power [44,45].

2.6. Molecular Dynamics Simulation

Molecular dynamics (MD) simulation is one of the most influential trends in computational drug design. MD implicates molecular mechanics concepts, Newton’s law of motion, and mathematical principles to simulate particulate movement, and gives a thorough analysis of their interaction forces while considering the physiological barriers that face the drug molecule inside the human body. Moreover, MD investigates the binding energy, unbinding energy, and conformational changes underlying interactions between the molecules and their receptors [46].
From MM/GBSA binding free energy results, the three compounds with the lowest binding energy were subjected to MD simulation to assess the stability of their interaction with CDK5 (PDB ID: 7VDQ). The academic version of the Desmond version 6.5 software package of D.E Shaw research was used to run the simulation process. First, the physiological simulation system was created by placing the nominated ligands with CDK5 3D structure in an orthorhombic box with the dimensions of 10 × 10 × 10 Ǻ, and transferable intermolecular potential three-point TIP3P was added to the box as a solvent model. Furthermore, the ions of Na+ and Cl− were added to maintain the electrostatic neutralization of the system and attain the physiologic concentration of 150 mM. The system energy was minimized using the energy parameters of the OPL3Se force field. To generate a state of equilibrium, an isothermal–isochoric (NVT) ensemble was used to maintain at 300 K accompanied by an isothermal–isobaric (NPT) ensemble to attain fixed values for pressure and temperature. For the sake of maintaining conditions of steady pressure and temperature, the Nose–Hoover chain thermostat and the Martyna–Tobias–Klein barostat were adopted. Finally, molecular dynamics simulation was run for 200 ns and the resultant trajectories were analyzed via Maestro’s Simulation Interaction Diagram panel and evaluated according to their root mean square deviation RMSD and root mean square fluctuation RMSF [47].

3. Results and Discussion

3.1. Molecular Docking and MM/GBSA Analysis

Molecular docking stands as an efficient and economically viable method within the realm of computational drug design. It serves the purpose of identifying and evaluating essential molecular interactions occurring between ligands and receptors [48,49]. In the present work, a collection of 47,450 compounds extracted from marine sources [29] underwent molecular docking analysis against CDK5, employing a high-throughput virtual screening (HTVS) mode. Subsequently, the top 1000 compounds were subjected to a second round of docking using the standard-precision (SP) mode, and the top 100 compounds from this selection were further evaluated using the extra-precision (XP) mode. The XP mode is known to be the most accurate and precise among all other modes, despite being computationally extensive. Hence, for the primary filtration of compounds, we relied on their XP docking scores. The scoring function that is built in with Glide utilizes both empirical and force-field-based factors for calculating the binding energy and filtering the best pose of docking [41,50]. The co-crystallized reference molecule, which is already reported as a potent and selective CDK5 inhibitor [23], presented a very high docking score of −13.516 kcal/mol, as summarized in Table 1. Comparably, the highest-ranked inhibitors displayed a scoring spectrum, ranging from −12.015 kcal/mol to as low as −13.863 kcal/mol, showcasing a diversity of binding affinities across the compounds studied.
It can be seen from the docking scores that there are variable binding affinities conveyed by molecular docking. However, it is evident that further refinements are necessary for these compounds before proceeding with additional processing or analysis. One example of refinement means is the use of MM/GBSA (molecular mechanics/generalized born surface area). Numerous investigations have consistently demonstrated that MM/GBSA generally surpasses the scoring functions employed in docking algorithms when it comes to accurately ranking ligands based on their binding affinity [51,52,53]. To achieve the goal of safely fine-tuning the docked compounds, we used the Prime module to calculate the MM/GBSA binding free energies. As reported in Table 1, the co-crystallized reference inhibitor had a binding energy of −74.93 kcal/mol and the rest of the top 10 compounds presented values ranging from −47.20 kcal/mol to −83.78 kcal/mol. Figure 1 summarizes the chemical structures of the top 10 compounds.
To proceed with more in-depth analysis, we chose the three most promising candidates: excoecariphenol B, excoecariphenol A, and zyzzyanone B. These compounds exhibit MM/GBSA values of −83.78 kcal/mol, −71.58 kcal/mol, and −71.05 kcal/mol, respectively.
According to the summary of 2D interactions observed in Figure 2 and detailed in Table 2, it is evident that the co-crystallized inhibitor effectively replicated critical hydrogen bonding interactions with specific amino acid residues, including Lys33, Gln51, Cys83, Asp86, and Asp144. Additionally, it established a noteworthy salt bridge interaction with Asp144 and a Pi-cation interaction with Tyr15.
Table 1. Docking scores of the top 10 ranking compounds with their MM/GBSA free binding energy results.
Table 1. Docking scores of the top 10 ranking compounds with their MM/GBSA free binding energy results.
Compound IDNameDocking ScoreXP G ScoreMM/GBSA
dG Bind
Taxonomy
CMNPD23126excoecariphenol B−13.343−13.343−83.78Excoecaria agallocha [54]
Co-crystallized reference2-[[7-[[2-fluoranyl-4-[3-(hydroxymethyl) pyrazol-1-yl] phenyl] amino]-1,6-naphthyridin-2-yl]-(1-methylpiperidin-4-yl) amino] ethanoic acid−13.516−13.516−74.93-
CMNPD23125excoecariphenol A−12.018−12.018−71.58Excoecaria agallocha [54]
CMNPD15113zyzzyanone B−12.791−12.791−71.05Zyzzya fuliginosa [55,56]
CMNPD15115zyzzyanone D−12.058−12.058−69.36Zyzzya fuliginosa [55,56]
CMNPD23074lamellarin A5−12.015−12.015−61.95Didemnum species [57]
CMNPD30137streptocarbazole E−13.863−13.863−60.50Streptomyces species [58]
CMNPD301353′-O-demethyl-4′-N-demethyl-4′-N-acetyl-4′-epi-staurosporine−12.547−12.547−56.15Streptomyces species [59]
CMNPD30138streptocarbazole C−12.190−12.190−53.82Streptomyces species [58]
CMNPD13175thalassiolin A−12.310−12.310−51.22Thalassia testudinum [60]
CMNPD29962-hydroxygarvin B−12.449−12.449−47.20Garveia annulata [61]
Notably, excoecariphenol B demonstrated its potential by forming hydrogen bonds with Gln51, Cys83, and Asp144, contributing to its binding affinity. This compound is a flavan-3-ol that belongs to the structural superfamily of phenylpropanoids and polyketides, and it is known to be produced by the mangrove plant Excoecaria agallocha. The extracts of this plant (excoecariphenols A–D) have previously been reported to act against hepatitis C virus [54].
Conversely, the other flavan-3-ol, excoecariphenol A, exhibited hydrogen bonding interactions with Gln51 and Cys83, underscoring its capacity to engage with the receptor. Meanwhile, the third compound, zyzzyanone B, displayed a broader spectrum of interactions, including hydrogen bonds with Gln51, Glu81, Asp86, and Cys83. It also featured a salt bridge interaction with Asp86 and a Pi–Pi stacking interaction with Phe80, further emphasizing its binding versatility. Zyzzyanone B is a pyrroloindole that descends from the structural superfamily of organoheterocyclic compounds, synthesized by the sponge Zyzzya fuliginosa. Zyzzyanone B has previously been reported to possess cytotoxic activity and antioxidant effects [55,62].
Furthermore, a closer examination of Table 2 reveals that both the co-crystallized reference inhibitor and excoecariphenol B established hydrophobic contacts with the same set of hydrophobic amino acid residues, namely Ile10, Val18, Ala31, Val64, Phe80, Phe82, Cys83, and Leu133. Additionally, excoecariphenol A reflected these contacts, except for Tyr15 and Val18, while zyzzyanone B maintained similar interactions, except for Val64. These observations highlight the potential binding similarities and distinctions among these compounds, shedding light on their interactions with the receptor. It has been previously documented that the active site of CDK5 closely resembles that of CDK1 and CDK2, posing a notable challenge in achieving selectivity. These three kinases share an almost identical ATP binding pocket. The only difference is that the residues Cys83 and Ser84 in CDK5 are Leu and Ser in CDK1, and Leu and His in CDK2. Importantly, there have been no reported instances of compounds interacting with these distinguishing features thus far. Of interest here, the residue Cys83 is engaged with our compounds (Figure 3) in hydrogen bonding, and this is a unique pattern of binding as there are no previously reported compounds acting similarly.
Henceforth, we propose the possibility that the investigated compounds could potentially exhibit a distinctive selectivity for binding to CDK5. This selectivity is substantiated by their comparable binding energies and binding patterns when compared to the known co-crystallized inhibitor. This suggests the potential for these compounds to display a distinctive selectivity for CDK5 kinase activity.
In order to strengthen their suitability in interacting with CDK5, subsequent molecular dynamics simulation is performed.

3.2. ADME and Toxicity Analysis

Utilizing the QikProp tool within the Maestro software platform version 12.8, we conducted an assessment of the pharmacokinetic properties for the compounds that exhibited superior docking scores when compared to the co-crystallized inhibitor. This analytical step was undertaken with the primary objective of determining the druggability of the top-docked compounds and exploring their potential adherence to Lipinski’s rule of five criteria. As allocated by Lipinski’s rule, a molecule is deemed druggable when it complies with the following criteria: the molecular weight does not exceed 500 Da, has a maximum of 5 hydrogen bond donors (donorHB ≤ 5), has no more than 10 hydrogen bond acceptors (acceptorHB ≤ 10), and has a predicted octanol/water partition coefficient (QPlogPo/w) lower than 5 [63,64].
In Table 3, it can be noted that zyzzyanone B fulfills all the criteria that are set by Lipinski’s rule of five. On the other hand, excoecariphenol A and excoecariphenol B both show non-conformity to the rule by violating the prescribed limits of the number of hydrogen bond donors or acceptors. Moreover, the molecular weight of excoecariphenol B was 500.473 Da, which is slightly higher than the value set by the rule (500). Additionally, two parameters related to the potential for cellular membrane penetration were subject to scrutiny: QPlogBB, an indicative of blood–brain barrier permeability; and QPP Caco-2, reflecting cell membrane permeability. It is worth noting that all compounds under investigation demonstrated favorable cellular permeability characteristics while displaying restricted ability to traverse the blood–brain barrier.
Furthermore, a pivotal parameter, QPlog HERG, responsible for predicting IC50 values pertaining to the blockage of HERG K+ channels, as presented in Table 2, conveyed that none of the compounds exhibited indications of cardiotoxicity.
Within the framework of ProTox-II (Table 4), oral toxicity is characterized in terms of the median lethal dose (LD50), denoted as milligrams per kilogram of body weight for the test population. Notably, both excoecariphenol A and excoecariphenol B fall into class 5, exhibiting predicted LD50 values of 2500 mg/kg. Conversely, zyzzyanone B is categorized as class 4, possessing a calculated LD50 of 1000 mg/kg. These findings underscore the safety profile of these compounds when orally ingested. These toxicity findings are further potentiated by the fact that all three candidate compounds possess no hepatotoxicity, mutagenicity, cytotoxicity, carcinogenicity or immunogenicity.

3.3. Molecular Dynamics (MD) Simulations Analysis

Primary molecular docking provides initial intuitions into the nature of interaction between a small molecule and its target receptor, but it does not give an account to the flexibility of the protein–ligand complexes, i.e., it does not simulate the biological environment. Fortunately, the integration of molecular docking with highly comprehensive approaches like molecular dynamics simulation increases the simulation power and yields more realistic and comprehensive computational analysis [65]. MD simulation is introduced as a valuable tool for closely examining how macromolecular systems behave dynamically when exposed to small molecules. These simulations produce numerous trajectories, each offering a wealth of information about how proteins and ligands interact. These trajectories serve as a rich data source, enabling in-depth exploration and a deeper understanding of the intricate interactions within the protein–ligand system [66].
In this work, we maintained a molecular dynamics simulation for the top three compounds that were selectively filtered via docking scores and MM/GBSA free binding energy, specifically excoecariphenol B, excoecariphenol A, and zyzzyanone B. Additionally, the co-crystallized reference inhibitor was simulated and held as a control against which the other compounds were compared. The aforementioned simulations were conducted for 200 nanoseconds, since the CDK5 enzyme is a relatively large entity and is anticipated to fluctuate the most during a simulation run.
The analysis of MD is achieved using different metrics, for instance, the root means square deviation (RMSD), which is a fundamental technique employed in the examination of atomic positional variations within amino acid structures in the presence of a ligand. The RMSD essentially measures the average square root of the deviations from the mean distances observed within the protein–ligand complex. The inspection of molecular dynamics (MD) trajectories and the subsequent assessment of RMSD plots provide a valuable initial idea about the stability of simulated protein–ligand complexes, as previously indicated [67]. For a globular protein, maintaining stability in a simulated system typically entails fluctuations within a range of 1–3 Å.
The molecular size greatly impacts the RMSD value; hence, as the size increases, it becomes difficult to directly address the MD simulation with a specific acceptable value. Importantly, as long as the RMSD value reaches a steady state of fluctuation, it can be fairly judged as a converged system [67]. For our studied compounds, Figure 4 depicts the RMSD plots for the four simulation runs during the 200 nanoseconds. It can be seen that all the compounds exhibited similar influence on the carbon alpha atoms of the protein, with an average value of 4.764 ± 0.64, indicating the occurrence of protein folding. Nonetheless, the general observation is that the protein fluctuated at the first 75 nanoseconds, then it converged and returned to fluctuate around the average value toward the simulation end-point. The fit of excoecariphenol B, excoecariphenol A, zyzzyanone B, and the co-crystallized reference inhibitor on the protein indicated that they converged after 125 ns, 40 ns, 25 ns and 75 ns, respectively.
The diversity in the behavior of amino acids and the precise sites of modifications within the protein are pivotal factors in comprehending the functional dynamics of the protein in the context of molecular dynamics simulations. Through a thorough examination of the trajectory data derived from root mean square fluctuation (RMSF) analysis, we can discern the dynamic behavior of individual amino acids [68].
The protein RMSF average values were 1.014 Å, 1.211 Å, 1.211 Å, and 1.211 Å for excoecariphenol B, excoecariphenol A, zyzzyanone B, and the co-crystallized reference inhibitor, respectively. As seen in Figure 5, there were areas of sharp peaks, particularly with excoecariphenol A as it presented contacts with Thr 14 and Tyr15, but providentially, they were brief and distant form the active site. The particular region of residues with a slight spike in the RMSF (up to 5 A) maybe due to the higher conformational flexibility of Thr14 and Tyr15. Fortunately, those residues are not essential for CDK5 binding. Overall, the RMFF patterns suggest that the simulated protein–ligand complexes are somewhat stable and with minimal fluctuations of the key amino acid residues.
The provided data in Figure 6 outline the various protein–ligand interactions observed during molecular dynamic simulations for three different compounds: excoecariphenol B, excoecariphenol A, and zyzzyanone B, as well as a reference inhibitor.
For excoecariphenol B, hydrogen bonds are primarily formed with Asp86 (77%), indicating strong interactions, and to a lesser extent with Lys128 (30%), Asn131 (20%), and Gly146 (10%). Moreover, water bridges are observed with Lys33 (5%), Ser47 (8%), Glu51 (5%), Asp126 (10%), and Asp144 (8%), showing additional contributions to ligand–receptor interactions.
As for the reference inhibitor, direct hydrogen bonds are formed with several residues, including Arg34 (120%), indicating a strong interaction. Additionally, hydrogen bonds are established with Val30 (5%), Ala21 (22%), Leu32 (10%), Lys33 (15%), and Lys75 (45%). Water bridges, a type of interaction considered significant in ligand binding, are notable with Phe19 (40%), Ala21 (20%), Val30 (20%), Leu32 (30%), Arg34 (20%), and Lys75 (15%). Moreover, hydrophobic contacts, which are pivotal for stabilizing ligand–receptor complexes, are present with Ile10 (35%), Tyr15 (60%), Kys20 (30%), Ile29 (15%), Lys75 (20%), and Phe82 (20%). An ionic interaction with Lys75 (5%) further contributes to the overall binding.
The compound excoecariphenol A is engaged in direct hydrogen bonds with Lys33 (10%), Ser47 (10%), Gln130 (110%), and Asp144 (60%), reflecting substantial interactions. In addition, water bridges are noted with Lys33 (50%), Arg36 (5%), Val44 (5%), Glu51 (8%), Cys83 (35%), Asp84 (20%), Asp86 (20%), Gln130 (5%), and Asp144 (40%), indicating their involvement in the ligand binding process. Notably, Gln130 is involved in both direct hydrogen bonds and water bridges.
For zyzzyanone B, direct hydrogen bonds are significantly established with Glu81 (140%) and Cys83 (20%), underlining strong binding. Alongside that, water bridges are prominent with Leu32 (45%), Val44 (5%), Ser47 (5%), Glu51 (7%), Val64 (45%), and Asp144 (80%). Moreover, hydrophobic contacts are present with Phe80 (100%), Phe82 (5%), and Ala31 (25%), highlighting their role in stabilization.
In comparison, excoecariphenol B forms strong hydrogen bonds with Asp86, while the reference inhibitor demonstrates a robust interaction with Arg34. Excoecariphenol A and zyzzyanone B also establish noteworthy hydrogen bonds primarily with Gln130 and Glu81, respectively. In addition, all compounds engage in water bridge interactions. Excoecariphenol A and zyzzyanone B exhibit prominent water bridge interactions with various residues. In contrast, excoecariphenol B and the reference inhibitor demonstrate water bridges with fewer residues. Conversely, hydrophobic interactions are not evident for all compounds. The reference inhibitor displayed significant hydrophobic contacts, especially with Tyr15 and Phe80. Zyzzyanone B forms a substantial hydrophobic interaction with Phe80.
In summary, each compound exhibits unique interaction patterns with the receptor, suggesting varying modes of binding. Excoecariphenol B and the reference inhibitor prioritize hydrogen bonding, while excoecariphenol A and zyzzyanone B rely on a combination of hydrogen bonding, water bridges, and hydrophobic contacts. These distinctive interaction profiles may have implications for the compounds’ selectivity and potency as ligands, emphasizing the importance of understanding these interactions in drug design and development.
Furthermore, in order to measure the total energy after MD production, post-MD MM/GBSA calculations were performed. The co-crystallized reference and the best three compounds’ trajectories underwent MM/GBSA calculations. Due to the substantial computational expense, only 10 frames were chosen from a total of 2000 frames. As can be noted in Figure 7, the reference showed the lowest energy value, with a mean value of −37.12 ± 6.212 kcal/mol. On the other hand, excoecariphenol B, excoecariphenol A, and zyzzyanone B presented energy averages of −21.42 ± 9.155 kcal/mol, −23.51 ± 2.453 kcal/mol, and −35.45 ± 4.042 kcal/mol, respectively. Notably, zyzzyanone B is the lowest of the candidate compounds, implying its superior binding.

4. Conclusions

CDK inhibitors hold immense promise in halting tumor development, progression, and metastasis, mediated through diverse pathways such as CDK5–FAK, PPARγ, PP2A, Hippo, and Wnt/β-catenin, among others. Despite growing demand, the clinical landscape has yet to witness the emergence of a selective CDK5 inhibitor. This research sought to address this gap by employing a range of computational techniques. Our study involved the screening of a substantial library of 47,450 marine natural products employing a multi-stage Glide docking process (HTVS, SP, and XP). Subsequently, the top 10 docked complexes were subjected to careful refinement via MM/GBSA analysis. Among these, three candidate compounds, specifically excoecariphenol B, excoecariphenol A, and zyzzyanone B, exhibited promising attributes. Through the extensive 200 nanoseconds of molecular dynamics simulations, we scrutinized the stability of these compounds, further enhancing our confidence in their potential as CDK5 inhibitors. Taken together, our in silico methodologies suggest that these marine-based compounds require further experimental investigation to ascertain their selectivity and therapeutic efficacy as novel CDK5 inhibitors. This research endeavors to contribute to the ongoing quest for innovative anti-cancer therapies, with a particular focus on CDK5 inhibition.

Author Contributions

Conceptualization, A.A.A.; methodology, A.A.A., M.A.A. and T.H.S.; software, A.A.A.; validation, Y.S.A., S.A.S. and N.M.A.; formal analysis A.F.E., R.Y.A. and E.A.A.; investigation, Y.S.A., S.A.S., N.M.A. and S.G.A.M.; resources, S.G.A.M., G.A.M., S.R.M.I., H.G.A.H., A.F.E., R.Y.A. and E.A.A.; writing—original draft preparation, A.A.A, T.H.S., M.A.A. and S.R.M.I.; writing—review and editing, S.G.A.M., G.A.M., S.R.M.I. and H.G.A.H.; visualization, A.A.A.; project administration, W.O., A.A. and A.E.S.; funding acquisition, W.O., A.A. and A.E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by funding from the Prince Sattam bin Abdulaziz University. Project number: (PSAU/2023/R/1444).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data presented in this study are available on public databases mentioned in the article, and the derived data are available on request from the corresponding author. The data are not publicly available due to confidentiality agreement with research.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ATPAdenosine triphosphate
CDK5Cyclin-dependent kinase 5
CHDCoronary heart disease
CMNPDComprehensive marine natural product database
FDAFood and Drug Administration
HTVSHigh-throughput virtual screening
MDMolecular dynamics
MM/GBSAMolecular mechanics/generalized born surface area
NPTIsothermal–isobaric
OPLSOptimized potentials for liquid simulations
RMSDRoot means square deviation
RMSFRoot mean square fluctuation
TIP3PTransferable interaction potential
PDBProtein data bank
SPStandard precision
XPExtra precision

References

  1. Mattiuzzi, C.; Lippi, G. Current Cancer Epidemiology. J. Epidemiol. Glob. Health 2019, 9, 217–222. [Google Scholar] [CrossRef] [PubMed]
  2. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, B.; Tan, V.B.C.; Lim, K.M.; Tay, T.E. The Activation and Inhibition of Cyclin-Dependent Kinase-5 by Phosphorylation. Biochemistry 2007, 46, 10841–10851. [Google Scholar] [CrossRef] [PubMed]
  4. Hellmich, M.R.; Pant, H.C.; Wada, E.; Battey, J.F. Neuronal Cdc2-like Kinase: A Cdc2-Related Protein Kinase with Predominantly Neuronal Expression. Proc. Natl. Acad. Sci. USA 1992, 89, 10867–10871. [Google Scholar] [CrossRef] [PubMed]
  5. Enserink, J.M.; Chymkowitch, P. Cell Cycle-Dependent Transcription: The Cyclin Dependent Kinase Cdk1 Is a Direct Regulator of Basal Transcription Machineries. Int. J. Mol. Sci. 2022, 23, 1293. [Google Scholar] [CrossRef]
  6. Futatsugi, A.; Utreras, E.; Rudrabhatla, P.; Jaffe, H.; Pant, H.C.; Kulkarni, A.B. Cyclin-Dependent Kinase 5 Regulates E2F Transcription Factor through Phosphorylation of Rb Protein in Neurons. Cell Cycle 2012, 11, 1603–1610. [Google Scholar] [CrossRef]
  7. Gonzalez, L.; Domingo-Muelas, A.; Duart-Abadia, P.; Nuñez, M.; Mikolcevic, P.; Llonch, E.; Cubillos-Rojas, M.; Cánovas, B.; Forrow, S.M.A.; Morante-Redolat, J.M.; et al. The Atypical CDK Activator RingoA/Spy1 Regulates Exit from Quiescence in Neural Stem Cells. iScience 2023, 26, 106202. [Google Scholar] [CrossRef]
  8. Otyepka, M.; Bártová, I.; Kříž, Z.; Koča, J. Different Mechanisms of CDK5 and CDK2 Activation as Revealed by CDK5/P25 and CDK2/Cyclin A Dynamics. J. Biol. Chem. 2006, 281, 7271–7281. [Google Scholar] [CrossRef]
  9. Tarricone, C.; Dhavan, R.; Peng, J.; Areces, L.B.; Tsai, L.-H.; Musacchio, A. Structure and Regulation of the CDK5-P25nck5a Complex. Mol. Cell 2001, 8, 657–669. [Google Scholar] [CrossRef]
  10. Semenov, I.; Akyuz, C.; Roginskaya, V.; Chauhan, D.; Corey, S.J. Growth Inhibition and Apoptosis of Myeloma Cells by the CDK Inhibitor Flavopiridol. Leuk. Res. 2002, 26, 271–280. [Google Scholar] [CrossRef]
  11. Bei, Y.; Cheng, N.; Chen, T.; Shu, Y.; Yang, Y.; Yang, N.; Zhou, X.; Liu, B.; Wei, J.; Liu, Q.; et al. CDK5 Inhibition Abrogates TNBC Stem-Cell Property and Enhances Anti-PD-1 Therapy. Adv. Sci. 2020, 7, 2001417. [Google Scholar] [CrossRef] [PubMed]
  12. De, S.; Holvey-Bates, E.G.; Mahen, K.; Willard, B.; Stark, G.R. The Ubiquitin E3 Ligase FBXO22 Degrades PD-L1 and Sensitizes Cancer Cells to DNA Damage. Proc. Natl. Acad. Sci. USA 2021, 118, e2112674118. [Google Scholar] [CrossRef] [PubMed]
  13. Wu, J.; Tan, Z.; Li, H.; Lin, M.; Jiang, Y.; Liang, L.; Ma, Q.; Gou, J.; Ning, L.; Li, X.; et al. Melatonin Reduces Proliferation and Promotes Apoptosis of Bladder Cancer Cells by Suppressing O-GlcNAcylation of Cyclin-Dependent-like Kinase 5. J. Pineal Res. 2021, 71, e12765. [Google Scholar] [CrossRef]
  14. Feldmann, G.; Mishra, A.; Hong, S.M.; Bisht, S.; Strock, C.J.; Ball, D.W.; Goggins, M.; Maitra, A.; Nelkin, B.D. Inhibiting the Cyclin-Dependent Kinase CDK5 Blocks Pancreatic Cancer Formation and Progression through the Suppression of Ras-Ral Signaling. Cancer Res. 2010, 70, 4460–4469. [Google Scholar] [CrossRef]
  15. Huang, J.; Chen, P.; Liu, K.; Liu, J.; Zhou, B.; Wu, R.; Peng, Q.; Liu, Z.X.; Li, C.; Kroemer, G.; et al. CDK1/2/5 Inhibition Overcomes IFNG-Mediated Adaptive Immune Resistance in Pancreatic Cancer. Gut 2021, 70, 890–899. [Google Scholar] [CrossRef]
  16. Lindqvist, J.; Imanishi, S.Y.; Torvaldson, E.; Malinen, M.; Remes, M.; Örn, F.; Palvimo, J.J.; Eriksson, J.E. Cyclin-Dependent Kinase 5 Acts as a Critical Determinant of AKT-Dependent Proliferation and Regulates Differential Gene Expression by the Androgen Receptor in Prostate Cancer Cells. Mol. Biol. Cell 2015, 26, 1971–1984. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, X.; Liu, P.; Zhang, C.; Chiewchengchol, D.; Zhao, F.; Yu, H.; Li, J.; Kambara, H.; Luo, K.Y.; Venkataraman, A.; et al. Positive Regulation of Interleukin-1β Bioactivity by Physiological ROS-Mediated Cysteine S-Glutathionylation. Cell Rep. 2017, 20, 224–235. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, X.X.; Xie, F.F.; Zhu, X.J.; Lin, F.; Pan, S.S.; Gong, L.H.; Qiu, J.G.; Zhang, W.J.; Jiang, Q.W.; Mei, X.L.; et al. Cyclin-Dependent Kinase Inhibitor Dinaciclib Potently Synergizes with Cisplatin in Preclinical Models of Ovarian Cancer. Oncotarget 2015, 6, 14926–14939. [Google Scholar] [CrossRef]
  19. Ehrlich, S.M.; Liebl, J.; Ardelt, M.A.; Lehr, T.; De Toni, E.N.; Mayr, D.; Brandl, L.; Kirchner, T.; Zahler, S.; Gerbes, A.L.; et al. Targeting Cyclin Dependent Kinase 5 in Hepatocellular Carcinoma—A Novel Therapeutic Approach. J. Hepatol. 2015, 63, 102–113. [Google Scholar] [CrossRef]
  20. Pozo, K.; Castro-Rivera, E.; Tan, C.; Plattner, F.; Schwach, G.; Siegl, V.; Meyer, D.; Guo, A.; Gundara, J.; Mettlach, G.; et al. The Role of Cdk5 in Neuroendocrine Thyroid Cancer. Cancer Cell 2013, 24, 499–511. [Google Scholar] [CrossRef]
  21. Sun, S.-S.; Zhou, X.; Huang, Y.-Y.; Kong, L.-P.; Mei, M.; Guo, W.-Y.; Zhao, M.-H.; Ren, Y.; Shen, Q.; Zhang, L. Targeting STAT3/MiR-21 Axis Inhibits Epithelial-Mesenchymal Transition via Regulating CDK5 in Head and Neck Squamous Cell Carcinoma. Mol. Cancer 2015, 14, 213. [Google Scholar] [CrossRef] [PubMed]
  22. Tang, W.; Lin, C.; Yu, Q.; Zhang, D.; Liu, Y.; Zhang, L.; Zhou, Z.; Zhang, J.; Ouyang, L. Novel Medicinal Chemistry Strategies Targeting CDK5 for Drug Discovery. J. Med. Chem. 2023, 66, 7140–7161. [Google Scholar] [CrossRef] [PubMed]
  23. Daniels, M.H.; Malojcic, G.; Clugston, S.L.; Williams, B.; Coeffet-Le Gal, M.; Pan-Zhou, X.R.; Venkatachalan, S.; Harmange, J.C.; Ledeboer, M. Discovery and Optimization of Highly Selective Inhibitors of CDK5. J. Med. Chem. 2022, 65, 3575–3596. [Google Scholar] [CrossRef]
  24. Senderowicz, A.M.; Headlee, D.; Stinson, S.F.; Lush, R.M.; Kalil, N.; Villalba, L.; Hill, K.; Steinberg, S.M.; Figg, W.D.; Tompkins, A.; et al. Phase I Trial of Continuous Infusion Flavopiridol, a Novel Cyclin-Dependent Kinase Inhibitor, in Patients with Refractory Neoplasms. J. Clin. Oncol. 1998, 16, 2986–2999. [Google Scholar] [CrossRef]
  25. Tan, A.R.; Yang, X.; Berman, A.; Zhai, S.; Sparreboom, A.; Parr, A.L.; Chow, C.; Brahim, J.S.; Steinberg, S.M.; Figg, W.D.; et al. Phase I Trial of the Cyclin-Dependent Kinase Inhibitor Flavopiridol in Combination with Docetaxel in Patients with Metastatic Breast Cancer. Clin. Cancer Res. 2004, 10, 5038–5047. [Google Scholar] [CrossRef]
  26. Shiradkar, M.R.; Akula, K.C.; Dasari, V.; Baru, V.; Chiningiri, B.; Gandhi, S.; Kaur, R. Clubbed Thiazoles by MAOS: A Novel Approach to Cyclin-Dependent Kinase 5/P25 Inhibitors as a Potential Treatment for Alzheimer’s Disease. Bioorg. Med. Chem. 2007, 15, 2601–2610. [Google Scholar] [CrossRef]
  27. Lenjisa, J.L.; Tadesse, S.; Khair, N.Z.; Kumarasiri, M.; Yu, M.; Albrecht, H.; Milne, R.; Wang, S. CDK5 in Oncology: Recent Advances and Future Prospects. Future Med. Chem. 2017, 9, 1939–1962. [Google Scholar] [CrossRef]
  28. Tekman, S.; Öner, N. A Pentose Identified in the Carbohydrate Group of Spongin. Nature 1963, 200, 77–78. [Google Scholar] [CrossRef]
  29. Lyu, C.; Chen, T.; Qiang, B.; Liu, N.; Wang, H.; Zhang, L.; Liu, Z. CMNPD: A Comprehensive Marine Natural Products Database towards Facilitating Drug Discovery from the Ocean. Nucleic Acids Res. 2021, 49, D509–D515. [Google Scholar] [CrossRef]
  30. Parate, S.; Kumar, V.; Hong, J.C.; Lee, K.W. Investigation of Marine-Derived Natural Products as Raf Kinase Inhibitory Protein (Rkip)-Binding Ligands. Mar. Drugs 2021, 19, 581. [Google Scholar] [CrossRef]
  31. Parate, S.; Kumar, V.; Lee, G.; Rampogu, S.; Hong, J.C.; Lee, K.W. Marine-Derived Natural Products as ATP-Competitive MTOR Kinase Inhibitors for Cancer Therapeutics. Pharmaceuticals 2021, 14, 282. [Google Scholar] [CrossRef]
  32. Veselovsky, A.; Ivanov, A. Strategy of Computer-Aided Drug Design. Curr. Drug Target Infectious Disord. 2003, 3, 33–40. [Google Scholar] [CrossRef]
  33. Saur, I.M.L.; Panstruga, R.; Schulze-Lefert, P. NOD-like Receptor-Mediated Plant Immunity: From Structure to Cell Death. Nat. Rev. Immunol. 2021, 21, 305–318. [Google Scholar] [CrossRef] [PubMed]
  34. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and Ligand Preparation: Parameters, Protocols, and Influence on Virtual Screening Enrichments. J. Comput. Aided. Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  35. Jacobson, M.P.; Friesner, R.A.; Xiang, Z.; Honig, B. On the Role of the Crystal Environment in Determining Protein Side-Chain Conformations. J. Mol. Biol. 2002, 320, 597–608. [Google Scholar] [CrossRef] [PubMed]
  36. Shoaib, T.H.; Abdelmoniem, N.; Mukhtar, R.M.; Alqhtani, A.T.; Alalawi, A.L.; Alawaji, R.; Althubyani, M.S.; Mohamed, S.G.A.; Mohamed, G.A.; Ibrahim, S.R.M.; et al. Molecular Docking and Molecular Dynamics Studies Reveal the Anticancer Potential of Medicinal-Plant-Derived Lignans as MDM2-P53 Interaction Inhibitors. Molecules 2023, 28, 6665. [Google Scholar] [CrossRef] [PubMed]
  37. Sherman, W.; Day, T.; Jacobson, M.P.; Friesner, R.A.; Farid, R. Novel Procedure for Modeling Ligand/Receptor Induced Fit Effects. J. Med. Chem. 2006, 49, 534–553. [Google Scholar] [CrossRef] [PubMed]
  38. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef] [PubMed]
  39. Dhasmana, A.; Raza, S.; Jahan, R.; Lohani, M.; Arif, J.M. High-Throughput Virtual Screening (HTVS) of Natural Compounds and Exploration of Their Biomolecular Mechanisms: An In Silico Approach. In New Look to Phytomedicine: Advancements in Herbal Products as Novel Drug Leads; Elsevier: Amsterdam, The Netherlands, 2018; pp. 523–548. ISBN 9780128146200. [Google Scholar]
  40. Repasky, M.P.; Murphy, R.B.; Banks, J.L.; Greenwood, J.R.; Tubert-Brohman, I.; Bhat, S.; Friesner, R.A. Docking Performance of the Glide Program as Evaluated on the Astex and DUD Datasets: A Complete Set of Glide SP Results and Selected Results for a New Scoring Function Integrating WaterMap and Glide. J. Comput. Aided. Mol. Des. 2012, 26, 787–799. [Google Scholar] [CrossRef]
  41. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [PubMed]
  42. Beard, H.; Cholleti, A.; Pearlman, D.; Sherman, W.; Loving, K.A. Applying Physics-Based Scoring to Calculate Free Energies of Binding for Single Amino Acid Mutations in Protein-Protein Complexes. PLoS ONE 2013, 8, e82849. [Google Scholar] [CrossRef]
  43. Sun, D.; Gao, W.; Hu, H.; Zhou, S. Why 90% of Clinical Drug Development Fails and How to Improve It? Acta Pharm. Sin. B 2022, 12, 3049–3062. [Google Scholar] [CrossRef] [PubMed]
  44. Ioakimidis, L.; Thoukydidis, L.; Mirza, A.; Naeem, S.; Reynisson, J. Benchmarking the Reliability of QikProp. Correlation between Experimental and Predicted Values. QSAR Comb. Sci. 2008, 27, 445–456. [Google Scholar] [CrossRef]
  45. Banerjee, P.; Eckert, A.O.; Schrey, A.K.; Preissner, R. ProTox-II: A Webserver for the Prediction of Toxicity of Chemicals. Nucleic Acids Res. 2018, 46, W257–W263. [Google Scholar] [CrossRef]
  46. Singh, S.; Bani Baker, Q.; Singh, D.B. Molecular Docking and Molecular Dynamics Simulation. In Bioinformatics: Methods and Applications; Elsevier: Amsterdam, The Netherlands, 2021; pp. 291–304. ISBN 9780323897754. [Google Scholar]
  47. Mukhtar, R.M.; Abdelmoniem, N.; Elrufaie, H.A.; Edris, A.; Ghaboosh, H.; Mahgoub, M.A.; Garelnabi, E.A.E.; Osman, W.; Sherif, A.E.; Ashour, A.; et al. Unlocking the Potential of Approved Drugs for the Allosteric Inhibition of Tropomyosin-Receptor Kinase A Using Molecular Docking and Molecular Dynamics Studies. Front. Chem. 2023, 11, 1205724. [Google Scholar] [CrossRef]
  48. Meng, X.-Y.; Zhang, H.-X.; Mezei, M.; Cui, M. Molecular Docking: A Powerful Approach for Structure-Based Drug Discovery. Curr. Comput. Aided-Drug Des. 2012, 7, 146–157. [Google Scholar] [CrossRef]
  49. Romano, T. Kroemer Structure-Based Drug Design: Docking and Scoring. Curr. Protein Pept. Sci. 2007, 8, 312–328. [Google Scholar] [CrossRef]
  50. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 2. Enrichment Factors in Database Screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef]
  51. Sun, H.; Li, Y.; Shen, M.; Tian, S.; Xu, L.; Pan, P.; Guan, Y.; Hou, T. Assessing the Performance of MM/PBSA and MM/GBSA Methods. 5. Improved Docking Performance Using High Solute Dielectric Constant MM/GBSA and MM/PBSA Rescoring. Phys. Chem. Chem. Phys. 2014, 16, 22035–22045. [Google Scholar] [CrossRef]
  52. Ahinko, M.; Niinivehmas, S.; Jokinen, E.; Pentikäinen, O.T. Suitability of MMGBSA for the Selection of Correct Ligand Binding Modes from Docking Results. Chem. Biol. Drug Des. 2019, 93, 522–538. [Google Scholar] [CrossRef] [PubMed]
  53. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the Performance of the Molecular Mechanics/Poisson Boltzmann Surface Area and Molecular Mechanics/Generalized Born Surface Area Methods. II. the Accuracy of Ranking Poses Generated from Docking. J. Comput. Chem. 2011, 32, 866–877. [Google Scholar] [CrossRef] [PubMed]
  54. Li, Y.; Yu, S.; Liu, D.; Proksch, P.; Lin, W. Inhibitory Effects of Polyphenols toward HCV from the Mangrove Plant Excoecaria agallocha L. Bioorganic Med. Chem. Lett. 2012, 22, 1099–1102. [Google Scholar] [CrossRef] [PubMed]
  55. Utkina, N.K.; Makarchenko, A.E.; Denisenko, V.A. Zyzzyanones B-D, Dipyrroloquinones from the Marine Sponge Zyzzya Fuliginosa. J. Nat. Prod. 2005, 68, 1424–1427. [Google Scholar] [CrossRef] [PubMed]
  56. Nijampatnam, B.; Dutta, S.; Velu, S.E. Recent Advances in Isolation, Synthesis, and Evaluation of Bioactivities of Bispyrroloquinone Alkaloids of Marine Origin. Chin. J. Nat. Med. 2015, 13, 561–577. [Google Scholar] [CrossRef] [PubMed]
  57. Plisson, F.; Huang, X.C.; Zhang, H.; Khalil, Z.; Capon, R.J. Lamellarins as Inhibitors of P-Glycoprotein-Mediated Multidrug Resistance in a Human Colon Cancer Cell Line. Chem. An Asian J. 2012, 7, 1616–1623. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, J.N.; Zhang, H.J.; Li, J.Q.; Ding, W.J.; Ma, Z.J. Bioactive Indolocarbazoles from the Marine-Derived Streptomyces Sp. DT-A61. J. Nat. Prod. 2018, 81, 949–956. [Google Scholar] [CrossRef]
  59. Xiao, F.; Li, H.; Xu, M.; Li, T.; Wang, J.; Sun, C.; Hong, K.; Li, W. Staurosporine Derivatives Generated by Pathway Engineering in a Heterologous Host and Their Cytotoxic Selectivity. J. Nat. Prod. 2018, 81, 1745–1751. [Google Scholar] [CrossRef]
  60. Rowley, D.C.; Hansen, M.S.T.; Rhodes, D.; Sotriffer, C.A.; Ni, H.; McCammon, J.A.; Bushman, F.D.; Fenical, W. Thalassiolins A-C: New Marine-Derived Inhibitors of HIV CDNA Integrase. Bioorganic Med. Chem. 2002, 10, 3619–3625. [Google Scholar] [CrossRef]
  61. Pereira, A.; Vottero, E.; Roberge, M.; Mauk, A.G.; Andersen, R.J. Indoleamine 2,3-Dioxygenase Inhibitors from the Northeastern Pacific Marine Hydroid Garveia Annulata. J. Nat. Prod. 2006, 69, 1496–1499. [Google Scholar] [CrossRef]
  62. Alonso, E.; Alvariño, R.; Leirós, M.; Tabudravu, J.N.; Feussner, K.; Dam, M.A.; Rateb, M.E.; Jaspars, M.; Botana, L.M. Evaluation of the Antioxidant Activity of the Marine Pyrroloiminoquinone Makaluvamines. Mar. Drugs 2016, 14, 197. [Google Scholar] [CrossRef]
  63. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef] [PubMed]
  64. Jorgensen, W.L.; Duffy, E.M. Prediction of Drug Solubility from Structure. Adv. Drug Deliv. Rev. 2002, 54, 355–366. [Google Scholar] [CrossRef] [PubMed]
  65. Santos, L.H.S.; Ferreira, R.S.; Caffarena, E.R. Integrating Molecular Docking and Molecular Dynamics Simulations. In Methods in Molecular Biology; Humana: New York, NY, USA, 2019; Volume 2053, pp. 13–34. ISBN 9781493997527. [Google Scholar]
  66. Shukla, R.; Tripathi, T. Molecular Dynamics Simulation of Protein and Protein–Ligand Complexes. In Computer-Aided Drug Design; Springer: Singapore, 2020; pp. 133–161. ISBN 9789811568152. [Google Scholar]
  67. Sargsyan, K.; Grauffel, C.; Lim, C. How Molecular Size Impacts RMSD Applications in Molecular Dynamics Simulations. J. Chem. Theory Comput. 2017, 13, 1518–1524. [Google Scholar] [CrossRef] [PubMed]
  68. Shoaib, T.H.; Ibraheem, W.; Abdelrahman, M.; Osman, W.; Sherif, A.E.; Ashour, A.; Ibrahim, S.R.M.; Ghazawi, K.F.; Miski, S.F.; Almadani, S.A.; et al. Exploring the Potential of Approved Drugs for Triple-Negative Breast Cancer Treatment by Targeting Casein Kinase 2: Insights from Computational Studies. PLoS ONE 2023, 18, e0289887. [Google Scholar] [CrossRef]
Figure 1. Chemical structures of the top 10 docked compounds with their CMNPD identification codes.
Figure 1. Chemical structures of the top 10 docked compounds with their CMNPD identification codes.
Metabolites 13 01090 g001
Figure 2. The 2D interactions of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). The bottom legends identify the types of interaction and bonds with their corresponding color codes.
Figure 2. The 2D interactions of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). The bottom legends identify the types of interaction and bonds with their corresponding color codes.
Metabolites 13 01090 g002aMetabolites 13 01090 g002b
Figure 3. The 3D interactions of the three docked inhibitors (in green), overlayed as they are depicted, forming hydrogen bonds (red dashed lines) with Cys83 of the CDK5 enzyme (PDB ID: 7VDQ).
Figure 3. The 3D interactions of the three docked inhibitors (in green), overlayed as they are depicted, forming hydrogen bonds (red dashed lines) with Cys83 of the CDK5 enzyme (PDB ID: 7VDQ).
Metabolites 13 01090 g003
Figure 4. RMSD plotting of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ).
Figure 4. RMSD plotting of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ).
Metabolites 13 01090 g004
Figure 5. RMSF plotting of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). Green lines represent loop regions.
Figure 5. RMSF plotting of the top inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). Green lines represent loop regions.
Metabolites 13 01090 g005
Figure 6. Protein–ligand contacts for the top three inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). The top legend depicts the color codes for the types of contacts that are formed throughout the simulation plotting of the top inhibitors and the enzyme.
Figure 6. Protein–ligand contacts for the top three inhibitors and the co-crystallized reference inhibitor with the CDK5 enzyme (PDB ID: 7VDQ). The top legend depicts the color codes for the types of contacts that are formed throughout the simulation plotting of the top inhibitors and the enzyme.
Metabolites 13 01090 g006
Figure 7. Post-MD MM/GBSA curves for the co-crystallized reference inhibitor and the top-three simulated compounds with the CDK5 target.
Figure 7. Post-MD MM/GBSA curves for the co-crystallized reference inhibitor and the top-three simulated compounds with the CDK5 target.
Metabolites 13 01090 g007
Table 2. Summary of the types of interactions that are formed by the docked inhibitors.
Table 2. Summary of the types of interactions that are formed by the docked inhibitors.
CompoundHydrogen BondsHydrophobic InteractionsOther Interactions
ResidueDistance (Å)
Excoecariphenol BGln51
Cys83
Asp144
1.93
1.70
2.39
Ile10, Tyr15, Val18, Ala31, Val64, Phe80, Phe82, Cys83, and Leu133-
Excoecariphenol AGln51
Cys83
1.83
1.68
Ile10, Ala31, Val64, Phe80, Phe82, Cys83, and Leu133-
Zyzzyanone BGln51
Glu81
Asp86
Cys83
1.85
2.25
2.11
1.94, 1.67
Ile10, Tyr15, Val18, Ala31, Phe80, Phe82, Cys83, and Leu133salt bridge-Asp86 Pi-Pi stacking-Phe80
ReferenceLys33
Gln51
Cys83
Asp86
Asp144
2.76
1.87
1.97, 1.84
2.08
1.98, 2.39
Ile10, Tyr15, Val18, Ala31, Val64, Phe80, Phe82, Cys83, and Leu133salt bridge-Asp144
Pi-cation-Tyr15
Table 3. Predicted ADME parameters for the top three compounds.
Table 3. Predicted ADME parameters for the top three compounds.
HBD aHBA bQPlog Po/w cQPlog S dQPlog
HERG e
QPP Caco fQPlog
BB g
Mwt hRule Of Five i
Excoecariphenol B1015.2−2.033−2.489−5.0251.377−4.213500.4733
Excoecariphenol A914.45−1.406−2.741−5.3054.115−3.69484.4742
Zyzzyanone B46.251.076−3.043−6.27522.892−1.397335.3620
Standard values≤5≤10−2.0–6.5−6.5˗0.5Below −5>25 poor
<500 great
−3–1.2>5000˗4
Note: (a) HBD (hydrogen bond donor). (b) HBA (hydrogen bond acceptor). (c) Predicted LogP (partition coefficient in octanol/water). (d) Predicted aqueous solubility. (e) Predicted IC50 for HERG K+ blockade. (f) Predicted Caco-2 cell permeability. (g) Predicted blood–brain partition coefficient. (h) Molecular weight. (i) Number of Lipinski’s rule of five violations.
Table 4. Predicted toxicity parameters for the top three compounds.
Table 4. Predicted toxicity parameters for the top three compounds.
CompoundOral ToxicityOrgan ToxicityToxicity Endpoints Prediction
Toxicity
Class
Predicted LD50
(mg/kg)
HepatotoxicityMutagenicityCyto
Toxicity
CarcinogenicityImmunotoxicity
Excoecariphenol B52500inactiveinactiveinactiveinactiveinactive
Excoecariphenol A52500inactiveinactiveinactiveinactiveinactive
Zyzzyanone B41000inactiveinactiveinactiveinactiveinactive
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shoaib, T.H.; Almogaddam, M.A.; Andijani, Y.S.; Saib, S.A.; Almaghrabi, N.M.; Elyas, A.F.; Azzouni, R.Y.; Awad, E.A.; Mohamed, S.G.A.; Mohamed, G.A.; et al. Marine-Derived Compounds for CDK5 Inhibition in Cancer: Integrating Multi-Stage Virtual Screening, MM/GBSA Analysis and Molecular Dynamics Investigations. Metabolites 2023, 13, 1090. https://doi.org/10.3390/metabo13101090

AMA Style

Shoaib TH, Almogaddam MA, Andijani YS, Saib SA, Almaghrabi NM, Elyas AF, Azzouni RY, Awad EA, Mohamed SGA, Mohamed GA, et al. Marine-Derived Compounds for CDK5 Inhibition in Cancer: Integrating Multi-Stage Virtual Screening, MM/GBSA Analysis and Molecular Dynamics Investigations. Metabolites. 2023; 13(10):1090. https://doi.org/10.3390/metabo13101090

Chicago/Turabian Style

Shoaib, Tagyedeen H., Mohammed A. Almogaddam, Yusra Saleh Andijani, Samaher Ahmad Saib, Najwa Mahmoud Almaghrabi, Abdulaziz Fahad Elyas, Rahmah Yasin Azzouni, Ehda Ahmad Awad, Shaimaa G. A. Mohamed, Gamal A. Mohamed, and et al. 2023. "Marine-Derived Compounds for CDK5 Inhibition in Cancer: Integrating Multi-Stage Virtual Screening, MM/GBSA Analysis and Molecular Dynamics Investigations" Metabolites 13, no. 10: 1090. https://doi.org/10.3390/metabo13101090

APA Style

Shoaib, T. H., Almogaddam, M. A., Andijani, Y. S., Saib, S. A., Almaghrabi, N. M., Elyas, A. F., Azzouni, R. Y., Awad, E. A., Mohamed, S. G. A., Mohamed, G. A., Ibrahim, S. R. M., Hussein, H. G. A., Osman, W., Ashour, A., Sherif, A. E., & Alzain, A. A. (2023). Marine-Derived Compounds for CDK5 Inhibition in Cancer: Integrating Multi-Stage Virtual Screening, MM/GBSA Analysis and Molecular Dynamics Investigations. Metabolites, 13(10), 1090. https://doi.org/10.3390/metabo13101090

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop