Discovery of Potent Neuraminidase Inhibitors Using a Combination of Pharmacophore-Based Virtual Screening and Molecular Simulation Approach

Identification of NS2B-NS3 Protease Inhibitors for Therapeutic Application in ZIKV Infection: A Pharmacophore-Based High-Throughput Virtual Screening and MD Simulations Approaches

Zika virus (ZIKV) pandemic and its implication in congenital malformations and severe neurological disorders had created serious threats to global health. ZIKV is a mosquito-borne flavivirus which spread rapidly and infect a large number of people in a shorter time-span. Due to the lack of effective therapeutics, this had become paramount urgency to discover effective drug molecules to encounter the viral infection. Various anti-ZIKV drug discovery efforts during the past several years had been unsuccessful to develop an effective cure. The NS2B-NS3 protein was reported as an attractive therapeutic target for inhibiting viral proliferation, due to its central role in viral replication and maturation of non-structural viral proteins. Therefore, the current in silico drug exploration aimed to identify the novel inhibitors of Zika NS2B-NS3 protease by implementing an e-pharmacophore-based high-throughput virtual screening. A 3D e-pharmacophore model was generated based on the five-featured (ADPRR) pharmacophore hypothesis. Subsequently, the predicted model is further subjected to the high-throughput virtual screening to reveal top hit molecules from the various small molecule databases. Initial hits were examined in terms of binding free energies and ADME properties to identify the candidate hit exhibiting a favourable pharmacokinetic profile. Eventually, molecular dynamic (MD) simulations studies were conducted to evaluate the binding stability of the hit molecule inside the receptor cavity. The findings of the in silico analysis manifested affirmative evidence for three hit molecules with -64.28, -55.15 and -50.16 kcal/mol binding free energies, as potent inhibitors of Zika NS2B-NS3 protease. Hence, these molecules holds the promising potential to serve as a prospective candidates to design effective drugs against ZIKV and related viral infections.

Discovery of novel thiophene derivatives as potent neuraminidase inhibitors

Neuraminidase (NA) is an important target for the treatment of influenza. In this study, a new lead NA inhibitor, 4 (ZINC01121127), was discovered by pharmacophore-based virtual screening and molecular dynamic (MD) simulation. Some novel NA inhibitors containing thiophene ring were synthesized by optimizing the skeleton of the lead compound 4. Compound 4b had the most potent inhibitory activity against NA (IC₅₀ = 0.03 μM), which was better than the positive control oseltamivir carboxylate (IC₅₀ = 0.06 μM). 4b (EC₅₀ = 1.59 μM) also exhibits excellent antiviral activity against A/chicken/Hubei/327/2004 (H5N1-DW), which is superior to the reference drug OSC (EC₅₀ = 5.97 μM). Molecular docking study shows that the thiophene moiety plays an essential role in compound 4b, which can bind well to the active site of NA. The good activity of 4b may be also ascribed to the extending of quinoline ring into the 150-cavity. The results of this study may provide an insightful help for the development of new NA inhibitors.

Inhibition of α-glucosidase and α-amylase by herbal compounds for the treatment of type 2 diabetes: A validation of in silico reverse docking with in vitro enzyme assays

BACKGROUND

α-Amylase and α-glucosidase are important therapeutic targets for the management of type 2 diabetes mellitus. The inhibition of these enzymes decreases postprandial hyperglycemia. In the present study, compounds found in commercially available herbs and spices were tested for their ability to inhibit α-amylase and α-glucosidase. These compounds were acetyleugenol, apigenin, cinnamic acid, eriodictyol, myrcene, piperine, and rosmarinic acid.

METHODS

The enzyme inhibitory nature of the compounds was evaluated using in silico docking analysis with Maestro software and was further confirmed by in vitro α-amylase and α-glucosidase biochemical assays.

RESULTS

The relationships between the in silico and in vitro results were well correlated; a more negative docking score was associated with a higher in vitro inhibitory activity. There was no significant (P > .05) difference between the inhibition constant (K_i ) value of acarbose, a widely prescribed α-glucosidase and α-amylase inhibitor, and those of apigenin, eriodictyol, and piperine. For α-amylase, there was no significant (P > .05) difference between the K_i value of acarbose and those of apigenin, cinnamic acid, and rosmarinic acid. The effect of the herbal compounds on cell viability was assessed with the sulforhodamine B (SRB) assay in C2C12 and HepG2 cells. Acetyleugenol, cinnamic acid, myrcene, piperine, and rosmarinic acid had similar (P > .05) IC₅₀ values to acarbose.

CONCLUSIONS

