Integrated Computational Approach for Virtual Hit Identification against Ebola Viral Proteins VP35 and VP40

Identification of Prospective Ebola Virus VP35 and VP40 Protein Inhibitors from Myxobacterial Natural Products

The Ebola virus (EBOV) is a lethal pathogen causing hemorrhagic fever syndrome which remains a global health challenge. In the EBOV, two multifunctional proteins, VP35 and VP40, have significant roles in replication, virion assembly, and budding from the cell and have been identified as druggable targets. In this study, we employed in silico methods comprising molecular docking, molecular dynamic simulations, and pharmacological properties to identify prospective drugs for inhibiting VP35 and VP40 proteins from the myxobacterial bioactive natural product repertoire. Cystobactamid 934-2, Cystobactamid 919-1, and Cittilin A bound firmly to VP35. Meanwhile, 2-Hydroxysorangiadenosine, Enhypyrazinone B, and Sorangiadenosine showed strong binding to the matrix protein VP40. Molecular dynamic simulations revealed that, among these compounds, Cystobactamid 919-1 and 2-Hydroxysorangiadenosine had stable interactions with their respective targets. Similarly, molecular mechanics Poisson-Boltzmann surface area (MMPBSA) calculations indicated close-fitting receptor binding with VP35 or VP40. These two compounds also exhibited good pharmacological properties. In conclusion, we identified Cystobactamid 919-1 and 2-Hydroxysorangiadenosine as potential ligands for EBOV that target VP35 and VP40 proteins. These findings signify an essential step in vitro and in vivo to validate their potential for EBOV inhibition.

Identification of novel inhibitors against VP40 protein of Marburg virus by integrating molecular modeling and dynamics approaches

Marburg virus (MV) is a highly etiological agent of haemorrhagic fever in humans and has spread across the world. Its outbreaks caused a 23-90% human death rate. However, there are currently no authorized preventive or curative measures yet. VP40 is the MV matrix protein, which builds protein shell underneath the viral envelope and confers hallmark filamentous. VP40 alone is able to induce assembly and budding of filamentous virus-like particles (VLPs), which resemble authentic virions. As a result, this research is credited with clarifying the function of VP40 and leading to the discovery of new therapeutic targets effective in combating MV disease (MVD). Virtual screening, molecular docking and molecular dynamics (MD) simulation were used to find the putative active chemicals based on a 3D pharmacophore model of the protein's active site cavity. Initially, andrographidine-C, a potent inhibitor was selected for the development of the pharmacophore model. Later, a library of 30,000 compounds along with the andrographidine-C was docked against VP40 protein. Three best hits including avanafil, diuvaretin and macrourone were subjected to further MD simulation analysis, as these compounds had better binding affinities as compared to andrographidine-C. Furthermore, throughout the 100 ns simulations, the back bone of VP40 protein in presence of avanafil, diuvaretin and macrourone remained stable which was further validated by MM-PBSA analysis. Additionally, all of these compounds depict maximum drug-like properties. The predicted drugs based on the ligand, avanafil, diuvaretin and macrourone could be exploited and developed as an alternative or complementary therapy for the treatment of MVD.Communicated by Ramaswamy H. Sarma.

Structural and Dynamical Basis of VP35-RBD Inhibition by Marine Fungi Compounds to Combat Marburg Virus Infection

The Marburg virus (MBV), a deadly pathogen, poses a serious threat to world health due to the lack of effective treatments, calling for an immediate search for targeted and efficient treatments. In this study, we focused on compounds originating from marine fungi in order to identify possible inhibitory compounds against the Marburg virus (MBV) VP35-RNA binding domain (VP35-RBD) using a computational approach. We started with a virtual screening procedure using the Lipinski filter as a guide. Based on their docking scores, 42 potential candidates were found. Four of these compounds-CMNPD17596, CMNPD22144, CMNPD25994, and CMNPD17598-as well as myricetin, the control compound, were chosen for re-docking analysis. Re-docking revealed that these particular compounds had a higher affinity for MBV VP35-RBD in comparison to the control. Analyzing the chemical interactions revealed unique binding properties for every compound, identified by a range of Pi-cation interactions and hydrogen bond types. We were able to learn more about the dynamic behaviors and stability of the protein-ligand complexes through a 200-nanosecond molecular dynamics simulation, as demonstrated by the compounds' consistent RMSD and RMSF values. The multidimensional nature of the data was clarified by the application of principal component analysis, which suggested stable conformations in the complexes with little modification. Further insight into the energy profiles and stability states of these complexes was also obtained by an examination of the free energy landscape. Our findings underscore the effectiveness of computational strategies in identifying and analyzing potential inhibitors for MBV VP35-RBD, offering promising paths for further experimental investigations and possible therapeutic development against the MBV.

Identification of potential inhibitors targeting Ebola virus VP35 protein: a computational strategy

