Design,synthesis and biological evaluation of novel zanamivir derivatives as potent neuraminidase inhibitors

Neuraminidase (NA) is an important antiviral drug target. Zanamivir is one of the most potent NA inhibitors. In this paper, a series of zanamivir derivatives as potential NA inhibitors were studied by combination of molecular modeling techniques including 3D-QSAR, molecular docking, and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q²?=?0.728 and r²?=?0.988, and the best CoMSIA (comparative molecular similarity indices analysis) model has q²?=?0.750 and r²?=?0.981, respectively. The built 3D-QSAR models show significant statistical quality and excellent predictive ability. Seven new NA inhibitors were designed and predicted. 20?ns of MD simulations were carried out and their binding free energies were calculated. Two designed compounds were selected to be synthesized and biologically evaluated by NA inhibition and virus inhibition assays. One compound (IC₅₀?=?0.670?µM, SI?>?149) exhibits excellent antiviral activity against A/WSN/33 H1N1, which is superior to the reference drug zanamivir (IC₅₀?=?0.873?µM, SI?>?115). The theoretical and experimental results may provide reference for development of new anti-influenza drugs.     

Design,synthesis and biological evaluation of novel oseltamivir derivatives as potent neuraminidase inhibitors

Neuraminidase (NA) is one of the particular potential targets for novel antiviral therapy. In this work, a series of neuraminidase inhibitors with the cyclohexene scaffold were studied based upon the combination of 3D-QSAR, molecular docking, and molecular dynamics techniques. The results indicate that the built 3D-QSAR models yield reliable statistical information: the correlation coefficient (r²) and cross-validation coefficient (q²) of CoMFA (comparative molecular field analysis) are 0.992 and 0.819; the r² and q² of CoMSIA (comparative molecular similarity analysis) are 0.992 and 0.863, respectively. Molecular docking and MD simulations were conducted to confirm the detailed binding mode of enzyme-inhibitor system. The new NA inhibitors had been designed, synthesized, and their inhibitory activities against group-1 neuraminidase were determined. One agent displayed excellent neuraminidase inhibition, with IC₅₀ value of 39.6?μM against NA, while IC₅₀ value for oseltamivir is 61.1?μM. This compound may be further investigated for the treatment of infection by the new type influenza virus.     

Design of the influenza virus inhibitors targeting the PA endonuclease using 3D-QSAR modeling,side-chain hopping,and docking

With the emergence of drug resistance and the structural determination of the PA N-terminal domain (PA_N), influenza endonucleases have become an attractive target for antiviral therapies for influenza infection. Here, we combined 3D-QSAR with side-chain hopping and molecular docking to produce novel structures as endonuclease inhibitors. First, a new molecular library was generated with side-chain hopping on an existing template molecule, L-742001, using an in-house fragment library that targets bivalent-cation-binding proteins. Then, the best 3D-QSAR model (AAAHR.500), with q² = 0.76 and r² = 0.97 from phase modeling, was constructed from 23 endonuclease inhibitors and validated with 17 test compounds. The AAAHR.500 model was then used to select effective candidates from the new molecular library. Combining 3D-QSAR with docking using Glide and Autodock, 13 compounds were considered the most likely candidate inhibitors. Docking studies showed that the binding modes of these compounds were consistent with the crystal structures of known inhibitors. These compounds could serve as potential endonuclease inhibitors for further biological activity tests.     

Pharmacophore modelling and atom-based 3D-QSAR studies on N-methyl pyrimidones as HIV-1 integrase inhibitors

Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to N-methyl pyrimidones as HIV-1 integrase inhibitors. Based on the ligand-based pharmacophore model, we got 5-point pharmacophore model AADDR, with two hydrogen bond acceptors (A), two hydrogen bond donors (D) and one aromatic ring (R). The generated pharmacophore-based alignment was used to derive a predictive atom-based 3D-QSAR model for the training set (r²?=?0.92, SD?=?0.16, F?=?84.8, N?=?40) and for test set (Q²?=?0.71, RMSE?=?0.06, Pearson R?=?0.90, N?=?10). From these results, AADDR pharmacophore feature was selected as best common pharmacophore hypothesis, and atom-based 3D-QSAR results also support the outcome by means of favourable and unfavourable regions of hydrophobic and electron-withdrawing groups for the most potent compound 30. These results can be useful for further design of new and potent HIV-1 IN inhibitors.     

