Novel HCV NS5B polymerase inhibitors derived from 4-(1',1'-dioxo-1',4'-dihydro-1'lambda(6)-benzo[1',2',4']thiadiazin-3'-yl)-5-hydroxy-2H-pyridazin-3-ones. Part 3: Further optimization of the 2-, 6-, and 7'-substituents and initial pharmacokinetic assessments

In-silico Investigations for the Identification of Novel Inhibitors Targeting Hepatitis C Virus RNA-dependent RNA Polymerase

BACKGROUND

Hepatitis C is an inflammatory condition of the liver caused by the hepatitis C virus, exhibiting acute and chronic manifestations with severity ranging from mild to severe and lifelong illnesses leading to liver cirrhosis and cancer. According to the World Health Organization's global estimates, a population of about 58 million have chronic hepatitis C virus infection, with around 1.5 million new infections occurring every year.

OBJECTIVE

The present study aimed to identify novel molecules targeting the Hepatitis C viral RNA Dependent RNA polymerases, which play a crucial role in genome replication, mRNA synthesis, etc. Methods: Structure-based virtual screening of chemical libraries of small molecules was done using AutoDock/Vina. The top-ranking pose for every ligand was complexed with the protein and used for further protein-ligand interaction analysis using the Protein-ligand interaction Profiler. Molecules from virtual screening were further assessed using the pkCSM web server. The proteinligand interactions were further subjected to molecular dynamics simulation studies to establish dynamic stability.

RESULTS

Molecular docking-based virtual screening of the database of small molecules, followed by screening based on pharmacokinetic and toxicity parameters, yielded eight probable RNA Dependent RNA polymerase inhibitors. The docking scores for the proposed candidates ranged from - 8.04 to -9.10 kcal/mol. The potential stability of the ligands bound to the target protein was demonstrated by molecular dynamics simulation studies.

CONCLUSION

Data from exhaustive computational studies proposed eight molecules as potential anti-viral candidates, targeting Hepatitis C viral RNA Dependent RNA polymerases, which can be further evaluated for their biological potential.

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

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.

A study of the allosteric inhibition of HCV RNA-dependent RNA polymerase and implementing virtual screening for the selection of promising dual-site inhibitors with low resistance potential

Structure-based pharmacophores were generated and validated using the bioactive conformations of different co-crystallized enzyme-inhibitor complexes for allosteric palm-1 and thumb-2 inhibitors of NS5B. Two pharmacophore models were obtained, one for palm-1 inhibitors with sensitivity = 0.929 and specificity = 0.983, and the other for thumb-2 inhibitors with sensitivity = 1 and specificity = 0.979. In addition, a quantitative structure activity relationship (QSAR) models were developed based on using the values of different scoring functions as descriptors predicting the activity on both allosteric binding sites (palm-1 and thumb-2). QSAR studies revealed good predictive and statistically significant two descriptor models (r²= .837, r²_adjusted = .792 and r²_prediction = .688 for palm-1 model and r²= .927, r²_adjusted = .908 and r²_prediction = .779 for thumb-2 model). External validation for the QSAR models assured their prediction power with r²_ext= .72 and .89 for palm-1 and thumb-2, respectively. Different docking protocols were examined for their validity to predict the correct binding poses of inhibitors inside their respective binding sites. Virtual screening was carried out on ZINC database using the generated pharmacophores, the selected valid docking algorithms and QSAR models to find compounds that could theoretically bind to both sites simultaneously.

Quantitative structure and bioactivity relationship study on HCV NS5B polymerase inhibitors

