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

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

Protein-Ligand Docking in the Machine-Learning Era

Molecular docking plays a significant role in early-stage drug discovery, from structure-based virtual screening (VS) to hit-to-lead optimization, and its capability and predictive power is critically dependent on the protein-ligand scoring function. In this review, we give a broad overview of recent scoring function development, as well as the docking-based applications in drug discovery. We outline the strategies and resources available for structure-based VS and discuss the assessment and development of classical and machine learning protein-ligand scoring functions. In particular, we highlight the recent progress of machine learning scoring function ranging from descriptor-based models to deep learning approaches. We also discuss the general workflow and docking protocols of structure-based VS, such as structure preparation, binding site detection, docking strategies, and post-docking filter/re-scoring, as well as a case study on the large-scale docking-based VS test on the LIT-PCBA data set.

Syntheses, in vitro, and in silico studies of rhodanine-based schiff bases as potential α-amylase inhibitors and radicals (DPPH and ABTS) scavengers

A two-step reaction method was used to synthesize a series of rhodanine-based Schiff bases (2-33) that were characterized using spectroscopic techniques. All compounds were assessed for α-amylase inhibitory and radical scavenging (DPPH and ABTS) activities. In comparison to the standard acarbose (IC₅₀ = 9.08 ± 0.07 µM), all compounds demonstrated good to moderate α-amylase inhibitory activity (IC₅₀ = 10.91 ± 0.08-61.89 ± 0.102 µM). Compounds also demonstrated significantly higher DPPH (IC₅₀ = 10.33 ± 0.02-96.65 ± 0.03 µM) and ABTS (IC₅₀ = 12.01 ± 0.12-97.47 ± 0.13 µM) radical scavenging activities than ascorbic acid (DPPH, IC₅₀ = 15.08 ± 0.03 µM; ABTS, IC₅₀ = 16.09 ± 0.17 µM). The limited structure-activity relationship (SAR) suggests that the position and nature of the substituted groups on the phenyl ring have a vital role in varying inhibitory potential. Among the series, compounds with an electron-withdrawing group at the para position showed the highest potency. Kinetic studies revealed that the compounds followed a competitive mode of inhibition. Molecular docking results are found to agree with experimental findings, showing that compounds reside in the active pocket due to the main rhodanine moiety.

Novel Elongator Protein 2 Inhibitors Mitigating Tumor Necrosis Factor-α Induced Osteogenic Differentiation Inhibition

Tumor necrosis factor-α is a common cytokine that increases in inflammatory processes, slows the differentiation of bone formation, and induces osteodystrophy in the long-term inflammatory microenvironment. Our previous study confirmed that the Elongation protein 2 (ELP2) plays a significant role in osteogenesis and osteogenic differentiation, which is considered a drug discovery target in diseases related to bone formation and differentiation. In this study, we applied an in silico virtual screening method to select molecules that bind to the ELP2 protein from a chemical drug molecule library and obtained 95 candidates. Then, we included 11 candidates by observing the docking patterns and the noncovalent bonds. The binding affinity of the ELP2 protein with the candidate compounds was examined by SPR analysis, and 5 out of 11 compounds performed good binding affinity to the mouse ELP2 protein. After in vitro cell differentiation assay, candidates 2# and 5# were shown to reduce differentiation inhibition after tumor necrosis factor-α stimulation, allowing further optimization and development for potential clinical treatment of inflammation-mediated orthopedic diseases.

The Fellowship of Privileged Scaffolds-One Structure to Inhibit Them All

Over the past few years, the application of privileged structure has emerged as a powerful approach to the discovery of new biologically active molecules. Privileged structures are molecular scaffolds with binding properties to the range of different biological targets. Moreover, privileged structures typically exhibit good drug-like properties, thus assuring more drug-like properties of modified compound. Our main objective is to discuss the privileged structures used for the development of antiviral agents.

StackHCV: a web-based integrative machine-learning framework for large-scale identification of hepatitis C virus NS5B inhibitors

