In silico drug repurposing and lipid bilayer molecular dynamics puzzled out potential breast cancer resistance protein (BCRP/ABCG2) inhibitors

Multidrug resistance (MDR) is a fundamental reason for the fiasco of carcinoma chemotherapy. A wide variety of anticarcinoma drugs are expelled from neoplasm cells through the ATP-binding cassette (ABC) transporter superfamily, rendering the neoplasm cells resistant to treatment. The ATP-binding cassette transporter G2 (ABCG2, gene symbol BCRP) is an ABC efflux transporter that plays a key function in MDR to antineoplastic therapies. For these reasons, the identification of medicaments as BCRP inhibitors could assist in discovering better curative approaches for breast cancer therapy. Because of the deficiency of prospective BCRP inhibitors, the SuperDRUG2 database was virtually screened for inhibitor activity towards the BCRP transporter using molecular docking computations. The most potent drug candidates were then characterized utilizing molecular dynamics (MD) simulations. Furthermore, molecular mechanics-generalized Born surface area (MM-GBSA) binding affinities of the most potent drug candidates were estimated. Based on the MM-GBSA binding affinities throughout 150 ns MD simulations, three drugs-namely zotarolimus (SD002595), temsirolimus (SD003393), and glecaprevir (SD006009)-revealed greater binding affinities towards BCRP transporter compared to the co-crystallized BWQ ligand with ΔGbinding values of -86.6 ± 5.6, -79.5 ± 8.0, -75.8 ± 4.6 and -59.5 ± 4.1 kcal/mol, respectively. The steadiness of these promising drugs bound with BCRP transporter was examined utilizing their structural and energetical analyses throughout a 150 ns MD simulation. To imitate the physiological environment, 150 ns MD simulations for the identified drugs bound with BCRP transporter were conducted in the 1-palmitoyl-2-oleoyl-phosphatidylcholine lipid bilayer. These findings identify zotarolimus, temsirolimus and glecaprevir as auspicious anti-MDR drug leads that warrant further experimental assays.Communicated by Ramaswamy H. Sarma.

Pyronaridine as a Bromodomain-Containing Protein 4-N-Terminal Bromodomain (BRD4-BD1) Inhibitor: In Silico Database Mining, Molecular Docking, and Molecular Dynamics Simulation

BRD4 (bromodomain-containing protein 4) is an epigenetic reader that realizes histone proteins and promotes the transcription of genes linked to cancer progression and non-cancer diseases such as acute heart failure and severe inflammation. The highly conserved N-terminal bromodomain (BD1) recognizes acylated lysine residues to organize the expression of genes. As such, BD1 is essential for disrupting BRD4 interactions and is a promising target for cancer treatment. To identify new BD1 inhibitors, a SuperDRUG2 database that contains more than 4600 pharmaceutical compounds was screened using in silico techniques. The efficiency of the AutoDock Vina1.1.2 software to anticipate inhibitor-BRD4-BD1 binding poses was first evaluated based on the co-crystallized R6S ligand in complex with BRD4-BD1. From database screening, the most promising BRD4-BD1 inhibitors were subsequently submitted to molecular dynamics (MD) simulations integrated with an MM-GBSA approach. MM-GBSA computations indicated promising BD1 binding with a benzonaphthyridine derivative, pyronaridine (SD003509), with an energy prediction (ΔGbinding) of -42.7 kcal/mol in comparison with -41.5 kcal/mol for a positive control inhibitor (R6S). Pharmacokinetic properties predicted oral bioavailability for both ligands, while post-dynamic analyses of the BRD4-BD1 binding pocket demonstrated greater stability for pyronaridine. These results confirm that in silico studies can provide insight into novel protein-ligand regulators, specifically that pyronaridine is a potential cancer drug candidate.

SuperNatural inhibitors to reverse multidrug resistance emerged by ABCB1 transporter: Database mining, lipid-mediated molecular dynamics, and pharmacokinetics study

