Multi-Phase In Silico Discovery of Potential SARS-CoV-2 RNA-Dependent RNA Polymerase Inhibitors among 3009 Clinical and FDA-Approved Related Drugs

Repurposing FDA-approved drugs for COVID-19: targeting the main protease through multi-phase in silico approach

BACKGROUND

The COVID-19 pandemic has created an urgent need for effective therapeutic agents. The SARS-CoV-2 Main Protease (Mpro) plays a crucial role in viral replication and immune evasion, making it a key target for drug development. While several studies have explored Mpro inhibition, identifying FDA-approved drugs with potential efficacy remains a critical research focus.

PURPOSE

This study aims to identify FDA-approved drugs that could inhibit SARS-CoV-2 Mpro. Using computational screening, we seek compounds that share structural similarities with a known co-crystallized ligand (PRD_002214) and exhibit strong binding affinity to the enzyme, providing viable candidates for COVID-19 treatment.

RESEARCH DESIGN

A systematic in silico approach was used, screening 3009 FDA-approved drugs. The initial screening focused on structural similarity to PRD_002214 (PDB ID: 6LU7), followed by molecular docking studies to predict binding affinity. Promising compounds were further analyzed through molecular dynamics (MD) simulations to evaluate their stability and interactions with Mpro over 100 ns.

STUDY SAMPLE

Of the 3009 FDA-approved drugs screened, 74 were selected for initial evaluation. After refinement, 28 compounds underwent docking analysis, with eight showing strong binding potential to Mpro.

ANALYSIS

Molecular docking assessed the binding affinity and interaction of the selected compounds with Mpro. MD simulations were conducted on the top compound, Atazanavir, to study its dynamic interactions. MM-GBSA, PLIP, and PCAT analyses were used to validate binding affinity and interactions.

RESULTS

Eight compounds, including Carfilzomib, Atazanavir, Darunavir, and others, exhibited promising binding affinities. Among them, Atazanavir showed the highest binding strength and was selected for further MD simulation studies. These simulations revealed that Atazanavir forms stable interactions with Mpro, demonstrating favorable binding and dynamic stability. The binding affinity was further confirmed through MM-GBSA, PLIP, and PCAT analyses, supporting Atazanavir's potential as an effective Mpro inhibitor.

CONCLUSIONS

In silico results suggest that Atazanavir is a promising candidate for targeting SARS-CoV-2 Mpro, with strong binding affinity and dynamic stability. These findings support its potential as a lead compound for further preclinical and clinical testing, though in vitro and in vivo validation are needed to confirm its therapeutic efficacy against COVID-19.

New apoptotic anti-triple-negative breast cancer theobromine derivative inhibiting EGFRWT and EGFRT790M: in silico and in vitro evaluation

A new theobromine-derived EGFR inhibitor (2-(3,7-Dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-1-yl)-N-(2,6-dimethylphenyl)acetamide) has been developed that has the essential structural characteristics to interact with EGFR's pocket. The designed compound is 2,6-di ortho methylphenyl)acetamide derivative of the well-known alkaloid, theobromine, (T-1-DOMPA). Firstly, deep DFT studies have been conducted to study the optimized chemical structure, molecular orbital and chemical reactivity analysis of T-1-DOMPA. Then, T-1-DOMPA's anticancer potentialities were estimated first through a structure-based computational approach. Utilizing molecular docking, molecular dynamics, MD, simulations over 100 ns, MM-PBSA and PLIP studies, T-1-DOMPA bonded to and inhibited the EGFR protein effectively. Subsequently, the ADMET profiles of T-1-DOMPA were computed before preparation, and its drug-likeness was anticipated. Therefore, T-1-DOMPA was prepared for the purposes of scrutinizing both the design and the results obtained in silico. The in vitro potential of T-1-DOMPA against triple-negative breast cancer cell lines, MDA- MB-231, was very promising with an IC50 value of1.8 µM, comparable to the reference drug (0.9 µM), and a much higher selectivity index of 2.6. Interestingly, T-1-DOMPA inhibited three other cancer cell lines (CaCO-2, HepG-2, and A549) with IC50 values of 1.98, 2.53, and 2.39 µM exhibiting selectivity index values of 2,4, 1.9, and 2, respectively. Additionally, T-1-DOMPA prevented effectively the MDA-MB-231cell line's healing and migration abilities. Also, T-1-DOMPA's abilities to induce apoptosis were confirmed by acridine orange/ethidium bromide (AO/EB) staining assay. Finally, T-1-DOMPA caused an up-regulation of the gene expression of the apoptotic gene, Caspase-3, in the treated MDA-MB-231cell.

Computer aided drug discovery (CADD) of a thieno[2,3-d]pyrimidine derivative as a new EGFR inhibitor targeting the ribose pocket

