Identification of potential inhibitors for HCV NS5b of genotype 4a by combining dynamic simulation,protein–ligand interaction fingerprint, 3D pharmacophore,docking and 3D QSAR

Abstract

HCV NS5B polymerase has been one of the most attractive targets for developing new drugs for HCV infection and many drugs were successfully developed, but all of them were designed for targeting Hepatitis C Virus genotype 1 (HCV GT1). Hepatitis C virus genotype 4a (HCV GT4a) dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS5B polymerase of GT4a using homology modeling, protein–ligand interaction fingerprint (PLIF), docking, pharmacophore, and 3D CoMFA quantitative structure activity relationship (QSAR). Firstly, a high-quality 3D model of HCV NS5B polymerase of GT4a was constructed using crystal structure of HCV NS5B polymerase of GT1 (PDB ID: 3hkw) as a template. Then, both the model and the template were simulated to compare conformational stability. PLIF was generated using five crystal structures of HCV NS5B (PDB ID: 4mia, 4mib, 4mk9, 4mka, and 4mkb), which revealed the most important residues and their interactions with the co-crystalized ligands. After that, a 3D pharmacophore model was developed from the generated PLIF data and then used as a screening filter for 17000328 drug-like zinc database compounds. 900 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. Finally, a 3D CoMFA QSAR was developed using 42 compounds as a training and 19 compounds as a test set. The 3D CoMFA QSAR was used to design and screen some potential inhibitors, these compounds were further evaluated by the docking stage. The highest ranked five hits from docking result (compounds (p1–p4) and compound q1) were selected for further analysis.

Communicated by Ramaswamy H. Sarma     

Discovery of hit molecules targeting allosteric site of hepatitis C virus NS5B polymerase

Abstract

Nonstructural protein 5B (NS5B), the RNA-dependent RNA polymerase of Hepatitis C Virus (HCV), plays a key role in viral amplification and is an attractive and most explored target for discovery of new therapeutic agents for Hepatitis C. Though safe and effective, NS5B inhibitors were launched in 2013 (Sovaldi) and 2014 (Harvoni, Viekira Pak), the high price tags of these medications limit their use among poor people in developing countries. Hence, still there exists a need for cost-effective and short duration anti-HCV agents especially those targeting niche patient population who were non-respondent to earlier therapies or with comorbid conditions. The present study describes the discovery of novel non-nucleoside (NNI) inhibitors of NS5B using a series of rational drug design techniques such as virtual screening, scaffold matching and molecular docking. 2D and 3D structure based virtual screening technique identified 300 hit compounds. Top 20 hits were screened out from identified hits using molecular docking technique. Four molecules, that are representative of 20 hits were evaluated for binding affinity under in vitro conditions using surface plasmon resonance-based assay and the results emphasized that compound with CoCoCo ID: 412075 could exhibit good binding response toward NS5B and could be a potential candidate as NS5B inhibitor.

Communicated by Ramaswamy H. Sarma     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Searching for anthranilic acid-based thumb pocket 2 HCV NS5B polymerase inhibitors through a combination of molecular docking, 3D-QSAR and virtual screening

A combination of the following computational methods: (i) molecular docking, (ii) 3-D Quantitative Structure Activity Relationship Comparative Molecular Field Analysis (3D-QSAR CoMFA), (iii) similarity search and (iv) virtual screening using PubChem database was applied to identify new anthranilic acid-based inhibitors of hepatitis C virus (HCV) replication. A number of known inhibitors were initially docked into the “Thumb Pocket 2” allosteric site of the crystal structure of the enzyme HCV RNA-dependent RNA polymerase (NS5B GT1b). Then, the CoMFA fields were generated through a receptor-based alignment of docking poses to build a validated and stable 3D-QSAR CoMFA model. The proposed model can be first utilized to get insight into the molecular features that promote bioactivity, and then within a virtual screening procedure, it can be used to estimate the activity of novel potential bioactive compounds prior to their synthesis and biological tests.     

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Abstract

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer’s test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski’s rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375).

Communicated by Ramaswamy H. Sarma     

Structure-based virtual screening,synthesis and SAR of novel inhibitors of hepatitis C virus NS5B polymerase