Several QSAR (quantitative structure-activity relationship) models for predicting the inhibitory activity of 333 hepatitis C virus (HCV) NS5B polymerase inhibitors were developed. All the inhibitors are HCV polymerase non-nucleoside analogue inhibitors (NNIs) fitting into the pocket of the NNI III binding site. For each molecule, global descriptors and 2D property autocorrelation descriptors were calculated from the program ADRIANA.Code. Pearson correlation analysis was used to select the significant descriptors for building models. The whole dataset was split into a training set and a test set randomly or using a Kohonen's self-organizing map (SOM). Then, the inhibitory activity of 333 HCV NS5B polymerase inhibitors was predicted using multilinear regression (MLR) analysis and support vector machine (SVM) method, respectively. For the test set of the best model (Model 2B), correlation coefficient of 0.91 was achieved. Some molecular descriptors, such as molecular complexity (Complexity), the number of hydrogen bonding donors (HDon) and the solubility of the molecule in water (log S) were found to be very important factors which determined the bioactivity of the HCV NS5B inhibitors. Some other molecular properties such as electrostatic and charge properties also played important roles in the interaction between the ligand and the protein. The selected molecular descriptors were further confirmed by analysing the interaction between two representative inhibitors and the polymerase in their crystal structures.

Combined 3D-QSAR, molecular docking, molecular dynamics simulation, and binding free energy calculation studies on the 5-hydroxy-2H-pyridazin-3-one derivatives as HCV NS5B polymerase inhibitors

The NS5B RNA-dependent RNA polymerase (RdRP) is a promising therapeutic target for developing novel anti-hepatitis C virus (HCV) drugs. In this work, a combined molecular modeling study was performed on a series of 193 5-hydroxy-2H-pyridazin-3-one derivatives as inhibitors of HCV NS5B Polymerase. The best 3D-QSAR models, including CoMFA and CoMSIA, are based on receptor (or docking). Furthermore, a 40-ns molecular dynamics (MD) simulation and binding free energy calculations using docked structures of NS5B with ten compounds, which have diverse structures and pIC50 values, were employed to determine the detailed binding process and to compare the binding modes of the inhibitors with different activities. On one side, the stability and rationality of molecular docking and 3D-QSAR results were validated by MD simulation. The binding free energies calculated by the MM-PBSA method gave a good correlation with the experimental biological activity. On the other side, by analyzing some differences between the molecular docking and the MD simulation results, we can find that the MD simulation could also remedy the defects of molecular docking. The analyses of the combined molecular modeling results have identified that Tyr448, Ser556, and Asp318 are the key amino acid residues in the NS5B binding pocket. The results from this study can provide some insights into the development of novel potent NS5B inhibitors.

Design and Development of NS5B Polymerase Non‐nucleoside Inhibitors for the Treatment of Hepatitis C Virus Infection

The hepatitis C virus (HCV) infects an estimated 130–170 million people worldwide and is associated with life‐threatening liver diseases. The recent introduction of the first two HCV direct‐acting antivirals (DAAs) as a complement to the interferon/ribavirin standard of care has provided patients with improved outcomes. Still, 25–30% of subjects infected with genotype 1 HCV do not respond adequately to treatment owing to the emergence of resistant virus and many suffer from severe side effects. A paradigm shift towards the development of interferon‐free combinations of DAAs with complementary modes of action is currently taking place. Virally encoded proteins and enzymes have become the target of HCV drug discovery efforts and several promising new agents are currently being evaluated in the clinic for treatment of chronic HCV infection. The NS5B RNA‐dependent RNA polymerase is responsible for replication of viral RNA and plays a pivotal role in the virus life cycle. NS5B is undoubtedly the most druggable HCV target and is susceptible to several classes of allosteric inhibitors that bind to four distinct sites on the enzyme. This chapter describes successful strategies that have led to the discovery of HCV NS5B antivirals. It is divided according to allosteric sites and describes how each of the known families of inhibitors was discovered, characterized and optimized to provide clinical candidates. When available, the strategies adopted by medicinal chemists to optimize initial leads and address challenges and liabilities encountered on the path to candidate selection are described, along with reported clinical outcomes.

Inhibitors for the hepatitis C virus RNA polymerase explored by SAR with advanced machine learning methods

Hepatitis C virus (HCV) is a global health challenge, affecting approximately 200 million people worldwide. In this study we developed SAR models with advanced machine learning classifiers Random Forest and k Nearest Neighbor Simulated Annealing for 679 small molecules with measured inhibition activity for NS5B genotype 1b. The activity was expressed as a binary value (active/inactive), where actives were considered molecules with IC50 ≤0.95 μM. We applied our SAR models to various drug-like databases and identified novel chemical scaffolds for NS5B inhibitors. Subsequent in vitro antiviral assays suggested a new activity for an existing prodrug, Candesartan cilexetil, which is currently used to treat hypertension and heart failure but has not been previously tested for anti-HCV activity. We also identified NS5B inhibitors with two novel non-nucleoside chemical motifs.