An effective approach to reverse multidrug resistance (MDR) is P-glycoprotein (P-gp, ABCB1) transport inhibition. To identify such molecular regulators, the SuperNatural II database, which comprises > 326,000 compounds, was virtually screened for ABCB1 transporter inhibitors. The Lipinski rule was utilized to initially screen the SuperNatural II database, identifying 128,126 compounds. Those natural compounds were docked against the ABCB1 transporter, and those with docking scores less than zosuquidar (ZQU) inhibitor were subjected to molecular dynamics (MD) simulations. Based on MM-GBA binding energy (ΔGbinding) estimations, UMHSN00009999 and UMHSN00097206 demonstrated ΔGbinding values of -68.3 and -64.1 kcal/mol, respectively, compared to ZQU with a ΔGbinding value of -49.8 kcal/mol. For an investigation of stability, structural and energetic analyses for UMHSN00009999- and UMHSN00097206-ABCB1 complexes were performed and proved the high steadiness of these complexes throughout 100 ns MD simulations. Pharmacokinetic properties of the identified compounds were also predicted. To mimic the physiological conditions, MD simulations in POPC membrane surroundings were applied to the UMHSN00009999- and UMHSN00097206-ABCB1 complexes. These results demonstrated that UMHSN00009999 and UMHSN00097206 are promising ABCB1 inhibitors for reversing MDR in cancer and warrant additional in-vitro/in-vivo studies.

In-Silico Mining of the Toxins Database (T3DB) towards Hunting Prospective Candidates as ABCB1 Inhibitors: Integrated Molecular Docking and Lipid Bilayer-Enhanced Molecular Dynamics Study

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 ΔG_binding 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.

Potential drug candidates as P-glycoprotein inhibitors to reverse multidrug resistance in cancer: an in silico drug discovery study

The failure of chemotherapy in the treatment of carcinoma is mainly due to the development of multidrug resistance (MDR), which is largely caused by the overexpression of P-glycoprotein (P-gp/ABCB1/MDR1). Until recently, the 3D structure of the P-gp transporter has not been experimentally resolved, which restricted the discovery of prospective P-gp inhibitors utilizing in silico techniques. In this study, the binding energies of 512 drug candidates in clinical or investigational stages were assessed as potential P-gp inhibitors employing in silico methods. On the basis of the available experimental data, the performance of the AutoDock4.2.6 software to predict the drug-P-gp binding mode was initially validated. Molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics-generalized Born surface area (MM-GBSA) binding energy computations were subsequently conducted to screen the investigated drug candidates. Based on the current results, five promising drug candidates, namely valspodar, dactinomycin, elbasvir, temsirolimus, and sirolimus, showed promising binding energies against P-gp transporter with ΔGbinding values of -126.7, -112.1, -111.9, -102.9, and -101.4 kcal/mol, respectively. The post-MD analyses revealed the energetical and structural stabilities of the identified drug candidates in complex with the P-gp transporter. Furthermore, in order to mimic the physiological conditions, the potent drugs complexed with the P-gp were subjected to 100 ns MD simulations in an explicit membrane-water environment. The pharmacokinetic properties of the identified drugs were predicted and demonstrated good ADMET characteristics. Overall, these results indicated that valspodar, dactinomycin, elbasvir, temsirolimus, and sirolimus hold promise as prospective P-gp inhibitors and warrant further invitro/invivo investigations.