Fast and accurate identification of inhibitors with potency against HCV NS5B polymerase is currently a challenging task. As conventional experimental methods is the gold standard method for the design and development of new HCV inhibitors, they often require costly investment of time and resources. In this study, we develop a novel machine learning-based meta-predictor (termed StackHCV) for accurate and large-scale identification of HCV inhibitors. Unlike the existing method, which is based on single-feature-based approach, we first constructed a pool of various baseline models by employing a wide range of heterogeneous molecular fingerprints with five popular machine learning algorithms (k-nearest neighbor, multi-layer perceptron, partial least squares, random forest and support vectors machine). Secondly, we integrated these baseline models in order to develop the final meta-based model by means of the stacking strategy. Extensive benchmarking experiments showed that StackHCV achieved a more accurate and stable performance as compared to its constituent baseline models on the training dataset and also outperformed the existing predictor on the independent test dataset. To facilitate the high-throughput identification of HCV inhibitors, we built a web server that can be freely accessed at http://camt.pythonanywhere.com/StackHCV . It is expected that StackHCV could be a useful tool for fast and precise identification of potential drugs against HCV NS5B particularly for liver cancer therapy and other clinical applications.

The Importance of Rhodanine Scaffold in Medicinal Chemistry: A Comprehensive Overview

After the clinical use of epalrestat that contains a rhodanine ring, in type II diabetes mellitus and diabetic complications, rhodanin-based compounds have become an important class of heterocyclic in the field of medicinal chemistry. Various modifications to the rhodanine ring have led to a broad spectrum of biological activity of these compounds. Synthesis of rhodanine derivatives, depended on advenced throughput scanning hits, frequently causes potent and selective modulators of targeted enzymes or receptors, which apply their pharmacological activities through different mechanisms of action. Rhodanine-based compounds will likely stay a privileged scaffold in drug discovery because of different probability of chemical modifications of the rhodanine ring. We have, therefore reviewed their biological activities and structure activity relationship.

Diaryl Ether: A Privileged Scaffold for Drug and Agrochemical Discovery

Diaryl ether (DE) is a functional scaffold existing widely both in natural products (NPs) and synthetic organic compounds. Statistically, DE is the second most popular and enduring scaffold within the numerous medicinal chemistry and agrochemical reports. Given its unique physicochemical properties and potential biological activities, DE nucleus is recognized as a fundamental element of medicinal and agrochemical agents aimed at different biological targets. Its drug-like derivatives have been extensively synthesized with interesting biological features including anticancer, anti-inflammatory, antiviral, antibacterial, antimalarial, herbicidal, fungicidal, insecticidal, and so on. In this review, we highlight the medicinal and agrochemical versatility of the DE motif according to the published information in the past decade and comprehensively give a summary of the target recognition, structure-activity relationship (SAR), and mechanism of action of its analogues. It is expected that this profile may provide valuable guidance for the discovery of new active ingredients both in drug and pesticide research.

HCVpred: A web server for predicting the bioactivity of hepatitis C virus NS5B inhibitors

