Selective binding of antiinfluenza drugs and their analogues to 'open' and 'closed' conformations of H5N1 neuraminidase

Designing new drug candidates as inhibitors against wild and mutant type neuraminidases: molecular docking, molecular dynamics and binding free energy calculations

Influenza virus is the cause of the death of millions of people with about 3-4 pandemics every hundred years in history. It also turns into a seasonal disease, bringing about approximately 5-15% of the population to be infected and 290,000-650,000 people to die every year. These numbers reveal that it is necessary to be on the alert to work towards influenza in order to protect public health. There are FDA-approved antiviral drugs such as oseltamivir and zanamivir recommended by the World Center for Disease Prevention. However, after the recent outbreaks such as bird flu and swine flu, increasing studies have shown that the flu virus has gained resistance to these drugs. So, there is an urgent need to find new drugs effective against this virus. This study aims to investigate new drug candidates targeting neuraminidase (NA) for the treatment of influenza by using computer aided drug design approaches. They involve virtual scanning, de novo design, rational design, docking, MD, MMGB/PBSA. The investigation includes H1N1, H5N1, H2N2 and H3N2 neuraminidase proteins and their mutant variants possessing resistance to FDA-approved drugs. Virtual screening consists of approximately 30 thousand molecules while de novo and rational designs produced over a hundred molecules. These approaches produced three lead molecules with binding energies for both non-mutant (-34.84, -59.99 and -60.66 kcal/mol) and mutant (-40.40, -58.93, -76.19 kcal/mol) H2N2 NA calculated by MM-PBSA compared with those of oseltamivir -25.64 and -18.40 respectively. The results offer new drug candidates against influenza infection.Communicated by Ramaswamy H. Sarma.

Five Novel Non-Sialic Acid-Like Scaffolds Inhibit In Vitro H1N1 and H5N2 Neuraminidase Activity of Influenza a Virus

Neuraminidase (NA) of influenza viruses enables the virus to access the cell membrane. It degrades the sialic acid contained in extracellular mucin. Later, it is responsible for releasing newly formed virions from the membrane of infected cells. Both processes become key functions within the viral cycle. Therefore, it is a therapeutic target for research of the new antiviral agents. Structure–activity relationships studies have revealed which are the important functional groups for the receptor–ligand interaction. Influenza virus type A NA activity was inhibited by five scaffolds without structural resemblance to sialic acid. Intending small organic compound repositioning along with drug repurposing, this study combined in silico simulations of ligand docking into the known binding site of NA, along with in vitro bioassays. The five proposed scaffolds are N-acetylphenylalanylmethionine, propanoic 3-[(2,5-dimethylphenyl) carbamoyl]-2-(piperazin-1-yl) acid, 3-(propylaminosulfonyl)-4-chlorobenzoic acid, ascorbic acid (vitamin C), and 4-(dipropylsulfamoyl) benzoic acid (probenecid). Their half maximal inhibitory concentration (IC₅₀) was determined through fluorometry. An acidic reagent 2′-O-(4-methylumbelliferyl)-α-dN-acetylneuraminic acid (MUNANA) was used as substrate for viruses of human influenza H1N1 or avian influenza H5N2. Inhibition was observed in millimolar ranges in a concentration-dependent manner. The IC₅₀ values of the five proposed scaffolds ranged from 6.4 to 73 mM. The values reflect a significant affinity difference with respect to the reference drug zanamivir (p < 0.001). Two compounds (N-acetyl dipeptide and 4-substituted benzoic acid) clearly showed competitive mechanisms, whereas ascorbic acid reflected non-competitive kinetics. The five small organic molecules constitute five different scaffolds with moderate NA affinities. They are proposed as lead compounds for developing new NA inhibitors which are not analogous to sialic acid.

Impact of tetramerization on the ligand recognition of N1 influenza neuraminidase via MMGBSA approach

Influenza virus neuraminidase (NA) is a homotetrameric surface protein that, in contrast to other non-influenza NAs, requires a quaternary assembly to exhibit enzymatic activity, suggesting that the oligomeric state significantly impacts the active site of influenza NA. Nevertheless, most structure-based drug design studies have been reported by employing the monomeric state in the closed or open-loop due to the computational cost of employing the tetrameric NA. In this work, we present MD simulations coupled to the MMGBSA approach of avian N1 type NA in its monomeric and tetrameric closed and open-loop state both with and without the inhibitor oseltamivir and its natural substrate, sialic acid. Structural and energetic analyses revealed that the tetrameric state impacts flexibility as well as the map of interactions participating in stabilizing the protein-ligand complexes with respect to the monomeric state. It was observed that the tetrameric state exerts dissimilar effects in binding affinity, characteristic of positive and negative cooperativity for oseltamivir and sialic acid, respectively. Based on our results, to perform a confident structure-based drug design, as well as to evaluate the impact of key mutations through MD simulations, it is important to consider the tetrameric state closed-loop state.