Chetomin, a SARS-CoV-2 3C-like Protease (3CLpro) Inhibitor: In Silico Screening, Enzyme Docking, Molecular Dynamics and Pharmacokinetics Analysis

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to over 6 million deaths. The 3C-like protease (3CLpro) enzyme of the SARS-CoV-2 virus is an attractive druggable target for exploring therapeutic drug candidates to combat COVID-19 due to its key function in viral replication. Marine natural products (MNPs) have attracted considerable attention as alternative sources of antiviral drug candidates. In looking for potential 3CLpro inhibitors, the MNP database (>14,000 molecules) was virtually screened against 3CLpro with the assistance of molecular docking computations. The performance of AutoDock and OEDocking software in anticipating the ligand-3CLpro binding mode was first validated according to the available experimental data. Based on the docking scores, the most potent MNPs were further subjected to molecular dynamics (MD) simulations, and the binding affinities of those molecules were computed using the MM-GBSA approach. According to MM-GBSA//200 ns MD simulations, chetomin (UMHMNP1403367) exhibited a higher binding affinity against 3CLpro than XF7, with ΔGbinding values of −55.5 and −43.7 kcal/mol, respectively. The steadiness and tightness of chetomin with 3CLpro were evaluated, revealing the high stabilization of chetomin (UMHMNP1403367) inside the binding pocket of 3CLpro throughout 200 ns MD simulations. The physicochemical and pharmacokinetic features of chetomin were also predicted, and the oral bioavailability of chetomin was demonstrated. Furthermore, the potentiality of chetomin analogues −namely, chetomin A-D− as 3CLpro inhibitors was investigated. These results warrant further in vivo and in vitro assays of chetomin (UMHMNP1403367) as a promising anti-COVID-19 drug candidate.

Naturally occurring plant-based anticancerous candidates as prospective ABCG2 inhibitors: an in silico drug discovery study

Abstract

ATP-binding cassette transporter G2 (ABCG2) is an efflux transporter related to the clinical multidrug resistance (MDR) phenomenon. Identifying ABCG2 inhibitors could help discover extraordinary curative strategies for carcinoma remediation. Hitherto, there is no medication drug inhibiting ABCG2 transporter, notwithstanding that a considerable number of drugs have been submitted to clinical-trial and investigational phases. In the search for unprecedented chemical compounds that could inhibit the ABCG2 transporter, an in silico screening was conducted on the Naturally Occurring Plant-based Anticancer Compound-Activity-Target (NPACT) database containing 1574 compounds. Inhibitor-ABCG2 binding affinities were estimated based on molecular docking and molecular minimization (MM) calculations and compared to a co-crystallized inhibitor (BWQ) acting as a reference inhibitor. Molecular dynamics (MD) simulations pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations were further executed for compounds with MM-GBSA//MM binding energies lower than BWQ (calc. − 60.5 kcal/mol). NPACT00968 and NPACT01545 demonstrated auspicious inhibitory activities according to binding affinities (ΔGbinding) over the 100 ns MD simulations that were nearly one and a half folds compared to BWQ (− 100.4, − 94.7, and − 62.9 kcal/mol, respectively). Throughout the 100 ns MD simulations, structural and energetical analyses unveiled outstanding stability of the ABCG2 transporter when bound with NPACT00968 and NPACT01545. In silico calculations hold a promise for those two inhibitors as drug candidates of ABCG2 transporter and emphasize that further in vitro and in vivo experiments are guaranteed.

Graphical abstract

In Silico Targeting Human Multidrug Transporter ABCG2 in Breast Cancer: Database Screening, Molecular Docking, and Molecular Dynamics Study

ABCG2 is a substantial member of the ABC transporter superfamily that plays a significant role in multidrug resistance in cancer. Until recently, the 3D structure of ABCG2 has not been resolved, which resulted in the limitation of developing potential ABCG2 inhibitors using structure-based drug discovery. Herein, eMolecules, ChEMBL, and ChEBI databases, containing >25 million compounds, were virtually screened against the ABCG2 transporter in homodimer form. Performance of AutoDock4.2.6 software to predict inhibitor-ABCG2 binding mode and affinity were validated on the basis of available experimental data. The explored databases were filtered based on docking scores. The most potent hits with binding affinities higher than that of experimental bound ligand (MZ29) were then selected and subjected to molecular mechanics minimization, followed by binding energy calculation using molecular mechanics-generalized Born surface area (MM-GBSA) approach. Furthermore, molecular dynamics simulations for 50 ns, followed by MM-GBSA binding energy calculations, were performed for the promising compounds, unveiling eight potential inhibitors with binding affinities <-55.8 kcal/mol. Structural and energetic analyses demonstrated the stability of the eight identified inhibitors over the 50 ns MD simulation. This research sheds light on the potentiality of the identified ABCG2 inhibitors as a therapeutic approach to overcome multidrug resistance cancer therapy.