Depending on the pharmacophoric characteristics of EGFR inhibitors, a new thieno[2,3-d]pyrimidine derivative has been developed. Firstly, the potential inhibitory effect of the designed compound against EGFR has been proven by docking experiments that showed correct binding modes and excellent binding energies of -98.44 and -88.00 kcal/mol, against EGFR wild-type and mutant type, respectively. Furthermore, MD simulations studies confirmed the precise energetic, conformational, and dynamic alterations that occurred after binding to EGFR. The correct binding was also confirmed by essential dynamics studies. To further investigate the general drug-like properties of the developed candidate, in silico ADME and toxicity studies have also been carried out. The thieno[2,3-d]pyrimidine derivative was synthesized following the earlier promising findings. Fascinatingly, the synthesized compound (4) showed promising inhibitory effects against EGFRWT and EGFRT790M with IC50 values of 25.8 and 182.3 nM, respectively. Also, it exhibited anticancer potentialities against A549 and MCF-7cell lines with IC50 values of 13.06 and 20.13 µM, respectively. Interestingly, these strong activities were combined with selectivity indices of 2.8 and 1.8 against the two cancer cell lines, respectively. Further investigations indicated the ability of compound 4 to arrest the cancer cells' growth at the G2/M phase and to increase early and late apoptosis percentages from 2.52% and 2.80 to 17.99% and 16.72%, respectively. Additionally, it was observed that compound 4 markedly increased the levels of caspase-3 and caspase-9 by 4 and 3-fold compared to the control cells. Moreover, it up-regulated the level of BAX by 3-fold and down-regulated the level of Bcl-2 by 3-fold affording a BAX/Bcl-2 ratio of 9.Communicated by Ramaswamy H. Sarma.

In silico and in vitro evaluation of the anti-virulence potential of patuletin, a natural methoxy flavone, against Pseudomonas aeruginosa

This study aimed to investigate the potential of patuletin, a rare natural flavonoid, as a virulence and LasR inhibitor against Pseudomonas aeruginosa. Various computational studies were utilized to explore the binding of Patuletin and LasR at a molecular level. Molecular docking revealed that Patuletin strongly interacted with the active pocket of LasR, with a high binding affinity value of -20.96 kcal/mol. Further molecular dynamics simulations, molecular mechanics generalized Born surface area (MM/GBSA), protein-ligand interaction profile (PLIP), and essential dynamics analyses confirmed the stability of the patuletin-LasR complex, and no significant structural changes were observed in the LasR protein upon binding. Key amino acids involved in binding were identified, along with a free energy value of -26.9 kcal/mol. In vitro assays were performed to assess patuletin's effects on P. aeruginosa. At a sub-inhibitory concentration (1/4 MIC), patuletin significantly reduced biofilm formation by 48% and 42%, decreased pyocyanin production by 24% and 14%, and decreased proteolytic activities by 42% and 20% in P. aeruginosa isolate ATCC 27853 (PA27853) and P. aeruginosa clinical isolate (PA1), respectively. In summary, this study demonstrated that patuletin effectively inhibited LasR activity in silico and attenuated virulence factors in vitro, including biofilm formation, pyocyanin production, and proteolytic activity. These findings suggest that patuletin holds promise as a potential therapeutic agent in combination with antibiotics to combat antibiotic-tolerant P. aeruginosa infections.

Isolation and in silico investigation of cannflavins from Cannabis sativa leaves as potential anti-SARS-CoV-2 agents targeting the Papain-Like Protease

This study aimed to isolate and identify three prenylflavonoids (cannflavin A, B, and C) from Cannabis sativa leaves using different chromatographic techniques. The potential of the isolated compounds against SARS-CoV-2 was suggested through several in silico analysis. Structural similarity studies against nine co-crystallized ligands of SARS-CoV-2's proteins indicated the similarities of the isolated cannflavins with the SARS-CoV-2 Papain-Like Protease (PLP) ligand, Y95. Then, flexible allignment study confirmed this similarity. Docking experiments showed successful binding of all cannflavins within the active pocket of PLP, with energies comparable to Y95. Among them, cannflavin A demonstrated the most similar binding mode, while cannflavin C exhibited the best energy. Molecular dynamics (MD) simulations and MM-GPSA confirmed the accurate binding of cannflavin A to the PLP. In silico ADMET studies indicated favourable drug-like properties for all three compounds, suggesting their potential as anti-SARS-CoV-2 agents. Further In vitro and In vivo investigations are necessary to validate these findings and establish their efficacy and safety profiles.

Computer-assisted drug discovery of potential natural inhibitors of the SARS-CoV-2 RNA-dependent RNA polymerase through a multi-phase in silico approach

BACKGROUND

The COVID-19 pandemic has led to significant loss of life and economic disruption worldwide. Currently, there are limited effective treatments available for this disease. SARS-CoV-2 RNA-dependent RNA polymerase (SARS-CoV-2 RdRp) has been identified as a potential target for drug development against COVID-19. Natural products have been shown to possess antiviral properties, making them a promising source for developing drugs against SARS-CoV-2.

OBJECTIVES

The objective of this study is to identify the most effective natural inhibitors of SARS-CoV-2 RdRp among a set of 4924 African natural products using a multi-phase in silico approach.

