Prediction of mutations engineered by randomness in H5N1 hemagglutinins of influenza A virus

This is the continuation of our studies on the prediction of mutation engineered by randomness in proteins from influenza A virus. In our previous studies, we have demonstrated that randomness plays a role in engineering mutations because the measures of randomness in protein are different before and after mutations. Thus we built a cause-mutation relationship to count the mutation engineered by randomness, and conducted several concept-initiated studies to predict the mutations in proteins from influenza A virus, which demonstrated the possibility of prediction of mutations along this line of thought. On the other hand, these concept-initiated studies indicate the directions forwards the enhancement of predictability, of which we need to use the neural network instead of logistic regression that was used in those concept-initiated studies to enhance the predictability. In this proof-of-concept study, we attempt to apply the neural network to modeling the cause-mutation relationship to predict the possible mutation positions, and then we use the amino acid mutating probability to predict the would-be-mutated amino acids at predicted positions. The results confirm the possibility of use of internal cause-mutation relationship with neural network model to predict the mutation positions and use of amino acid mutating probability to predict the would-be-mutated amino acids.     

Prediction of Mutations Initiated by Internal Power in H3N2 Hemagglutinins of Influenza A Virus from North America

Probably the best way to predict mutations is to find the cause for mutations, by which the cause–mutation relationship can be built. However, many causes which have resulted in mutations in the past might not leave any trace due to the changes in environments. As well, the current proteins may not be sensitive to the causes, which led to mutations in the past, because of evolution. Thus we might have recorded many mutations, but few of their corresponding causes, and it would be difficult to establish the one-to-one cause–mutation relationship. However, the internal power engineering mutations within a protein would exist, of which randomness should play an important role. Since 1999, we have developed three methods to quantify the randomness within a protein by which we can build a cause–mutation relationship because we can classify the occurrence and non-occurrence of mutation as unity and zero, and transfer this relationship into the classification problem, which can be solved using logistic regression. Recently, we used the logistic regression to predict the mutation positions in H5N1 hemagglutinins from influenza A virus, and applied the amino-acid mutating probability to predict the would-be-mutated amino acids at predicted positions as the concept-initiated study. However, we still need to conduct many proof-of-concept studies to test whether this cause–mutation relationship is independent of protein subtypes, whether the logistic regression is powerful enough, etc. In this study, we attempted to use the logistic regression to predict the mutation positions in H3N2 hemagglutinins of influenza A virus from North America to answer the questions in the proof-of-concept stage.     

Prediction of mutation trend in hemagglutinins and neuraminidases from influenza A viruses by means of cross-impact analysis

In this study, we use the cross-impact analysis to define the relationship among impact, mutation, and outbreak of bird flu. Then we use the distribution rank, which is developed by us over last several years, to quantify the mutations from amino acid sequences of 134 hemagglutinins and 97 neuraminidases. With the help of Bayesian equation, we calculate the probability of occurring of mutation in H5, H6, and H9 hemagglutinins, and N1 and N2 neuraminidases. Finally, we estimate the probability of occurring of mutation with different intensities of an impact. Although we have no means to predict an impact, which is severe enough to lead to the mutations in hemagglutinins and neuraminidases resulting in the outbreak of bird flu, we can in principle monitor the changes in distribution rank along the time course, and predict the trend of mutations, even to predict the degree of outbreak of bird flu.     

Prediction of mutations in H5N1 hemagglutinins from influenza A virus

In this study, we determine the mutation relation among 333 H5N1 hemagglutinins of influenza A viruses according to their amino acid and RNA codon sequences. Then, we calculate seven probabilistic numbers, which have been developed by us since 1999, for each amino acid in these hemagglutinins. With the seven numeric numbers as independents and the probability of occurrence of mutation at each hemagglutinin position as dependent, we use the logistic regression to model 967 missense point mutations from 333 hemagglutinins to get the population estimates. Thereafter, we predict the future mutation positions in H5N1 hemagglutinin. Finally, we use the translation probabilities between RNA codons and mutated amino acids to predict the would-be-mutated amino acids in H5N1 hemagglutinin.     

Improvement of model for prediction of hemagglutinin mutations in H5N1 influenza viruses with distinguishing of arginine, leucine and serine

In a continuation of our attempt to predict mutations in proteins from influenza A virus, this study attempted to answer the question of whether distinguishing between arginine, leucine and serine can improve the predictability as these residues are governed by different probabilistic mechanism translating from RNA codons to amino acids. In this study, we made the prediction based on the mutation relation among 299 H5N1 hemagglutinins of influenza A virus. Then, we compared the results based on the distinguishing of arginine, leucine and serine with the results without distinguishing of arginine, leucine and serine. The results show that the prediction together with distinguishing between arginine, leucine and serine is better than prediction without distinguishing between these residues.     

