Novel HCV NS5B polymerase inhibitors derived from 4-(1',1'-dioxo-1',4'-dihydro-1'lambda6-benzo[1',2',4']thiadiazin-3'-yl)-5-hydroxy-2H-pyridazin-3-ones. Part 1: exploration of 7'-substitution of benzothiadiazine

An Efficient Modern Strategy to Screen Drug Candidates Targeting RdRp of SARS-CoV-2 With Potentially High Selectivity and Specificity

Desired drug candidates should have both a high potential binding chance and high specificity. Recently, many drug screening strategies have been developed to screen compounds with high possible binding chances or high binding affinity. However, there is still no good solution to detect whether those selected compounds possess high specificity. Here, we developed a reverse DFCNN (Dense Fully Connected Neural Network) and a reverse docking protocol to check a given compound’s ability to bind diversified targets and estimate its specificity with homemade formulas. We used the RNA-dependent RNA polymerase (RdRp) target as a proof-of-concept example to identify drug candidates with high selectivity and high specificity. We first used a previously developed hybrid screening method to find drug candidates from an 8888-size compound database. The hybrid screening method takes advantage of the deep learning-based method, traditional molecular docking, molecular dynamics simulation, and binding free energy calculated by metadynamics, which should be powerful in selecting high binding affinity candidates. Also, we integrated the reverse DFCNN and reversed docking against a diversified 102 proteins to the pipeline for assessing the specificity of those selected candidates, and finally got compounds that have both predicted selectivity and specificity. Among the eight selected candidates, Platycodin D and Tubeimoside III were confirmed to effectively inhibit SARS-CoV-2 replication in vitro with EC₅₀ values of 619.5 and 265.5 nM, respectively. Our study discovered that Tubeimoside III could inhibit SARS-CoV-2 replication potently for the first time. Furthermore, the underlying mechanisms of Platycodin D and Tubeimoside III inhibiting SARS-CoV-2 are highly possible by blocking the RdRp cavity according to our screening procedure. In addition, the careful analysis predicted common critical residues involved in the binding with active inhibitors Platycodin D and Tubeimoside III, Azithromycin, and Pralatrexate, which hopefully promote the development of non-covalent binding inhibitors against RdRp.

An Integrated Deep Learning and Molecular Dynamics Simulation-Based Screening Pipeline Identifies Inhibitors of a New Cancer Drug Target TIPE2

The TIPE2 (tumor necrosis factor-alpha-induced protein 8-like 2) protein is a major regulator of cancer and inflammatory diseases. The availability of its sequence and structure, as well as the critical amino acids involved in its ligand binding, provides insights into its function and helps greatly identify novel drug candidates against TIPE2 protein. With the current advances in deep learning and molecular dynamics simulation-based drug screening, large-scale exploration of inhibitory candidates for TIPE2 becomes possible. In this work, we apply deep learning-based methods to perform a preliminary screening against TIPE2 over several commercially available compound datasets. Then, we carried a fine screening by molecular dynamics simulations, followed by metadynamics simulations. Finally, four compounds were selected for experimental validation from 64 candidates obtained from the screening. With surprising accuracy, three compounds out of four can bind to TIPE2. Among them, UM-164 exhibited the strongest binding affinity of 4.97 µM and was able to interfere with the binding of TIPE2 and PIP2 according to competitive bio-layer interferometry (BLI), which indicates that UM-164 is a potential inhibitor against TIPE2 function. The work demonstrates the feasibility of incorporating deep learning and MD simulation in virtual drug screening and provides high potential inhibitors against TIPE2 for drug development.

A novel virtual screening procedure identifies Pralatrexate as inhibitor of SARS-CoV-2 RdRp and it reduces viral replication in vitro

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus poses serious threats to the global public health and leads to worldwide crisis. No effective drug or vaccine is readily available. The viral RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target. A hybrid drug screening procedure was proposed and applied to identify potential drug candidates targeting RdRp from 1906 approved drugs. Among the four selected market available drug candidates, Pralatrexate and Azithromycin were confirmed to effectively inhibit SARS-CoV-2 replication in vitro with EC50 values of 0.008μM and 9.453 μM, respectively. For the first time, our study discovered that Pralatrexate is able to potently inhibit SARS-CoV-2 replication with a stronger inhibitory activity than Remdesivir within the same experimental conditions. The paper demonstrates the feasibility of fast and accurate anti-viral drug screening for inhibitors of SARS-CoV-2 and provides potential therapeutic agents against COVID-19.

Application of carbamyl in structural optimization

Carbamyl is considered a privileged structure in medicinal chemistry. It has a wide range of biological activities such as antimicrobial, anticancer, anti-epilepsy, for which the best evidence is a number of marketed carbamyl-containing drugs. Carbamyl is formed of primary amine and carbonyl moieties that act as hydrogen bond donors and hydrogen acceptors with residues of targets respectively, which are benefit for improving pharmacological activities. In other cases, the introduced carbamyl improves drug-like properties including oral bioavailability. In this review, we introduce the carbamyl-containing drugs and the application of carbamyl in structural optimization as a result of enhancing activities or/and drug-like properties.

Discovering protein-ligand chalcogen bonding in the protein data bank using endocyclic sulfur-containing heterocycles as ligand search subsets