Ebola virus (EBOV) poses a severe threat as a highly infectious pathogen, causing devastating hemorrhagic fever in both humans and animals. The EBOV virus VP35 protein plays a crucial role in viral replication and exhibits the ability to suppress the host interferon cascade, leading to immune system depletion. As a potential drug target, VP35 protein inhibition holds promise for combating EBOV. To discover new drug candidates, we employed a computer-aided drug design approach, focusing on compounds capable of inhibiting VP35 protein replication. In this connection, a pharmacophore model was generated using molecular interactions between the VP35 protein and its inhibitor. ZINC and Cambridge database were screened using validated pharmacophore model. Further the compounds were filtered based on Lipinski's rule of five and subjected to MD simulation and relative binding free energy calculation. Six compounds manifest a significant docking score and strong binding interaction towards VP35 protein. MD simulations further confirmed the remarkable stability of these six complexes. Relative binding free energy calculations also showed significant ΔG value in the range of -132.3 and -49.3 kcal/mol. This study paves the way for further optimization of these compounds as potential inhibitors of VP35, facilitating subsequent experimental in vitro studies.Communicated by Ramaswamy H. Sarma.

Folate Conjugated Polyethylene Glycol Probe for Tumor-Targeted Drug Delivery of 5-Fluorouracil

A targeted delivery system is primarily intended to carry a potent anticancer drug to specific tumor sites within the bodily tissues. In the present study, a carrier system has been designed using folic acid (FA), bis-amine polyethylene glycol (PEG), and an anticancer drug, 5-fluorouracil (5-FU). FA and PEG were joined via an amide bond, and the resulting FA-PEG-NH₂ was coupled to 5-FU producing folate-polyethylene glycol conjugated 5-fluorouracil (FA-PEG-5-FU). Spectroscopic techniques (UV-Vis, ¹HNMR, FTIR, and HPLC) were used for the characterization of products. Prodrug (FA-PEG-5-FU) was analyzed for drug release profile (in vitro) up to 10 days and compared to a standard anticancer drug (5-FU). Folate conjugate was also analyzed to study its folate receptors (FR) mediated transport and in vitro cytotoxicity assays using HeLa cancer cells/Vero cells, respectively, and antitumor activity in tumor-bearing mice models. Folate conjugate showed steady drug release patterns and improved uptake in the HeLa cancer cells than Vero cells. Folate conjugate treated mice group showed smaller tumor volumes; specifically after the 15th day post-treatment, tumor sizes were decreased significantly compared to the standard drug group (5-FU). Molecular docking findings demonstrated importance of Trp138, Trp140, and Lys136 in the stabilization of flexible loop flanking the active site. The folic acid conjugated probe has shown the potential of targeted drug delivery and sustained release of anticancer drug to tumor lesions with intact antitumor efficacy.

In Silico Screening of Marine Compounds as an Emerging and Promising Approach against Estrogen Receptor Alpha-Positive Breast Cancer

Presently, the majority of breast tumors are estrogen receptor (ER) positive. Breast cancer (BC) is defined by uncontrolled cell proliferation (CP) in breast tissue. BCs are caused by the overexpression of genes that promote CP in breast cells. The discovery of effective inhibitors is an excellent chemopreventive method. Our in silico approach analysis offers a pharmacoinformatics methodology for identifying lead molecules targeting cochaperone HSP90 and the epidermal growth factor receptors (EGFR) and human epidermal growth factor receptor 2 (HER2)/neu receptor. BC has been associated with the high expression of these targets. The use of drug-likeness filters aided in determining the therapeutic properties of possible lead compounds. In this study, docking-based virtual screening (VS) was performed. Database of about 450 cancer marine compounds was used. The X-ray-assisted structure of ERα with 4-OHT (PDB code: 3ERT) was chosen for 4-OHT. A docking-based virtual screening was performed on the dataset supplied using the molecular operating environment (MOE) dock application. The binding energy (BE) and explanation of the protein inhibitor interaction (PII) are crucial findings for future both in terms of dry or wet lab research. The GBVI/WAS binding-free energy assessment (in kcal/mol) scores were used to grade the compounds. Compounds with a BE of less than -9.500 kcal/mol were deemed to be the most effective inhibitors. For further analysis, the top seven structurally diverse scaffolds were selected. Seven marine compounds exhibited the best docking score, which validates them to be potent anti-BC compounds. These compounds' bioactive potential and prospective drug-likeness profile make them promising leads for further experimental research.

Computational Study on Potential Novel Anti-Ebola Virus Protein VP35 Natural Compounds

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/β production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of −8.2, −8.2, −8.1 and −8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC₅₀), found using a random forest classifier, ranged from 3.35 to 11.99 μM, while the Ki values ranged from 0.97 to 1.37 μM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of −8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein–ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Fragment-based in silico design of SARS-CoV-2 main protease inhibitors

3CLpro is essential for SARS-CoV-2 replication and infection; its inhibition using small molecules is a potential therapeutic strategy. In this study, a comprehensive crystallography-guided fragment-based drug discovery approach was employed to design new inhibitors for SARS-CoV-2 3CLpro. All small molecules co-crystallized with SARS-CoV-2 3CLpro with structures deposited in the Protein Data Bank were used as inputs. Fragments sitting in the binding pocket (87) were grouped into eight geographical types. They were interactively coupled using various synthetically reasonable linkers to generate larger molecules with divalent binding modes taking advantage of two different fragments' interactions. In total, 1,251 compounds were proposed, and 7,158 stereoisomers were screened using Glide (standard precision and extra precision), AutoDock Vina, and Prime MMGBSA. The top 22 hits having conformations approaching the linear combination of their constituent fragments were selected for MD simulation on Desmond. MD simulation suggested 15 of these did adopt conformations very close to their constituent pieces with far higher binding affinity than either constituent domain alone. These structures could provide a starting point for the further design of SARS-CoV-2 3CLpro inhibitors with improved binding, and structures are provided.