METHODS

The study utilized remdesivir (RTP), the co-crystallized ligand of RdRp, as a starting point to select compounds that have the most similar chemical structures among the examined set of compounds. Molecular fingerprints and structure similarity studies were carried out in the first part of the study. The second part of the study included molecular docking against SARS-CoV-2 RdRp (PDB ID: 7BV2) and Molecular Dynamics (MD) simulations including the calculation of RMSD, RMSF, Rg, SASA, hydrogen bonding, and PLIP. Moreover, the calculations of Molecular mechanics with generalised Born and surface area solvation (MM-GBSA) Lennard-Jones and Columbic electrostatic interaction energies have been conducted. Additionally, in silico ADMET and toxicity studies were performed to examine the drug likeness degrees of the selected compounds.

RESULTS

Eight compounds were identified as the most effective natural inhibitors of SARS-CoV-2 RdRp. These compounds are kaempferol 3-galactoside, kaempferol 3-O-β-D-glucopyranoside, mangiferin methyl ether, luteolin 7-O-β-D-glucopyranoside, quercetin-O-β-D-3-glucopyranoside, 1-methoxy-3-indolylmethyl glucosinolate, naringenin, and asphodelin A 4'-O-β-D-glucopyranoside.

CONCLUSION

The results of this study provide valuable information for the development of natural product-based drugs against COVID-19. However, the elected compounds should be further studied in vitro and in vivo to confirm their efficacy in treating COVID-19.

Synthesis, biological evaluation and computer-aided discovery of new thiazolidine-2,4-dione derivatives as potential antitumor VEGFR-2 inhibitors

In this study, novel VEGFR-2-targeting thiazolidine-2,4-dione derivatives with potential anticancer properties were designed and synthesized. The ability of the designed derivatives to inhibit VEGFR-2 and stop the growth of three different cancer cell types (HT-29, A-549, and HCT-116) was examined in vitro. The IC50 value of compound 15, 0.081 μM, demonstrated the best anti-VEGFR-2 potency. Additionally, compound 15 showed remarkable anti-proliferative activities against the tested cancer cell lines, with IC50 values ranging from 13.56 to 17.8 μM. Additional flow cytometric investigations showed that compound 15 increased apoptosis in HT-29 cancer cells (from 3.1% to 31.4%) arresting their growth in the S phase. Furthermore, compound 15's apoptosis induction in the same cell line was confirmed by increasing the levels of BAX (4.8-fold) and decreasing Bcl-2 (2.8-fold). Also, compound 15 noticeably increased caspase-8 and caspase-9 levels by 1.7 and 3.2-fold, respectively. Computational methods were used to perform molecular analysis of the VEGFR-2-15 complex. Molecular dynamics simulations and molecular docking were utilized to analyze the complex's kinetic and structural characteristics. Protein-ligand interaction profiler analysis (PLIP) determined the 3D interactions and binding conformation of the VEGFR-2-15 complex. DFT analyses also provided insights into the 3D geometry, reactivity, and electronic characteristics of compound 15. Computational ADMET and toxicity experiments were conducted to determine the potential of the synthesized compounds for therapeutic development. The study's findings suggest that compound 15 might be an effective anticancer lead compound and could guide future attempts to develop new drugs.

A novel thieno[2,3-d]pyrimidine derivative inhibiting vascular endothelial growth factor receptor-2: A story of computer-aided drug discovery

Following the pharmacophoric features of vascular endothelial growth factor receptor 2 (VEGFR-2) inhibitors, a novel thieno[2,3-d]pyrimidine derivative has been designed and its activity against VEGFR-2 has been demonstrated by molecular docking studies that showed an accurate binding mode and an excellent binding energy. Furthermore, the recorded binding was confirmed by a series of molecular dynamics simulation studies, which also revealed precise energetic, conformational, and dynamic changes. Additionally, molecular mechanics with generalized Born and surface area solvation and polymer-induced liquid precursors studies were conducted and verified the results of the MD simulations. Next, in silico absorption, distribution, metabolism, excretion, and toxicity studies have also been conducted to examine the general drug-like nature of the designed candidate. According to the previous results, the thieno[2,3-d]pyrimidine derivative was synthesized. Fascinatingly, it inhibited VEGFR-2 (IC50  = 68.13 nM) and demonstrated strong inhibitory activity toward human liver (HepG2), and prostate (PC3) cell lines with IC50 values of 6.60 and 11.25 µM, respectively. As well, it was safe and showed a high selectivity index against normal cell lines (WI-38). Finally, the thieno[2,3-d]pyrimidine derivative arrested the growth of the HepG2 cells at the G2/M phase inducing both early and late apoptosis. These results were further confirmed through the ability of the thieno[2,3-d]pyrimidine derivative to induce significant changes in the apoptotic genes levels of caspase-3, caspase-9, Bcl-2 associated X-protein, and B-cell lymphoma 2.

In silico, in vitro VEGFR-2 inhibition, and anticancer activity of a 3-(hydrazonomethyl)naphthalene-2-ol derivative