Hepatitis C virus (HCV) is one of the major causes of liver disease affecting an estimated 170 million people culminating in 300,000 deaths from cirrhosis or liver cancer. NS5B is one of three potential therapeutic targets against HCV (i.e., the other two being NS3/4A and NS5A) that is central to viral replication. In this study, we developed a classification structure-activity relationship (CSAR) model for identifying substructures giving rise to anti-HCV activities among a set of 578 non-redundant compounds. NS5B inhibitors were described by a set of 12 fingerprint descriptors and predictive models were constructed from 100 independent data splits using the random forest algorithm. The modelability (MODI index) of the data set was determined to be robust with a value of 0.88 exceeding established threshold of 0.65. The predictive performance was deduced by the accuracy, sensitivity, specificity, and Matthews correlation coefficient, which was found to be statistically robust (i.e., the former three parameters afforded values in excess of 0.8 while the latter statistical parameter provided a value >0.7). An in-depth analysis of the top 20 important descriptors revealed that aromatic ring and alkyl side chains are important for NS5B inhibition. Finally, the predictive model is deployed as a publicly accessible HCVpred web server (available at http://codes.bio/hcvpred/) that would allow users to predict the biological activity as being active or inactive against HCV NS5B. Thus, the knowledge and web server presented herein can be used in the design of more potent and specific drugs against the HCV NS5B.

Spectroscopic and In Silico DNA Binding Studies on the Interaction of Some New N-Substituted Rhodanines with Calf-thymus DNA: In Vitro Anticancer Activities

BACKGROUND

In this era of science, cancer is a black dot on the face of humankind. Consequently, the search of promising anticancer agents continues.

AIMS

Here we designed and synthesized new N-substituted rhodanines (RD1-7), evaluated their multispectroscopic interaction with calf thymus DNA, in silico and anticancer studies against MDA-MB-231cancer cell line.

METHODS

By MTT assay rhodanine RD1 was found to be the most potent with IC50 value of 72.61 μM. In addition, DNA binding studies (UV-vis and fluorescence) revealed strong binding affinity of RD1-7 with DNA (Kb in the range of 1.5-7.4 × 105 M-1). Moreover, molecular docking study, experimental DNA binding and anticancer studies are all well agreed to each other.

RESULTS

It was observed that H-bonding and hydrophobic attractions were responsible for stability of DNAcompound adducts. Besides, the reported rhodanines (RD1-7) were found as minor groove binders of DNA. Concisely, RD1-7 indicated promising pharmacological properties and hence, shows auspicious future for the development of novel anticancer agents.

CONCLUSION

The reported rhodanines showed excellent anticancer properties. Therefore, the described rhodanines may be used as potential anticancer agents in the future.

The antiviral Sofosbuvir against mucormycosis: an in silico perspective

Aim: Mucormycosis (zygomycosis) is a rare fungal infection that affects humans (40–100% mortality). Rhizopus oryzae is the primary fungus responsible for 70% of mucormycosis cases. RNA-dependent RNA polymerase (RdRp) is a vital enzyme accountable for the RNA polymerization process in different organisms, including R. oryzae. Blocking this enzyme has been previously reported as a successful strategy to eradicate viral infections. Materials & methods: AutoDock Vina is utilized for the calculation of binding affinities of Sofosbuvir, Ribavirin and uridine triphosphate nucleotide to the fungal RdRp model built by homology modeling (no solved structures available). Results: Sofosbuvir shows excellent binding affinity to the fungal RdRp in silico. Conclusion: In this study, R. oryzae RdRp is suggested to be a possible protein target against the nucleotide inhibitor, Sofosbuvir.

Anti-hepatitis-C virus activity and QSAR study of certain thiazolidinone and thiazolotriazine derivatives as potential NS5B polymerase inhibitors

The present study reports on evaluation of anti-HCV activity and QSAR of certain arylidenethiazolidinone derivatives as potential inhibitors of HCV-NS5B polymerase. The pursued compounds involving, 5-aryliden-3-arylacetamidothiazolidin-2,4-diones 4-6(a-f), 5-arylidine-2-(N-arylacetamido)-iminothiazolidin-4-one (10) and their rigid counterparts 5-arylidinethiazolotriazines 13-15(a-f), were synthesized and their structures confirmed by spectral and elemental analyses. The results of NS5B polymerase inhibition assay revealed compound 4e, as the most active inhibitor (IC₅₀ = 0.035 μM), which is four folds greater than that of the reference agent, VCH-759, (IC₅₀ = 0.14 μM). Meanwhile, compounds 4b, 4c, 5a, and 5c, and 13b, 14e and 15c displayed equipotency to 2 folds higher activity than VCH-759 (IC₅₀ values: 0.085, 0.14, 0.14, 0.10, 0.12, 0.09 and 0.07 μM, respectively). Assessment of the anti-HCV activity (GT1a) using human hepatoma cell line (Huh-7.5) illustrates superior activity of 4e (EC₅₀ = 3.80 μM) relative to VCH-759 (EC₅₀ = 5.29 μM). Cytotoxicity evaluation on, Transformed normal cell lines (Human Liver Epithelial-2, THLE-2 and Proximal Tubular Epithelial, RPTEC/TERT1), demonstrate enhanced safety profile of 4e (CC₅₀ = 102.77, 161.37 μM, respectively) compared to VCH-759 (CC₅₀ = 61.83, 81.28 μM, respectively). Molecular docking of the synthesized derivatives to NS5B polymerase allosteric site (PDB: 2HWH) showed similar binding modes to that of the co-crystallized ligand. Moreover, QSAR models were established for the studied thiazolidinones and thiazolotriazines to investigate the molecular characteristics contributing to the observed NS5B polymerase inhibition activity. The obtained results inspire further investigations of thiazolidinones and thiazolotriazine aiming at affording more potent, safe and orally active non-nucleoside NS5B polymerase inhibitors as anti-HCV drug candidates.

Synthesis and Cytotoxic Activity against K562 and MCF7 Cell Lines of SomeN-(5-Arylidene-4-oxo-2-thioxothiazolidin-3-yl)-2-((4-oxo-3-phenyl-3,4-dihydroquinazoline-2-yl)thio)acetamide Compounds

Ethyl 2-((4-oxo-3-phenyl-3,4-dihydroquinazolin-2-yl)thio)acetate (3) which was synthesized starting from anthranilic acid (1) via 2-thioxo-3-phenylquinazolin-4(3H)-one (2) reacted with hydrazine hydrate to afford 2-((4-oxo-3-phenyl-3,4-dihydroquinazolin-2-yl)thio)acetohydrazide (4). Reaction of (4) with thiocarbonyl-bis-thioglycolic acid gave a new compound nameN-(4-oxo-2-thioxothiazolidin-3-yl)-2-((4-oxo-3-phenyl-3,4-dihydroquinazolin-2-yl)thio)acetamide (5). Knoevenagel condensation of (5) with appropriate aldehydes gave fourteen (Z)-N-(5-arylidene-4-oxo-2-thioxothiazolidin-3-yl)-2-((4-oxo-3-phenyl-3,4-dihydroquinazolin-2-yl)thio)acetamide compounds (6a–o) with moderate yield. The chemical structure of the compounds was elucidated on the basis of IR,1H-NMR,13C-NMR, and HR-MS spectral data. The 5-arylidene-2-thioxothiazolidinone compounds exhibited mild-to-moderate cytotoxic activity against both K562 (chronic myelogenous leukemia) cells and MCF7 (breast cancer) cells.

Highly efficient microwave synthesis of rhodanine and 2-thiohydantoin derivatives and determination of relationships between their chemical structures and antibacterial activity

Here we report studies on the synthesis of 12 new heterocyclic derivatives that differ in three structural motifs and the simultaneous evaluation of the impact of these three variables on the biological properties. The examined compounds are based on rhodanine and 2-thiohydantoin cores equipped with hydrogen or carboxymethyl substituents at the N-3 position and linked to a triphenylamine moiety through 1,4-phenylene, 1,4-naphthalenylene and 1,9-anthracenylene spacers at the C-5 position of the heterocycles. All the compounds were synthesized very quickly, selectively and in high yields according to the developed microwave-assisted Knoevenagel condensation protocol, and they were characterized thoroughly with NMR, FT-IR and ESI-HRMS techniques. The derivatives were tested for their activity against selected strains of Gram-positive and Gram-negative bacteria and yeast. Two compounds showed good activity against Gram-positive bacteria, and all of them showed low cytotoxicity against three cell lines of the human immune system. Based on membrane permeability assays it was demonstrated that the active compounds do not penetrate the cell membrane, and thus they must act on the bacterial cell surface. Finally, we proved that the evaluated structure modifications had a synergistic effect and the simultaneous presence of a 1,4-phenylene spacer and carboxymethyl group at N-3 caused the highest boost in antimicrobial activity.

An efficient microwave-assisted synthesis of rhodanine and 2-thiohydantoin derivatives, and the correlation between their chemical structure and biological properties is reported.

Evolution of efficacious pangenotypic hepatitis C virus therapies

Hepatitis C compromises the quality of life of more than 350 million individuals worldwide. Over the last decade, therapeutic regimens for treating hepatitis C virus (HCV) infections have undergone rapid advancements. Initially, structure-based drug design was used to develop molecules that inhibit viral enzymes. Subsequently, establishment of cell-based replicon systems enabled investigations into various stages of HCV life cycle including its entry, replication, translation, and assembly, as well as role of host proteins. Collectively, these approaches have facilitated identification of important molecules that are deemed essential for HCV life cycle. The expanded set of putative virus and host-encoded targets has brought us one step closer to developing robust strategies for efficacious, pangenotypic, and well-tolerated medicines against HCV. Herein, we provide an overview of the development of various classes of virus and host-directed therapies that are currently in use along with others that are undergoing clinical evaluation.

Structure activity relationship studies on rhodanines and derived enethiol inhibitors of metallo-β-lactamases

Metallo-β-lactamases (MBLs) enable bacterial resistance to almost all classes of β-lactam antibiotics. We report studies on enethiol containing MBL inhibitors, which were prepared by rhodanine hydrolysis. The enethiols inhibit MBLs from different subclasses. Crystallographic analyses reveal that the enethiol sulphur displaces the di-Zn(II) ion bridging 'hydrolytic' water. In some, but not all, cases biophysical analyses provide evidence that rhodanine/enethiol inhibition involves formation of a ternary MBL enethiol rhodanine complex. The results demonstrate how low molecular weight active site Zn(II) chelating compounds can inhibit a range of clinically relevant MBLs and provide additional evidence for the potential of rhodanines to be hydrolysed to potent inhibitors of MBL protein fold and, maybe, other metallo-enzymes, perhaps contributing to the complex biological effects of rhodanines. The results imply that any medicinal chemistry studies employing rhodanines (and related scaffolds) as inhibitors should as a matter of course include testing of their hydrolysis products.

Quantum molecular modeling of hepatitis C virus inhibition through non-structural protein 5B polymerase receptor binding of C5-arylidene rhodanines

We have carried out high-level quantum chemical computations followed by molecular docking studies on a set of 17C₅-arylidene rhodanine isomers to provide insights into the binding modes with different reported binding pockets of the nonstructural protein 5B (NS5B) polymerase that contribute to the hepatitis C virus (HCV) inhibition. We optimized the multi-target profile of the selected rhodanine analogs to investigate potential non-nucleotide inhibitors (NNIs) by quantum chemical optimization of the 18 isomers followed by docking with quantum chemically optimized structures of each isomer with NS5B polymerase at multiple binding pockets. The binding affinities of the PP-I, PP-II and TP-II pockets of NS5B polymerase were analyzed for all the 17 isomers of 2-[(5Z)-5-(2,4-dichlorobenzylidene)-4-oxo-2-thioxo-1,3-thiazolidin-3-yl]-3-phenylpropanoic acid. On the basis of binding propensity at the different pockets and inhibitor constants, we ranked these isomers as potential candidates for the HCV inhibition. We have identified four isomers as promising NNIs of NS5B polymerase with comparable binding and inhibition to the standard (1,3) dichloro substituted isomer that exhibits in vitro activity and several other isomers as candidates in a "multi-targeted drug" approach.

5-Ene-4-thiazolidinones - An efficient tool in medicinal chemistry

The presented review is an attempt to summarize a huge volume of data on 5-ene-4-thiazolidinones being a widely studied class of small molecules used in modern organic and medicinal chemistry. The manuscript covers approaches to the synthesis of 5-ene-4-thiazolidinone derivatives: modification of the C5 position of the basic core; synthesis of the target compounds in the one-pot or multistage reactions or transformation of other related heterocycles. The most prominent pharmacological profiles of 5-ene derivatives of different 4-thiazolidinone subtypes belonging to hit-, lead-compounds, drug-candidates and drugs as well as the most studied targets have been discussed. Currently target compounds (especially 5-en-rhodanines) are assigned as frequent hitters or pan-assay interference compounds (PAINS) within high-throughput screening campaigns. Nevertheless, the crucial impact of the presence/nature of C5 substituent (namely 5-ene) on the pharmacological effects of 5-ene-4-thiazolidinones was confirmed by the numerous listed findings from the original articles. The main directions for active 5-ene-4-thiazolidinones optimization have been shown: i) complication of the fragment in the C5 position; ii) introduction of the substituents in the N3 position (especially fragments with carboxylic group or its derivatives); iii) annealing in complex heterocyclic systems; iv) combination with other pharmacologically attractive fragments within hybrid pharmacophore approach. Moreover, the utilization of 5-ene-4-thiazolidinones in the synthesis of complex compounds with potent pharmacological application is described. The chemical transformations cover mainly the reactions which involve the exocyclic double bond in C5 position of the main core and correspond to the abovementioned direction of the 5-ene-4-thiazolidinone modification.

