Oseltamivir resistance depends on the position 273 amino acid of neuraminidase of the type A influenza virus (H1N1), circulating in human population

The structure of neuraminidase of the type A influenza virus (H1N1) spreading in the human population was analyzed. The obtained results indicate a significant correlation between the oseltamivir sensitivity and the nature of the amino acid localized not only to neuraminidase position 274, but also to position 273 of this protein. Phenylalanine at position 273 in neuraminidase indicates a higher propensity to influenza virus mutation H274Y, leading to the appearance of resistant strains. It is suggested that the mutation at position 273 may be one of the characteristics allowing type A influenza virus to be ascribed to a pandemic or a seasonal type.     

Molecular distribution of amino acid substitutions on neuraminidase from the 2009 (H1N1) human influenza pandemic virus

The pandemic influenza AH1N1 (2009) caused an outbreak of human infection that spread to the world. Neuraminidase (NA) is an antigenic surface glycoprotein, which is essential to the influenza infection process, and is the target of anti-flu drugs oseltamivir and zanamivir. Currently, NA inhibitors are the pillar pharmacological strategy against seasonal and global influenza. Although mutations observed after NA-inhibitor treatment are characterized by changes in conserved amino acids of the enzyme catalytic site, it is possible that specific amino acid substitutions (AASs) distant from the active site such as H274Y, could confer oseltamivir or zanamivir resistance. To better understand the molecular distribution pattern of NA AASs, we analyzed NA AASs from all available reported pandemic AH1N1 NA sequences, including those reported from America, Africa, Asia, Europe, Oceania, and specifically from Mexico. The molecular distributions of the AASs were obtained at the secondary structure domain level for both the active and catalytic sites, and compared between geographic regions. Our results showed that NA AASs from America, Asia, Europe, Oceania and Mexico followed similar molecular distribution patterns. The compiled data of this study showed that highly conserved amino acids from the NA active site and catalytic site are indeed being affected by mutations. The reported NA AASs follow a similar molecular distribution pattern worldwide. Although most AASs are distributed distantly from the active site, this study shows the emergence of mutations affecting the previously conserved active and catalytic site. A significant number of unique AASs were reported simultaneously on different continents.     

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.     

Another look at the molecular mechanism of the resistance of H5N1 influenza A virus neuraminidase (NA) to oseltamivir (OTV)

In the context of a recent pandemic threat by the worldwide spread of H5N1 avian influenza, the high resistance of H5N1 virus to the most widely used commercial drug, oseltamivir (Tamiflu), is currently an important research topic. Herein, molecular bases of the mechanism of H5N1 NA resistance to oseltamivir were elucidated using a computational approach in a systematic fashion. Using the crystal structure of the complex of H5N1 NA with OTV (PDB ID: 2hu0) as the starting point, the question, how mutations at His274 by both smaller side chain (Gly, Ser, Asn, Gln) and larger side chain (Phe, Tyr) residues influence the sensitivity of N1 to oseltamivir, was addressed and correlated with the experimental data. The smaller side chain residue mutations of His274 resulted in slightly enhanced or unchanged NA sensitivity to OTV, while His274Phe and His274Tyr reduced the susceptibility of OTV to N1. In contrast to the binding free energies, the net charges of Glu276 and Arg224, making charge-charge interactions with Glu276, were established to be more sensitive to detecting subtle conformational differences induced at the key residue Glu276 by the His274X mutations. This study provides deeper insights into the possibility of developing viable drug-resistant mutants.     

Molecular dynamics simulation of oseltamivir resistance in neuraminidase of avian influenza H5N1 virus

The outbreak of avian influenza virus H5N1 has raised a global concern because of its high virulence and mutation rate. Although two classes of antiviral drugs, M2 ion channel protein inhibitors and neuraminidase inhibitors, are expected to be important in controlling the early stages of a potential pandemic. Different strains of influenza viruses have differing degrees of resistance against the antivirals. In order to analyze the detailed information on the viral resistance, molecular dynamics simulations were carried out for the neuraminidase (NA) complex with oseltamivir. The carboxylate of Glu276 of H252Y NA faces toward the O-ethyl-propyl group of oesltamivir, Glu276 of wild-type NA adopts a conformation pointing away from the oesltamivir. τ2 and τ3 torsional angles fluctuation of the oesltamivir are relatively high for the H252Y mutant NA complex. In addition, there are fewer hydrogen bonds between the oesltamivir and H252Y mutation NA. The results show that H252Y mutation NA has high resistance against the drug.     