Bacterial Neuraminidase Rescues Influenza Virus Replication from Inhibition by a Neuraminidase Inhibitor

Influenza virus neuraminidase (NA) cleaves terminal sialic acid residues on oligosaccharide chains that are receptors for virus binding, thus playing an important role in the release of virions from infected cells to promote the spread of cell-to-cell infection. In addition, NA plays a role at the initial stage of viral infection in the respiratory tract by degrading hemagglutination inhibitors in body fluid which competitively inhibit receptor binding of the virus. Current first line anti-influenza drugs are viral NA-specific inhibitors, which do not inhibit bacterial neuraminidases. Since neuraminidase producing bacteria have been isolated from oral and upper respiratory commensal bacterial flora, we posited that bacterial neuraminidases could decrease the antiviral effectiveness of NA inhibitor drugs in respiratory organs when viral NA is inhibited. Using in vitro models of infection, we aimed to clarify the effects of bacterial neuraminidases on influenza virus infection in the presence of the NA inhibitor drug zanamivir. We found that zanamivir reduced progeny virus yield to less than 2% of that in its absence, however the yield was restored almost entirely by the exogenous addition of bacterial neuraminidase from Streptococcus pneumoniae. Furthermore, cell-to-cell infection was severely inhibited by zanamivir but restored by the addition of bacterial neuraminidase. Next we examined the effects of bacterial neuraminidase on hemagglutination inhibition and infectivity neutralization activities of human saliva in the presence of zanamivir. We found that the drug enhanced both inhibitory activities of saliva, while the addition of bacterial neuraminidase diminished this enhancement. Altogether, our results showed that bacterial neuraminidases functioned as the predominant NA when viral NA was inhibited to promote the spread of infection and to inactivate the neutralization activity of saliva. We propose that neuraminidase from bacterial flora in patients may reduce the efficacy of NA inhibitor drugs during influenza virus infection. (295 words).     

Neuraminidase pharmacophore model derived from diverse classes of inhibitors

A three-dimensional pharmacophore model was developed based on 22 currently available inhibitors, which were carefully selected with great diversity in both molecular structure and bioactivity, for discovering new potent neuraminidase (NA) inhibitors to fight against avian influenza virus. The best hypothesis (Hypo1), consisting of five features, namely, one positive ionizable group, one negative ionizable group, one hydrophobic point, and two hydrogen-bond donors, has a correlation coefficient of 0.902, a root mean square deviation of 1.392, and a cost difference of 72.88, suggesting that a highly predictive pharmacophore model was successfully obtained. The application of the model shows great success in predicting the activities of 88 known NA inhibitors in our test set with a correlation coefficient of 0.818 with a cross-validation of 98% confidence level. Accordingly, our model should be reliable in identifying structurally diverse compounds with desired biological activity.     

Design and development of topoisomerase inhibitors using molecular modelling studies

Topoisomerase inhibitors are used as anticancer and antibacterial agents. A series of novel 2,4,6-tri-substituted pyridine derivatives reported as topoisomerase inhibitors were used for quantitative structure–activity relationship (QSAR) study. In order to understand the structural requirement of these topoisomerase inhibitors, a ligand-based pharmacophore and atom-based 3D-QSAR model have been developed. A five-point pharmacophore with one hydrophobic group (H4), four aromatic rings (R5, R6, R7 and R8) was obtained. The pharmacophore hypothesis yielded a 3D-QSAR model with good partial least-square (PLS) statistic results. The training set correlation is characterized by PLS factors (r² = 0.7892, SD = 0.2948, F = 49.9, P = 1.379). The test set correlation is characterized by PLS factors (q² = 0.7776, root mean squared error = 0.2764, Pearson R = 0.8926). The docking study revealed the binding orientations of these inhibitors at active site amino acid residues of topoisomerases enzyme. The results of pharmacophore hypothesis and 3D-QSAR provided the detail structural insights as well as highlighted the important binding features of novel 2,4,6-tri-substituted pyridine derivatives and can be developed as potent topoisomerase inhibitors.

Figure

Open in a separate windowKey structural requirement for topoisomerase activity     