Recent developments with rhodanine as a scaffold for drug discovery

INTRODUCTION

Rhodanines, as one of the 4-thiazolidinones subtypes, are recognized as privileged heterocycles in medicinal chemistry. The main achievements include the development of drug-like molecules with numerous biological activities as well as approved drugs. Among rhodanines, 5-ene-rhodanines are of special interest, and are often claimed as pan assay interference compounds due to Michael acceptor functionality. Areas covered: Herein, the synthetic protocols for rhodanines and their transformation are reviewed. Biological activity is briefly discussed as well as biotargets, mode of actions and optimization directions. Furthermore, the utilization of 5-ene-rhodanines in Michael additions are discussed while both pro and contra arguments have been outlined within medicinal chemistry application. Expert opinion: Rhodanines remain privileged heterocycles in drug discovery. They are accessible building blocks for optimization and transformation into related heterocycles, simplified analogues and fused heterocycles with a thiazolidine framework. Michael acceptor functionality, as well as the thesis about low selectivity towards biotargets of rhodanines, must be confirmed experimentally and it cannot be based on just the presence of conjugated α,β-unsaturated carbonyl. Moreover, the positive aspects of Michael acceptors must be considered as well as their multitarget properties. New criteria for target affinity must be found. In conclusion, rhodanines are generally not problematic per se.