Blue Biotechnology: Computational Screening of Sarcophyton Cembranoid Diterpenes for SARS-CoV-2 Main Protease Inhibition

The coronavirus pandemic has affected more than 150 million people, while over 3.25 million people have died from the coronavirus disease 2019 (COVID-19). As there are no established therapies for COVID-19 treatment, drugs that inhibit viral replication are a promising target; specifically, the main protease (Mpro) that process CoV-encoded polyproteins serves as an Achilles heel for assembly of replication-transcription machinery as well as down-stream viral replication. In the search for potential antiviral drugs that target Mpro, a series of cembranoid diterpenes from the biologically active soft-coral genus Sarcophyton have been examined as SARS-CoV-2 Mpro inhibitors. Over 360 metabolites from the genus were screened using molecular docking calculations. Promising diterpenes were further characterized by molecular dynamics (MD) simulations based on molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. According to in silico calculations, five cembranoid diterpenes manifested adequate binding affinities as Mpro inhibitors with ΔGbinding < -33.0 kcal/mol. Binding energy and structural analyses of the most potent Sarcophyton inhibitor, bislatumlide A (340), was compared to darunavir, an HIV protease inhibitor that has been recently subjected to clinical-trial as an anti-COVID-19 drug. In silico analysis indicates that 340 has a higher binding affinity against Mpro than darunavir with ΔGbinding values of -43.8 and -34.8 kcal/mol, respectively throughout 100 ns MD simulations. Drug-likeness calculations revealed robust bioavailability and protein-protein interactions were identified for 340; biochemical signaling genes included ACE, MAPK14 and ESR1 as identified based on a STRING database. Pathway enrichment analysis combined with reactome mining revealed that 340 has the capability to re-modulate the p38 MAPK pathway hijacked by SARS-CoV-2 and antagonize injurious effects. These findings justify further in vivo and in vitro testing of 340 as an antiviral agent against SARS-CoV-2.

Prospective Drug Candidates as Human Multidrug Transporter ABCG2 Inhibitors: an In Silico Drug Discovery Study

Breast cancer resistance protein (ABCG2) is a human ATP-binding cassette (ABC) that plays a paramount role in multidrug resistance (MDR) in cancer therapy. The discovery of ABCG2 inhibitors could assist in designing unprecedented therapeutic strategies for cancer treatment. There is as yet no approved drug targeting ABCG2, although a large number of drug candidates have been clinically investigated. In this work, binding affinities of 181 drug candidates in clinical-trial or investigational stages as ABCG2 inhibitors were inspected using in silico techniques. Based on available experimental data, the performance of AutoDock4.2.6 software was first validated to predict the inhibitor-ABCG2 binding mode and affinity. Combined molecular docking calculations and molecular dynamics (MD) simulations, followed by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations, were then performed to filter out the studied drug candidates. From the estimated docking scores and MM-GBSA binding energies, six auspicious drug candidates-namely, pibrentasvir, venetoclax, ledipasvir, avatrombopag, cobicistat, and revefenacin-exhibited auspicious binding energies with value < -70.0 kcal/mol. Interestingly, pibrentasvir, venetoclax, and ledipasvir were observed to show even higher binding affinities with the ABCG2 transporter with binding energies of < -80.0 kcal/mol over long MD simulations of 100 ns. The stabilities of these three promising candidates in complex with ABCG2 transporter were demonstrated by their energetics and structural analyses throughout the 100 ns MD simulations. The current study throws new light on pibrentasvir, venetoclax, and ledipasvir as curative options for multidrug resistant cancers by inhibiting ABCG2 transporter.

In Silico Mining of Terpenes from Red-Sea Invertebrates for SARS-CoV-2 Main Protease (Mpro) Inhibitors

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.

