Integrating computational methods guided the discovery of phytochemicals as potential Pin1 inhibitors for cancer: pharmacophore modeling, molecular docking, MM-GBSA calculations and molecular dynamics studies

Pin1 is a pivotal player in interactions with a diverse array of phosphorylated proteins closely linked to critical processes such as carcinogenesis and tumor suppression. Its axial role in cancer initiation and progression, coupled with its overexpression and activation in various cancers render it a potential candidate for the development of targeted therapeutics. While several known Pin1 inhibitors possess favorable enzymatic profiles, their cellular efficacy often falls short. Consequently, the pursuit of novel Pin1 inhibitors has gained considerable attention in the field of medicinal chemistry. In this study, we employed the Phase tool from Schrödinger to construct a structure-based pharmacophore model. Subsequently, 449,008 natural products (NPs) from the SN3 database underwent screening to identify compounds sharing pharmacophoric features with the native ligand. This resulted in 650 compounds, which then underwent molecular docking and binding free energy calculations. Among them, SN0021307, SN0449787 and SN0079231 showed better docking scores with values of -9.891, -7.579 and -7.097 kcal/mol, respectively than the reference compound (-6.064 kcal/mol). Also, SN0021307, SN0449787 and SN0079231 exhibited lower free binding energies (-57.12, -49.81 and -46.05 kcal/mol, respectively) than the reference ligand (-37.75 kcal/mol). Based on these studies, SN0021307, SN0449787, and SN0079231 showed better binding affinity that the reference compound. Further the validation of these findings, molecular dynamics simulations confirmed the stability of the ligand-receptor complex for 100 ns with RMSD ranging from 0.6 to 1.8 Å. Based on these promising results, these three phytochemicals emerge as promising lead compounds warranting comprehensive biological screening in future investigations. These compounds hold great potential for further exploration regarding their efficacy and safety as Pin1 inhibitors, which could usher in new avenues for combating cancer.

Design of new Mcl-1 inhibitors for cancer using fragments hybridization, molecular docking, and molecular dynamics studies

Apoptosis is a critical process that regulates cell survival and death and plays an essential role in cancer development. The Bcl-2 protein family, including myeloid leukemia 1 (Mcl-1), is a key regulator of the intrinsic apoptosis pathway, and its overexpression in many human cancers has prompted efforts to develop Mcl-1 inhibitors as potential anticancer agents. In this study, we aimed to design new Mcl-1 inhibitors using various computational techniques. First, we used the Mcl-1 receptor-ligand complex to build an e-pharmacophore hypothesis and screened a library of 567,000 fragments from the Enamine database. We obtained 410 fragments and used them to design 92,384 novel compounds, which we then docked into the Mcl-1 binding cavity using HTVS, SP, and XP docking modes of Glide. To assess their suitability as drug candidates, we conducted MM-GBSA calculations and ADME prediction, leading to the identification of 10 compounds with excellent binding affinity and favorable pharmacokinetic properties. To further investigate the interaction strength, we performed molecular dynamics simulations on the top three Mcl-1 receptor-ligand complexes to study their interaction stability. Overall, our findings suggest that these compounds have promising potential as anticancer agents, pending further experimental validation such as Mcl-1 apoptosis Assay. By combining experimental methods with various in silico approaches, these techniques prove to be invaluable for identifying novel drug candidates with distinct therapeutic applications using fragment-based drug design. This methodology has the potential to expedite the drug discovery process while also reducing its costs.Communicated by Ramaswamy H. Sarma.

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

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.

Molecular Docking and Molecular Dynamics Studies Reveal the Anticancer Potential of Medicinal-Plant-Derived Lignans as MDM2-P53 Interaction Inhibitors

The interaction between the tumor suppressor protein p53 and its negative regulator, the MDM2 oncogenic protein, has gained significant attention in cancer drug discovery. In this study, 120 lignans reported from Ferula sinkiangensis and Justicia procumbens were assessed for docking simulations on the active pocket of the MDM2 crystal structure bound to Nutlin-3a. The docking analysis identified nine compounds with higher docking scores than the co-crystallized reference. Subsequent AMDET profiling revealed satisfactory pharmacokinetic and safety parameters for these natural products. Three compounds, namely, justin A, 6-hydroxy justicidin A, and 6'-hydroxy justicidin B, were selected for further investigation due to their strong binding affinities of -7.526 kcal/mol, -7.438 kcal/mol, and -7.240 kcal/mol, respectively, which surpassed the binding affinity of the reference inhibitor Nutlin-3a (-6.830 kcal/mol). To assess the stability and reliability of the binding of the candidate hits, a molecular dynamics simulation was performed over a duration of 100 ns. Remarkably, the thorough analysis demonstrated that all the hits exhibited stable molecular dynamics profiles. Based on their effective binding to MDM2, favorable pharmacokinetic properties, and molecular dynamics behavior, these compounds represent a promising starting point for further refinement. Nevertheless, it is essential to synthesize the suggested compounds and evaluate their activity through in vitro and in vivo experiments.