Effect of an asparagine-to-serine mutation at position 294 in neuraminidase on the pathogenicity of highly pathogenic H5N1 influenza A virus

Like the histidine-to-tyrosine substitution at position 274 in neuraminidase (NA H274Y), an asparagine-to-serine mutation at position 294 in this protein (NA N294S) confers oseltamivir resistance to highly pathogenic H5N1 influenza A viruses. However, unlike viruses with the NA H274Y mutation, the properties of viruses possessing NA N294S are not well understood. Here, we assessed the effect of the NA N294S substitution on the replication and pathogenicity of human H5N1 viruses and on the efficacy of the NA inhibitors oseltamivir and zanamivir in mouse and ferret models. Although NA N294S-possessing H5N1 viruses were attenuated in mice and ferrets compared to their oseltamivir-sensitive counterparts, one of the infected ferrets died from systemic infection, demonstrating the potential lethality in ferrets of oseltamivir-resistant H5N1 viruses with the NA N294S substitution. The efficacy of oseltamivir, but not that of zanamivir, against an NA N294S-possessing virus was substantially impaired both in ferrets and in vitro. These results demonstrate the considerable pathogenicity of NA N294S substitution-possessing H5N1 viruses and underscore the importance of monitoring the emergence of the NA N294S mutation in circulating H5N1 viruses.     

Receptor recognition mechanism of human influenza A H1N1 (1918), avian influenza A H5N1 (2004), and pandemic H1N1 (2009) neuraminidase

Influenza A neuraminidase (NA) is a target for anti-influenza drugs. The function of this enzyme is to cleave a glycosidic linkage of a host cell receptor that links sialic acid (Sia) to galactose (Gal), to allow the virus to leave an infected cell and propagate. The receptor is an oligosaccharide on the host cell surface. There are two types of oligosaccharide receptor; the first, which is found mainly on avian epithelial cell surfaces, links Sia with Gal by an α2,3 glycosidic linkage; in the second, found mainly on human epithelial cell surfaces, linkage is via an α2,6 linkage. Some researchers believe that NAs from different viruses show selectivity for each type of linkage, but there is limited information available to confirm this hypothesis. To see if the linkage type is more specific to any particular NA, a number of NA-receptor complexes of human influenza A H1N1 (1918), avian influenza A H5N1 (2004), and a pandemic strain of H1N1 (2009) were constructed using homology modeling and molecular dynamics simulation. The results show that the two types of receptor analogues bound to NAs use different mechanisms. Moreover, it was found that a residue unique to avian virus NA is responsible for the recognition of the Siaα2,3Gal receptor, and a residue unique to human virus NA is responsible for the recognition of Siaα2,6Gal. We believe that this finding could explain how NAs of different virus origins always possess some unique residues.     

Computational study on new natural polycyclic compounds of H1N1 influenza virus neuraminidase

A new strain of influenza A (H1N1) virus is a major cause of morbidity and mortality around the world. The neuraminidase of the influenza virus has been the most potential target for the anti-influenza drugs such as oseltamivir and zanamivir. However, the emergence of drug-resistant variants of these drugs makes a pressing need for the development of new neuraminidase inhibitors for controlling illness and transmission. Here a 3D structure model of H1N1 avian influenza virus neuraminidase type 1 (N1) was constructed based on the structure of the template H5N1 avian influenza virus N1. Upon application of virtual screening technique for N1 inhibitors, two novel compounds (ZINC database ID: ZINC02128091, ZINC02098378) were found as the most favorable interaction energy with N1. Docking results showed that the compounds bound not only in the active pocket, but also in a new hydrophobic cave which contains Arg368, Trp399, Ile427, Pro431 and Lys432 of N1. Our result suggested that both of the screened compounds containing the hydrophobic group bring a strong conjugation effect with Arg293, Arg368 Lys432 of N1 by pi-pi interaction. However, the control inhibitors zanamivir and oseltamivir do not have this effect. The details of N1-compound binding structure obtained will be valuable for the development of a new anti-influenza virus agent.     

On the structure-based design of novel inhibitors of H5N1 influenza A virus neuraminidase (NA)