The significance of naturally occurring neuraminidase quasispecies of H5N1 avian influenza virus on resistance to oseltamivir: a point of concern

Viral adaptability and survival arise due to the presence of quasispecies populations that are able to escape the immune response or produce drug-resistant variants. However, the presence of H5N1 virus with natural mutations acquired without any drug selection pressure poses a great threat. Cloacal samples collected from the 2004-2005 epidemics in Thailand from Asian open-billed storks revealed one major and several minor quasispecies populations with mutations on the oseltamivir (OTV)-binding site of the neuraminidase gene (NA) without prior exposure to a drug. Therefore, this study investigated the binding between the NA-containing novel mutations and OTV drug using molecular dynamic simulations and plaque inhibition assay. The results revealed that the mutant populations, S236F mutant, S236F/C278Y mutant, A250V/V266A/P271H/G285S mutant and C278Y mutant, had a lower binding affinity with OTV as compared with the WT virus due to rearrangement of amino acid residues and increased flexibility in the 150-loop. This result was further emphasized through the IC50 values obtained for the major population and WT virus, 104.74 nM and 18.30 nM, respectively. Taken together, these data suggest that H5N1 viruses isolated from wild birds have already acquired OTV-resistant point mutations without any exposure to a drug.

Discovery of Influenza A virus neuraminidase inhibitors using support vector machine and Naïve Bayesian models

Neuraminidase (NA) is a critical enzyme in the life cycle of influenza virus, which is known as a successful paradigm in the design of anti-influenza agents. However, to date there are no classification models for the virtual screening of NA inhibitors. In this work, we built support vector machine and Naïve Bayesian models of NA inhibitors and non-inhibitors, with different ratios of active-to-inactive compounds in the training set and different molecular descriptors. Four models with sensitivity or Matthews correlation coefficients greater than 0.9 were chosen to predict the NA inhibitory activities of 15,600 compounds in our in-house database. We combined the results of four optimal models and selected 60 representative compounds to assess their NA inhibitory profiles in vitro. Nine NA inhibitors were identified, five of which were oseltamivir derivatives with large C-5 substituents exhibiting potent inhibition against H1N1 NA with IC50 values in the range of 12.9-185.0 nM, and against H3N2 NA with IC50 values between 18.9 and 366.1 nM. The other four active compounds belonged to novel scaffolds, with IC50 values ranging 39.5-63.8 μM against H1N1 NA and 44.5-114.1 μM against H3N2 NA. This is the first time that classification models of NA inhibitors and non-inhibitors are built and their prediction results validated experimentally using in vitro assays.

Molecular modeling studies demonstrate key mutations that could affect the ligand recognition by influenza AH1N1 neuraminidase

The goal of this study was to identify neuraminidase (NA) residue mutants from human influenza AH1N1 using sequences from 1918 to 2012. Multiple alignment studies of complete NA sequences (5732) were performed. Subsequently, the crystallographic structure of the 1918 influenza (PDB ID: 3BEQ-A) was used as a wild-type structure and three-dimensional (3-D) template for homology modeling of the mutated selected NA sequences. The 3-D mutated NAs were refined using molecular dynamics (MD) simulations (50 ns). The refined 3-D models were used to perform docking studies using oseltamivir. Multiple sequence alignment studies showed seven representative mutations (A232V, K262R, V263I, T264V, S367L, S369N, and S369K). MD simulations applied to 3-D NAs showed that each NA had different active-site shapes according to structural surface visualization and docking results. Moreover, Cartesian principal component analyses (cPCA) show structural differences among these NA structures caused by mutations. These theoretical results suggest that the selected mutations that are located outside of the active site of NA could affect oseltamivir recognition and could be associated with resistance to oseltamivir.

Calculating binding free energies for protein-carbohydrate complexes

A variety of computational techniques may be applied to compute theoretical binding free energies for protein-carbohydrate complexes. Elucidation of the intermolecular interactions, as well as the thermodynamic effects, that contribute to the relative strength of receptor binding can shed light on biomolecular recognition, and the resulting initiation or inhibition of a biological process. Three types of free energy methods are discussed here, including MM-PB/GBSA, thermodynamic integration, and a non-equilibrium alternative utilizing SMD. Throughout this chapter, the well-known concanavalin A lectin is employed as a model system to demonstrate the application of these methods to the special case of carbohydrate binding.

Synthesis of acylguanidine zanamivir derivatives as neuraminidase inhibitors and the evaluation of their bio-activities

