The binding properties of the H5N1 influenza virus neuraminidase as inferred from molecular modeling

The avian influenza H5N1 virus has emerged as an important pathogen, causing severe disease in humans and posing a pandemic threat. Substrate specificity is crucial for the virus to obtain the ability to spread from avian to human. Therefore, an investigation of the binding properties of ligands at the molecular level is important for understanding the catalytic mechanism of the avian influenza virus neuraminidase and for designing novel and specific inhibitors of H5N1 neuraminidase. Based on the available crystal structure of H5N1, we have characterized the binding properties between sialic acid, methyl 3’sialyllactoside, methyl 6’sialyllactoside and the H5N1 influenza virus neuraminidase using molecular docking and molecular dynamics simulations. Obtained molecular dynamics trajectories were analyzed in terms of ligand conformations, N1-ligand interactions, and in terms of loop flexibility. It was found that in the N1-SA complex the sialic acid ring undergoes a transition from the B _2,5 to the ² C ₅ conformation. However, in the N1-3SL and N1-6SL complexes sialic acid remained in the distorted boat conformation. The obtained results indicate that 3SL has only weak interactions with the 150-loop, whereas the N1-6SL complex shows strong interactions. Most of the differences arise from the various conformations around the glycosidic linkage, between the sialic acid and galactose, which facilitate the above interactions of 6SL with the enzyme, and as a consequence the interactions between the 150- and 430- loops. This finding suggests that the altered flexibility of loops in and around the active site is one of the reasons why the avian N1 preferentially cleaves sialic acid from α-(2-3)-Gal glycoconjugates over α-(2-6)-Gal. These molecular modeling results are consistent with available experimental results on the specificity of N1.     

Analogs of zanamivir with modified C4-substituents as the inhibitors against the group-1 neuraminidases of influenza viruses

Unlike the group-2 neuraminidase, the group-1 neuraminidase of influenza virus possesses a flexible loop (the 150-loop) and a cavity (the 150-cavity) adjacent to the active site, and renders a conformational change from the ‘open’ form to the ‘closed’ form on binding with substrate (sialo-glycoprotein) or inhibitor (e.g., zanamivir). Zanamivir derivative 8a having an extended (piperazinocarbonyl)propyl substituent at the internal N-position of the guanidino group is designed as a possible inhibitor on the basis of computer docking to the open form of N1 subtype neuraminidase. Indeed, compound 8a exhibits strong neuraminidase inhibition and good anti-influenza activity against H1N1 virus with IC₅₀ = 2.15 μM and EC₅₀ = 0.77 μM, respectively. This study may provide a clue to future design of better group-1 neuraminidase inhibitors.     

Design,synthesis and biological evaluation of substituted flavones and aurones as potential anti-influenza agents

We designed a series of substituted flavones and aurones as non-competitive H1N1 neuraminidase (NA) inhibitors and anti-influenza agents. The molecular docking studies showed that the designed flavones and aurones occupied 150-cavity and 430-cavity of H1N1-NA. We then synthesized these compounds and evaluated these for cytotoxicity, reduction in H1N1 virus yield, H1N1-NA inhibition and kinetics of inhibition. The virus yield reduction assay and H1N1-NA inhibition assay demonstrated that the compound 1f (4-methoxyflavone) had the lowest EC₅₀ of 9.36 nM and IC₅₀ of 8.74 μM respectively. Moreover, kinetic studies illustrated that compounds 1f and 2f had non-competitive inhibition mechanism.     

Virtual screening pipeline and ligand modelling for H5N1 neuraminidase

The H5N1 virus neuraminidase structure was solved in two different conformations depending on the inhibitor concentration. In the absence of oseltamivir or at a low concentration, the neuraminidase structure assumes an open form that closes at a high oseltamivir concentration due to the shift of the so-called 150-loop near the active site. Although the close conformation is similar to all the other structurally known neuraminidase types, it doesn’t appear to be the most likely physiological condition for N1.To investigate the specific ligand binding properties of the open form, we screened by docking simulation, a large dataset of ligands and compared the results with closed form. The virtual screening procedure was implemented in a docking pipeline that also performs a step-by-step, target specific, filtering approach for data reduction. The selected ligands display binding ability involving multiple sites of interaction including the active site and an adjacent cavity made available by the 150-loop shift. Two ligands are especially interesting and are proposed as substituents to design oseltamivir derivatives specifically suited for the open conformation.     