The structure-based design of novel H5N1 neuraminidase inhibitors is currently a research topic of vital importance owing to both a recent pandemic threat by the worldwide spread of H5N1 avian influenza and the high resistance of H5N1 virus to the most widely used commercial drug, oseltamivir-OTV (Tamiflu). A specific criterion used in this work for determining fully acceptable conformations of potential inhibitors is a previous experimental proposal of exploiting potential benefits for drug design offered by the ‘150-cavity’ adjacent to the NA active site. Using the crystal structure of H5N1 NA (PDB ID: 2hty) as the starting point, in a set of 54 inhibitors previously proposed by modifying the side chains of oseltamivir, 4 inhibitors were identified using two different computational strategies (ArgusLab4.0.1, FlexX-E3.0.1) both to lower the binding free energy (BFE) of oseltamivir and to have partially acceptable conformations. These 4 oseltamivr structure-based analogues were found to adopt the most promising conformations by identifying the guanidinium side chain of Arg156 as a prospective partner for making polar contacts, but none of the modified 4-amino groups of oseltamivir in the 4 favorable conformations was found to make polar contacts with the guanidinium side chain of Arg156. Hence, the structures of two additional inhibitors were designed and shown to further lower the binding free energy of OTV relative to the previous 54 inhibitors. These two novel structures clearly suggest that it may be possible for a new substituent to be developed by functional modifications at position of the 4-amino group of oseltamivir in order to make polar contacts with the guanidinium side chain of Arg156, and thereby enhance the binding of a more potent inhibitor. Several standpoints of vital importance for designing novel structures of potentially more effective H5N1 NA inhibitors are established.     

Functional association between influenza A (H1N1) virus and human

Influenza A (H1N1) virus is a severe threat worldwide. It is important to gain a better understanding of the mechanism of the infection. In the paper, we established a computational framework to investigate the crosstalk between the virus and the host, by finding out the proteins that the virus is attacking. The targeted proteins were predicted by taking human proteins laid on the same GO functions or processes as the virus proteins. One hundred and one core proteins were identified. The results provide some knowledge of the possible biological processes and molecular interactions caused by the viral infection, including the host responses.     

Studies on neuraminidase from influenza virus A(H3N2) obtained by two procedures

• 1.1. Neuraminidase was obtained by (A) bromelain solubilization or (B) by treatment with N-lauroylsarcosine.
• 2.2. 5-N-acetyl-2-O(3-methoxyphenyl)-α-d-neuraminic acid, employed as substrate, avoids the interference produced by the thiobarbituric acid method, and is not interfered by the ampholytes.
• 3.3. Only about 20% of original enzyme activity was lost after electrofocusing. The sample from procedure A showed two peaks, corresponding to pis 4.4 and 5.6. The sample from procedure B, having a higher activity, showed only one peak at pI 4.4.
• 4.4. Samples A and B showed different K_m and hydrolysis rate with N-acetylneuraminyl-lactose and glycophorin A. It was not found significantly different with other substrates: α₁-acid glycoprotein, brain gangliosides. 5-N-acetyl-2-O-(3-methoxyphenyl)-α-d-neuraminic acid and 2'-(4-methyl umbelliferyl)-α-d-N-acetylneuraminic acid.

     

Decreased neuraminidase activity is important for the adaptation of H5N1 influenza virus to human airway epithelium

Highly pathogenic avian H5N1 influenza viruses remain a pandemic threat. Antiviral drugs such as neuraminidase (NA) inhibitors will be crucial for disease control in the event of a pandemic. Should drug-resistant H5N1 viruses develop, all defense strategies will be compromised. To determine the likelihood and mechanisms of emergence of NA inhibitor-resistant H5N1 variants in humans, we serially passaged two H5N1 viruses, A/Hong Kong/213/03 and A/Turkey/65-1242/06, in normal human bronchial epithelial (NHBE) cells in the presence of oseltamivir, zanamivir, or peramivir. To monitor the emergence of changes associated with the adaptation of H5N1 viruses to humans, we passaged the strains in the absence of drugs. Under pressure of each NA inhibitor, A/Turkey/65-1242/06 developed mutations in the hemagglutinin (HA) (H28R and P194L/T215I) and NA (E119A) proteins that reduced virus binding to α2,3-sialyl receptor and NA activity. Oseltamivir pressure selected a variant of A/Hong Kong/213/03 virus with HA P194S mutation that decreased viral binding to α2,6 receptor. Under peramivir pressure, A/Hong Kong/213/03 virus developed a novel NA mutation, R156K, that reduced binding to all three drugs, caused about 90% loss of NA activity, and compromised replication in NHBE cells. Both strains were eliminated in NHBE cells when they were cultivated in the absence of drugs. Here, we show for the first time that decreased NA activity mediated through NA inhibitors is essential for the adaptation of pandemic H5N1 influenza virus to humans. This ability of decreased NA activity to promote H5N1 infection underlines the necessity to optimize management strategies for a plausible H5N1 pandemic.     