Pharmacoinformatics and molecular dynamics simulation studies reveal potential covalent and FDA-approved inhibitors of SARS-CoV-2 main protease 3CLpro

The SARS-CoV-2 was confirmed to cause the global pandemic of coronavirus disease 2019 (COVID-19). The 3-chymotrypsin-like protease (3CLpro), an essential enzyme for viral replication, is a valid target to combat SARS-CoV and MERS-CoV. In this work, we present a structure-based study to identify potential covalent inhibitors containing a variety of chemical warheads. The targeted Asinex Focused Covalent (AFCL) library was screened based on different reaction types and potential covalent inhibitors were identified. In addition, we screened FDA-approved protease inhibitors to find candidates to be repurposed against SARS-CoV-2 3CLpro. A number of compounds with significant covalent docking scores were identified. These compounds were able to establish a covalent bond (C-S) with the reactive thiol group of Cys145 and to form favorable interactions with residues lining the substrate-binding site. Moreover, paritaprevir and simeprevir from FDA-approved protease inhibitors were identified as potential inhibitors of SARS-CoV-2 3CLpro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CLpro were characterized by molecular dynamics (MD) simulations. The identified compounds are potential inhibitors worthy of further development as COVID-19 drugs. Importantly, the identified FDA-approved anti-hepatitis-C virus (HCV) drugs paritaprevir and simeprevir could be ready for clinical trials to treat infected patients and help curb COVID-19. Communicated by Ramaswamy H. Sarma.

Antiparkinsonian activity of Cucurbita pepo seeds along with possible underlying mechanism

Cucurbita pepo is used as a vegetable in Pakistan and its seeds are also rich in tocopherol. Data showed the pivotal role of tocopherol in the treatment of Parkinson's disease (PD). The current study was designed to probe into the antiparkinson activity of methanolic extract of C. pepo (MECP) seeds in the haloperidol-induced Parkinson rat model. Behavioral studies showed improvement in motor functions. The increase in catalase, superoxide dismutase, glutathione levels whereas the decreases in the malondialdehyde and nitrite levels were noted in a dose-dependent manner. Acetylcholine-esterase (AchE) activity was increased. Molecular docking results revealed significant binding interaction of selected phytoconstituents within an active site of target protein AchE (PDB ID: 4EY7). Furthermore, α-synuclein was up regulated with down regulation of TNF-α and IL-1β in the qRT-PCR study. Subsequently, ADMET results on the basis of structure to activity predictions in terms of pharmacokinetics and toxicity estimations show that selected phytochemicals exhibited moderately acceptable properties. These properties add knowledge towards the structural features which could improve the bioavailability of selected phytochemicals before moving towards the initial phase of the drug development. Our integrated drug discovery scheme concluded that C. pepo seeds could ameliorate symptoms of PD and may prove a lead remedy for the treatment of PD.

Molecular Docking Studies Reveal Rhein from rhubarb (Rheum rhabarbarum) as a Putative Inhibitor of ATP-binding Cassette Super-family G member 2

BACKGROUND

ATP-binding cassette Super-family G member 2 protein is an active ATPbinding cassette transporter with the potential to combat cancer stem cells.

OBJECTIVE

Due to the lack of potential ATP-binding cassette Super-family G member 2 inhibitors, we screened natural inhibitors, which could be a safe source to control multidrug resistance by blocking the regulation of ATP-binding cassette Super-family G member 2 protein.

METHODS

Three-dimensional structure of ATP-binding cassette Super-family G member 2 protein downloaded from the protein databank and chemical structures of 166 selected compounds of the training dataset were retrieved from PubChem. Drug-likeness and docking analysis was conducted to shortlist the dataset for pharmacophore generation. LigandScout 4.1.5 used for pharmacophorebased screening of Zbc library of ZINC database and Autodock Vina were utilized for molecular docking against the predicted active pocket of the target protein to evaluate the potential association of protein and ligands. The physiochemical properties of novel compounds were calculated by admetSAR respectively.

RESULTS

Through pharmacophore-based screening, ZINC4098704 (Rhein) was identified as a lead compound which demonstrates the least binding energy (-8.5) and the highest binding affinity with the target protein and showed optimal physiochemical profile. This compound is highly recommended for a laboratory test to confirm its activity as an ATP-binding cassette Super-family G member 2 inhibitors.

CONCLUSION

Our computer-based study systematically selected natural lead compounds, which could be effective in inhibiting ATP-binding cassette Super-family G member 2 and may help reverse the effect of multidrug resistance to increase the effectiveness of chemotherapy in cancer treatment.

Machine learning techniques applied to the drug design and discovery of new antivirals: a brief look over the past decade

Introduction: Drug design and discovery of new antivirals will always be extremely important in medicinal chemistry, taking into account known and new viral diseases that are yet to come. Although machine learning (ML) have shown to improve predictions on the biological potential of chemicals and accelerate the discovery of drugs over the past decade, new methods and their combinations have improved their performance and established promising perspectives regarding ML in the search for new antivirals.Areas covered: The authors consider some interesting areas that deal with different ML techniques applied to antivirals. Recent innovative studies on ML and antivirals were selected and analyzed in detail. Also, the authors provide a brief look at the past to the present to detect advances and bottlenecks in the area.Expert opinion: From classical ML techniques, it was possible to boost the searches for antivirals. However, from the emergence of new algorithms and the improvement in old approaches, promising results will be achieved every day, as we have observed in the case of SARS-CoV-2. Recent experience has shown that it is possible to use ML to discover new antiviral candidates from virtual screening and drug repurposing.