Structural characterization of the 1918 influenza virus H1N1 neuraminidase

Influenza virus neuraminidase (NA) plays a crucial role in facilitating the spread of newly synthesized virus in the host and is an important target for controlling disease progression. The NA crystal structure from the 1918 "Spanish flu" (A/Brevig Mission/1/18 H1N1) and that of its complex with zanamivir (Relenza) at 1.65-A and 1.45-A resolutions, respectively, corroborated the successful expression of correctly folded NA tetramers in a baculovirus expression system. An additional cavity adjacent to the substrate-binding site is observed in N1, compared to N2 and N9 NAs, including H5N1. This cavity arises from an open conformation of the 150 loop (Gly147 to Asp151) and appears to be conserved among group 1 NAs (N1, N4, N5, and N8). It closes upon zanamivir binding. Three calcium sites were identified, including a novel site that may be conserved in N1 and N4. Thus, these high-resolution structures, combined with our recombinant expression system, provide new opportunities to augment the limited arsenal of therapeutics against influenza.     

The inhibitory performance of flavonoid cyanidin-3-sambubiocide against H274Y mutation in H1N1 influenza virus

Oseltamivir (Tamiflu) is the most accepted antiviral drug that targets the neuraminidase (NA) protein to inhibit the viral release from the host cell. Few H1N1 influenza strains with the H274Y mutation creates drug resistance to oseltamivir. In this study, we report that flavonoid cyanidin-3-sambubiocide (C3S) compound acts as a potential inhibitor against H274Y mutation. The drug resistance mechanism and inhibitory activity of C3S and oseltamivir against wild-type (WT) and H274Y mutant-type (MT) have been studied and compared based on the results of molecular docking, molecular dynamics, and quantum chemical methods. Oseltamivir has been found less binding affinity with MT. C3S has more binding affinity with WT and MT proteins. From the dynamical study, the 150th loop of the MT protein has found more deformation than WT. A single H274Y mutation induces the conformational changes in the 150th loop which leads to produce more resistance to oseltamivir. The 150th cavity is more attractive target for C3S to stop the conformational changes in the MT, than 430th cavity of NA protein. The C3S is stabilized with MT by more number of hydrogen bonds than oseltamivir. The electrostatic interaction energy shows a stronger C3S binding with MT and this compound may be more effective against oseltamivir-resistant virus strains.     

Exploring the mechanism of zanamivir resistance in a neuraminidase mutant: a molecular dynamics study

It is critical to understand the molecular basis of the drug resistance of influenza viruses to efficiently treat this infectious disease. Recently, H1N1 strains of influenza A carrying a mutation of Q136K in neuraminidase were found. The new strain showed a strong Zanamivir neutralization effect. In this study, normal molecular dynamics simulations and metadynamics simulations were employed to explore the mechanism of Zanamivir resistance. The wild-type neuraminidase contained a 3(10) helix before the 150 loop, and there was interaction between the 150 and 430 loops. However, the helix and the interaction between the two loops were disturbed in the mutant protein due to interaction between K136 and nearby residues. Hydrogen-bond network analysis showed weakened interaction between the Zanamivir drug and E276/D151 on account of the electrostatic interaction between K136 and D151. Metadynamics simulations showed that the free energy landscape was different in the mutant than in the wild-type neuraminidase. Conformation with the global minimum of free energy for the mutant protein was different from the wild-type conformation. While the drug fit completely into the active site of the wild-type neuraminidase, it did not match the active site of the mutant variant. This study indicates that the altered hydrogen-bond network and the deformation of the 150 loop are the key factors in development of Zanamivir resistance. Furthermore, the Q136K mutation has a variable effect on conformation of different N1 variants, with conformation of the 1918 N1 variant being more profoundly affected than that of the other N1 variants studied in this paper. This observation warrants further experimental investigation.     

Investigation of flexibility of neuraminidase 150-loop using tamiflu derivatives in influenza A viruses H1N1 and H5N1