In agreement with the general features of VEGFR-2 inhibitors, a new naphthalene analog (compound 7) has been designed and synthesized. The inhibitory potential of compound 7 was indicated by the proper binding and the perfect energy of -21.10 kcal/mol compared to sorafenib (-21.22) in the molecular docking studies. Next, six MD simulation studies over 100 ns (RMSD, RMSF, SASA, RoG, hydrogen bonding, and distance between the center of mass) confirmed the accurate interaction of compound 7 with the catalytic pocket of VEGFR-2. Similarly, an MM-GBSA established proper binding showing an exact total binding energy of -36.95 ± 3.03 kcal/Mol. Additionally, the MM-GBSA experiment indicated the vital amino acids in the binding process. Types and number of interactions of compound 7 with catalytic pocket of VEGFR-2 were determined through Protein-Ligand Interaction Profiler (PLIP). As a new compound, the DFT was employed to optimize the molecular structure of compound 7. The DFT experiments also verified the interaction features of compound 7 with the VEGFR-2 active site. In silico ADMET experiments revealed the general drug-likeness of compound 7. Fascinatingly, the in vitro examinations were consistent with the in silico experiments as compound 7 inhibited the VEGFR-2 enzyme with an IC50 value of 37 nM. Captivatingly, compound 7 inhibited both MCF-7 and HCT 116 cancer cells exhibiting IC50 values of 10.56 and 7.07 µM exhibiting excellent selectivity indexes of 9.04 and 13.50, respectively.Communicated by Ramaswamy H. Sarma.

Discovery of new quinoline and isatine derivatives as potential VEGFR-2 inhibitors: design, synthesis, antiproliferative, docking and MD simulation studies

A new set of quinoline and isatine derivatives were synthesized as antiangiogenic VEGFR-2 inhibitors. On a biological level, the in vitro ability of the obtained candidates to inhibit VEGFR-2 was found to be strong with IC50 values in the range of 76.64-175.50 nM. To investigate the cytotoxicity and safety, all compounds were tested against a panel of four cancer cell lines (A549, Caco2, HepG2 and MDA) as well as two normal cell lines (Vero and WI-38). Interestingly, compound 12 exhibited noticeable cytotoxicity against A549, Caco2 and MDA with IC50 values of 5.40, 0.58 and 0.94 µM, respectively. These results were better and comparable to that of doxorubicin (0.70, 0.82 and 0.90 µM, respectively) with more than three folds higher selectivity index against the Caco2 cell lines. Compound 9 prevented the healing of the cancer cells at a low concentration. Also, the compound's potential to induce programmed cell death in Caco-2 was proved through the significant down regulating of the expression of Bcl2, Bcl-xl and Survivin in addition to the slight upregulation of the TGF-β gene. The cell cycle analysis indicated that compound 9 arrested the Caco-2 cells in the G2/M phase. Interestingly, the molecular docking studies against VEGFR-2 revealed the correct binding of the targeted compounds similar to sorafenib. Furthermore, MD experiments validated the binding of compound 12 with VEGFR-2 over 100 ns, as well as MM-PBSA analysis that confirmed the precise binding with optimum energy. Finally, ADMET analysis showed the general drug-likeness and confirmed the safety of the tested compounds.Communicated by Ramaswamy H. Sarma.