Exploring the potential of approved drugs for triple-negative breast cancer treatment by targeting casein kinase 2: Insights from computational studies

Triple-negative breast cancer (TNBC) is an aggressive malignancy that requires effective targeted drug therapy. In this study, we employed in silico methods to evaluate the efficacy of seven approved drugs against human ck2 alpha kinase, a significant modulator of TNBC metastasis and invasiveness. Molecular docking revealed that the co-crystallized reference inhibitor 108600 achieved a docking score of (-7.390 kcal/mol). Notably, among the seven approved drugs tested, sunitinib, bazedoxifene, and etravirine exhibited superior docking scores compared to the reference inhibitor. Specifically, their respective docking scores were -10.401, -7.937, and -7.743 kcal/mol. Further analysis using MM/GBSA demonstrated that these three top-ranked drugs possessed better binding energies than the reference ligand. Subsequent molecular dynamics simulations identified etravirine, an FDA-approved antiviral drug, as the only repurposed drug that demonstrated a stable and reliable binding mode with the human ck2 alpha protein, based on various analysis measures including RMSD, RMSF, and radius of gyration. Principal component analysis indicated that etravirine exhibited comparable stability of motion as a complex with human ck2 alpha protein, similar to the co-crystallized inhibitor. Additionally, Density functional theory (DFT) calculations were performed on a complex of etravirine and a representative gold atom positioned at different sites relative to the heteroatoms of etravirine. The results of the DFT calculations revealed low-energy complexes that could potentially serve as guides for experimental trials involving gold nanocarriers of etravirine, enhancing its delivery to malignant cells and introducing a new drug delivery route. Based on the results obtained in this research study, etravirine shows promise as a potential antitumor agent targeting TNBC, warranting further investigation through experimental and clinical assessments.

Modulation of NRF2/KEAP1-Mediated Oxidative Stress for Cancer Treatment by Natural Products Using Pharmacophore-Based Screening, Molecular Docking, and Molecular Dynamics Studies

Oxidative stress plays a significant role in the development of cancer. Inhibiting the protein-protein interaction (PPI) between Keap1 and Nrf2 offers a promising strategy to activate the Nrf2 antioxidant pathway, which is normally suppressed by the binding of Keap1 to Nrf2. This study aimed to identify natural compounds capable of targeting the kelch domain of KEAP1 using structure-based drug design methods. A pharmacophore model was constructed based on the KEAP1-inhibitor complex, leading to the selection of 6178 compounds that matched the model. Subsequently, docking and MM/GBSA analyses were conducted, resulting in the identification of 10 compounds with superior binding energies compared to the reference compound. From these, three compounds (ZINC000002123788, ZINC000002111341, and ZINC000002125904) were chosen for further investigation. Ligand-residue interaction analysis revealed specific interactions between these compounds and key residues, indicating their stability within the binding site. ADMET analysis confirmed that the selected compounds possessed desirable drug-like properties. Furthermore, molecular dynamics simulations were performed, demonstrating the stability of the ligand-protein complexes over a 100 ns duration. These findings underscore the potential of the selected natural compounds as agents targeting KEAP1 and provide valuable insights for future experimental studies.

Discovery of dual-target natural antimalarial agents against DHODH and PMT of Plasmodium falciparum: pharmacophore modelling, molecular docking, quantum mechanics, and molecular dynamics simulations

Malaria is a lethal disease that claims thousands of lives worldwide annually. The objective of this study was to identify new natural compounds that can target two P. falciparum enzymes; P. falciparum Dihydroorotate dehydrogenase (PfDHODH) and P. falciparum phosphoethanolamine methyltransferase (PfPMT). To accomplish this, e-pharmacophore modelling and molecular docking were employed against PfDHODH. Following this, 1201 natural compounds with docking scores of ≤ -7 kcal/mol were docked into the active site of the second enzyme PMT. The top nine compounds were subjected to further investigation using MM-GBSA free binding energy calculations and ADME analysis. The results revealed favourable free binding energy values better than the references, as well as acceptable pharmacokinetic properties. Compounds ZINC000013377887, ZINC000015113777, and ZINC000085595753 were scrutinized to assess their interaction stability with the PfDHODH enzyme, and chemical stability reactivity using molecular dynamics (MD) simulation and density functional theory (DFT) calculations. These findings indicate that the three natural compounds are potential candidates for dual PfDHODH and PfPMT inhibitors for malaria treatment.