Several of the herbal compounds studied could regulate postprandial hyperglycemia. Using herbal plants has several advantages including low cost, natural origin, and easy cultivation. These compounds can easily be consumed as teas or as herbs and spices to flavor food.

Synthesis and bio-evaluation of a novel selective butyrylcholinesterase inhibitor discovered through structure-based virtual screening

In recent years, butyrylcholinesterase (BChE) has gradually gained worldwide interests as a novel target for treating Alzheimer's disease (AD). Here, two pharmacophore models were generated using Schrödinger suite and used to virtually screen ChemDiv database, from which three hits were obtained. Among them, 2513-4169 displayed the highest inhibitory activity and selectivity against BChE (eeAChE IC₅₀ > 10 μM, eqBChE IC₅₀ = 3.73 ± 1.90 μM). Molecular dynamic (MD) simulation validated the binding pattern of 2513-4169 in BChE, and it could form a various of receptor-ligand interactions with adjacent residues. In vitro cytotoxicity assay proved the safety of 2513-4169 on diverse neural cell lines. Moreover, the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) assay performed on SH-SY5Y cells proved the neuroprotective effect of 2513-4169 against toxic Aβ_1-42. In vivo behavioral study further confirmed the great efficacy of 2513-4169 on reversing Aβ_1-42-induced cognitive impairment of mice and clearing the toxic Aβ_1-42 in brains. Moreover, 2513-4169 was proved to be able to cross blood-brain barrier (BBB) through a parallel artificial membrane permeation assay of BBB (PAMPA-BBB). Taken together, 2513-4169 is a promising lead compound for future optimization to discover anti-AD treating agents.

Discovery of Potential Flavonoid Inhibitors Against COVID-19 3CL Proteinase Based on Virtual Screening Strategy

The outbreak of 2019 novel coronavirus (COVID-19) has caused serious threat to public health. Discovery of new anti-COVID-19 drugs is urgently needed. Fortunately, the crystal structure of COVID-19 3CL proteinase was recently resolved. The proteinase has been identified as a promising target for drug discovery in this crisis. Here, a dataset including 2030 natural compounds was screened and refined based on the machine learning and molecular docking. The performance of six machine learning (ML) methods of predicting active coronavirus inhibitors had achieved satisfactory accuracy, especially, the AUC (Area Under ROC Curve) scores with fivefold cross-validation of Logistic Regression (LR) reached up to 0.976. Comprehensive ML prediction and molecular docking results accounted for the compound Rutin, which was approved by NMPA (National Medical Products Administration), exhibited the best AUC and the most promising binding affinity compared to other compounds. Therefore, Rutin might be a promising agent in anti-COVID-19 drugs development.

Density Functional Theory and Molecular Simulation Studies for Prioritizing Anaplastic Lymphoma Kinase Inhibitors

Targeting anaplastic lymphoma kinase (ALK) is one of the important treatment strategies for the treatment of non-small cell lung cancer (NSCLC). In the present perspective, multidimensional approaches were used for the identification of ALK inhibitors. Initially, an e-pharmacophore model was generated using the PHASE algorithm and was used as a 3D query to screen 468,200 molecules of ASINEX database. Prior to the screening process, the model was evaluated for its significance and the ability to differentiate actives from inactives, using enrichment analysis. Subsequently, the hierarchical docking protocol and binding free energy calculations were instigated using GLIDE algorithm and Prime module, respectively. Further, the pharmacokinetic/pharmacodynamics (PK/PD) properties and toxicities of the hit compounds were envisaged respectively using QikProp program, Osiris explorer, and Protox-II algorithm. These approaches retrieved two hits namely BAS 00137817 and BAS 00680055 with acceptable absorption, distribution, metabolism, excretion and toxicity (ADMET) properties and higher affinity towards ALK protein. Additionally, density functional theory calculations and molecular dynamics simulations were performed to validate the inhibitory activity of the lead compounds. It is noteworthy to mention that all the hits constitute of particular scaffolds which play a major role in the downregulation of some ALK-positive lung cancer pathways. We speculate that the outcomes of this research are of substantial prominence in the rational designing of novel and efficacious ALK inhibitors.

E-pharmacophore hypothesis strategy to discover potent inhibitor for influenza treatment