Molecular modeling and docking revealed superiority of IDX-184 as HCV polymerase inhibitor

Aim: IDX-184 is a nonstructural 5b nucleoside inhibitor (NI) that was under clinical trials against HCV. This work adopts a molecular modeling approach in order to study the interaction between IDX-184 and HCV polymerase from four different genotypes. Methods: Comparisons to the native nucleotide (Guanosine triphosphate) and other NIs were performed using interaction descriptors, calculated using semiempirical quantum mechanics and molecular docking. Results: IDX-184 shows potent binding to the active site of the polymerases. In addition, IDX-184 was better than Sofosbuvir and Ribavirin when docked into polymerase active site (even with experimentally solved structure). Conclusion: Analysis of the interaction descriptors and docking complexes suggests IDX-184 as a superior NI against the studied HCV subtypes.

Applications of computer-aided approaches in the development of hepatitis C antiviral agents

INTRODUCTION

Hepatitis C virus (HCV) is a global health problem that causes several chronic life-threatening liver diseases. The numbers of people affected by HCV are rising annually. Since 2011, the FDA has approved several anti-HCV drugs; while many other promising HCV drugs are currently in late clinical trials. Areas covered: This review discusses the applications of different computational approaches in HCV drug design. Expert opinion: Molecular docking and virtual screening approaches have emerged as a low-cost tool to screen large databases and identify potential small-molecule hits against HCV targets. Ligand-based approaches are useful for filtering-out compounds with rich physicochemical properties to inhibit HCV targets. Molecular dynamics (MD) remains a useful tool in optimizing the ligand-protein complexes and understand the ligand binding modes and drug resistance mechanisms in HCV. Despite their varied roles, the application of in-silico approaches in HCV drug design is still in its infancy. A more mature application should aim at modelling the whole HCV replicon in its active form and help to identify new effective druggable sites within the replicon system. With more technological advancements, the roles of computer-aided methods are only going to increase several folds in the development of next-generation HCV drugs.