Targeting Ebola virus replication through pharmaceutical intervention

Introduction. The consistent emergence/reemergence of filoviruses into a world that previously lacked an approved pharmaceutical intervention parallels an experience repeatedly played-out for most other emerging pathogenic zoonotic viruses. Investment to preemptively develop effective and low-cost prophylactic and therapeutic interventions against viruses that have high potential for emergence and societal impact should be a priority.Areas covered. Candidate drugs can be characterized into those that interfere with cellular processes required for Ebola virus (EBOV) replication (host-directed), and those that directly target virally encoded functions (direct-acting). We discuss strategies to identify pharmaceutical interventions for EBOV infections. PubMed/Web of Science databases were searched to establish a detailed catalog of these interventions.Expert opinion. Many drug candidates show promising in vitro inhibitory activity, but experience with EBOV shows the general lack of translation to in vivo efficacy for host-directed repurposed drugs. Better translation is seen for direct-acting antivirals, in particular monoclonal antibodies. The FDA-approved monoclonal antibody treatment, Inmazeb™ is a success story that could be improved in terms of impact on EBOV-associated disease and mortality, possibly by combination with other direct-acting agents targeting distinct aspects of the viral replication cycle. Costs need to be addressed given EBOV emergence primarily in under-resourced countries.

A Molecular Modeling Approach to Identify Potential Antileishmanial Compounds Against the Cell Division Cycle (cdc)-2-Related Kinase 12 (CRK12) Receptor of Leishmania donovani

The huge burden of leishmaniasis caused by the trypanosomatid protozoan parasite Leishmania is well known. This illness was included in the list of neglected tropical diseases targeted for elimination by the World Health Organization. However, the increasing evidence of resistance to existing antimonial drugs has made the eradication of the disease difficult to achieve, thus warranting the search for new drug targets. We report here studies that used computational methods to identify inhibitors of receptors from natural products. The cell division cycle-2-related kinase 12 (CRK12) receptor is a plausible drug target against Leishmania donovani. This study modelled the 3D molecular structure of the L. donovani CRK12 (LdCRK12) and screened for small molecules with potential inhibitory activity from African flora. An integrated library of 7722 African natural product-derived compounds and known inhibitors were screened against the LdCRK12 using AutoDock Vina after performing energy minimization with GROMACS 2018. Four natural products, namely sesamin (NANPDB1649), methyl ellagic acid (NANPDB1406), stylopine (NANPDB2581), and sennecicannabine (NANPDB6446) were found to be potential LdCRK12 inhibitory molecules. The molecular docking studies revealed two compounds NANPDB1406 and NANPDB2581 with binding affinities of -9.5 and -9.2 kcal/mol, respectively, against LdCRK12 which were higher than those of the known inhibitors and drugs, including GSK3186899, amphotericin B, miltefosine, and paromomycin. All the four compounds were predicted to have inhibitory constant (Ki) values ranging from 0.108 to 0.587 μM. NANPDB2581, NANPDB1649 and NANPDB1406 were also predicted as antileishmanial with Pa and Pi values of 0.415 and 0.043, 0.391 and 0.052, and 0.351 and 0.071, respectively. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) computations reinforced their good binding mechanisms. Most compounds were observed to bind in the ATP binding pocket of the kinase domain. Lys488 was predicted as a key residue critical for ligand binding in the ATP binding pocket of the LdCRK12. The molecules were pharmacologically profiled as druglike with inconsequential toxicity. The identified molecules have scaffolds that could form the backbone for fragment-based drug design of novel leishmanicides but warrant further studies to evaluate their therapeutic potential.

Discovery of human coronaviruses pan-papain-like protease inhibitors using computational approaches

The papain-like protease (PLpro) is vital for the replication of coronaviruses (CoVs), as well as for escaping innate-immune responses of the host. Hence, it has emerged as an attractive antiviral drug-target. In this study, computational approaches were employed, mainly the structure-based virtual screening coupled with all-atom molecular dynamics (MD) simulations to computationally identify specific inhibitors of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) PLpro, which can be further developed as potential pan-PLpro based broad-spectrum antiviral drugs. The sequence, structure, and functional conserveness of most deadly human CoVs PLpro were explored, and it was revealed that functionally important catalytic triad residues are well conserved among SARS-CoV, SARS-CoV-2, and middle east respiratory syndrome coronavirus (MERS-CoV). The subsequent screening of a focused protease inhibitors database composed of ∼7,000 compounds resulted in the identification of three candidate compounds, ADM_13083841, LMG_15521745, and SYN_15517940. These three compounds established conserved interactions which were further explored through MD simulations, free energy calculations, and residual energy contribution estimated by MM-PB(GB)SA method. All these compounds showed stable conformation and interacted well with the active residues of SARS-CoV-2 PLpro, and showed consistent interaction profile with SARS-CoV PLpro and MERS-CoV PLpro as well. Conclusively, the reported SARS-CoV-2 PLpro specific compounds could serve as seeds for developing potent pan-PLpro based broad-spectrum antiviral drugs against deadly human coronaviruses. Moreover, the presented information related to binding site residual energy contribution could lead to further optimization of these compounds.