Rutin and flavone analogs as prospective SARS-CoV-2 main protease inhibitors: In silico drug discovery study

Coronavirus disease 2019 (COVID-19) is a new pandemic characterized by quick spreading and illness of the respiratory system. To date, there is no specific therapy for Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2). Flavonoids, especially rutin, have attracted considerable interest as a prospective SARS-CoV-2 main protease (M^pro) inhibitor. In this study, a database containing 2017 flavone analogs was prepared and screened against SARS-CoV-2 M^pro using the molecular docking technique. According to the results, 371 flavone analogs exhibited good potency towards M^pro with docking scores less than −9.0 kcal/mol. Molecular dynamics (MD) simulations, followed by molecular mechanics-generalized Born surface area (MM/GBSA) binding energy calculations, were performed for the top potent analogs in complex with M^pro. Compared to rutin, PubChem-129-716-607 and PubChem-885-071-27 showed better binding affinities against SARS-CoV-2 M^pro over 150 ns MD course with ΔG_binding values of −69.0 and −68.1 kcal/mol, respectively. Structural and energetic analyses demonstrated high stability of the identified analogs inside the SARS-CoV-2 M^pro active site over 150 ns MD simulations. The oral bioavailabilities of probable SARS-CoV-2 M^pro inhibitors were underpinned using drug-likeness parameters. A comparison of the binding affinities demonstrated that the MM/GBSA binding energies of the identified flavone analogs were approximately three and two times less than those of lopinavir and baicalein, respectively. In conclusion, PubChem-129-716-607 and PubChem-885-071-27 are promising anti-COVID-19 drug candidates that warrant further clinical investigations.

In silico drug discovery of major metabolites from spices as SARS-CoV-2 main protease inhibitors

Coronavirus Disease 2019 (COVID-19) is an infectious illness caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), originally identified in Wuhan, China (December 2019) and has since expanded into a pandemic. Here, we investigate metabolites present in several common spices as possible inhibitors of COVID-19. Specifically, 32 compounds isolated from 14 cooking seasonings were examined as inhibitors for SARS-CoV-2 main protease (M^pro), which is required for viral multiplication. Using a drug discovery approach to identify possible antiviral leads, in silico molecular docking studies were performed. Docking calculations revealed a high potency of salvianolic acid A and curcumin as M^pro inhibitors with binding energies of -9.7 and -9.2 kcal/mol, respectively. Binding mode analysis demonstrated the ability of salvianolic acid A and curcumin to form nine and six hydrogen bonds, respectively with amino acids proximal to M^pro's active site. Stabilities and binding affinities of the two identified natural spices were calculated over 40 ns molecular dynamics simulations and compared to an antiviral protease inhibitor (lopinavir). Molecular mechanics-generalized Born surface area energy calculations revealed greater salvianolic acid A affinity for the enzyme over curcumin and lopinavir with energies of -44.8, -34.2 and -34.8 kcal/mol, respectively. Using a STRING database, protein-protein interactions were identified for salvianolic acid A included the biochemical signaling genes ACE, MAPK14 and ESR1; and for curcumin, EGFR and TNF. This study establishes salvianolic acid A as an in silico natural product inhibitor against the SARS-CoV-2 main protease and provides a promising inhibitor lead for in vitro enzyme testing.

Carotane sesquiterpenes from Ferula vesceritensis: in silico analysis as SARS-CoV-2 binding inhibitors

Two sesquiterpenes, 8α-anisate-dauc-4-ene-3,9-dione (webiol anisate) (1) and 10α-acetoxy-6α-benzoate-jaeschkeanadiol (2) as well as, ten known analogues (3-10), and two sesquiterpene coumarins (11-12) were isolated from an organic root extract of Ferula vesceritensis (Fam. Apiaceae). Chemical structures were elucidated based on IR, 1D- and 2D-NMR and HRMS, spectroscopic analyses. With molecular overlap observed between two protease inhibitors that are being examined as anti-COVID-19 drugs, and sesquiterpenes isolated here, metabolite molecular docking calculations were made using the main protease (Mpro), which is required for viral multiplication as well as RNA-dependent RNA polymerase (RdRp). In silico binding-inhibition analysis predicted that select F. vesceritensis sesquiterpenes can bind to these enzymes required for viral replication. Structures of the isolated constituents were also consistent with the chemo-systematic grouping of F. vesceritensis secondary metabolites with other Ferula species.