The surface protein of Influenza virus, Neuraminidase (NA), is believed to play a critical role in the release of new viral particle and thus spreads infection. It has been recognized as a valid drug target for anti-influenza therapy. Despite the number of available approved drugs for the influenza infection treatment, the emergence of resistant variants with novel mutations are the foremost challenges for the currently used NA inhibitors. Thus, the current investigation was carried out to ascertain potent inhibitors using computational strategies such as e-pharmacophore based virtual screening and docking approach. A three-dimensional e-pharmacophore hypothesis was generated based on the chemical features of complexes of the drugs and NA protein using PHASE module of Schrödinger suite. The generated hypothesis consisted of one hydrogen bond acceptor (A), two hydrogen bond donors (D), one negatively charged group (N) and one aromatic ring (R), ADDNR. The hypothesis was further evaluated for its integrity using enrichment analysis and used to filter out molecules with similar pharmacophoric features from approved, investigational and experimental subsets of DrugBank and ZINC database. In addition, ligand filtration was performed to curb down the molecules to an efficient collection of hit molecules by using Lipinski “rule of five and ADME analysis by using Qikprop module. Overall, the results from our analysis suggest that compound lisinopril and formoterol could serve as potent antiviral compounds for the treatment of influenza A virus infection. It is worth mentioning that the results correlate well with literature evidences.

Discovery of novel anaplastic lymphoma kinase inhibitors: Structure and energy-based pharmacophore strategy

The clinical outcomes in patients with non-small cell lung cancer have improved, as a result of anaplastic lymphoma kinase (ALK) inhibition. Therefore in the current study, substantial effort has been made to identify ALK inhibitors through systematic virtual screening experiment consisting of e-pharmacophore and pharmacophore perception techniques. Initially, a pharmacophore model (AAAHP.193) and an e-pharmacophore model (DDRRR) encompassing the whole dataset of 12 known ALK inhibitors were developed. The hypotheses could retrieve effective compounds from DrugBank database (8621 molecules), which were then subjected to molecular docking and ADME prediction. These approaches resulted in the identification of five hits, namely, nebivolol, HDY, D42, 796, and LZE having higher Glide docking scores and promising ADME properties with augmented CNS involvement. Moreover, molecular dynamics simulations were performed to validate the inhibitory activity of the hit compounds, and density functional theory calculations were carried out to scrutinize the chemical reactivity of the hits. Subsequent interaction and scaffold analysis identified prominent interactions of the hits with ALK kinase domain and scaffolds with anti-tumor activity against lung cancer cell lines. We strongly believe that the study provides an outlook for the sighting of novel and potent ALK inhibitors in the near future.

Exploring the Lead Compounds for Zika Virus NS2B-NS3 Protein: an e-Pharmacophore-Based Approach

The rapid spread of the Zika virus and its association with the abnormal brain development constitute a global health emergency. With a continuing spread of the mosquito vector, the exposure is expected to accelerate in the coming years. Despite number of efforts, there is still no proper vaccine or medicine to combat this virus. Of note, the NS2B-NS3 protein is proven to be the potential target for the Zika virus therapeutics. Hence, e-pharmacophore-based drug design strategy was employed to identify potent inhibitors of NS2B-NS3 protein from ASINEX database consisting of 467,802 molecules. A 3D e-pharmacophore model was generated using PHASE module of Schrödinger Suite. The generated model consists of one hydrogen bond acceptor (A), two hydrogen bond donors (D), and two aromatic rings (R), ADDRR. The model was further evaluated for its ability to screen actives using enrichment analysis. Subsequently, high-throughput virtual screening protocol was employed, and the resultant hit molecules were also examined for its binding free energies and ADME properties using Prime MM-GBSA and Qikprop module of Schrodinger packages, respectively. Finally, the screened hit molecule was subjected to molecular dynamics simulation to examine its stability. Overall, the results from our analysis suggest that compound BAS 19192837 could be a potent inhibitor for the NS2B-NS3 protein of the Zika virus. It is also noteworthy to mention that our results are in good agreement with literature evidences. We hope that this result is of immense importance in designing potential drug molecules to combat the spread of Zika virus in the near future.

Taxifolin as dual inhibitor of Mtb DNA gyrase and isoleucyl-tRNA synthetase: in silico molecular docking, dynamics simulation and in vitro assays

DNA gyrase and aminoacyl-tRNA synthetases are two essential bacterial enzymes involved in DNA replication, transcription and translation. Flavonoids are plant secondary metabolites with variable phenolic structures. In this study, eight flavonoids structurally similar to quercetin were selected and their ADMET properties were evaluated. Molecular docking and free energy calculations were carried out to examine the binding of these flavonoids to the ATP-binding site and editing domain of DNA gyrase and Isoleucyl-tRNA synthetase, respectively. Taxifolin was found out to be the top lead molecule in both the docking studies with a good number of interactions with the active site amino acids. Further, binding of taxifolin to the proteins was extensively studied using 50 ns molecular dynamics simulation. In vitro anti-tuberculosis activity of taxifolin was evaluated and compared with the standard drugs. Minimal inhibition concentration of taxifolin was found to be ≤ 12.5 μg/ml.