Hepatitis C virus (HCV) NS5B polymerase is a key target for the development of therapeutic agents aimed at the treatment of HCV infections. Here we report on the identification of novel allosteric inhibitors of HCV NS5B through a combination of structure-based virtual screening, synthesis and structure–activity relationship (SAR) optimization approach. Virtual screening of 260,000 compounds from the ChemBridge database against the tetracyclic indole inhibitor binding pocket of NS5B (allosteric pocket-1, AP-1), sequentially down-sized the library by 4 orders of magnitude to yield 23 candidates. In vitro evaluation of the NS5B inhibitory activity of the in-silico selected compounds resulted in 17% hit rate, identifying two novel chemotypes. Of these, compound 3, bearing the rhodanine scaffold, proved amenable for productive SAR exploration and synthetic modification. As a result, 25 derivatives that exhibited IC₅₀ values ranging from 7.7 to 68.0 μM were developed. Docking analysis of lead compound 28 within the tetracyclic indole- and benzylidene-binding allosteric pockets (AP-1 and AP-3, respectively) of NS5B revealed topological similarities between these two pockets. Compound 28, a novel rhodanine analog with NS5B inhibitory potency in the low micromolar level range may be a promising lead for future development of more potent NS5B inhibitors.     

Identification of potential tumour-associated carbonic anhydrase isozyme IX inhibitors: atom-based 3D-QSAR modelling,pharmacophore-based virtual screening and molecular docking studies

Abstract

Tumour hypoxia results in dramatic changes in the gene expression, proliferation and survival of tumour cells. The tumour cells shift towards anaerobic glycolysis which results in change of pH in their microenvironment. In response to this stress, over expression of carbonic anhydrase IX (CA IX) genes is observed in many solid tumours. So, selective inhibition of CA IX can be a promising target for anti-cancer drugs. In this work in silico tools like atom-based 3D-QSAR modelling, pharmacophore-based virtual screening and molecular docking were used to identify potential CA IX inhibitors. Based on the training set used in the QSAR model, twenty pharmacophore models were generated. Out of these, HHHR_1, AHHR_1, DHHHR_1, AHHHR_1 model was used to screen a database of 1,50,000 compounds retrieved from ZINC 15 database. R² and Q² was 0.9864 and 0.8799, respectively, for the developed QSAR model. 163 compounds showed a phase screen score above 2.4 in which ZINC02260669 was the highest ranked (screen score, 2.852058) compound in all the four models. Built QSAR model was used to predict the activity of all these 163 compounds and ZINC72370966 showed the highest predicted activity with pKi value of 7.649. These compounds were docked against CA IX (human) protein (PDB ID 5FL6) and molecular docking results showed favourable binding interactions for the best ten identified hits. This work gives design insights and some potential scaffolds which can be developed as CA IX inhibitors.

Communicated by Ramaswamy H. Sarma     

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with ?259.41, ?110.09, ?188.25, ?163.30, ?202.10, and ?158.79 kJ mol^?1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.     

Structure-based design, synthesis, and biological evaluation of 1,1-dioxoisothiazole and benzo[b]thiophene-1,1-dioxide derivatives as novel inhibitors of hepatitis C virus NS5B polymerase

A novel series of HCV NS5B polymerase inhibitors comprising 1,1-dioxoisothiazoles and benzo[b]thiophene-1,1-dioxides were designed, synthesized, and evaluated. SAR studies guided by structure-based design led to the identification of a number of potent NS5B inhibitors with nanomolar IC₅₀ values. The most potent compound exhibited IC₅₀ less than 10 nM against the genotype 1b HCV polymerase and EC₅₀ of 70 nM against a genotype 1b replicon in cell culture. The DMPK properties of selected compounds were also evaluated.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Structural and energetic understanding of novel natural inhibitors of Mycobacterium tuberculosis malate synthase