Classification of HCV NS5B polymerase inhibitors using support vector machine

Using a support vector machine (SVM), three classification models were built to predict whether a compound is an active or weakly active inhibitor based on a dataset of 386 hepatitis C virus (HCV) NS5B polymerase NNIs (non-nucleoside analogue inhibitors) fitting into the pocket of the NNI III binding site. For each molecule, global descriptors, 2D and 3D property autocorrelation descriptors were calculated from the program ADRIANA.Code. Three models were developed with the combination of different types of descriptors. Model 2 based on 16 global and 2D autocorrelation descriptors gave the highest prediction accuracy of 88.24% and MCC (Matthews correlation coefficient) of 0.789 on test set. Model 1 based on 13 global descriptors showed the highest prediction accuracy of 86.25% and MCC of 0.732 on external test set (including 80 compounds). Some molecular properties such as molecular shape descriptors (InertiaZ, InertiaX and Span), number of rotatable bonds (NRotBond), water solubility (LogS), and hydrogen bonding related descriptors performed important roles in the interactions between the ligand and NS5B polymerase.

Recent advances in drug discovery of benzothiadiazine and related analogs as HCV NS5B polymerase inhibitors

Hepatitis C virus (HCV) is a major health burden, with an estimated 170 million chronically infected individuals worldwide, and a leading cause of liver transplantation. Patients are at increased risk of developing liver cirrhosis, hepatocellular carcinoma and even liver failure. In the past two decades, several approaches have been adopted to inhibit non-structural viral proteins. The RNA-dependent RNA polymerase (NS5B) of HCV is one of the attractive validated targets for development of new drugs to block HCV infection. In this review, we report the recent progress made towards identifying and developing benzothiadiazines as HCV NS5B polymerase inhibitors. The substituted benzothiadiazine class was identified by HTS in 2002 as an NS5B inhibitor. Further optimization and modification of the core has improved the potency and pharmacokinetic properties of substituted benzothiadiazines. Research on palm site-binding benzothiadiazine analogs and related derivatives and analogs is discussed in this article.

Ligand-induced changes in hepatitis C virus NS5B polymerase structure

Hepatitis C virus (HCV) RNA-dependent RNA polymerase (NS5B) is required for viral replication. Crystal structures of the NS5B apoprotein show that the finger and thumb domains interact to encircle the active site, and that inhibitors defined by P495 resistance that bind to the thumb-finger interface displace the Δ1 finger loop and disrupt this structure. Since crystal structures may not reveal all of the conformations of a protein in solution we have developed an alternative method, using limited trypsin protease digestion, to investigate the impact of inhibitors as well as substrates on the movement of the Δ1 loop. This assay can be used to study NS5B under conditions that support enzymatic activity. In the absence of inhibitors, no specific region of NS5B was hypersensitive to trypsin, and no specific intermediate cleavage products were formed. Binding of P495-site inhibitors to NS5B induced specific trypsin hypersensitivity at lysine residues 50 and 51. Previously characterized inhibitors and mutant polymerases were used to link this specific trypsin hypersensitivity to movement of the Δ1 loop. Trypsin hypersensitivity identical to the inhibitor pattern was also induced by the binding of the RNA template. The addition of primer to the NS5B-template complex eliminated the hypersensitivity. The data are consistent with displacement of the Δ1 finger loop from the thumb by the binding of template, and reversal by the addition of primer or NTP. Our results complement inhibitor-enzyme co-crystal studies, and the assay provides a rapid and sensitive method to study dynamic changes in HCV NS5B polymerase conformation under conditions that support functional activity.

Non-nucleoside inhibitors of hepatitis C virus polymerase: current progress and future challenges