This study focuses on design, synthesis and in vitro evaluation of inhibitory potency of two series of sialylmimetic that target an exosite (“150-cavity”) adjacent to the active site of influenza neuraminidases from A/California/07/2009 (H1N1) pandemic strain and A/chicken/Nakorn-Patom/Thailand/CU-K2-2004 (H5N1). The structure-activity analysis as well as 3-D structure of the complex of parental compound with the pandemic neuraminidase p09N1 revealed high flexibility of the 150-cavity towards various modification of the neuraminidase inhibitors. Furthermore, our comparison of two methods for inhibition constant determination performed at slightly different pH values suggest that the experimental conditions of the measurement could dramatically influence the outcome of the analysis in the compound-dependent manner. Therefore, previously reported K_i values determined at non-physiological pH should be carefully scrutinized.     

Mapping the sequence mutations of the 2009 H1N1 influenza A virus neuraminidase relative to drug and antibody binding sites

 

In this work, we study the consequences of sequence variations of the "2009 H1N1" (swine or Mexican flu) influenza A virus strain neuraminidase for drug treatment and vaccination. We find that it is phylogenetically more closely related to European H1N1 swine flu and H5N1 avian flu rather than to the H1N1 counterparts in the Americas. Homology-based 3D structure modeling reveals that the novel mutations are preferentially located at the protein surface and do not interfere with the active site. The latter is the binding cavity for 3 currently used neuraminidase inhibitors: oseltamivir (Tamiflu^®), zanamivir (Relenza^®) and peramivir; thus, the drugs should remain effective for treatment. However, the antigenic regions of the neuraminidase relevant for vaccine development, serological typing and passive antibody treatment can differ from those of previous strains and already vary among patients.

Reviewers

This article was reviewed by Sandor Pongor and L. Aravind.     

Neuraminidase inhibitor R-125489--a promising drug for treating influenza virus: steered molecular dynamics approach

Two neuraminidase inhibitors, oseltamivir and zanamivir, are important drug treatments for influenza. Oseltamivir-resistant mutants of the influenza virus A/H1N1 and A/H5N1 have emerged, necessitating the development of new long-acting antiviral agents. One such agent is a new neuraminidase inhibitor R-125489 and its prodrug CS-8958. An atomic level understanding of the nature of this antiviral agents binding is still missing. We address this gap in our knowledge by applying steered molecular dynamics (SMD) simulations to different subtypes of seasonal and highly pathogenic influenza viruses. We show that, in agreement with experiments, R-125489 binds to neuraminidase more tightly than CS-8958. Based on results obtained by SMD and the molecular mechanics-Poisson–Boltzmann surface area method, we predict that R-125489 can be used to treat not only wild-type but also tamiflu-resistant N294S, H274Y variants of A/H5N1 virus as its binding affinity does not vary much across these systems. The high correlation level between theoretically determined rupture forces and experimental data on binding energies for the large number of systems studied here implies that SMD is a promising tool for drug design.     

Reassortment between Avian H5N1 and Human Influenza Viruses Is Mainly Restricted to the Matrix and Neuraminidase Gene Segments

Highly pathogenic avian influenza H5N1 viruses have devastated the poultry industry in many countries of the eastern hemisphere. Occasionally H5N1 viruses cross the species barrier and infect humans, sometimes with a severe clinical outcome. When this happens, there is a chance of reassortment between H5N1 and human influenza viruses. To assess the potential of H5N1 viruses to reassort with contemporary human influenza viruses (H1N1, H3N2 and pandemic H1N1), we used an in vitro selection method to generate reassortant viruses, that contained the H5 hemagglutinin gene, and that have a replication advantage in vitro. We found that the neuraminidase and matrix gene segments of human influenza viruses were preferentially selected by H5 viruses. However, these H5 reassortant viruses did not show a marked increase in replication in MDCK cells and human bronchial epithelial cells. In ferrets, inoculation with a mixture of H5N1-pandemic H1N1 reassortant viruses resulted in outgrowth of reassortant H5 viruses that had incorporated the neuraminidase and matrix gene segment of pandemic 2009 H1N1. This virus was not transmitted via aerosols or respiratory droplets to naïve recipient ferrets. Altogether, these data emphasize the potential of avian H5N1 viruses to reassort with contemporary human influenza viruses. The neuraminidase and matrix gene segments of human influenza viruses showed the highest genetic compatibility with HPAI H5N1 virus.     