The nature of amino acid substitution in antigenic drift of hemagglutinin N3 and neuraminidase N2 from the influenza virus

The nature of amino acid replacements in 16 drift variants of hemagglutinin H3 subtype and 5 drift variants of neuraminidase N2 subtype of the influenza A virus were studied. The dependences of relative replacement frequencies and relative quantities of frequent replacements upon differences of properties of substituted residues are plotted. In contrast to most of the known proteins, amino acid replacements in hemagglutinin and neuraminidase depend weakly on the physico-chemical parameters of amino acids. For the antigenic determinants studied the replacement frequencies were compared to those calculated according to two models: one for conservative replacements and the other for accidental mutation of the genetic code. The differences in the nature of amino acid replacements are found in four antigenic determinants of hemagglutinin. The replacements in experimentally selected proteins are shown to go beyond limitations of natural variants. The explanations of the reasons of low epidemicity of some strains and ineffective attempt to imitate the natural antigenic drift of viruses by using experimental selection are proposed. The causes of time-limited circulation of H3N2 influenza virus subtype are discussed.     

Chalcones as novel influenza A (H1N1) neuraminidase inhibitors from Glycyrrhiza inflata

The emergence of highly pathogenic influenza A virus strains, such as the new H1N1 swine influenza (novel influenza), represents a serious threat to global human health. During our course of an anti-influenza screening program on natural products, one new licochalcone G (1) and seven known (2-8) chalcones were isolated as active principles from the acetone extract of Glycyrrhiza inflata. Compounds 3 and 6 without prenyl group showed strong inhibitory effects on various neuraminidases from influenza viral strains, H1N1, H9N2, novel H1N1 (WT), and oseltamivir-resistant novel H1N1 (H274Y) expressed in 293T cells. In addition, the efficacy of oseltamivir with the presence of compound 3 (5 μM) was increased against H274Y neuraminidase. This evidence of synergistic effect makes this inhibitor to have a potential possibility for control of pandemic infection by oseltamivir-resistant influenza virus.     

An influenza A (H1N1) virus, closely related to swine influenza virus, responsible for a fatal case of human influenza.

In July 1991, an influenza A virus, designated A/Maryland/12/91 (A/MD), was isolated from the bronchial secretions of a 27-year-old animal caretaker. He had been admitted to the hospital with bilateral pneumonia and died of acute respiratory distress syndrome 13 days later. Antigenic analyses with postinfection ferret antisera and monoclonal antibodies to recent H1 swine hemagglutinins indicated that the hemagglutinin of this virus was antigenically related to, but distinguishable from, those of other influenza A (H1N1) viruses currently circulating in swine. Oligonucleotide mapping of total viral RNAs revealed differences between A/MD and other contemporary swine viruses. However, partial sequencing of each RNA segment of A/MD demonstrated that all segments were related to those of currently circulating swine viruses. Sequence analysis of the entire hemagglutinin, nucleoprotein, and matrix genes of A/MD revealed a high level of identity with other contemporary swine viruses. Our studies on A/MD emphasize that H1N1 viruses in pigs obviously continue to cross species barriers and infect humans.     

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.     

Computational studies of H5N1 influenza virus resistance to oseltamivir

Influenza A (H5N1) virus is one of the world's greatest pandemic threats. Neuraminidase (NA) inhibitors, oseltamivir and zanamivir, prevent the spread of influenza, but drug‐resistant viruses have reduced their effectiveness. Resistance depends on the binding properties of NA‐drug complexes. Key residue mutations within the active site of NA glycoproteins diminish binding, thereby resulting in drug resistance. We performed molecular simulations and calculations to characterize the mechanisms of H5N1 influenza virus resistance to oseltamivir and predict potential drug‐resistant mutations. We examined two resistant NA mutations, H274Y and N294S, and one non‐drug‐resistant mutation, E119G. Six‐nanosecond unrestrained molecular dynamic simulations with explicit solvent were performed using NA‐oseltamivir complexes containing either NA wild‐type H5N1 virus or a variant. MM_PBSA techniques were then used to rank the binding free energies of these complexes. Detailed analyses indicated that conformational change of E276 in the Pocket 1 region of NA is a key source of drug resistance in the H274Y mutant but not in the N294S mutant.