In December 2019, a COVID-19 epidemic was discovered in Wuhan, China, and since has disseminated around the world impacting human health for millions. Herein, in-silico drug discovery approaches have been utilized to identify potential natural products (NPs) as Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) main protease (M pro) inhibitors. The MolPort database that contains over 100,000 NPs was screened and filtered using molecular docking techniques. Based on calculated docking scores, the top 5,000 NPs/natural-like products (NLPs) were selected and subjected to molecular dynamics (MD) simulations followed by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. Combined 50 ns MD simulations and MM-GBSA calculations revealed nine potent NLPs with binding affinities (DG binding) > À48.0 kcal/mol. Interestingly, among the identified NLPs, four bis([1,3]dioxolo)pyran-5-carboxamide derivatives showed DG binding > À56.0 kcal/ mol, forming essential short hydrogen bonds with HIS163 and GLY143 amino acids via dioxolane oxygen atoms. Structural and energetic analyses over 50 ns MD simulation demonstrated NLP-M pro complex stability. Drug-likeness predictions revealed the prospects of the identified NLPs as potential drug candidates. The findings are expected to provide a novel contribution to the field of COVID-19 drug discovery.

No abstract

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent for the COVID-19 pandemic, which generated more than 1.82 million deaths in 2020 alone, in addition to 83.8 million infections. Currently, there is no antiviral medication to treat COVID-19. In the search for drug leads, marine-derived metabolites are reported here as prospective SARS-CoV-2 inhibitors. Two hundred and twenty-seven terpene natural products isolated from the biodiverse Red-Sea ecosystem were screened for inhibitor activity against the SARS-CoV-2 main protease (Mpro) using molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics/generalized Born surface area binding energy calculations. On the basis of in silico analyses, six terpenes demonstrated high potency as Mpro inhibitors with ΔGbinding ≤ −40.0 kcal/mol. The stability and binding affinity of the most potent metabolite, erylosides B, were compared to the human immunodeficiency virus protease inhibitor, lopinavir. Erylosides B showed greater binding affinity towards SARS-CoV-2 Mpro than lopinavir over 100 ns with ΔGbinding values of −51.9 vs. −33.6 kcal/mol, respectively. Protein–protein interactions indicate that erylosides B biochemical signaling shares gene components that mediate severe acute respiratory syndrome diseases, including the cytokine- and immune-signaling components BCL2L1, IL2, and PRKC. Pathway enrichment analysis and Boolean network modeling were performed towards a deep dissection and mining of the erylosides B target–function interactions. The current study identifies erylosides B as a promising anti-COVID-19 drug lead that warrants further in vitro and in vivo testing.

Penicillin-binding proteins (PBPs) catalyze the final stages for peptidoglycan cell-wall bio-synthesis. Mutations in the PBP2a subunit can attenuate β-lactam antibiotic activity, resulting in unimpeded cell-wall formation and methicillin-resistant Staphylococcus aureus (MRSA). A double mutation in PBP2a (i.e., N146K and E150K) is resistant to β-lactam inhibitors; however, (E)-3-(2-(4-cyanostyryl)-4-oxoquinazolin-3(4H)-yl) benzoic acid (QNZ), a heterocyclic antibiotic devoid of a β-lactam ring, interacts non-covalently with PBP2a allosteric site and inhibits PBP enzymatic activity. In the search for novel inhibitors that target this PBP2a allosteric site in acidic medium, an in silico screening was performed. Chemical databases including eMolecules, ChEMBL, and ChEBI were virtually screened for candidate inhibitors with a physicochemical similarity to QNZ. PBP2a binding affinities from the screening were calculated based on molecular docking with co-crystallized ligand QNZ serving as a reference. Molecular minimization calculations were performed for inhibitors with docking scores lower than QNZ (calc. −8.3 kcal/mol) followed by combined MD simulations and MM-GBSA binding energy calculations. Compounds eMol26313223 and eMol26314565 exhibited promising inhibitor activities based on binding affinities (ΔGbinding) that were twice that of QNZ (−38.5, −34.5, and −15.4 kcal/mol, respectively). Structural and energetic analyses over a 50 ns MD simulation revealed high stability for the inhibitors when complexed with the double mutated PBP2a. The pharmacokinetic properties of the two inhibitors were predicted using an in silico ADMET analysis. Calculated binding affinities hold promise for eMol26313223 and eMol26314565 as allosteric inhibitors of PBP2a in acidic medium and establish that further in vitro and in vivo inhibition experimentation is warranted.