Computational Approaches for the Discovery of Novel Hepatitis C Virus NS3/4A and NS5B Inhibitors

Nonstructural 5B (NS5B) polymerase and Nonstructural 3/4A (NS3/4A) protease have proven to be promising targets for the development of anti-HCV (Hepatitis C Virus) agents. The NS5B polymerase is of paramount importance in HCV viral replication; therefore, employing NS5B inhibitors was considered an effective way for the treatment of HCV. Identifying inhibitors against NS3/4A serine protease represents another attractive approach applied in anti-HCV drug discovery, which is evidenced by its crucial role of in the biogenesis of the viral replication activity. In this chapter, many different computational approaches including Quantitative Structure-Activity Relationship (QSAR) and virtual screening in anti-HCV drug discovery were considered and discussed in detail. Virtual Screening (VS) techniques, including ligand-based and structure-based, and QSAR have been utilized for the discovery of NS5B inhibitors. Moreover, using various in silico protocols and workflows, a number of studies have been conducted with an aim of identifying potential NS3/4A blockage agents.

5-Carba-pterocarpens: A new scaffold with anti-HCV activity

The synthesis of a series of 5-carba-pterocarpens derivatives involving the cyclization of α-aryl-α-tetralones is described. Several compounds demonstrated potent activity and selectivity in vitro against HCV replicon reporter cells. The best profile in Huh7/Rep-Feo1b replicon reporter cells was observed with 2h (EC50 = 5.5 μM/SI = 20), while 2e was the most active in Huh7.5-FGR-JC1-Rluc2A replicon reporter cells (EC50 = 1.5 μM/SI = 70). Hydroxy groups at A- and D-rings are essential for anti-HCV activity, and substitutions in the A-ring at positions 3 and 4 resulted in enhanced activity of the compounds.

Discovery of Novel Hepatitis C Virus NS5B Polymerase Inhibitors by Combining Random Forest, Multiple e-Pharmacophore Modeling and Docking