A series of acylguanidine-modified zanamivir analogs were synthesized and their inhibitory activities against the NAs of avian influenza viruses (H1N1 and H3N2) were evaluated. In particular, zanamivir derivative , with a hydrophobic naphthalene substituent, exhibits the best inhibitory activity against group-1 NA with an IC50 of 20 nM.

Mutation-induced loop opening and energetics for binding of tamiflu to influenza N8 neuraminidase

Tamiflu, also known as oseltamivir (OTV), binds to influenza A neuraminidase (H5N1) with very high affinity (0.32 nM). However, this inhibitor binds to other neuraminidases as well. In the present work, a systematic computational study is performed to investigate the mechanism underlying the binding of oseltamivir to N8 neuraminidase (NA) in "open" and "closed" conformations of the 150-loop through molecular dynamics simulations and the popular and well established molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation method. Whereas the closed conformation is stable for wild type N8, it transforms into the open conformation for the mutants Y252H, H274Y, and R292K, indicating that bound to oseltamivir these mutants are preferentially in the open conformation. Our calculations show that the binding of wild type oseltamivir to the closed conformation of N8 neuraminidase is energetically favored compared to the binding to the open conformation. We observe water mediated binding of oseltamivir to the N8 neuraminidase in both conformations which is not seen in the case of binding of the same drug to the H5N1 neuraminidase. The decomposition of the binding free energy reveals the mechanisms underlying the binding and changes in affinity due to mutations. Considering the mutant N8 variants in the open conformation adopted during the simulations, we observe a significant loss in the size of the total binding free energy for the N8(Y252H)-OTV, N8(H274Y)-OTV, and N8(R292K)-OTV complexes compared to N8(WT)-OTV, mainly due to the decrease in the size of the intermolecular electrostatic energy. For R292K, an unfavorable shift in the van der Waals interactions also contributes to the drug resistance. The mutations cause a significant expansion in the active site cavity, increasing its solvent accessible surface compared to the crystal structures of both the open and closed conformations. Our study underscores the need to consider dynamics in rationalizing the structure-function relationships of various antiviral inhibitor-NA complexes.

Long time scale GPU dynamics reveal the mechanism of drug resistance of the dual mutant I223R/H275Y neuraminidase from H1N1-2009 influenza virus

Multidrug resistance of the pandemic H1N1-2009 strain of influenza has been reported due to widespread treatment using the neuraminidase (NA) inhibitors, oseltamivir (Tamiflu), and zanamivir (Relenza). From clinical data, the single I223R (IR(1)) mutant of H1N1-2009 NA reduced efficacy of oseltamivir and zanamivir by 45 and 10 times, (1) respectively. More seriously, the efficacy of these two inhibitors against the double mutant I223R/H275Y (IRHY(2)) was significantly reduced by a factor of 12 374 and 21 times, respectively, compared to the wild-type.(2) This has led to the question of why the efficacy of the NA inhibitors is reduced by the occurrence of these mutations and, specifically, why the efficacy of oseltamivir against the double mutant IRHY was significantly reduced, to the point where oseltamivir has become an ineffective treatment. In this study, 1 μs of molecular dynamics (MD) simulations was performed to answer these questions. The simulations, run using graphical processors (GPUs), were used to investigate the effect of conformational change upon binding of the NA inhibitors oseltamivir and zanamivir in the wild-type and the IR and IRHY mutant strains. These long time scale dynamics simulations demonstrated that the mechanism of resistance of IRHY to oseltamivir was due to the loss of key hydrogen bonds between the inhibitor and residues in the 150-loop. This allowed NA to transition from a closed to an open conformation. Oseltamivir binds weakly with the open conformation of NA due to poor electrostatic interactions between the inhibitor and the active site. The results suggest that the efficacy of oseltamivir is reduced significantly because of conformational changes that lead to the open form of the 150-loop. This suggests that drug resistance could be overcome by increasing hydrogen bond interactions between NA inhibitors and residues in the 150-loop, with the aim of maintaining the closed conformation, or by designing inhibitors that can form a hydrogen bond to the mutant R223 residue, thereby preventing competition between R223 and R152.

Molecular-level simulation of pandemic influenza glycoproteins

Computational simulation of pandemic diseases provides important insight into many disease features that may benefit public health. This is especially true for the influenza virus, a continuing global pandemic threat. Molecular or atomic-level investigation of influenza has predominantly focused on the two major virus glycoproteins, neuraminidase (NA) and hemagglutinin (HA). In this chapter, we walk the readers through major considerations for studying pandemic influenza glycoproteins, from choosing the most useful choice of system(s) to avoiding common pitfalls in experimental design and execution. While a brief discussion of several potential simulation and docking techniques is presented, we emphasize molecular dynamics (MD) and Brownian dynamics (BD) simulation techniques and molecular docking, within the context of biologically outstanding questions in influenza research.