Fate of 130 hemagglutinins from different influenza A viruses

In this study, we use our probabilistic models to analyze 130 hemagglutinins from different influenza A virus in order to gain the insight into their fate. The results provide three lines of evidence regarding the H5, H6, and H9 hemagglutinins: (i) the H5 hemagglutinins are more sensitive to mutations, this is the current state of the H5, H6, and H9 hemagglutinins; (ii) the H5 hemagglutinins had experienced more mutations in the past, this is the history of the H5, H6, and H9 hemagglutinins; and (iii) the H6 hemagglutinins has a bigger potential towards future mutations, this is the future of the H5, H6, and H9 hemagglutinins. Furthermore, this study gives two clues on the mutation tendency that is a degeneration process and the species susceptibility that is the chickens and ducks.     

Determination of mutation trend in hemagglutinins by means of translation probability between RNA codons and mutated amino acids

In this study, we used the 183 translation probabilities between RNA codons and mutated amino acids to construct the theoretical distributions of mutated amino acids in hemagglutinins of influenza A virus. We then compared the actual distributions of mutated amino acids from 953 hemagglutinins with their theoretical ones. The results demonstrated that mutated amino acids generally follow the direction of the theoretical distributions governed by RNA codons. This, in turn, highlights the mutation trend of amino acids in hemagglutinins and provides a method for estimating possible mutations in a protein according to its theoretical distributions of mutated amino acids.     

Three sampling strategies to predict mutations in H5N1 hemagglutitins from influenza A virus

After several studies on prediction of mutation, we examine the effect of three sampling strategies, the sampling based on years, the sampling based on number of mutations, and the sampling based on the unpredictable portion of amino-acid pairs, on the prediction performance in H5N1 hemagglutinins. The results show that the sampling strategy does play an important role in prediction, which should be taken into account when predicting the next generation of mutations in proteins from influenza A virus.     

In silico analysis of drug-resistant mutant of neuraminidase (N294S) against oseltamivir

The recent H1N1 influenza pandemic has attracted worldwide attention due to the high infection rate. Oseltamivir is a new class of anti-viral agent approved for the treatment and prevention of influenza infections. The principal target for this drug is a virus surface glycoprotein, neuraminidase (NA), which facilitates the release of nascent virus and thus spreads infection. Until recently, only a low prevalence of neuraminidase inhibitor (NAI) resistance (<1 %) had been detected in circulating viruses. However, there have been reports of significant numbers of A (H1N1) influenza strains with a N294S neuraminidase mutation that was highly resistant to the NAI, oseltamivir. Hence, in the present study, we highlight the effect of point mutation-induced oseltamivir resistance in H1N1 subtype neuraminidases by molecular simulation approach. The docking analysis reveals that mutation (N294S) significantly affects the binding affinity of oseltamivir with mutant type NA. This is mainly due to the decrease in the flexibility of binding site residues and the difference in prevalence of hydrogen bonds in the wild and mutant structures. This study throws light on the possible effects of drug-resistant mutations on the large functionally important collective motions in biological systems.     

Hemagglutinin Stalk-Based Universal Vaccine Constructs Protect against Group 2 Influenza A Viruses

Current influenza virus vaccines contain H1N1 (phylogenetic group 1 hemagglutinin), H3N2 (phylogenetic group 2 hemagglutinin), and influenza B virus components. These vaccines induce good protection against closely matched strains by predominantly eliciting antibodies against the membrane distal globular head domain of their respective viral hemagglutinins. This domain, however, undergoes rapid antigenic drift, allowing the virus to escape neutralizing antibody responses. The membrane proximal stalk domain of the hemagglutinin is much more conserved compared to the head domain. In recent years, a growing collection of antibodies that neutralize a broad range of influenza virus strains and subtypes by binding to this domain has been isolated. Here, we demonstrate that a vaccination strategy based on the stalk domain of the H3 hemagglutinin (group 2) induces in mice broadly neutralizing anti-stalk antibodies that are highly cross-reactive to heterologous H3, H10, H14, H15, and H7 (derived from the novel Chinese H7N9 virus) hemagglutinins. Furthermore, we demonstrate that these antibodies confer broad protection against influenza viruses expressing various group 2 hemagglutinins, including an H7 subtype. Through passive transfer experiments, we show that the protection is mediated mainly by neutralizing antibodies against the stalk domain. Our data suggest that, in mice, a vaccine strategy based on the hemagglutinin stalk domain can protect against viruses expressing divergent group 2 hemagglutinins.     