Identification of novel human USP2 inhibitor and its putative role in treatment of COVID-19 by inhibiting SARS-CoV-2 papain-like (PLpro) protease

Human ubiquitin carboxyl-terminal hydrolase-2 (USP2) inhibitors, such as thiopurine analogs, have been reported to inhibit SARS-CoV papain-like proteases (PLpro). The PLpro have significant functional implications in the innate immune response during SARS-CoV-2 infection and considered an important antiviral target. Both proteases share strikingly similar USP fold with right-handed thumb-palm-fingers structural scaffold and conserved catalytic triad Cys-His-Asp/Asn. In this urgency situation of COVID-19 outbreak, there is a lack of in-vitro facilities readily available to test SARS-CoV-2 inhibitors in whole-cell assays. Therefore, we adopted an alternate route to identify potential USP2 inhibitor through integrated in-silico efforts. After an extensive virtual screening protocol, the best compounds were selected and tested. The compound Z93 showed significant IC50 value against Jurkat (9.67 μM) and MOTL-4 cells (11.8 μM). The binding mode of Z93 was extensively analyzed through molecular docking, followed by MD simulations, and molecular interactions were compared with SARS-CoV-2. The relative binding poses of Z93 fitted well in the binding site of both proteases and showed consensus π-π stacking and H-bond interactions with histidine and aspartate/asparagine residues of the catalytic triad. These results led us to speculate that compound Z93 might be the first potential chemical lead against SARS-CoV-2 PLpro, which warrants in-vitro evaluations.

A Putative Prophylactic Solution for COVID-19: Development of Novel Multiepitope Vaccine Candidate against SARS-COV-2 by Comprehensive Immunoinformatic and Molecular Modelling Approach

The outbreak of 2019-novel coronavirus (SARS-CoV-2) that causes severe respiratory infection (COVID-19) has spread in China, and the World Health Organization has declared it a pandemic. However, no approved drug or vaccines are available, and treatment is mainly supportive and through a few repurposed drugs. The urgency of the situation requires the development of SARS-CoV-2-based vaccines. Immunoinformatic and molecular modelling are time-efficient methods that are generally used to accelerate the discovery and design of the candidate peptides for vaccine development. In recent years, the use of multiepitope vaccines has proved to be a promising immunization strategy against viruses and pathogens, thus inducing more comprehensive protective immunity. The current study demonstrated a comprehensive in silico strategy to design stable multiepitope vaccine construct (MVC) from B-cell and T-cell epitopes of essential SARS-CoV-2 proteins with the help of adjuvants and linkers. The integrated molecular dynamics simulations analysis revealed the stability of MVC and its interaction with human Toll-like receptors (TLRs), which trigger an innate and adaptive immune response. Later, the in silico cloning in a known pET28a vector system also estimated the possibility of MVC expression in Escherichia coli. Despite that this study lacks validation of this vaccine construct in terms of its efficacy, the current integrated strategy encompasses the initial multiple epitope vaccine design concepts. After validation, this MVC can be present as a better prophylactic solution against COVID-19.

Structural elucidation of SARS-CoV-2 vital proteins: Computational methods reveal potential drug candidates against main protease, Nsp12 polymerase and Nsp13 helicase

Recently emerged SARS-CoV-2 caused a major outbreak of coronavirus disease 2019 (COVID-19) and instigated a widespread fear, threatening global health safety. To date, no licensed antiviral drugs or vaccines are available against COVID-19 although several clinical trials are under way to test possible therapies. During this urgent situation, computational drug discovery methods provide an alternative to tiresome high-throughput screening, particularly in the hit-to-lead-optimization stage. Identification of small molecules that specifically target viral replication apparatus has indicated the highest potential towards antiviral drug discovery. In this work, we present potential compounds that specifically target SARS-CoV-2 vital proteins, including the main protease, Nsp12 RNA polymerase and Nsp13 helicase. An integrative virtual screening and molecular dynamics simulations approach has facilitated the identification of potential binding modes and favourable molecular interaction profile of corresponding compounds. Moreover, the identification of structurally important binding site residues in conserved motifs located inside the active site highlights relative importance of ligand binding based on residual energy decomposition analysis. Although the current study lacks experimental validation, the structural information obtained from this computational study has paved way for the design of targeted inhibitors to combat COVID-19 outbreak.

Discovery of novel Hepatitis C virus inhibitor targeting multiple allosteric sites of NS5B polymerase

HCV is a viral infection posing a severe global threat when left untreated progress to end-stage liver disease, including cirrhosis and HCC. The NS5B polymerase of HCV is the most potent target that harbors four allosteric binding sites that could interfere with the HCV infection. We present the discovery of a novel synthetic compound that harbors the potential of NS5B polymerase inhibition. All eight compounds belonging to the benzothiazine family of heterocycles displayed no cellular cytotoxicity in HepG2 cells at nontoxic dose concentration (200 μM). Subsequently, among eight compounds of the series, merely compound 5b exhibited significant inhibition of the expression of the HCV NS5B gene as compared to DMSO control in semi-quantitative PCR. Based on our western blot result, 5b at the range of 50, 100 and 200 μM induced 20, 40, and 70% inhibition of NS5B protein respectively. To estimate the binding potential, 5b was docked at respective allosteric sites followed by molecular dynamics (MD) simulations for a period of 20 ns. In addition, binding free energy calculation by MM-GB/PBSA method revealed a conserved interaction profile of residues lining the allosteric sites in agreement with the reported NS5B co-crystallized inhibitors. The presented results provide important information about a novel compound 5b which may facilitate the the discovery of novel inhibitors that tends to target multiple sites on NS5B polymerase.