I222 Neuraminidase Mutations Further Reduce Oseltamivir Susceptibility of Indonesian Clade 2.1 Highly Pathogenic Avian Influenza A(H5N1) Viruses

We have tested the susceptibility to neuraminidase inhibitors of 155 clade 2.1 H5N1 viruses from Indonesia, isolated between 2006–2008 as well as 12 clade 1 isolates from Thailand and Cambodia from 2004–2007 using a fluorometric MUNANA-based enzyme inhibition assay. The Thailand and Cambodian clade 1 isolates tested here were all susceptible to oseltamivir and zanamivir, and sequence comparison indicated that reduced oseltamivir susceptibility we observed previously with clade 1 Cambodian isolates correlated with an S246G neuraminidase mutation. Eight Indonesian viruses (5%), all bearing I222 neuraminidase mutations, were identified as mild to extreme outliers for oseltamivir based on statistical analysis by box plots. IC₅₀s were from 50 to 500-fold higher than the reference clade 1 virus from Viet Nam, ranging from 43–75 nM for I222T/V mutants and from 268–349 nM for I222M mutants. All eight viruses were from different geographic locales; all I222M variants were from central Sumatra. None of the H5N1 isolates tested demonstrated reduced susceptibility to zanamivir (IC₅₀s all <5 nM). All I222 mutants showed loss of slow binding specifically for oseltamivir in an IC₅₀ kinetics assay. We identified four other Indonesian isolates with higher IC₅₀s which also demonstrated loss of slow binding, including one virus with an I117V mutation. There was a minimal effect on the binding of zanamivir and peramivir for all isolates tested. As H5N1 remains a potential pandemic threat, the incidence of mutations conferring reduced oseltamivir susceptibility is concerning and emphasizes the need for greater surveillance of drug susceptibility.     

In Silico Design of an Oseltamivir Derivative with Increased Affinity against Wild-Type and Mutant Variants of Neuraminidase and Hemagglutinin of Influenza A H1N1 Virus

Antiviral resistance has turned into a world concern nowadays. Influenza A H1N1 emerged as a problem at the world level due to the neuraminidase (NA) mutations. The NA mutants conferred resistance to oseltamivir and zanamivir. Several efforts were conducted to develop better anti-influenza A H1N1 drugs. Our research group combined in silico methods to create a compound derived from oseltamivir to be tested in vitro against influenza A H1N1. Here we show the results of a new compound derived from oseltamivir but with specific chemical modifications, with significant affinity either on NA (in silico and in vitro assays) or HA (in silico) from influenza A H1N1 strain. We include docking and molecular dynamics (MD) simulations of the oseltamivir derivative at the binding site onto NA and HA of influenza A H1N1. Additionally, the biological experimental results show that oseltamivir derivative decreases the lytic-plaque formation on viral susceptibility assays, and it does not show cytotoxicity. Finally, oseltamivir derivative assayed on viral NA showed a concentration-dependent inhibition behavior at nM, depicting a high affinity of the compound for the enzyme, corroborated with the MD simulations results, placing our designed oseltamivir derivative as a potential antiviral against influenza A H1N1.     

Analogue inhibitors by modifying oseltamivir based on the crystal neuraminidase structure for treating drug-resistant H5N1 virus

The worldwide spread of H5N1 avian influenza and the increasing reports about its resistance to the existing drugs have made a priority for the development of the new anti-influenza molecules. The crystal structure of H5N1 avian influenza neuraminidase reported recently by Russell et al. [R.J. Russell, L.F. Haire, D.J. Stevens, P.J. Collins, Y. P. Lin, G.M. Blackburn, A.J. Hay, S.J. Gamblin, J.J. Skehel, The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design, Nature 443 (2006) 45-49] have provided new opportunities for drug design in this regard. It is revealed through the structure that the active sites of the group-1 neuraminidases, which contain the N1 subtype, have a very different three-dimensional structure from those of group-2 neuraminidases. The key difference is in the 150-loop cavity adjacent to the conserved active site in neuraminidase. Based on these findings and by modifying oseltamivir, six analog inhibitors were proposed as candidates for developing inhibitors against H5N1 virus, particularly against the oseltamivir-resistant H5N1 virus strain.     