Resistance Mutation R292K Is Induced in Influenza A(H6N2) Virus by Exposure of Infected Mallards to Low Levels of Oseltamivir

Resistance to neuraminidase inhibitors (NAIs) is problematic as these drugs constitute the major treatment option for severe influenza. Extensive use of the NAI oseltamivir (Tamiflu®) results in up to 865 ng/L of its active metabolite oseltamivir carboxylate (OC) in river water. There one of the natural reservoirs of influenza A, dabbling ducks, can be exposed. We previously demonstrated that an influenza A(H1N1) virus in mallards (Anas platyrhynchos) exposed to 1 µg/L of OC developed oseltamivir resistance through the mutation H274Y (N2-numbering). In this study, we assessed the resistance development in an A(H6N2) virus, which belongs to the phylogenetic N2 group of neuraminidases with distinct functional and resistance characteristics. Mallards were infected with A(H6N2) while exposed to 120 ng/L, 1.2 µg/L or 12 µg/L of OC in their sole water source. After 4 days with 12 µg/L of OC exposure, the resistance mutation R292K emerged and then persisted. Drug sensitivity was decreased ≈13,000-fold for OC and ≈7.8-fold for zanamivir. Viral shedding was similar when comparing R292K and wild-type virus indicating sustained replication and transmission. Reduced neuraminidase activity and decrease in recovered virus after propagation in embryonated hen eggs was observed in R292K viruses. The initial, but not the later R292K isolates reverted to wild-type during egg-propagation, suggesting a stabilization of the mutation, possibly through additional mutations in the neuraminidase (D113N or D141N) or hemagglutinin (E216K). Our results indicate a risk for OC resistance development also in a N2 group influenza virus and that exposure to one NAI can result in a decreased sensitivity to other NAIs as well. If established in influenza viruses circulating among wild birds, the resistance could spread to humans via re-assortment or direct transmission. This could potentially cause an oseltamivir-resistant pandemic; a serious health concern as preparedness plans rely heavily on oseltamivir before vaccines can be mass-produced.     

一株人感染高致病性禽流感(H5N1)病毒基因组分子特征分析

【背景】H5N1禽流感病毒可以感染人类导致重症呼吸道感染,致死率高。【目的】研究我中心确认的一例人感染高致病性禽流感H5N1病毒A/Nanjing/1/2015的可能起源及基因组分子特征。【方法】对病人痰液样本中的H5N1病毒进行全基因组测序,使用CLC Genomics Workbench 9.0对序列进行拼接,使用BLAST和MEGA 5.22软件进行同源性比对和各片段分子特征分析。【结果】该株禽流感病毒属于H5亚型的2.3.2.1c家系,其8个片段均与江浙地区禽类中分离的病毒高度同源,未发现有明显的重配。分子特征显示,该病毒血凝素(Hemagglutinin,HA)蛋白裂解位点为PQRERRRR/G,受体结合位点呈现禽类受体特点,但出现D94N、S133A和T188I氨基酸置换增强了病毒对人类受体的亲和性。神经氨酸酶(Neuraminidase,NA)蛋白颈部在49-68位缺失20个氨基酸,非结构蛋白1 (Non-structure protein,NS1)存在P42S置换和80-84位氨基酸的缺失。其他蛋白中也存在多个增强病毒致病力和对人类细胞亲和力的氨基酸突变。对耐药位点分析发现存在对奥司他韦的耐药突变H_274Y,病毒对金刚烷胺仍旧敏感。【结论】人感染高致病性禽流感H5N1病毒A/Nanjing/1/2015属于2.3.2.1c家系,禽类来源,关键位点较保守,但仍出现了多个氨基酸的进化与变异使其更利于感染人类。H5N1禽流感病毒进化活跃,持续动态监测不能放松。     

Matrix protein 1: A comparative in silico study on different strains of influenza A H5N1 Virus

The importance of influenza viruses as worldwide infectious agents is well recognized. Specific mutations and evolution in influenza viruses is difficult to predict. We studied specific mutations in matrix protein 1 (M1) of H5N1 influenza A virus together with properties associated with it using prediction tools developed in Bioinformatics. Changes in hydrophobicity, polarity and secondary structure at the site of mutation were noticed and documented to gain insight towards its infection.     

Synthesis, cloning and determination of the primary structure of a full-size DNA copy of the neuraminidase gene from influenza virus type A subtype H1N1