Predicting the potency of hERG K<Superscript>+</Superscript> channel inhibition by combining 3D-QSAR pharmacophore and 2D-QSAR models

Blockade of the hERG K⁺ channel has been identified as the most important mechanism of QT interval prolongation and thus inducing cardiac risk. In this work, an ensemble of 3D-QSAR pharmacophore models was constructed to provide insight into the determinants of the interactions between the hERG K⁺ channel and channel inhibitors. To predict hERG inhibitory activities, the predicted values from the ensemble of models were averaged, and the results thus obtained showed that the predictive ability of the combined 3D-QSAR pharmacophore model was greater that those of the individual models. Also, using the same training and test sets, a 2D-QSAR model based on a heuristic machine-learning method was developed in order to analyze the physicochemical characters of hERG inhibitors. The models indicated that the inhibitors have certain key inhibitory features in common, including hydrophobicity, aromaticity, and flexibility. A final model was developed by combining the combined 3D-QSAR pharmacophore with the 2D-QSAR model, and this final model outperformed any other individual model, showing the highest predictive ability and the lowest deviation. This model can not only predict hERG inhibitory potency accurately, thus allowing fast cardiac safety evaluation, but it provides an effective tool for avoiding hERG inhibitory liability and thus enhanced cardiac risk in the design and optimization of new chemical entities.     

Structural insights of JAK2 inhibitors: pharmacophore modeling and ligand-based 3D-QSAR studies of pyrido-indole derivatives

Abstract

In this study we have performed pharmacophore modeling and built a 3D QSAR model for pyrido-indole derivatives as Janus Kinase 2 inhibitors. An efficient pharmacophore has been identified from a data set of 51 molecules and the identified pharmacophore hypothesis consisted of one hydrogen bond acceptor, two hydrogen bond donors and three aromatic rings, i.e. ADDRRR. A powerful 3D-QSAR model has also been constructed by employing Partial Least Square regression analysis with a regression coefficient of 0.97 (R²) and Q² of 0.95, and Pearson-R of 0.98.     

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

RWJ-270201 (BCX-1812): a novel neuraminidase inhibitor for influenza.

The influenza virus neuraminidase (NA) is important in the pathogenesis of infection and, thus, is an attractive target for agents used in the treatment and prophylaxis of influenza. This article describes preclinical and early clinical data related to RWJ-270201 (BCX-1812), a novel, orally active NA inhibitor that was rationally designed for having potent and selective activity against influenza A and B viruses. RWJ-270201 is a unique NA inhibitor with a cyclopentane ring structure and high selectivity for the influenza NA. RWJ-270201 has efficacy comparable to or better than earlier NA inhibitors against a wide range of influenza A and B isolates, including recently emerged and avian strains, both in vitro and in a lethal murine model of influenza. Based on the high selectivity and efficacy of RWJ-270201 against both type A and B influenza strains in preclinical studies as well as murine pharmacodynamic studies supporting the potential for once-daily administration, clinical trials were initiated in order to determine the tolerability and antiviral activity of RWJ-270201 in humans. To date, clinical studies have indicated that RWJ-270201 is well tolerated and has antiviral activity in human experimental influenza models when administered orally once daily.     

Pharmacophore modeling, 3D-QSAR,and in silico ADME prediction of N-pyridyl and pyrimidine benzamides as potent antiepileptic agents

Biological mechanism attributing mutations in KCNQ2/Q3 results in benign familial neonatal epilepsy (BFNE), a rare form of epilepsy and thus neglected. It offers a potential target for antiepileptic drug discovery. In the present work, a pharmacophore-based 3D-QSAR model was generated for a series of N-pyridyl and pyrimidine benzamides possessing KCNQ2/Q3 opening activity. The pharmacophore model generated contains one hydrogen bond donor (D), one hydrophobic (H), and two aromatic rings (R). They are the crucial molecular write-up detailing predicted binding efficacy of high affinity and low affinity ligands for KCNQ2/Q3 opening activity. Furthermore, it has been validated by using a biological correlation between pharmacophore hypothesis-based 3D-QSAR variables and functional fingerprints of openers responsible for the receptor binding and also by docking of these benzamides into the validated homology model. Excellent statistical computational tools of QSAR model such as good correlation coefficient (R^2?>^?0.80), higher F value (F?>?39), and excellent predictive power (Q²?>?0.7) with low standard deviation (SD <0.3) strongly suggest that the developed model could be used for prediction of antiepileptic activity of newer analogs. A preliminary pharmacokinetic profile of these derivatives was also performed on the basis of QikProp predictions.     