The NS5B polymerase is one of the most attractive targets for developing new drugs to block Hepatitis C virus (HCV) infection. We describe the discovery of novel potent HCV NS5B polymerase inhibitors by employing a virtual screening (VS) approach, which is based on random forest (RB-VS), e-pharmacophore (PB-VS), and docking (DB-VS) methods. In the RB-VS stage, after feature selection, a model with 16 descriptors was used. In the PB-VS stage, six energy-based pharmacophore (e-pharmacophore) models from different crystal structures of the NS5B polymerase with ligands binding at the palm I, thumb I and thumb II regions were used. In the DB-VS stage, the Glide SP and XP docking protocols with default parameters were employed. In the virtual screening approach, the RB-VS, PB-VS and DB-VS methods were applied in increasing order of complexity to screen the InterBioScreen database. From the final hits, we selected 5 compounds for further anti-HCV activity and cellular cytotoxicity assay. All 5 compounds were found to inhibit NS5B polymerase with IC50 values of 2.01-23.84 μM and displayed anti-HCV activities with EC50 values ranging from 1.61 to 21.88 μM, and all compounds displayed no cellular cytotoxicity (CC50 > 100 μM) except compound N2, which displayed weak cytotoxicity with a CC50 value of 51.3 μM. The hit compound N2 had the best antiviral activity against HCV, with a selective index of 32.1. The 5 hit compounds with new scaffolds could potentially serve as NS5B polymerase inhibitors through further optimization and development.

Recent advances on the synthesis of hepatitis C virus NS5B RNA-dependent RNA-polymerase inhibitors

Hepatitis C is a viral liver infection considered as the major cause of cirrhosis and hepatocellular carcinoma (HCC). The HCV NS5B polymerase, an RNA-dependent RNA polymerase, is essential for HCV replication, which is able to catalyze the synthesis of positive (genomic) and negative (template) strand HCV RNA, but has no functional equivalent in mammalian cells. Therefore, the NS5B polymerase has emerged as an attractive target for the development of specifically targeted antiviral therapy for HCV (DAA, for direct-acting antivirals). Recently, a growing number of compounds have been reported as the NS5B polymerase inhibitors, some of which especially have been licensed in clinical trials. This review describes recent advances on the synthesis of the NS5B polymerase inhibitors, focusing on the merits and demerits of their synthetic methods. In particular, inspiration from the synthesis and the future direction of the NS5B polymerase inhibitors are highlighted.

Multistage virtual screening and identification of novel HIV-1 protease inhibitors by integrating SVM, shape, pharmacophore and docking methods

The HIV-1 protease has proven to be a crucial component of the HIV replication machinery and a reliable target for anti-HIV drug discovery. In this study, we applied an optimized hierarchical multistage virtual screening method targeting HIV-1 protease. The method sequentially applied SVM (Support Vector Machine), shape similarity, pharmacophore modeling and molecular docking. Using a validation set (270 positives, 155,996 negatives), the multistage virtual screening method showed a high hit rate and high enrichment factor of 80.47% and 465.75, respectively. Furthermore, this approach was applied to screen the National Cancer Institute database (NCI), which contains 260,000 molecules. From the final hit list, 6 molecules were selected for further testing in an in vitro HIV-1 protease inhibitory assay, and 2 molecules (NSC111887 and NSC121217) showed inhibitory potency against HIV-1 protease, with IC50 values of 62 μM and 162 μM, respectively. With further chemical development, these 2 molecules could potentially serve as HIV-1 protease inhibitors.

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.

Synthesis and antibacterial activity of analogs of 5-arylidene-3-(4-methylcoumarin-7-yloxyacetylamino)-2-thioxo-1,3-thiazoli-din-4-one

In an effort to develop new antimicrobial agents, 3-(4-methylcoumarin-7-yloxyacetylamino)-2-thioxo-1,3-thiazolidin-4-one (4) was synthesized by reaction of thiocarbonylbisthioglycolic acid with ethyl (4-methyl-2-oxo-2H-chromen-7-yloxy)aceto- hydrazide (3), which was prepared in turn from 7-hydroxy-4-methylcoumarin (1). The condensation of compound 4 with different aromatic aldehydes afforded a series of 5-(arylidene)-3-(4-methylcoumarin-7-yloxyacetyl-amino)-2-thioxo-1,3-thiozolidin-4-one analogs 5a–h. The structures of these synthetic compounds were elucidated on the basis of IR, ¹H-NMR and ¹³C-NMR spectral data and ESI-MS spectrometric analysis. Compounds 5a–h were examined for their antibacterial activity against several strains of Gram-positive and Gram-negative bacteria.

2-Heteroarylimino-5-arylidene-4-thiazolidinones as a new class of non-nucleoside inhibitors of HCV NS5B polymerase