Persistent infection by Mycobacterium tuberculosis requires the glyoxylate shunt. This is a bypass to the tricarboxylic acid cycle in which isocitrate lyase (ICL) and malate synthase (MS) catalyze the net incorporation of carbon during mycobacterial growth on acetate or fatty acids as the primary carbon source. To identify a potential antitubercular compound, we performed a structure-based screening of natural compounds from the ZINC database (n = 1 67 740) against the M tuberculosis MS (MtbMS) structure. The ligands were screened against MtbMS, and 354 ligands were found to have better docking score. These compounds were assessed for Lipinski and absorption, distribution, metabolism, excretion, and toxicity prediction where 15 compounds were found to fit well for redocking studies. After refinement by molecular docking and drug-likeness analysis, four potential inhibitors (ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522) were identified. These four ligands with phenyl-diketo acid were further subjected to molecular dynamics simulation to compare the dynamics and stability of the protein structure after ligand binding. The binding energy analysis was calculated to determine the intermolecular interactions. Our results suggested that the four compounds had a binding free energy of −201.96, −242.02, −187.03, and −169.02 kJ·mol⁻¹, for compounds with IDs ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522, respectively. We concluded that two compounds (ZINC1483899 and ZINC1754310) displayed considerable structural and pharmacological properties and could be probable drug candidates to fight against M tuberculosis parasites.     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Design,synthesis, and structure-activity relationships of novel imidazo[4,5-c]pyridine derivatives as potent non-nucleoside inhibitors of hepatitis C virus NS5B

The hepatitis C virus (HCV) NS5B polymerase is an attractive target for the development of novel and selective inhibitors of HCV replication. In this paper, the design, synthesis, and preliminary SAR studies of novel inhibitors of HCV NS5B polymerase based on the structure of tegobuvir have been described. The efforts to optimize the antiviral potency and reduce the treatment side effects with respect to genotype 1b resulted in the discovery of compound 3, which exhibited an EC₅₀ of 1.163?nM and a CC₅₀ >200?nM in a cell-based HCV replicon system assay. Additionally, testing for inhibition of the hERG channel showed a marked improvement over tegobuvir and the pharmacokinetic properties of compound 3 indicated that it was worthy of further investigation as a non-nucleoside inhibitor of HCV NS5B polymerase.     

Targeting the NTPase site of Zika virus NS3 helicase for inhibitor discovery

Abstract

The major threats linked to Zika virus (ZIKV) are microcephaly, Guillain-Barre syndrome, and the ability to transfer through sexual transmission. Despite these threats, Zika specific FDA approved drugs or vaccines are not available as of yet. Additionally, the involvement of pregnant women makes the drug screening process lengthy and complicated in terms of safety and minimum toxicity of the molecules. Since NS3 helicase of ZIKV performs the critical function of unwinding double-stranded RNA during replication, it is considered as a promising drug target to block ZIKV replication. In the present study, we have exploited the NTPase site of ZIKV NS3 helicase for screening potential inhibitor compounds by molecular docking, and molecular dynamics (MD) simulation approaches. NS3 helicase hydrolyzes the ATP to use its energy for unwinding RNA. We have chosen twenty natural compounds from ZINC library with known antiviral properties and a helicase focused library (HFL) of small molecules from Life Chemicals compounds. After going through docking, the top hit molecules from ZINC and HFL library were further analysed by MD simulations to find out stable binding poses. Finally, we have reported the molecules with potential of binding at NTPase pocket of ZIKV NS3 helicase, which could be further tested on virus through in vitro experiments to check their efficacy.

Communicated by Ramaswamy H. Sarma     

Non-nucleoside inhibitors of HCV polymerase NS5B. Part 2: Synthesis and structure–activity relationships of benzothiazine-substituted quinolinediones

A new series of benzothiazine-substituted quinolinediones were evaluated as inhibitors of HCV polymerase NS5B. SAR studies on this series revealed a methyl sulfonamide group as a high affinity feature. Analogues with this group showed submicromolar potencies in the HCV cell based replicon assay. Pharmacokinetic and toxicology studies were also performed on a selected compound (34) to evaluate in vivo properties of this new class of inhibitors of HCV NS5B polymerase.     

Design and evaluation of novel tetracyclic benzofurans as palm site allosteric inhibitors of HCV NS5B polymerase

Hepatitis C virus (HCV) NS5B polymerase is a prime target for the development of direct-acting antiviral drugs for the treatment of chronic HCV infection. Several novel and potent HCV NS5B non-nucleoside inhibitors with unique tetracyclic bezonfuran-based structures were prepared and evaluated. Similar to clinical developmental compound MK-8876, N-linked (compounds 1 and 2) and C-linked (compounds 3 and 4) tetracyclic structures maintained broad spectrum anti-replicon potency profiles and demonstrated moderate to excellent oral bioavailability and pharmacokinetic parameters across the three preclinical animal species. To better understand the importance of tetracyclic structures related to pan genotypic potency profiles especially against clinically relevant GT1a variants, the teracycles with different ring size were prepared and in vitro evaluations suggested compounds with six number ring have better overall potency profiles.