Design,synthesis, and biological activity of thiazole derivatives as novel influenza neuraminidase inhibitors

A series of novel influenza neuraminidase (NA) inhibitors based on thiazole core were synthesized and evaluated for their ability to inhibit NA of influenza A virus (H₃N₂). All compounds were synthesized in good yields starting from commercially available 2-amino-4-thiazole-acetic ester using a suitable synthetic strategy. These compounds showed moderate inhibitory activity against influenza A NA. The most potent compound of this series is compound 4d (IC₅₀?=?3.43 μM), which is about 20-fold less potent than oseltamivir, and could be used to design novel influenza NA inhibitors that exhibit increased activity based on thiazole ring.     

Inhibition of Influenza A Virus Infection In Vitro by Peptides Designed In Silico

Influenza A viruses are enveloped, segmented negative single-stranded RNA viruses, capable of causing severe human respiratory infections. Currently, only two types of drugs are used to treat influenza A infections, the M2 H⁺ ion channel blockers (amantadine and rimantadine) and the neuraminidase inhibitors (NAI) (oseltamivir and zanamivir). Moreover, the emergence of drug-resistant influenza A virus strains has emphasized the need to develop new antiviral agents to complement or replace the existing drugs. Influenza A virus has on the surface a glycoprotein named hemagglutinin (HA) which due to its important role in the initial stage of infection: receptor binding and fusion activities of viral and endosomal membranes, is a potential target for new antiviral drugs. In this work we designed nine peptides using several bioinformatics tools. These peptides were derived from the HA1 and HA2 subunits of influenza A HA with the aim to inhibit influenza A virus infection. The peptides were synthetized and their antiviral activity was tested in vitro against several influenza A viral strains: Puerto Rico/916/34 (H1N1), (H1N1)pdm09, swine (H1N1) and avian (H5N2). We found these peptides were able to inhibit the influenza A viral strains tested, without showing any cytotoxic effect. By docking studies we found evidence that all the peptides were capable to bind to the viral HA, principally to important regions on the viral HA stalk, thus could prevent the HA conformational changes required to carry out its membranes fusion activity.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Efficacy of the New Neuraminidase Inhibitor CS-8958 against H5N1 Influenza Viruses

Currently, two neuraminidase (NA) inhibitors, oseltamivir and zanamivir, which must be administrated twice daily for 5 days for maximum therapeutic effect, are licensed for the treatment of influenza. However, oseltamivir-resistant mutants of seasonal H1N1 and highly pathogenic H5N1 avian influenza A viruses have emerged. Therefore, alternative antiviral agents are needed. Recently, a new neuraminidase inhibitor, R-125489, and its prodrug, CS-8958, have been developed. CS-8958 functions as a long-acting NA inhibitor in vivo (mice) and is efficacious against seasonal influenza strains following a single intranasal dose. Here, we tested the efficacy of this compound against H5N1 influenza viruses, which have spread across several continents and caused epidemics with high morbidity and mortality. We demonstrated that R-125489 interferes with the NA activity of H5N1 viruses, including oseltamivir-resistant and different clade strains. A single dose of CS-8958 (1,500 µg/kg) given to mice 2 h post-infection with H5N1 influenza viruses produced a higher survival rate than did continuous five-day administration of oseltamivir (50 mg/kg twice daily). Virus titers in lungs and brain were substantially lower in infected mice treated with a single dose of CS-8958 than in those treated with the five-day course of oseltamivir. CS-8958 was also highly efficacious against highly pathogenic H5N1 influenza virus and oseltamivir-resistant variants. A single dose of CS-8958 given seven days prior to virus infection also protected mice against H5N1 virus lethal infection. To evaluate the improved efficacy of CS-8958 over oseltamivir, the binding stability of R-125489 to various subtypes of influenza virus was assessed and compared with that of other NA inhibitors. We found that R-125489 bound to NA more tightly than did any other NA inhibitor tested. Our results indicate that CS-8958 is highly effective for the treatment and prophylaxis of infection with H5N1 influenza viruses, including oseltamivir-resistant mutants.     