Hepatitis C virus (HCV) NS5B polymerase is an important and attractive target for the development of anti-HCV drugs. Here we report on the design, synthesis and evaluation of twenty-four novel allosteric inhibitors bearing the 4-thiazolidinone scaffold as inhibitors of HCV NS5B polymerase. Eleven compounds tested were found to inhibit HCV NS5B with IC₅₀ values ranging between 19.8 and 64.9 μM. Compound 24 was the most active of this series with an IC₅₀ of 5.6 μM. A number of these derivatives further exhibited strong inhibition against HCV 1b and 2a genotypes in cell based antiviral assays. Molecular docking analysis predicted that the thiazolidinone derivatives bind to the NS5B thumb pocket-II (TP-II). Our results suggest that further optimization of the thiazolidinone scaffold may be possible to yield new derivatives with improved enzyme- and cell-based activity.

Molecular modeling comparison of the performance of NS5b polymerase inhibitor (PSI-7977) on prevalent HCV genotypes

The current available treatment for hepatitis C virus (HCV)-the causative of liver cirrhosis and development of liver cancer-is a dual therapy using modified interferon and ribavirin. While this regimen increases the sustained viral response rate up to 40-80 % in different genotypes, unfortunately, it is poorly tolerated by patients. PSI-7977, a prodrug for PSI-7409, is a Non-Structural 5b (NS5b) polymerase nucleoside inhibitor that is currently in phase III clinical trials. The activated PSI-7977 is a direct acting antiviral (DAA) drug that acts on NS5b polymerase of HCV through a coordination bond with the two Mg(+2) present at the GDD active site motif. The present work utilizes a molecular modeling approach for studying the interaction between the activated PSI-7977 and the 12 amino acids constituting a 5 Å region surrounding the GDD active triad motif for HCV genotypes 1a, 2b, 3b and 4a. The analysis of the interaction parameters suggests that PSI-7977 is probably a better DAA drug for HCV genotypes 1a and 3b rather than genotypes 2b and 4a.

Design and synthesis of L- and D-phenylalanine derived rhodanines with novel C5-arylidenes as inhibitors of HCV NS5B polymerase

Hepatitis C virus (HCV) NS5B polymerase is a key target for anti-HCV therapeutics development. Herein, we report the synthesis and in vitro evaluation of anti-NS5B polymerase activity of a molecular hybrid of our previously reported lead compounds 1 (IC50=7.7 μM) and 2 (IC50=10.6 μM) as represented by hybrid compound 27 (IC50=6.7 μM). We have explored the optimal substituents on the terminal phenyl ring of the 3-phenoxybenzylidene moiety in 27, by generating a set of six analogs. This resulted in the identification of compound 34 with an IC50 of 2.6 μM. To probe the role of stereochemistry towards the observed biological activity, we synthesized and evaluated the D-isomers 41 (IC50=19.3 μM) and 45 (IC50=5.4 μM) as enantiomers of the l-isomers 27 and 34, respectively. The binding site of compounds 32 and 34 was mapped to palm pocket-I (PP-I) of NS5B. The docking models of 34 and 45 within the PP-I of NS5B were investigated to envisage the molecular mechanism of inhibition.

Discovery of new scaffolds for rational design of HCV NS5B polymerase inhibitors

Hepatitis C virus (HCV) NS5B polymerase is a key target for the development of anti-HCV drugs. Here we report on the identification of novel allosteric inhibitors of HCV NS5B through a combination of structure-based virtual screening and in vitro NS5B inhibition assays. One hundred and sixty thousand compounds from the Otava database were virtually screened against the thiazolone inhibitor binding site on NS5B (thumb pocket-2, TP-2), resulting in a sequential down-sizing of the library by 2.7 orders of magnitude to yield 59 NS5B non-nucleoside inhibitor (NNI) candidates. In vitro evaluation of the NS5B inhibitory activity of the 59 selected compounds resulted in a 14% hit rate, yielding 8 novel structural scaffolds. Of these, compound 1 bearing a 4-hydrazinoquinazoline scaffold was the most active (IC(50) = 16.0 μM). The binding site of all 8 NNIs was mapped to TP-2 of NS5B as inferred by a decrease in their inhibition potency against the M423T NS5B mutant, employed as a screen for TP-2 site binders. At 100 μM concentration, none of the eight compounds exhibited any cytotoxicity, and all except compound 8 exhibited between 40 and 60% inhibition of intracellular NS5B polymerase activity in BHK-NS5B-FRLuc reporter cells. These inhibitor scaffolds will form the basis for future optimization and development of more potent NS5B inhibitors.