The current standard of care for hepatitis C virus (HCV) infection is a combination of PEGylated interferon and ribavirin, which offer limited efficacy and significant side effects. Novel HCV-specific inhibitors, including those directed at the viral polymerase, have become the focus of HCV drug-discovery efforts in the past decade. In addition to the active site targeted by traditional nucleoside inhibitors, at least four different allosteric-binding sites have been reported for the HCV polymerase, which offer ample opportunities for small-molecule inhibitors. In this review, we summarize the recent progress in the discovery of non-nucleoside HCV polymerase inhibitors with a focus on novel chemical matters, their clinical efficacy, safety and potential for combination therapy.

Discovery of tricyclic 5,6-dihydro-1H-pyridin-2-ones as novel, potent, and orally bioavailable inhibitors of HCV NS5B polymerase

A novel series of non-nucleoside small molecules containing a tricyclic dihydropyridinone structural motif was identified as potent HCV NS5B polymerase inhibitors. Driven by structure-based design and building on our previous efforts in related series of molecules, we undertook extensive SAR studies, in which we identified a number of metabolically stable and very potent compounds in genotype 1a and 1b replicon assays. This work culminated in the discovery of several inhibitors, which combined potent in vitro antiviral activity against both 1a and 1b genotypes, metabolic stability, good oral bioavailability, and high C(12) (PO)/EC(50) ratios.

5,5'- and 6,6'-dialkyl-5,6-dihydro-1H-pyridin-2-ones as potent inhibitors of HCV NS5B polymerase

The discovery of 5,5'- and 6,6'-dialkyl-5,6-dihydro-1H-pyridin-2-ones as potent inhibitors of the HCV RNA-dependent RNA polymerase (NS5B) is described. Several of these agents also display potent antiviral activity in cell culture experiments (EC50 <0.10 microM). In vitro DMPK data for selected compounds as well as crystal structures of representative inhibitors complexed with the NS5B protein are also disclosed.

5,6-Dihydro-1H-pyridin-2-ones as potent inhibitors of HCV NS5B polymerase

5,6-Dihydro-1H-pyridin-2-one analogs were discovered as a novel class of inhibitors of genotype 1 HCV NS5B polymerase. Among these, compound 4ad displayed potent inhibitory activities in biochemical and replicon assays (IC(50) (1b)<10nM; IC(50) (1a)<25nM, EC(50) (1b)=16nM), good in vitro DMPK properties, as well as moderate oral bioavailability in monkeys (F=24%).

Novel HCV NS5B polymerase inhibitors derived from 4-(1',1'-dioxo-1',4'-dihydro-1'lambda(6)-benzo[1',2',4']thiadiazin-3'-yl)-5-hydroxy-2H-pyridazin-3-ones. Part 5: Exploration of pyridazinones containing 6-amino-substituents

The synthesis of 4-(1',1'-dioxo-1',4'-dihydro-1'lambda(6)-benzo[1',2',4']thiadiazin-3'-yl)-5-hydroxy-2H-pyridazin-3-ones bearing 6-amino substituents as potent inhibitors of the HCV RNA-dependent RNA polymerase (NS5B) is described. Several of these agents also display potent antiviral activity in cell culture experiments (EC(50)<0.10 microM). In vitro DMPK data (microsome t(1/2), Caco-2 P(app)) for many of the compounds are also disclosed, and a crystal structure of a representative inhibitor complexed with the NS5B protein is discussed.

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(50) values. The most potent compound exhibited IC(50) less than 10nM against the genotype 1b HCV polymerase and EC(50) of 70 nM against a genotype 1b replicon in cell culture. The DMPK properties of selected compounds were also evaluated.

Pyrrolo[1,2-b]pyridazin-2-ones as potent inhibitors of HCV NS5B polymerase

Pyrrolo[1,2-b]pyridazin-2-one analogs were discovered as a novel class of inhibitors of genotype 1 HCV NS5B polymerase. Structure-based design led to the discovery of compound 3 k, which displayed potent inhibitory activities in biochemical and replicon assays (IC(50) (1b)<10nM; EC(50) (1b)=12 nM) as well as good stability towards human liver microsomes (HLM t(1/2)>60 min).