Complete nucleotide sequence of the cloned full-length DNA copy of the influenza virus A (H1N1) neuraminidase gene has been determined. The predicted amino acid sequence is compared with sequences of neuraminidases from other influenza virus strains. A section of the neuraminidase is found to be homologous to the chicken lysozyme catalytic centre.     

Transmission of Influenza Virus in a Mammalian Host Is Increased by PB2 Amino Acids 627K or 627E/701N

Since 2003, more than 380 cases of H5N1 influenza virus infection of humans have been reported. Although the resultant disease in these cases was often severe or fatal, transmission of avian influenza viruses between humans is rare. The precise nature of the barrier blocking human-to-human spread is unknown. It is clear, however, that efficient human-to-human transmission of an antigenically novel influenza virus would result in a pandemic. Influenza viruses with changes at amino acids 627 or 701 of the PB2 protein have been isolated from human cases of highly pathogenic H5 and H7 avian influenza. Herein, we have used the guinea pig model to test the contributions of PB2 627 and 701 to mammalian transmission. To this end, viruses carrying mutations at these positions were generated in the A/Panama/2007/99 (H3N2) and A/Viet Nam/1203/04 (H5N1) backgrounds. In the context of either rPan99 or rVN1203, mutation of lysine 627 to the avian consensus residue glutamic acid was found to decrease transmission. Introduction of an asparagine at position 701, in conjunction with the K627E mutation, resulted in a phenotype more similar to that of the parental strains, suggesting that this residue can compensate for the lack of 627K in terms of increasing transmission in mammals. Thus, our data show that PB2 amino acids 627 and 701 are determinants of mammalian inter-host transmission in diverse virus backgrounds.     

Mutations in H5N1 Influenza Virus Hemagglutinin that Confer Binding to Human Tracheal Airway Epithelium

The emergence in 2009 of a swine-origin H1N1 influenza virus as the first pandemic of the 21st Century is a timely reminder of the international public health impact of influenza viruses, even those associated with mild disease. The widespread distribution of highly pathogenic H5N1 influenza virus in the avian population has spawned concern that it may give rise to a human influenza pandemic. The mortality rate associated with occasional human infection by H5N1 virus approximates 60%, suggesting that an H5N1 pandemic would be devastating to global health and economy. To date, the H5N1 virus has not acquired the propensity to transmit efficiently between humans. The reasons behind this are unclear, especially given the high mutation rate associated with influenza virus replication. Here we used a panel of recombinant H5 hemagglutinin (HA) variants to demonstrate the potential for H5 HA to bind human airway epithelium, the predominant target tissue for influenza virus infection and spread. While parental H5 HA exhibited limited binding to human tracheal epithelium, introduction of selected mutations converted the binding profile to that of a current human influenza strain HA. Strikingly, these amino-acid changes required multiple simultaneous mutations in the genomes of naturally occurring H5 isolates. Moreover, H5 HAs bearing intermediate sequences failed to bind airway tissues and likely represent mutations that are an evolutionary “dead end.” We conclude that, although genetic changes that adapt H5 to human airways can be demonstrated, they may not readily arise during natural virus replication. This genetic barrier limits the likelihood that current H5 viruses will originate a human pandemic.     

Persistent host markers in pandemic and H5N1 influenza viruses

Avian influenza viruses have adapted to human hosts, causing pandemics in humans. The key host-specific amino acid mutations required for an avian influenza virus to function in humans are unknown. Through multiple-sequence alignment and statistical testing of each aligned amino acid, we identified markers that discriminate human influenza viruses from avian influenza viruses. We applied strict thresholds to select only markers which are highly preserved in human influenza virus isolates over time. We found that a subset of these persistent host markers exist in all human pandemic influenza virus sequences from 1918, 1957, and 1968, while others are acquired as the virus becomes a seasonal influenza virus. We also show that human H5N1 influenza viruses are significantly more likely to contain the amino acid predominant in human strains for a few persistent host markers than avian H5N1 influenza viruses. This sporadic enrichment of amino acids present in human-hosted viruses may indicate that some H5N1 viruses have made modest adaptations to their new hosts in the recent past. The markers reported here should be useful in monitoring potential pandemic influenza viruses.     

Structural differences among hemagglutinins of influenza A virus subtypes are reflected in their antigenic architecture: analysis of H9 escape mutants