In Silico Studies against Viral Sexually Transmitted Diseases

Sexually Transmitted Diseases (STDs) refer to a variety of clinical syndromes and infections caused by pathogens that can be acquired and transmitted through sexual activity. Among STDs widely reported in the literature, viral sexual diseases have been increasing in a number of cases globally. This emphasizes the need for prevention and treatment. Among the methods widely used in drug planning are Computer-Aided Drug Design (CADD) studies and molecular docking which have the objective of investigating molecular interactions between two molecules to better understand the three -dimensional structural characteristics of the compounds. This review will discuss molecular docking studies applied to viral STDs, such as Ebola virus, Herpes virus and HIV, and reveal promising new drug candidates with high levels of specificity to their respective targets.

Perspectives towards antiviral drug discovery against Ebola virus

Ebola virus disease (EVD), caused by Ebola viruses, resulted in more than 11 500 deaths according to a recent 2018 WHO report. With mortality rates up to 90%, it is nowadays one of the most deadly infectious diseases. However, no Food and Drug Administration‐approved Ebola drugs or vaccines are available yet with the mainstay of therapy being supportive care. The high fatality rate and absence of effective treatment or vaccination make Ebola virus a category‐A biothreat pathogen. Fortunately, a series of investigational countermeasures have been developed to control and prevent this global threat. This review summarizes the recent therapeutic advances and ongoing research progress from research and development to clinical trials in the development of small‐molecule antiviral drugs, small‐interference RNA molecules, phosphorodiamidate morpholino oligomers, full‐length monoclonal antibodies, and vaccines. Moreover, difficulties are highlighted in the search for effective countermeasures against EVD with additional focus on the interplay between available in silico prediction methods and their evidenced potential in antiviral drug discovery.

Pharmacoinformatics-based identification of potential bioactive compounds against Ebola virus protein VP24

BACKGROUND

The impact of Ebola virus disease (EVD) is devastating with concomitant high fatalities. Currently, various drugs and vaccines are at different stages of development, corroborating the need to identify new therapeutic molecules. The VP24 protein of the Ebola virus (EBOV) plays a key role in the pathology and replication of the EVD. The VP24 protein interferes with the host immune response to viral infections and promotes nucleocapsid formation, thus making it a viable drug target. This study sought to identify putative lead compounds from the African flora with potential to inhibit the activity of the EBOV VP24 protein using pharmacoinformatics and molecular docking.

METHODS

An integrated library of 7675 natural products originating from Africa obtained from the AfroDB and NANPDB databases, as well as known inhibitors were screened against VP24 (PDB ID: 4M0Q) utilising AutoDock Vina after energy minimization using GROMACS. The top 19 compounds were physicochemically and pharmacologically profiled using ADMET Predictor™, SwissADME and DataWarrior. The mechanisms of binding between the molecules and EBOV VP24 were characterised using LigPlot+. The performance of the molecular docking was evaluated by generating a receiver operating characteristic (ROC) by screening known inhibitors and decoys against EBOV VP24. The prediction of activity spectra for substances (PASS) and machine learning-based Open Bayesian models were used to predict the anti-viral and anti-Ebola activity of the molecules, respectively.

RESULTS

Four natural products, namely, ZINC000095486070, ZINC000003594643, ZINC000095486008 and sarcophine were found to be potential EBOV VP24-inhibitiory molecules. The molecular docking results showed that ZINC000095486070 had high binding affinity of -9.7 kcal/mol with EBOV VP24, which was greater than those of the known VP24-inhibitors used as standards in the study including Ouabain, Nilotinib, Clomiphene, Torimefene, Miglustat and BCX4430. The area under the curve of the generated ROC for evaluating the performance of the molecular docking was 0.77, which was considered acceptable. The predicted promising molecules were also validated using induced-fit docking with the receptor using Schrödinger and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. The molecules had better binding mechanisms and were pharmacologically profiled to have plausible efficacies, negligible toxicity as well as suitable for designing anti-Ebola scaffolds. ZINC000095486008 and sarcophine (NANPDB135) were predicted to possess anti-viral activity, while ZINC000095486070 and ZINC000003594643 to be anti-Ebola compounds.

CONCLUSION

The identified compounds are potential inhibitors worthy of further development as EBOV biotherapeutic agents. The scaffolds of the compounds could also serve as building blocks for designing novel Ebola inhibitors.