In-silico drug repurposing and molecular dynamics puzzled out potential SARS-CoV-2 main protease inhibitors

Herein, the DrugBank database which contains 10,036 approved and investigational drugs was explored deeply for potential drugs that target SARS-CoV-2 main protease (M^pro). Filtration process of the database was conducted using three levels of accuracy for molecular docking calculations. The top 35 drugs with docking scores > −11.0 kcal/mol were then subjected to 10 ns molecular dynamics (MD) simulations followed by molecular mechanics–generalized Born surface area (MM-GBSA) binding energy calculations. The results showed that DB02388 and Cobicistat (DB09065) exhibited potential binding affinities towards M^pro over 100 ns MD simulations, with binding energy values of −49.67 and −46.60 kcal/mol, respectively. Binding energy and structural analyses demonstrated the higher stability of DB02388 over Cobicistat. The potency of DB02388 and Cobicistat is attributed to their abilities to form several hydrogen bonds with the essential amino acids inside the active site of M^pro. Compared to DB02388 and Cobicistat, Darunavir showed a much lower binding affinity of −34.83 kcal/mol. The present study highlights the potentiality of DB02388 and Cobicistat as anti-COVID-19 drugs for clinical trials.

Communicated by Ramaswamy H. Sarma.

Natural-like products as potential SARS-CoV-2 Mpro inhibitors: in-silico drug discovery

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 (ΔG_binding) > -48.0 kcal/mol. Interestingly, among the identified NLPs, four bis([1,3]dioxolo)pyran-5-carboxamide derivatives showed ΔG_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.Communicated by Ramaswamy H. Sarma.

Repurposing of FDA-approved antivirals, antibiotics, anthelmintics, antioxidants, and cell protectives against SARS-CoV-2 papain-like protease

SARS-CoV-2 or Coronavirus disease 19 (COVID-19) is a rapidly spreading, highly contagious, and sometimes fatal disease for which drug discovery and vaccine development are critical. SARS-CoV-2 papain-like protease (PL^pro) was used to virtually screen 1697 clinical FDA-approved drugs. Among the top results expected to bind with SARS-CoV-2 PL^pro strongly were three cell protectives and antioxidants (NAD+, quercitrin, and oxiglutatione), three antivirals (ritonavir, moroxydine, and zanamivir), two antimicrobials (doripenem and sulfaguanidine), two anticancer drugs, three benzimidazole anthelmintics, one antacid (famotidine), three anti-hypertensive ACE receptor blockers (candesartan, losartan, and valsartan) and other miscellaneous systemically or topically acting drugs. The binding patterns of these drugs were superior to the previously identified SARS CoV PL^pro inhibitor, 6-mercaptopurine (6-MP), suggesting a potential for repurposing these drugs to treat COVID-19. The objective of drug repurposing is the rapid relocation of safe and approved drugs by bypassing the lengthy pharmacokinetic, toxicity, and preclinical phases. The ten drugs with the highest estimated docking scores with favorable pharmacokinetics were subjected to molecular dynamics (MD) simulations followed by molecular mechanics/generalized Born surface area (MM/GBSA) binding energy calculations. Phenformin, quercetin, and ritonavir all demonstrated prospective binding affinities for COVID-19 PL^pro over 50 ns MD simulations, with binding energy values of −56.6, −40.9, and −37.6 kcal/mol, respectively. Energetic and structural analyses showed phenformin was more stable than quercetin and ritonavir. The list of the drugs provided herein constitutes a primer for clinical application in COVID-19 patients and guidance for further antiviral studies.

Communicated by Ramaswamy H. Sarma

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling, molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors -namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics-generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ -34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.