Structural and functional analysis of laninamivir and its octanoate prodrug reveals group specific mechanisms for influenza NA inhibition

The 2009 H1N1 influenza pandemic (pH1N1) led to record sales of neuraminidase (NA) inhibitors, which has contributed significantly to the recent increase in oseltamivir-resistant viruses. Therefore, development and careful evaluation of novel NA inhibitors is of great interest. Recently, a highly potent NA inhibitor, laninamivir, has been approved for use in Japan. Laninamivir is effective using a single inhaled dose via its octanoate prodrug (CS-8958) and has been demonstrated to be effective against oseltamivir-resistant NA in vitro. However, effectiveness of laninamivir octanoate prodrug against oseltamivir-resistant influenza infection in adults has not been demonstrated. NA is classified into 2 groups based upon phylogenetic analysis and it is becoming clear that each group has some distinct structural features. Recently, we found that pH1N1 N1 NA (p09N1) is an atypical group 1 NA with some group 2-like features in its active site (lack of a 150-cavity). Furthermore, it has been reported that certain oseltamivir-resistant substitutions in the NA active site are group 1 specific. In order to comprehensively evaluate the effectiveness of laninamivir, we utilized recombinant N5 (typical group 1), p09N1 (atypical group 1) and N2 from the 1957 pandemic H2N2 (p57N2) (typical group 2) to carry out in vitro inhibition assays. We found that laninamivir and its octanoate prodrug display group specific preferences to different influenza NAs and provide the structural basis of their specific action based upon their novel complex crystal structures. Our results indicate that laninamivir and zanamivir are more effective against group 1 NA with a 150-cavity than group 2 NA with no 150-cavity. Furthermore, we have found that the laninamivir octanoate prodrug has a unique binding mode in p09N1 that is different from that of group 2 p57N2, but with some similarities to NA-oseltamivir binding, which provides additional insight into group specific differences of oseltamivir binding and resistance.

The influenza neuraminidase (NA) enzyme is the most successful drug target against the seasonal and pandemic flu. The 2009 H1N1 flu pandemic led to record sales of the NA inhibitors oseltamivir (Tamiflu) and zanamivir (Relenza). Recently, a new drug, laninamivir (Inavir), has been approved for use in Japan can also be administered effectively using a single dose via its octanoate prodrug (CS-8958), however its effectiveness against oseltamivir-resistant influenza infection has not been demonstrated in clinical studies. In this study we comprehensively evaluate the effectiveness of laninamivir and its prodrug using NA from different groups with different active site features. We expressed and purified a group 2 NA from the 1957 pandemic H2N2, an atypical group 1 NA from the 2009 H1N1 pandemic and a group 1 NA from avian H12N5. NA inhibition was assayed and NAs were further crystallized with each inhibitor to determine the structural basis of their action. We found that laninamivir inhibition is highly potent for each NA, however binding and inhibition of laninamivir and its prodrug showed group specific preferences. Our results provide the structural and functional basis of NA inhibition using classical and novel inhibitors, with NAs from multiple serotypes with different properties.

Study of Tamiflu sensitivity to variants of A/H5N1 virus using different force fields

An accurate estimation of binding free energy of a ligand to receptor ΔG(bind) is one of the most important problems in drug design. The success of solution of this problem is expected to depend on force fields used for modeling a ligand-receptor complex. In this paper, we consider the impact of four main force fields, AMBER99SB, CHARMM27, GROMOS96 43a1, and OPLS-AA/L, on the binding affinity of Oseltamivir carboxylate to the wild-type and Y252H, N294S, and H274Y mutants of glycoprotein neuraminidase from the pandemic A/H5N1 virus. Having used the molecular mechanic-Poisson-Boltzmann surface area method, we have shown that ΔG(bind), obtained by AMBER99SB, OPLS-AA/L, and CHARMM27, shows the high correlation with the available experimental data. They correctly capture the binding ranking Y252H → WT → N294S → H274Y observed in experiments (Collins, P. J. et al. Nature 2008, 453, 1258). In terms of absolute values of binding scores, results obtained by AMBER99SB are in the nearest range with experiments, while OPLS-AA/L, which is applied to study binding of Oseltamivir to the influenza virus for the first time, gives rather big negative values for ΔG(bind). GROMOS96 43a1 provides a lower correlation as it supports Oseltamivir to be more resistant to N294S than H274Y. Our study suggests that force fields have pronounced influence on theoretical estimations of binding free energy of a ligand to receptor. The effect of all-atom models on dynamics of the binding pocket as well as on the hydrogen-bond network between Oseltamivir and receptors is studied in detail. The hydrogen network, obtained by GROMOS, is weakest among four studied force fields.