The chalcogen bond, the noncovalent, electrostatic attraction between covalently bonded atoms in group 16 and Lewis bases, is present in protein-ligand interactions based on X-ray structures deposited in the Protein Data Bank (PDB). Discovering protein-ligand chalcogen bonding in the PDB employed a strategy that focused on searching the database for protein complexes of five-membered, heterocyclic ligands containing endocyclic sulfur with endo electron-withdrawing groups (isothiazoles; thiazoles; 1,2,3-, 1,2.4-, 1,2,5-, 1,3,4-thiadiazoles) and thiophenes with exo electron-withdrawing groups, e.g., 2-chloro, 2-bromo, 2-amino, 2-alkylthio. Out of 930 ligands investigated, 33 or 3.5% have protein-ligand S---O interactions of which 31 are chalcogen bonds and two appear to be S---HO hydrogen bonds. The bond angles for some of the chalcogen bonds found in the PDB are less than 90°, and an electrostatic model is proposed to explain this phenomenon.

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.

In vitro activity and resistance profile of dasabuvir, a nonnucleoside hepatitis C virus polymerase inhibitor

Dasabuvir (ABT-333) is a nonnucleoside inhibitor of the RNA-dependent RNA polymerase encoded by the hepatitis C virus (HCV) NS5B gene. Dasabuvir inhibited recombinant NS5B polymerases derived from HCV genotype 1a and 1b clinical isolates, with 50% inhibitory concentration (IC50) values between 2.2 and 10.7 nM, and was at least 7,000-fold selective for the inhibition of HCV genotype 1 polymerases over human/mammalian polymerases. In the HCV subgenomic replicon system, dasabuvir inhibited genotype 1a (strain H77) and 1b (strain Con1) replicons with 50% effective concentration (EC50) values of 7.7 and 1.8 nM, respectively, with a 13-fold decrease in inhibitory activity in the presence of 40% human plasma. This level of activity was retained against a panel of chimeric subgenomic replicons that contained HCV NS5B genes from 22 genotype 1 clinical isolates from treatment-naive patients, with EC50s ranging between 0.15 and 8.57 nM. Maintenance of replicon-containing cells in medium containing dasabuvir at concentrations 10-fold or 100-fold greater than the EC50 resulted in selection of resistant replicon clones. Sequencing of the NS5B coding regions from these clones revealed the presence of variants, including C316Y, M414T, Y448C, Y448H, and S556G, that are consistent with binding to the palm I site of HCV polymerase. Consequently, dasabuvir retained full activity against replicons known to confer resistance to other polymerase inhibitors, including the S282T variant in the nucleoside binding site and the M423T, P495A, P495S, and V499A single variants in the thumb domain. The use of dasabuvir in combination with inhibitors targeting HCV NS3/NS4A protease (ABT-450 with ritonavir) and NS5A (ombitasvir) is in development for the treatment of HCV genotype 1 infections.

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.

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.

Efficient in silico assay of inhibitors of hepatitis C Virus RNA-dependent RNA polymerase by structure-based virtual screening and in vitro evaluation

To identify a new protective or therapeutic intervention for hepatitis C virus (HCV) infection, we performed efficient structure-based virtual screening to identify novel inhibitory agents for HCV. To this end, we selected NS5B, an RNA-dependent RNA polymerase (RdRp), as the target for the treatment of HCV infection. To decipher the dockable nature of various RdRp X-ray crystals, we docked the crystal ligand (inhibitor) to the crystal receptor (enzyme). The accuracy of regeneration of the crystal pose indicates the amenability of the RdRp binding pocket for structure-based virtual screening. We also utilized a consensus scoring scheme to reduce false positives, thereby ensuring efficient virtual screening. In this study, each molecule that ranked in the top 1% among all screening molecules gained 1 consensus point in a scoring function. Thus, after virtual screening of 57,177 chemicals from the Maybridge Screening collection, 14 molecules gained 8 points across 11 scoring functions. One of them, an isoxazole, showed significant dose-dependent inhibition of HCV RdRp activity and replication. In this study, we have developed a structure-based virtual screening method using HCV RdRp for efficient identification of novel inhibitors.

Synthetic utility of ammonium salts in a Cu-catalyzed three-component reaction as a facile coupling partner

Ammonium salts were found to be a convenient and inexpensive reagent in the Cu-catalyzed three-component reaction with terminal alkynes and sulfonyl or phosphoryl azides leading to N-unprotected amidines. Thus obtained amidines bearing 2-bromobenzenesulfonyl moiety were efficiently cyclized by the Cu-catalyzed intramolecular N-arylation to give an important pharmacophore skeleton of 2H-1,2,4-benzothiadiazine 1,1-dioxides. Conveniently, two tandem catalytic procedures could be readily operated in one pot.

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).

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

5-Hydroxy-3(2H)-pyridazinone derivatives were investigated as inhibitors of genotype 1 HCV NS5B polymerase. Lead optimization led to the discovery of compound 3a, which displayed potent inhibitory activities in biochemical and replicon assays [IC(50) (1b)<10nM; IC(50) (1a)=22 nM; EC(50) (1b)=5nM], good stability toward human liver microsomes (HLM t(1/2)>60 min), and high ratios of liver to plasma concentrations 12h after a single oral administration to rats.