New2-((2-(2,4-dinitrophenyl)hydrazineeylidene) derivatives: design, synthesis, in silico, and in vitro anticancer studies

A series of novel hydrazone compounds have been synthesized by the condensation of hydrazines and different substituted salicylaldehydes at a molar ratio of 1:1 in one step reaction and characterized by FT-IR, ESI-MS, 1H NMR, and single crystal x-ray diffraction. The crystal structure of the compound shows a trans configuration around the C = N bond and triclinic system with P -1/-p 1. Synthesized compounds were screened for cytotoxicity activities against A375 (melanoma), HT-29 (Colon), and A549 (lung) cancer cell lines. Among them, compound 2 exhibited the highest cytotoxic effect against the A375 cell line (IC50 = 0.30 µM) and HT-29 cell line (1.68 µM), compared to those of apatinib as a reference standard drug (0.28, 1.49 µM, respectively). The cytocompatibility assay on the L929 normal cell line and the hemolysis assay on human RBC were used to validate the non-toxic action. From DFT calculation, the various parameters such as HOMO-LUMO energies, Hirshfeld, and MEP have been studied. Furthermore, in silico molecular docking with three receptors was studied. Among four compounds, compound 2 has the lowest binding energy against cyclin dependent kinase (ΔGb = -9.3 kcal/mol). In addition to this, molecular dynamics (MD) simulation was also performed. Based on this study, these novel hydrazones can be considered a promising anticancer agent due to their potent cytotoxicity activities and computational analysis.Communicated by Ramaswamy H. Sarma.

Anticancer derivative of the natural alkaloid, theobromine, inhibiting EGFR protein: Computer-aided drug discovery approach

A new semisynthetic derivative of the natural alkaloid, theobromine, has been designed as a lead antiangiogenic compound targeting the EGFR protein. The designed compound is an (m-tolyl)acetamide theobromine derivative, (T-1-MTA). Molecular Docking studies have shown a great potential for T-1-MTA to bind to EGFR. MD studies (100 ns) verified the proposed binding. By MM-GBSA analysis, the exact binding with optimal energy of T-1-MTA was also identified. Then, DFT calculations were performed to identify the stability, reactivity, electrostatic potential, and total electron density of T-1-MTA. Furthermore, ADMET analysis indicated the T-1-MTA's general likeness and safety. Accordingly, T-1-MTA has been synthesized to be examined in vitro. Intriguingly, T-1-MTA inhibited the EGFR protein with an IC50 value of 22.89 nM and demonstrated cytotoxic activities against the two cancer cell lines, A549, and HCT-116, with IC50 values of 22.49, and 24.97 μM, respectively. Interestingly, T-1-MTA's IC50 against the normal cell lines, WI-38, was very high (55.14 μM) indicating high selectivity degrees of 2.4 and 2.2, respectively. Furthermore, the flow cytometry analysis of A549 treated with T-1-MTA showed significantly increased ratios of early apoptosis (from 0.07% to 21.24%) as well as late apoptosis (from 0.73% to 37.97%).

Design and synthesis of thiazolidine-2,4-diones hybrids with 1,2-dihydroquinolones and 2-oxindoles as potential VEGFR-2 inhibitors: in-vitro anticancer evaluation and in-silico studies

A thiazolidine-2,4-dione nucleus was molecularly hybridised with the effective antitumor moieties; 2-oxo-1,2-dihydroquinoline and 2-oxoindoline to obtain new hybrids with potential activity against VEGFR-2. The cytotoxic effects of the synthesised derivatives against Caco-2, HepG-2, and MDA-MB-231 cell lines were investigated. Compound 12a was found to be the most potent candidate against the investigated cell lines with IC₅₀ values of 2, 10, and 40 µM, respectively. Furthermore, the synthesised derivatives were tested in vitro for their VEGFR-2 inhibitory activity showing strong inhibition. Moreover, an in vitro viability study against Vero non-cancerous cell line was investigated and the results reflected a high safety profile of all tested compounds. Compound 12a was further investigated for its apoptotic behaviour by assessing the gene expression of four genes (Bcl2, Bcl-xl, TGF, and Survivin). Molecular dynamic simulations authenticated the high affinity, accurate binding, and perfect dynamics of compound 12a against VEGFR-2.

Benzoxazole derivatives as new VEGFR-2 inhibitors and apoptosis inducers: design, synthesis, in silico studies, and antiproliferative evaluation

In this study, a set of novel benzoxazole derivatives were designed, synthesised, and biologically evaluated as potential VEGFR-2 inhibitors. Five compounds (12d, 12f, 12i, 12l, and 13a) displayed high growth inhibitory activities against HepG2 and MCF-7 cell lines and were further investigated for their VEGFR-2 inhibitory activities. The most potent anti-proliferative member 12 l (IC₅₀ = 10.50 μM and 15.21 μM against HepG2 and MCF-7, respectively) had the most promising VEGFR-2 inhibitory activity (IC₅₀ = 97.38 nM). A further biological evaluation revealed that compound 12l could arrest the HepG2 cell growth mainly at the Pre-G1 and G1 phases. Furthermore, compound 12l could induce apoptosis in HepG2 cells by 35.13%. likely, compound 12l exhibited a significant elevation in caspase-3 level (2.98-fold) and BAX (3.40-fold), and a significant reduction in Bcl-2 level (2.12-fold). Finally, docking studies indicated that 12l exhibited interactions with the key amino acids in a similar way to sorafenib.

The Discovery of Potential SARS-CoV-2 Natural Inhibitors among 4924 African Metabolites Targeting the Papain-like Protease: A Multi-Phase In Silico Approach

Four compounds, hippacine, 4,2'-dihydroxy-4'-methoxychalcone, 2',5'-dihydroxy-4-methoxychalcone, and wighteone, were selected from 4924 African natural metabolites as potential inhibitors against SARS-CoV-2 papain-like protease (PLpro, PDB ID: 3E9S). A multi-phased in silico approach was employed to select the most similar metabolites to the co-crystallized ligand (TTT) of the PLpro through molecular fingerprints and structural similarity studies. Followingly, to examine the binding of the selected metabolites with the PLpro (molecular docking. Further, to confirm this binding through molecular dynamics simulations. Finally, in silico ADMET and toxicity studies were carried out to prefer the most convenient compounds and their drug-likeness. The obtained results could be a weapon in the battle against COVID-19 via more in vitro and in vivo studies.

Design, synthesis, anti-proliferative evaluation, docking, and MD simulations studies of new thiazolidine-2,4-diones targeting VEGFR-2 and apoptosis pathway

We report herein, the design and synthesis of thiazolidine-2,4-diones derivatives as new inhibitors for VEGFR-2. The designed members were assessed for their in vitro anticancer activity against four cancer cell lines; A549, Caco-2, HepG-2 and MDA-MB-231. Compound 14a showed the most potent effects against Caco-2, and HepG-2 cell lines (IC50 = of 1.5 and 31.5 μM, respectively). Next, the in vitro VEGFR-2 inhibitory activity, safety profiles and selectivity indices were examined for all the synthesized members against the normal Vero cell line. Compound 14a (the safest member against Caco-2 cell line) was further investigated for its ability to inhibit Caco-2 cells migration and healing. Moreover, the apoptotic induction of compound 14a against Caco-2 cell line was investigated by assessing against four apoptotic genes (Bcl2, Bcl-xl, TGF, and Survivin). The results revealed that compound 14a can exert apoptosis through significant reduction of Bcl2, Survivin, and TGF gene expression levels. Finally, deep computational studies including molecular docking, ADMET, toxicity studies, and MD simulation were carried out. Also, the DFT calculations were performed and discussed, and the results confirmed the inhibitory reactivity of 14a against VEGFR-2. Compound 14a is expected to be used as a potential lead in the development of new VEGFR-2 inhibitors with increased potency.

New Anticancer Theobromine Derivative Targeting EGFRWT and EGFRT790M: Design, Semi-Synthesis, In Silico, and In Vitro Anticancer Studies

Based on the pharmacophoric features of EGFR inhibitors, a new semisynthetic theobromine-derived compound was designed to interact with the catalytic pocket of EGFR. Molecular docking against wild (EGFR^WT; PDB: 4HJO) and mutant (EGFR^T790M; PDB: 3W2O) types of EGFR-TK indicated that the designed theobromine derivative had the potential to bind to that pocket as an antiangiogenic inhibitor. The MD and MM-GBSA experiments identified the exact binding with optimum energy and dynamics. Additionally, the DFT calculations studied electrostatic potential, stability, and total electron density of the designed theobromine derivative. Both in silico ADMET and toxicity analyses demonstrated its general likeness and safety. We synthesized the designed theobromine derivative (compound XI) which showed an IC₅₀ value of 17.23 nM for EGFR inhibition besides IC₅₀ values of 21.99 and 22.02 µM for its cytotoxicity against A549 and HCT-116 cell lines, respectively. Interestingly, compound XI expressed a weak cytotoxic potential against the healthy W138 cell line (IC₅₀ = 49.44 µM, 1.6 times safer than erlotinib), exhibiting the high selectivity index of 2.2. Compound XI arrested the growth of A549 at the G2/M stage and increased the incidence of apoptosis.

Discovery of Potential SARS-CoV-2 Papain-like Protease Natural Inhibitors Employing a Multi-Phase In Silico Approach

As an extension of our research against COVID-19, a multiphase in silico approach was applied in the selection of the three most common inhibitors (Glycyrrhizoflavone (76), Arctigenin (94), and Thiangazole (298)) against papain-like protease, PLpro (PDB ID: 4OW0), among 310 metabolites of natural origin. All compounds of the exam set were reported as antivirals. The structural similarity between the examined compound set and S88, the co-crystallized ligand of PLpro, was examined through structural similarity and fingerprint studies. The two experiments pointed to Brevicollin (28), Cryptopleurine (41), Columbamine (46), Palmatine (47), Glycyrrhizoflavone (76), Licochalcone A (87), Arctigenin (94), Termilignan (98), Anolignan B (99), 4,5-dihydroxy-6″-deoxybromotopsentin (192), Dercitin (193), Tryptanthrin (200), 6-Cyano-5-methoxy-12-methylindolo [2, 3A] carbazole (211), Thiangazole (298), and Phenoxan (300). The binding ability against PLpro was screened through molecular docking, disclosing the favorable binding modes of six metabolites. ADMET studies expected molecules 28, 76, 94, 200, and 298 as the most favorable metabolites. Then, molecules 76, 94, and 298 were chosen through in silico toxicity studies. Finally, DFT studies were carried out on glycyrrhizoflavone (76) and indicated a high level of similarity in the molecular orbital analysis. The obtained data can be used in further in vitro and in vivo studies to examine and confirm the inhibitory effect of the filtered metabolites against PLpro and SARS-CoV-2.

Discovery of Some Heterocyclic Molecules as Bone Morphogenetic Protein 2 (BMP-2)-Inducible Kinase Inhibitors: Virtual Screening, ADME Properties, and Molecular Docking Simulations

Bone morphogenetic proteins (BMPs) are growth factors that have a vital role in the production of bone, cartilage, ligaments, and tendons. Tumors’ upregulation of bone morphogenetic proteins (BMPs) and their receptors are key features of cancer progression. Regulation of the BMP kinase system is a new promising strategy for the development of anti-cancer drugs. In this work, based on a careful literature study, a library of benzothiophene and benzofuran derivatives was subjected to different computational techniques to study the effect of chemical structure changes on the ability of these two scaffolds to target BMP-2 inducible kinase, and to reach promising candidates with proposed activity against BMP-2 inducible kinase. The results of screening against Lipinski’s and Veber’s Rules produced twenty-one outside eighty-four compounds having drug-like molecular nature. Computational ADMET studies favored ten compounds (11, 26, 27, 29, 30, 31, 34, 35, 65, and 72) with good pharmacokinetic profile. Computational toxicity studies excluded compound 34 to elect nine compounds for molecular docking studies which displayed eight compounds (26, 27, 29, 30, 31, 35, 65, and 72) as promising BMP-2 inducible kinase inhibitors. The nine fascinating compounds will be subjected to extensive screening against serine/threonine kinases to explore their potential against these critical proteins. These promising candidates based on benzothiophene and benzofuran scaffolds deserve further clinical investigation as BMP-2 kinase inhibitors for the treatment of cancer.

Structure-Based Virtual Screening, Docking, ADMET, Molecular Dynamics, and MM-PBSA Calculations for the Discovery of Potential Natural SARS-CoV-2 Helicase Inhibitors from the Traditional Chinese Medicine

Continuing our antecedent work against COVID-19, a set of 5956 compounds of traditional Chinese medicine have been virtually screened for their potential against SARS-CoV-2 helicase (PDB ID: 5RMM). Initially, a fingerprint study with VXG, the ligand of the target enzyme, disclosed the similarity of 187 compounds. Then, a molecular similarity study declared the most similar 40 compounds. Subsequently, molecular docking studies were carried out to examine the binding modes and energies. Then, the most appropriate 26 compounds were subjected to in silico ADMET and toxicity studies to select the most convenient inhibitors to be: (1R,2S)-ephedrine (57), (1R,2S)-norephedrine (59), 2-(4-(pyrrolidin-1-yl)phenyl)acetic acid (84), 1-phenylpropane-1,2-dione (195), 2-methoxycinnamic acid (246), 2-methoxybenzoic acid (364), (R)-2-((R)-5-oxopyrrolidin-3-yl)-2-phenylacetic acid (405), (Z)-6-(3-hydroxy-4-methoxystyryl)-4-methoxy-2H-pyran-2-one (533), 8-chloro-2-(2-phenylethyl)-5,6,7-trihydroxy-5,6,7,8-tetrahydrochromone (637), 3-((1R,2S)-2-(dimethylamino)-1-hydroxypropyl)phenol (818), (R)-2-ethyl-4-(1-hydroxy-2-(methylamino)ethyl)phenol (5159), and (R)-2-((1S,2S,5S)-2-benzyl-5-hydroxy-4-methylcyclohex-3-en-1-yl)propane-1,2-diol (5168). Among the selected 12 compounds, the metabolites, compound 533 showed the best docking scores. Interestingly, the MD simulation studies for compound 533, the one with the highest docking score, over 100 ns showed its correct binding to SARS-CoV-2 helicase with low energy and optimum dynamics. Finally, MM-PBSA studies showed that 533 bonded favorably to SARS-CoV-2 helicase with a free energy value of −83 kJ/mol. Further, the free energy decomposition study determined the essential amino acid residues that contributed favorably to the binding process. The obtained results give a huge hope to find a cure for COVID-19 through further in vitro and in vivo studies for the selected compounds.

Modified Benzoxazole-Based VEGFR-2 Inhibitors and Apoptosis Inducers: Design, Synthesis, and Anti-Proliferative Evaluation

This work is one of our efforts to discover potent anticancer agents. We modified the most promising derivative of our previous work concerned with the development of VEGFR-2 inhibitor candidates. Thirteen new compounds based on benzoxazole moiety were synthesized and evaluated against three human cancer cell lines, namely, breast cancer (MCF-7), colorectal carcinoma (HCT116), and hepatocellular carcinoma (HepG2). The synthesized compounds were also evaluated against VEGFR-2 kinase activity. The biological testing fallouts showed that compound 8d was more potent than standard sorafenib. Such compound showed IC50 values of 3.43, 2.79, and 2.43 µM against the aforementioned cancer cell lines, respectively, compared to IC50 values of 4.21, 5.30, and 3.40 µM reported for sorafenib. Compound 8d also was found to exert exceptional VEGFR-2 inhibition activity with an IC50 value of 0.0554 μM compared to sorafenib (0.0782 μM). In addition, compound 8h revealed excellent cytotoxic effects with IC50 values of 3.53, 2.94, and 2.76 µM against experienced cell lines, respectively. Furthermore, compounds 8a and 8e were found to inhibit VEGFR-2 kinase activity with IC50 values of 0.0579 and 0.0741 μM, exceeding that of sorafenib. Compound 8d showed a significant apoptotic effect and arrested the HepG2 cells at the pre-G1 phase. In addition, it exerted a significant inhibition for TNF-α (90.54%) and of IL-6 (92.19%) compared to dexamethasone (93.15%). The molecular docking studies showed that the binding pattern of the new compounds to VEGFR-2 kinase was similar to that of sorafenib.

A Multistage In Silico Study of Natural Potential Inhibitors Targeting SARS-CoV-2 Main Protease

Among a group of 310 natural antiviral natural metabolites, our team identified three compounds as the most potent natural inhibitors against the SARS-CoV-2 main protease (PDB ID: 5R84), Mpro. The identified compounds are sattazolin and caprolactin A and B. A validated multistage in silico study was conducted using several techniques. First, the molecular structures of the selected metabolites were compared with that of GWS, the co-crystallized ligand of Mpro, in a structural similarity study. The aim of this study was to determine the thirty most similar metabolites (10%) that may bind to the Mpro similar to GWS. Then, molecular docking against Mpro and pharmacophore studies led to the choice of five metabolites that exhibited good binding modes against the Mpro and good fit values against the generated pharmacophore model. Among them, three metabolites were chosen according to ADMET studies. The most promising Mpro inhibitor was determined by toxicity and DFT studies to be caprolactin A (292). Finally, molecular dynamics (MD) simulation studies were performed for caprolactin A to confirm the obtained results and understand the thermodynamic characteristics of the binding. It is hoped that the accomplished results could represent a positive step in the battle against COVID-19 through further in vitro and in vivo studies on the selected compounds.

The Computational Preventive Potential of the Rare Flavonoid, Patuletin, Isolated from Tagetes patula, against SARS-CoV-2

The rare flavonoid, patuletin, was isolated from the flowers of Tagetes patula growing in Egypt. The rarity of the isolated compound inspired us to scrutinize its preventive effect against COVID-19 utilizing a multi-step computational approach. Firstly, a structural similarity study was carried out against nine ligands of nine SARS-CoV-2 proteins. The results showed a large structural similarity between patuletin and F86, the ligand of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp). Then, a 3D-Flexible alignment study of patuletin and F86 verified the proposed similarity. To determine the binding opportunity, patuletin was docked against the RdRp showing a correct binding inside its active pocket with an energy of −20 kcal/mol that was comparable to that of F86 (−23 kcal/mol). Following, several MD simulations as well as MM-PBSA studies authenticated the accurate binding of patuletin in the RdRp via the correct dynamic and energetic behaviors over 100 ns. Additionally, in silico ADMET studies showed the general safety and drug-likeness of patuletin.

The Assessment of Anticancer and VEGFR-2 Inhibitory Activities of a New 1H-Indole Derivative: In Silico and In Vitro Approaches

Corresponding to the reported features of anti-VEGFR-2-approved compounds, a new 1H-indole derivative (compound 7) was designed. The inhibitory potential of the designed compound was revealed via a molecular docking study that showed the appropriate binding. Then, MD simulation (six studies) over a period of 100 ns was performed to confirm the precise binding and optimum energy. Additionally, MM-GBSA reaffirmed the perfect binding, exhibiting a total precise energy of −40.38 Kcal/Mol. The MM-GBSA experiments named the essential amino acids in the protein–ligand interaction, employing the binding energy decomposition and revealing the diversity of interactions of compound 7 inside the VEGFR-2 enzyme. As compound 7 is new, DFT experiments were utilized for molecular structure optimization. Additionally, the DFT results validated the coherent interaction of compound 7 with the VEGFR-2 enzyme. A good value of drug-likeness of compound 7 was acknowledged via in silico ADMET studies. Interestingly, the experimental in vitro prohibitory potential of compound 7 was better than that of sorafenib, demonstrating an IC50 value of 25 nM. Notably, the strong inhibitory effects of compound 10 against two cancer cell lines (MCF-7 and HCT 116) were established with IC50 values of 12.93 and 11.52 μM, disclosing high selectivity indexes of 6.7 and 7.5, respectively.

Multi-Step In Silico Discovery of Natural Drugs against COVID-19 Targeting Main Protease

In continuation of our antecedent work against COVID-19, three natural compounds, namely, Luteoside C (130), Kahalalide E (184), and Streptovaricin B (278) were determined as the most promising SARS-CoV-2 main protease (Mpro) inhibitors among 310 naturally originated antiviral compounds. This was performed via a multi-step in silico method. At first, a molecular structure similarity study was done with PRD_002214, the co-crystallized ligand of Mpro (PDB ID: 6LU7), and favored thirty compounds. Subsequently, the fingerprint study performed with respect to PRD_002214 resulted in the election of sixteen compounds (7, 128, 130, 156, 157, 158, 180, 184, 203, 204, 210, 237, 264, 276, 277, and 278). Then, results of molecular docking versus Mpro PDB ID: 6LU7 favored eight compounds (128, 130, 156, 180, 184, 203, 204, and 278) based on their binding affinities. Then, in silico toxicity studies were performed for the promising compounds and revealed that all of them have good toxicity profiles. Finally, molecular dynamic (MD) simulation experiments were carried out for compounds 130, 184, and 278, which exhibited the best binding modes against Mpro. MD tests revealed that luteoside C (130) has the greatest potential to inhibit SARS-CoV-2 main protease.