A monoclonal antibody targeting a highly conserved epitope in influenza B neuraminidase provides protection against drug resistant strains

All influenza viral neuraminidases (NA) of both type A and B viruses have only one universally conserved sequence located between amino acids 222–230. A monoclonal antibody against this region has been previously reported to provide broad inhibition against all nine subtypes of influenza A NA; yet its inhibitory effect against influenza B viral NA remained unknown. Here, we report that the monoclonal antibody provides a broad inhibition against various strains of influenza B viruses of both Victoria and Yamagata genetic lineage. Moreover, the growth and NA enzymatic activity of two drug resistant influenza B strains (E117D and D197E) are also inhibited by the antibody even though these two mutations are conformationally proximal to the universal epitope. Collectively, these data suggest that this unique, highly-conserved linear sequence in viral NA is exposed sufficiently to allow access by inhibitory antibody during the course of infection; it could represent a potential target for antiviral agents and vaccine-induced immune responses against diverse strains of type B influenza virus.     

Influenza A virus neuraminidase limits viral superinfection

Enveloped viruses use multiple mechanisms to inhibit infection of a target cell by more than one virion. These mechanisms may be of particular importance for the evolution of segmented viruses, because superinfection exclusion may limit the frequency of reassortment of viral genes. Here, we show that cellular expression of influenza A virus neuraminidase (NA), but not hemagglutinin (HA) or the M2 proton pump, inhibits entry of HA-pseudotyped retroviruses. Cells infected with H1N1 or H3N2 influenza A virus were similarly refractory to HA-mediated infection and to superinfection with a second influenza A virus. Both HA-mediated entry and viral superinfection were rescued by the neuraminidase inhibitors oseltamivir carboxylate and zanamivir. These inhibitors also prevented the removal of alpha-2,3- and alpha-2,6-linked sialic acid observed in cells expressing NA or infected with influenza A viruses. Our data indicate that NA alone among viral proteins limits influenza A virus superinfection.     

A strategy for theoretical binding constant,Ki, calculations for neuraminidase aromatic inhibitors designed on the basis of the active site structure of influenza virus neuraminidase

Neuraminidase (NA) is one of the two major surface antigens of influenza virus. It plays an indispensable role in the release and spread of progeny virus particles during infection. NA inhibitors reduce virus infection in animals. To improve the clinical efficacy of NA inhibitors, we have begun the design of non-carbohydrate inhibitors based on the active site structure of NA. The approach is an iterative process of ligand modeling and electrostatic calculations followed by chemical synthesis of compounds, biological testing, and NA-inhibitor complex structure determination by X-ray crystallography. A strategy has been developed to calculate K_i for newly designed inhibitors. The calculations using the DelPhi program were performed for carbohydrate inhibitors and three preliminary benzoic acid inhibitors of neuraminidase (BANA) that have been synthesized and shown to bind to the active site of NA in the crystal structure. The calculated K_is of these inhibitors have an enlightening agreement with their in vitro biological activities. This demonstrates that the calculations produce informative results on the affinity of modeled inhibitors. GRID maps were also calculated and several pockets were identified for accepting possible new ligands. The calculated K_is for newly designed ligands suggest that these potential compounds will have high inhibitory activities. © 1995 Wiley-Liss, Inc.     

Identification of novel urease inhibitors: pharmacophore modeling,virtual screening and molecular docking studies

Abstract

Pharmacophore modeling and atom-based three-dimensional quantitative structure–activity relationship (3D-QSAR) have been developed on N-acylglycino- and hippurohydroxamic acid derivatives, which are known potential inhibitors of urease. This is followed by virtual screening and ADMET (absorption, distribution, metabolism, excretion and toxicity) studies on a large library of known drugs in order to get lead molecules as Helicobacter pylori urease inhibitors. A suitable three-featured pharmacophore model comprising one H-bond acceptor and two H-bond donor features (ADD.10) has been found to be the best QSAR model. An external library of compounds (～3000 molecules), pre-filtered using Lipinski’s rule of five, has been further screened using the pharmacophore model ADD.10. By analyzing the fitness of the hits with respect to the pharmacophore model and their binding interaction inside the urease active site, four molecules have been predicted to be extremely good urease inhibitors. Two of these have significant potential and should be taken up for further drug-designing process.