In silico structural elucidation of RNA-dependent RNA polymerase towards the identification of potential Crimean-Congo Hemorrhagic Fever Virus inhibitors

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (-ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Inhibition of Oncogenic Kinases: An In Vitro Validated Computational Approach Identified Potential Multi-Target Anticancer Compounds

Tumorigenesis in humans is a multistep progression that imitates genetic changes leading to cell transformation and malignancy. Oncogenic kinases play a central role in cancer progression, rendering them putative targets for the design of anti-cancer drugs. The presented work aims to identify the potential multi-target inhibitors of oncogenic receptor tyrosine kinases (RTKs) and serine/threonine kinases (STKs). For this, chemoinformatics and structure-based virtual screening approaches were combined with an in vitro validation of lead hits on both cancerous and non-cancerous cell lines. A total of 16 different kinase structures were screened against ~739,000 prefiltered compounds using diversity selection, after which the top hits were filtered for promising pharmacokinetic properties. This led to the identification of 12 and 9 compounds against RTKs and STKs, respectively. Molecular dynamics (MD) simulations were carried out to better comprehend the stability of the predicted hit kinase-compound complexes. Two top-ranked compounds against each kinase class were tested in vitro for cytotoxicity, with compound F34 showing the most promising inhibitory activity in HeLa, HepG2, and Vero cell lines with IC₅₀ values of 145.46 μM, 175.48 μM, and 130.52 μM, respectively. Additional docking of F34 against various RTKs was carried out to support potential multi-target inhibition. Together with reliable MD simulations, these results suggest the promising potential of identified multi-target STK and RTK scaffolds for further kinase-specific anti-cancer drug development toward combinatorial therapies.

Ebola Virus Bayesian Machine Learning Models Enable New in Vitro Leads

We have previously described the first Bayesian machine learning models from FDA-approved drug screens, for identifying compounds active against the Ebola virus (EBOV). These models led to the identification of three active molecules in vitro: tilorone, pyronaridine, and quinacrine. A follow-up study demonstrated that one of these compounds, tilorone, has 100% in vivo efficacy in mice infected with mouse-adapted EBOV at 30 mg/kg/day intraperitoneal. This suggested that we can learn from the published data on EBOV inhibition and use it to select new compounds for testing that are active in vivo. We used these previously built Bayesian machine learning EBOV models alongside our chemical insights for the selection of 12 molecules, absent from the training set, to test for in vitro EBOV inhibition. Nine molecules were directly selected using the model, and eight of these molecules possessed a promising in vitro activity (EC50 < 15 μM). Three further compounds were selected for an in vitro evaluation because they were antimalarials, and compounds of this class like pyronaridine and quinacrine have previously been shown to inhibit EBOV. We identified the antimalarial drug arterolane (IC50 = 4.53 μM) and the anticancer clinical candidate lucanthone (IC50 = 3.27 μM) as novel compounds that have EBOV inhibitory activity in HeLa cells and generally lack cytotoxicity. This work provides further validation for using machine learning and medicinal chemistry expertize to prioritize compounds for testing in vitro prior to more costly in vivo tests. These studies provide further corroboration of this strategy and suggest that it can likely be applied to other pathogens in the future.

Molecular dynamics and combined docking studies for the identification of Zaire ebola virus inhibitors

Ebola virus (EBOV) is a lethal human pathogen with a risk of global spread of its zoonotic infections, and Ebolavirus Zaire specifically has the highest fatality rate amongst other species. There is a need for continuous effort towards having therapies, as a single licensed treatment to neutralize the EBOV is yet to come into reality. This present study virtually screened the MCULE database containing almost 36 million compounds against the structure of a Zaire Ebola viral protein (VP) 35 and a consensus scoring of both MCULE and CLCDDW docking programs remarked five compounds as potential hits. These compounds, with binding energies ranging from -7.9 to -8.9 kcal/mol, were assessed for predictions of their physicochemical and bioactivity properties, as well as absorption, distribution, metabolism, excretion, and toxicity (ADMET) criteria. The results of the 50 ns molecular dynamics simulations showed the presence of dynamic stability between ligand and protein complexes, and the structures remained significantly unchanged at the ligand-binding site throughout the simulation period. Both docking analysis and molecular dynamics simulation studies suggested strong binding affinity towards the receptor cavity and these selected compounds as potential inhibitors against the Zaire Ebola VP 35. With respect to inhibition constant values, bioavailability radar and other physicochemical properties, compound A (MCULE-1018045960-0-1) appeared to be the most promising hit compound. However, the ligand efficiency and ligand efficiency scale need improvement during optimization, and also validation via in vitro and in vivo studies are necessary to finally make a lead compound in treating Ebola virus diseases. Communicated by Ramaswamy H. Sarma.

High Throughput and Computational Repurposing for Neglected Diseases

Purpose

Neglected tropical diseases (NTDs) represent are a heterogeneous group of communicable diseases that are found within the poorest populations of the world. There are 23 NTDs that have been prioritized by the World Health Organization, which are endemic in 149 countries and affect more than 1.4 billion people, costing these developing economies billions of dollars annually. The NTDs result from four different causative pathogens: protozoa, bacteria, helminth and virus. The majority of the diseases lack effective treatments. Therefore, new therapeutics for NTDs are desperately needed.

Methods

We describe various high throughput screening and computational approaches that have been performed in recent years. We have collated the molecules identified in these studies and calculated molecular properties.

Results

Numerous global repurposing efforts have yielded some promising compounds for various neglected tropical diseases. These compounds when analyzed as one would expect appear drug-like. Several large datasets are also now in the public domain and this enables machine learning models to be constructed that then facilitate the discovery of new molecules for these pathogens.

Conclusions

In the space of a few years many groups have either performed experimental or computational repurposing high throughput screens against neglected diseases. These have identified compounds which in many cases are already approved drugs. Such approaches perhaps offer a more efficient way to develop treatments which are generally not a focus for global pharmaceutical companies because of the economics or the lack of a viable market. Other diseases could perhaps benefit from these repurposing approaches.

Electronic supplementary material

The online version of this article (10.1007/s11095-018-2558-3) contains supplementary material, which is available to authorized users.

Antiviral Agents Against Ebola Virus Infection: Repositioning Old Drugs and Finding Novel Small Molecules

Ebola virus (EBOV) causes a deadly hemorrhagic syndrome in humans with mortality rate up to 90%. First reported in Zaire in 1976, EBOV outbreaks showed a fluctuating trend during time and fora long period it was considered a tragic disease confined to the isolated regions of the African continent where the EBOV fear was perpetuated among the poor communities. The extreme severity of the recent 2014-16 EBOV outbreak in terms of fatality rate and rapid spread out of Africa led to the understanding that EBOV is a global health risk and highlights the necessity to find countermeasures against it. In the recent years, several small molecules have been shown to display in vitro and in vivo efficacy against EBOV and some of them have advanced into clinical trials. In addition, also existing drugs have been tested for their anti-EBOV activity and were shown to be promising candidates. However, despite the constant effort addressed to identify anti-EBOV therapeutics, no approved drugs are available against EBOV yet. In this chapter, we describe the main EBOV life cycle steps, providing a detailed picture of the druggable viral and host targets that have been explored so far by different technologies. We then summarize the small molecules, nucleic acid oligomers, and antibody-based therapies reported to have an effect either in in silico, or in biochemical and cell-based assays or in animal models and clinical trials, listing them according to their demonstrated or putative mechanism of action.

In silico and in vivo characterization of cabralealactone, solasodin and salvadorin in a rat model: potential anti-inflammatory agents

Background

The present study investigates the hepato- and DNA-protective effects of standardized extracts of Cleome brachycarpa (cabralealactone), Solanum incanum (solasodin), and Salvadora oleioides (salvadorin) in rats.

Materials and methods

Hepatotoxicity was induced with intraperitoneal injection of carbon tetrachloride (CCl₄) (1 mL/kg b.wt.) once a week for 12 weeks. The hepato- and DNA protective effects of the extracts in different combinations were compared with that of a standard drug Clavazin (200 mg/kg b.wt.). Tissue alanine aminotransferase, alpha-fetoprotein, tumor necrosis factor alpha (TNF-α), isoprostanes-2α, malondialdehyde, and 8-hydroxydeoxyguanosine, the significant hallmarks of oxidative stress, were studied.

Results

Histopathological findings of the liver sections from the rat group which received CCl₄+cabralealactone, solasodin, and salvadorin demonstrated improved centrilobular hepatocyte regeneration with moderate areas of congestion and infiltration comparable with Clavazin. For in silico study, the identified compounds were subjected to molecular docking with cyclooxygenase-2 and TNF-α followed by a molecular dynamics study, which indicated their potential as anti-inflammatory agents.

Conclusion

Cabralealactone, solasodin, and salvadorin confer some hepatoprotective and DNA-damage protective effects against CCl₄-induced toxicity. They successfully restored the normal architecture of hepatocytes and have the potential to be used as inhibitor to main culprits, that is, cyclooxygenase-2 and TNF-α. They can combat oxidative stress and liver injuries both as mono and combinational therapies. However, combination therapy has more ameliorating effects.

A Systematic Review of Computational Drug Discovery, Development, and Repurposing for Ebola Virus Disease Treatment

Ebola virus disease (EVD) is a deadly global public health threat, with no currently approved treatments. Traditional drug discovery and development is too expensive and inefficient to react quickly to the threat. We review published research studies that utilize computational approaches to find or develop drugs that target the Ebola virus and synthesize its results. A variety of hypothesized and/or novel treatments are reported to have potential anti-Ebola activity. Approaches that utilize multi-targeting/polypharmacology have the most promise in treating EVD.

Structure-based virtual screening and molecular docking for the identification of potential multi-targeted inhibitors against breast cancer

Breast cancer is characterized by an uncontrolled growth of cells in breast tissue. Genes that foster cell growth in breast cells are overexpressed, giving rise to breast tumors. The identification of effective inhibitors represents a rational chemopreventive strategy. The current in silico study provides a pharmacoinformatic approach for the identification of active compounds against a co-chaperone HSP90 and the human epidermal growth factor receptors EGFR and HER2/neu receptor. The elevated levels of expression of these target proteins have been documented in breast cancer. The utilization of drug-likeness filters helped to evaluate the pharmacological activity of potential lead compounds. Those fulfilling this criterion were subjected to energy minimization for 1000 steepest descent steps at a root means square gradient of 0.02 with an Amber ff12SB force field. Based on molecular docking results and binding interaction analysis, this study represents five chemical compounds (S-258282355, S-258012947, S-259417539, S-258002927, and S-259411474) that indicate high binding energies that range between -8.7 to -10.3 kcal/mol. With high cytochrome P inhibitory promiscuity activity, these multi-targeted potential hits portray not only good physiochemical interactions but also an excellent profile of absorption, distribution, metabolism, excretion, and toxicity, which hypothesizes that these compounds can be developed as anticancer drugs in the near future.