Molecular Dynamics Simulations Suggest that Electrostatic Funnel Directs Binding of Tamiflu to Influenza N1 Neuraminidases

Oseltamivir (Tamiflu) is currently the frontline antiviral drug employed to fight the flu virus in infected individuals by inhibiting neuraminidase, a flu protein responsible for the release of newly synthesized virions. However, oseltamivir resistance has become a critical problem due to rapid mutation of the flu virus. Unfortunately, how mutations actually confer drug resistance is not well understood. In this study, we employ molecular dynamics (MD) and steered molecular dynamics (SMD) simulations, as well as graphics processing unit (GPU)-accelerated electrostatic mapping, to uncover the mechanism behind point mutation induced oseltamivir-resistance in both H5N1 “avian” and H1N1pdm “swine” flu N1-subtype neuraminidases. The simulations reveal an electrostatic binding funnel that plays a key role in directing oseltamivir into and out of its binding site on N1 neuraminidase. The binding pathway for oseltamivir suggests how mutations disrupt drug binding and how new drugs may circumvent the resistance mechanisms.     

Design,synthesis, and evaluation of carboxyl-modified oseltamivir derivatives with improved lipophilicity as neuraminidase inhibitors

In this study, a series of carboxyl-modified oseltamivir analogs with improved lipophilicity were designed and synthesized, and their inhibitory activities against neuraminidase from influenza A virus H5N1 subtype were evaluated. The results demonstrated that compound 5m exhibited potent inhibitory activity (IC₅₀?=?1.30?±?0.23?μM), and it targeted the recently discovered 430-cavity. Compound 5m (Log D?=??0.12) is more lipophilic than oseltamivir carboxylate (Log D?=??1.69) at pH 7.4, which is potentially propitious to improved membrane permeability and oral drug absorption. Meanwhile, 5m showed high stability in human liver microsomes. The findings of this study can be valuable in identifying neuraminidase inhibitors with optimal lipophilicity and in the exploration of 430-cavity.     

Changes in the Length of the Neuraminidase Stalk Region Impact H7N9 Virulence in Mice

The neuraminidase stalk of the newly emerged H7N9 influenza virus possesses a 5-amino-acid deletion. This study focuses on characterizing the biological functions of H7N9 with varied neuraminidase stalk lengths. Results indicate that the 5-amino-acid deletion had no impact on virus infectivity or replication in vitro or in vivo compared to that of a virus with a full-length stalk, but enhanced virulence in mice was observed for H7N9 encoding a 19- to 20-amino-acid deletion, suggesting that N9 stalk length impacts virulence in mammals, as N1 stalk length does.     

Heterosubtype Neutralizing Responses to Influenza A (H5N1) Viruses Are Mediated by Antibodies to Virus Haemagglutinin

Background

It is increasingly clear that influenza A infection induces cross-subtype neutralizing antibodies that may potentially confer protection against zoonotic infections. It is unclear whether this is mediated by antibodies to the neuraminidase (NA) or haemagglutinin (HA). We use pseudoviral particles (H5pp) coated with H5 haemagglutinin but not N1 neuraminidase to address this question. In this study, we investigate whether cross-neutralizing antibodies in persons unexposed to H5N1 is reactive to the H5 haemagglutinin.

Methodology/Principal Findings

We measured H5-neutralization antibody titers pre- and post-vaccination using the H5N1 micro-neutralization test (MN) and H5pp tests in subjects given seasonal vaccines and in selected sera from European elderly volunteers in a H5N1 vaccine trial who had detectable pre-vaccination H5N1 MN antibody titers. We found detectable (titer ≥20) H5N1 neutralizing antibodies in a minority of pre-seasonal vaccine sera and evidence of a serological response to H5N1 in others after seasonal influenza vaccination. There was excellent correlation in the antibody titers between the H5N1 MN and H5pp tests. Similar correlations were found between MN and H5pp in the pre-vaccine sera from the cohort of H5N1 vaccine trial recipients.