The P-glycoprotein (P-gp/ABCB1) is responsible for a xenobiotic efflux pump that shackles intracellular drug accumulation. Additionally, it is included in the dud of considerable antiviral and anticancer chemotherapies because of the multidrug resistance (MDR) phenomenon. In the search for prospective anticancer drugs that inhibit the ABCB1 transporter, the Natural Product Activity and Species Source (NPASS) database, containing >35,000 molecules, was explored for identifying ABCB1 inhibitors. The performance of AutoDock4.2.6 software to anticipate ABCB1 docking score and pose was first assessed according to available experimental data. The docking scores of the NPASS molecules were predicted against the ABCB1 transporter. Molecular dynamics (MD) simulations were conducted for molecules with docking scores lower than taxol, a reference inhibitor, pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. On the basis of MM-GBSA calculations, five compounds revealed promising binding affinities as ABCB1 inhibitors with ΔGbinding < −105.0 kcal/mol. The binding affinity and stability of the identified inhibitors were compared to the chemotherapeutic agent. Structural and energetical analyses unveiled great steadiness of the investigated inhibitors within the ABCB1 active site throughout 100 ns MD simulations. Conclusively, these findings point out that NPC104372, NPC475164, NPC2313, NPC197736, and NPC477344 hold guarantees as potential ABCB1 drug candidates and warrant further in vitro/in vivo tests.

Sirtuin 2 (SIRT2) is a member of the sirtuin protein family, which includes lysine deacylases that are NAD+-dependent and organize several biological processes. Different forms of cancer have been associated with dysregulation of SIRT2 activity. Hence, identifying potent inhibitors for SIRT2 has piqued considerable attention in the drug discovery community. In the current study, the Natural Products Atlas (NPAtlas) database was mined to hunt potential SIRT2 inhibitors utilizing in silico techniques. Initially, the performance of the employed docking protocol to anticipate ligand-SIRT2 binding mode was assessed according to the accessible experimental data. Based on the predicted docking scores, the most promising NPAtlas molecules were selected and submitted to molecular dynamics (MD) simulations, followed by binding energy computations. Based on the MM-GBSA binding energy estimations over a 200 ns MD course, three NPAtlas compounds, namely NPA009578, NPA006805, and NPA001884, were identified with better ΔGbinding towards SIRT2 protein than the native ligand (SirReal2) with values of − 59.9, − 57.4, − 53.5, and − 49.7 kcal/mol, respectively. On the basis of structural and energetic assessments, the identified NPAtlas compounds were confirmed to be steady over a 200 ns MD course. The drug-likeness and pharmacokinetic characteristics of the identified NPAtlas molecules were anticipated, and robust bioavailability was predicted. Conclusively, the current results propose potent inhibitors for SIRT2 deserving more in vitro/in vivo investigation.

No abstract

Multidrug resistance (MDR) is one of the most problematic issues in chemotherapeutic carcinoma therapy. The ABCB1 transporter, a drug efflux pump overexpressed in cancer cells, has been thoroughly investigated for its association with MDR. Thus, discovering ABCB1 inhibitors can reverse the MDR in cancer cells. In the current work, a molecular docking technique was utilized for hunting the most prospective ABCB1 inhibitors from the Toxin and Toxin-Target Database (T3DB). Based on the docking computations, the most promising T3DB compounds complexed with the ABCB1 transporter were subjected to molecular dynamics (MD) simulations over 100 ns. Utilizing the MM-GBSA approach, the corresponding binding affinities were computed. Compared to ZQU (calc. −49.8 kcal/mol), Emamectin B1a (T3D1043), Emamectin B1b (T3D1044), Vincristine (T3D4016), Vinblastine (T3D4017), and Vindesine (T3D2479) complexed with ABCB1 transporter demonstrated outstanding binding affinities with ΔGbinding values of −93.0, −92.6, −93.8, −92.2, and −90.8 kcal/mol, respectively. The structural and energetic investigations confirmed the constancy of the identified T3DB compounds complexed with the ABCB1 transporter during the 100 ns MD course. To mimic the physiological conditions, MD simulations were conducted for those identified inhibitors complexed with ABCB1 transporter in the presence of a POPC membrane. These findings revealed that Emamectin B1a, Emamectin B1b, Vincristine, Vinblastine, and Vindesine are promising ABCB1 inhibitors that can reverse the MDR. Therefore, subjecting those compounds to further in-vitro and in-vivo investigations is worthwhile.