We used a panel of monoclonal antibodies to H9 hemagglutinin to select 18 escape mutants of mouse-adapted influenza A/Swine/Hong Kong/9/98 (H9N2) virus. Cross-reactions of the mutants with the antibodies and the sequencing of hemagglutinin genes revealed two minimally overlapping epitopes. We mapped the amino acid changes to two areas of the recently reported three-dimensional structure of A/Swine/Hong Kong/9/98 hemagglutinin. The grouping of the antigenically relevant amino acid positions in H9 hemagglutinin differs from the pattern observed in H3 and H5 hemagglutinins. Several positions in site B of H3 hemagglutinin are distributed in two sites of H9 hemagglutinin. Unlike any subtype analyzed so far, H9 hemagglutinin does not contain an antigenic site corresponding to site A in H3 hemagglutinin. Positions 145 and 193 (H3 numbering), which in H3 hemagglutinin belong to sites A and B, respectively, are within one site in H9 hemagglutinin. This finding is consistent with the peculiarity of the three-dimensional structure of the H9 molecule, that is, the absence from H9 hemagglutinin of the lateral loop that forms site A in H3 and the equivalent site in H5 hemagglutinins. The escape mutants analyzed displayed phenotypic variations, including decreased virulence for mice and changes in affinity for sialyl substrates. Our results demonstrate a correlation between intersubtype differences in three-dimensional structure and variations among subtypes in the distribution of antigenic areas. Our findings also suggest that covariation and pleiotropic effects of antibody-selected mutations may be important in the evolution of H9 influenza virus, a possible causative agent of a future pandemic.     

禽流感病毒H5N1全基因组克隆与分析

为明确广东地区分离的一株禽流感病毒H5N1的遗传背景,建立流感病毒反向遗传的平台。对该株禽流感病毒进行了空斑纯化与组织细胞培养,检测其在MDCK细胞中的增殖特性;利用H5N1病毒通用引物,通过RT-PCR对该病毒全基因组的8条片段进行全长克隆及测序分析;将H5N1的8条全长基因组片段分别插入反向遗传通用载体中,构建禽流感病毒H5N1的感染性克隆。结果表明,该H5N1毒株在MDCK细胞中可不依赖胰酶进行有效增殖与复制,可使MDCK细胞出现典型细胞病变,具有高致病性禽流感病毒的细胞增殖特征。RT-PCR克隆得到该H5N1毒株的PB2、PB1、PA、HA、NP、NA、M和NS八条全长片段,经测序分析确认该毒株的基因序列,其内部编码序列出现多处突变,其中HA连接肽为多个连续碱性氨基酸,表明该毒株可不依赖胰酶进行有效复制,与细胞培养结果一致,未出现抗药性的遗传突变。PCR与测序证明,插入H5N1八个全长基因组片段的载体序列完全正确,表明成功构建了该毒株的感染性克隆。为明确病毒遗传信息,建立流感病毒反向遗传的平台,为进一步研究禽流感病毒相关疫苗提供了研究基础。     

Impact of the H275Y and I223V Mutations in the Neuraminidase of the 2009 Pandemic Influenza Virus In Vitro and Evaluating Experimental Reproducibility

The 2009 pandemic H1N1 (H1N1pdm09) influenza virus is naturally susceptible to neuraminidase (NA) inhibitors, but mutations in the NA protein can cause oseltamivir resistance. The H275Y and I223V amino acid substitutions in the NA of the H1N1pdm09 influenza strain have been separately observed in patients exhibiting oseltamivir-resistance. Here, we apply mathematical modelling techniques to compare the fitness of the wild-type H1N1pdm09 strain relative to each of these two mutants. We find that both the H275Y and I223V mutations in the H1N1pdm09 background significantly lengthen the duration of the eclipse phase (by 2.5 h and 3.6 h, respectively), consistent with these NA mutations delaying the release of viral progeny from newly infected cells. Cells infected by H1N1pdm09 virus carrying the I223V mutation display a disadvantageous, shorter infectious lifespan (17 h shorter) than those infected with the wild-type or MUT-H275Y strains. In terms of compensating traits, the H275Y mutation in the H1N1pdm09 background results in increased virus infectiousness, as we reported previously, whereas the I223V exhibits none, leaving it overall less fit than both its wild-type counterpart and the MUT-H275Y strain. Using computer simulated competition experiments, we determine that in the presence of oseltamivir at doses even below standard therapy, both the MUT-H275Y and MUT-I223V dominate their wild-type counterpart in all aspects, and the MUT-H275Y outcompetes the MUT-I223V. The H275Y mutation should therefore be more commonly observed than the I223V mutation in circulating H1N1pdm09 strains, assuming both mutations have a similar impact or no significant impact on between-host transmission. We also show that mathematical modelling offers a relatively inexpensive and reliable means to quantify inter-experimental variability and assess the reproducibility of results.