Conclusions/Significance

Heterosubtype neutralizing antibody to H5N1 in healthy volunteers unexposed to H5N1 is mediated by cross-reaction to the H5 haemagglutinin.     

Characterization of Pharmacologically Active Anti-Peptide Antibodies Directed Against the First and Second Extracellular Loops of the Serotonin 5-HT1A Receptor

Abstract: The immunological properties and the functional role of the first (loop I) and second (loop II) extracellular loops of the human serotonin 5-HT_1A receptor were studied with three populations of anti-peptide antibodies: Ab-1 (loop I; sequence Y-Q-V-L-N-K-W-T-L-G-Q-V-T-C-D-L; residues 96–111), Ab-2 (loop II; sequence G-W-R-T-P-E-D-R-S-D-P-D-A-C-T-I-S-K-D-H-G; residues 173–193), and Ab-12 (produced against loop I but cross-reacting with loop II). Chemical modification of peptide amino acid residues revealed the importance of the polyanionic stretch near the N-terminal domain of loop II for Ab-2 antibody binding and the role of the cysteine residues in both loops for the binding of Ab-1 and Ab-12 antibodies. Antibodies Ab-2 and Ab-12 recognized only the nonglycosylated form of the receptor (42 kDa) on immunoblots with transfected HeLa cells expressing the human 5-HT_1A receptor but recognized the glycosylated forms (55 and 65 kDa) of rat 5-HT_1A receptor from hippocampus membranes. The Ab-1 antibodies recognized no protein band from any cell type studied. Preincubation of transfected HeLa cell membranes with Ab-2 antibodies revealed two affinity binding sites of the 5-HT_1A receptor (K_DH = 0.54 ± 0.09 nM and K_DL = 13.74 ± 4.9 nM) for the agonist 8-hydroxy-2-(di-n-[³H]propylamino)tetralin ([³H]8-OH-DPAT) binding, but Ab-1 and Ab-12 revealed only one site (K_D of ≈2.5 nM). In contrast to the Ab-2 antibodies, Ab-1 and Ab-12 antibodies decreased the B_max of the [³H]8-OH-DPAT binding to 42 and 31%, respectively. These findings suggest that there are at least two epitopes on the extracellular loops: one inducing a high-affinity state for agonist binding and the other interfering with the accessibility of the ligand binding pocket.     

Structural characterization of the heme‐based oxygen sensor,AfGcHK, its interactions with the cognate response regulator,and their combined mechanism of action in a bacterial two‐component signaling system

The oxygen sensor histidine kinase AfGcHK from the bacterium Anaeromyxobacter sp. Fw 109‐5 forms a two‐component signal transduction system together with its cognate response regulator (RR). The binding of oxygen to the heme iron of its N‐terminal sensor domain causes the C‐terminal kinase domain of AfGcHK to autophosphorylate at His183 and then transfer this phosphate to Asp52 or Asp169 of the RR protein. Analytical ultracentrifugation revealed that AfGcHK and the RR protein form a complex with 2:1 stoichiometry. Hydrogen‐deuterium exchange coupled to mass spectrometry (HDX‐MS) suggested that the most flexible part of the whole AfGcHK protein is a loop that connects the two domains and that the heme distal side of AfGcHK, which is responsible for oxygen binding, is the only flexible part of the sensor domain. HDX‐MS studies on the AfGcHK:RR complex also showed that the N‐side of the H9 helix in the dimerization domain of the AfGcHK kinase domain interacts with the helix H1 and the β‐strand B2 area of the RR protein's Rec1 domain, and that the C‐side of the H8 helix region in the dimerization domain of the AfGcHK protein interacts mostly with the helix H5 and β‐strand B6 area of the Rec1 domain. The Rec1 domain containing the phosphorylable Asp52 of the RR protein probably has a significantly higher affinity for AfGcHK than the Rec2 domain. We speculate that phosphorylation at Asp52 changes the overall structure of RR such that the Rec2 area containing the second phosphorylation site (Asp169) can also interact with AfGcHK. Proteins 2016; 84:1375–1389. © 2016 Wiley Periodicals, Inc.