Accumulation of amino acid substitutions promotes irreversible structural changes in the hemagglutinin of human influenza AH3 virus during evolution

In order to clarify the effect of an accumulation of amino acid substitutions on the hemadsorption character of the influenza AH3 virus hemagglutinin (HA) protein, we introduced single-point amino acid changes into the HA1 domain of the HA proteins of influenza viruses isolated in 1968 (A/Aichi/2/68) and 1997 (A/Sydney/5/97) by using PCR-based random mutation or site-directed mutagenesis. These substitutions were classified as positive or negative according to their effects on the hemadsorption activity. The rate of positive substitutions was about 50% for both strains. Of 44 amino acid changes that were identical in the two strains with regard to both the substituted amino acids and their positions in the HA1 domain, 22% of the changes that were positive in A/Aichi/2/68 were negative in A/Sydney/5/97 and 27% of the changes that were negative in A/Aichi/2/68 were positive in A/Sydney/5/97. A similar discordance rate was also seen for the antigenic sites. These results suggest that the accumulation of amino acid substitutions in the HA protein during evolution promoted irreversible structural changes and therefore that antigenic changes in the H3HA protein may not be limited.     

Accumulation pattern of amino acid substitutions in protein evolution

Summary A simple method for the evolutionary analysis of amino acid sequence data is presented and used to examine whether the number of variable sites (NVS) of a protein is constant during its evolution. The NVSs for hemoglobin and for mitochondrial cytochrome c are each found to be almost constant, and the ratio between the NVSs is close to the ratio between the unit evolutionary periods. This indicates that the substitution rate per variable site is almost uniform for these proteins, as the neutral theory claims. An advantage of the present analysis is that it can be done without knowledge of paleontological divergence times and can be extended to bacterial proteins such as bacterial c-type cytochromes. It is suggested that the NVS of cytochrome c has been almost constant even over the long period (ca. 3.0 billion years) of bacterial evolution but that at least two different substitution rates are necessary to describe the accumulated changes in the sequence. This two clock interpretation is consistent with fossil evidence for the appearance times of photosynthetic bacteria and eukaryotes.     

A single amino acid substitution in 1918 influenza virus hemagglutinin changes receptor binding specificity

The receptor binding specificity of influenza viruses may be important for host restriction of human and avian viruses. Here, we show that the hemagglutinin (HA) of the virus that caused the 1918 influenza pandemic has strain-specific differences in its receptor binding specificity. The A/South Carolina/1/18 HA preferentially binds the alpha2,6 sialic acid (human) cellular receptor, whereas the A/New York/1/18 HA, which differs by only one amino acid, binds both the alpha2,6 and the alpha2,3 sialic acid (avian) cellular receptors. Compared to the conserved consensus sequence in the receptor binding site of avian HAs, only a single amino acid at position 190 was changed in the A/New York/1/18 HA. Mutation of this single amino acid back to the avian consensus resulted in a preference for the avian receptor.     

Changes in the amino acid sequence of hemagglutinin during sequential adaptation of human influenza virus A to replication in mice

The primary structure of hemagglutinin (HA) gene of Influenza virus A/USSR/90/77 (H1N1) variants after 3 and 11 passages has been determined. In the HA1 coding region of mice-adapted virus (11 passages) there are two amino acid substitutions: Thr 89----Ala and Asn 127----Asp. At the first stage of adaptation (3-rd passage) only a single mutation was detected: Asn 127----Asp. The adaptation is accompanied by the loss of specific carbohydrate attachment sites adjacent to the receptor-binding site located at HA1 subunit with a concomitant variation in antigenicity.     

Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus

It has been proposed that antigenic evolution of hemagglutinin 1 (HA1) for H3N2 human influenza A virus was punctuated. In the population genetic analysis, however, it was controversial whether positive selection operated on HA1 in a punctuated manner for the branches of the phylogenetic tree where transitions to new antigenic clusters occurred (C branches), or continuously. In the molecular evolutionary analysis, positive selection was detected for the trunk (T) branches but the relationship between antigenic evolution and positive selection was unclear. Here molecular evolutionary analysis was conducted to examine natural selection operating on HA1 of H3N2 human influenza A virus by dividing HA1 into epitopes A-E and other sites, as well as dividing the phylogenetic tree into the C branches overlapping with the T branches (C-T branches), those not overlapping with the T branches (C-NT branches), the T branches not overlapping with the C branches (NC-T branches), and other branches (NC-NT branches). Positive selection was detected for C, T, and NC-T branches, whereas evolution for the NC-NT branches appeared to be mainly neutral. Positive selection appeared to have operated throughout the trunk, which covered the entire time period of the phylogenetic tree, suggesting that positive selection operated continuously on HA1 during evolution of H3N2 human influenza A virus.     

Nature of amino acid substitutions in homologous proteins during evolution

Replacement substitutions of mitochondrial cytochrome $\text{c}$ and α- and β-chains of haemoglobin have been studied by considering the structural similarity among amino acid residues at the secondary and tertiary structural levels. Secondary structural similarity explains ～70% while tertiary structural similarity explains ～50% of observed replacements for most of the cases. These structural similarities could not account for all the replacement substitutions. The study was extended to consider the composition of codons, and the chemical nature and polarity of the replacing and replaced residues. These also could not individually account for all the affected replacements. In general, no property of amino acid residues is conserved for substitutions occurring at any single position during evolution of proteins.     

Evidence for N-glycan shielding of antigenic sites during evolution of human influenza A virus hemagglutinin

After the emergence of influenza A viruses in the human population, the number of N-glycosylation sites (NGS) in the globular head region of hemagglutinin (HA) has increased continuously for several decades. It has been speculated that the addition of NGS to the globular head region of HA has conferred selective advantages to the virus by preventing the binding of antibodies (Ab) to antigenic sites (AS). Here, the effect of N-glycosylation on the binding of Ab to AS in human influenza A virus subtype H3N2 (A/H3N2) was examined by inferring natural selection at AS and other sites (NAS) that are located close to and distantly from the NGS in the three-dimensional structure of HA through a comparison of the rates of synonymous (d(S)) and nonsynonymous (d(N)) substitutions. When positions 63, 122, 126, 133, 144, and 246 in the globular head region of HA were non-NGS, the d(N)/d(S) was >1 and positive selection was detected at the AS located near these positions. However, the d(N)/d(S) value decreased and the evidence of positive selection disappeared when these positions became NGS. In contrast, d(N)/d(S) at the AS distantly located from the positions mentioned above and at the NAS of any location were generally <1 and did not decrease when these positions changed from non-NGS to NGS. These results suggest that the attachment of N-glycans to the NGS in the globular head region of HA prevented the binding of Ab to AS in the evolutionary history of human A/H3N2 virus.     

Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.     

Compensatory evolution of net-charge in influenza A virus hemagglutinin

The propagation of influenza A virus depends on the balance between the activities of hemagglutinin (HA) for binding to host cells and neuraminidase (NA) for releasing from infected cells (HA-NA balance). Since the host cell membrane and the sialic acid receptor are negatively charged, the amino acid substitutions increasing (charge+) and decreasing (charge-) the positive charge of HA subunit 1 (HA1) enhance and reduce, respectively, the binding avidity and affinity. The positive charge of HA1 in human influenza A virus bearing subtype H3N2 (A/H3N2 virus) was observed to have increased during evolution, but the evolutionary mechanism for this observation was unclear because this may disrupt the HA-NA balance. Here we show, from the phylogenetic analysis of HA for human A/H3N2 and A/H1N1 viruses, that the relative frequencies of charge+ and charge- substitutions were elevated on the branches where the number of N-glycosylation sites (NGS) increased and decreased, respectively, compared to those where the number of NGS did not change. On the latter branches, the net-charge of HA1 appeared to have been largely maintained to preserve its structure and function. Since the charge+ and charge- substitutions in HA1 have opposite effects to the gain and loss of NGS on the binding and release of the virus, the net-charge of HA1 may have evolved to compensate for the effect of the gain and loss of NGS, probably through changing the avidity. Apparently, the relative frequency of charge- substitutions in HA1 of A/H3N2 virus was elevated after the introduction of oseltamivir, and that of charge+ substitutions in HA1 of A/H1N1 virus was elevated after the spread of oseltamivir resistance. These observations may also be explained by the compensatory effect of the net-charge in HA1 on the NA activity for keeping the HA-NA balance.     

Single amino acid substitutions in the hemagglutinin can alter the host range and receptor binding properties of H1 strains of influenza A virus.

We have previously characterized an influenza A (H1N1) virus which has host-dependent growth and receptor binding properties and have shown that a mutation which removes an oligosaccharide from the tip of the hemagglutinin (HA) by changing Asn-129 to Asp permits this virus to grow to high titer in MDBK cells, (C. M. Deom, A. J. Caton, and I. T. Schulze, Proc. Natl. Acad. Sci. USA 83:3771-3775, 1986). We have now isolated monoclonal antibodies specific for the mutant HA and have used escape mutants to identify alterations in HA sequence which reduce virus yields from MDBK cells without reducing those from chicken embryo fibroblasts. Two types of escape mutants which grow equally well in chicken embryo fibroblasts were obtained. Those with the parent phenotype contain Asn at residue 129 and are glycosylated at that site. Those with the mutant phenotype are unchanged at residue 129 but have a Gly to Glu substitution at residue 158, which is close to residue 129 on the HA1 subunit. Binding assays with neoglycoproteins containing N-acetylneuraminic acid in either alpha 2,3 or alpha 2,6 linkage to galactose showed that the MDBK-synthesized oligosaccharides at Asn-129 reduce binding to both of these receptors, leaving the HA's preference for alpha 2,6 linkages unchanged. Glu at residue 158 greatly reduces binding to both receptors without reducing virus yields from MDBK cells. We conclude that changes in the receptor binding properties of the HA can result either from direct alteration of the HA protein by host cell glycosylation or from mutations in the HA gene and that these changes generate heterogeneity that can contribute to the survival of influenza A virus populations in nature.     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

Volume and polarity changes accompanied by amino acid substitutions in protein evolution.

We evaluated the volume and polarity changes accompanied by amino acid substitutions along branches of the phylogenetic trees of cytochrome c, myoglobin and hemoglobin alpha and beta chains. In most cases the volume changes accompanied by the substitutions were found to be much larger than the volume of cavities existing in the interior of X-ray-analysed proteins. This implies that the interior of the proteins is very flexible and the necessary space for a larger amino acid residue substitution can be provided by adjusting nearby structures. Also, the volume and polarity changes are not particularly dependent on whether the substituted site is located in the exterior or interior of the proteins. This result supports the concept of the covarions by Fitch and Markowitz, when combined with the known fact that the exterior sites are more variable than the interior ones during protein evolution.     

Reciprocal amino acid substitutions in the evolution of homologous peptides

Amino acid sequences of peptides are often inferred from their amino acid compositions by comparison with homologous peptides of known sequence. The probabilities are considered that by such an approach errors are made due to the occurrence of balanced double changes, i.e. reciprocal substitutions, between two homologous peptides of identical compositions. Formulae are derived for the calculation of these probabilities, depending on peptide length and evolutionary distance. However, such calculations requiring too much computer time, the probabilities for reciprocal substitutions are estimated by simulation of evolutionary changes in peptides. It can be concluded from the resulting data that for many purposes the possible errors in amino acid sequences partially inferred from amino acid compositions are acceptably small.     

Two types of amino acid substitutions in protein evolution

Summary The frequency of amino acid substitutions, relative to the frequency expected by chance, decreases linearly with the increase in physico-chemical differences between amino acid pairs involved in a substitution. This correlation does not apply to abnormal human hemoglobins. Since abnormal hemoglobins mostly reflect the process of mutation rather than selection, the correlation manifest during protein evolution between substitution frequency and physico-chemical difference in amino acids can be attributed to natural selection. Outside of abnormal proteins, the correlation also does not apply to certain regions of proteins characterized by rapid rates of substitution. In these cases again, except for the largest physico-chemical differences between amino acid pairs, the substitution frequencies seem to be independent of the physico-chemical parameters. The limination of the substituents involving the largest physicochemical differences can once more be attributed to natural selection. For smaller physico-chemical differences, natural selection, if it is operating in the polypeptide regions, must be based on parameters other than those examined.     

A single amino acid change in the cytoplasmic domain alters the polarized delivery of influenza virus hemagglutinin

In the polarized kidney cell line MDCK, the influenza virus hemagglutinin (HA) has been well characterized as a model for apically sorted membrane glycoproteins. Previous work from our laboratory has shown that a single amino acid change in the cytoplasmic sequence of HA converts it from a protein that is excluded from coated pits to one that is efficiently internalized. Using trypsin or antibodies to mark protein on the surface, we have shown in MDCK cells that HA containing this mutation is no longer transported to the apical surface but instead is delivered directly to the basolateral plasma membrane. We propose that a cytoplasmic feature similar to an endocytosis signal can cause exclusive basolateral delivery.     

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.     

Evolution of the hemagglutinin gene of human influenza A virus H3 subtype

An evolutional tree of human influenza viruses of the H3N2-subtype is suggested on the basis of combined published primary structures of the hemagglutinin HA1-subunit. Possible differences between natural and sequenced structures are discussed. A tendency to reversions in the course of antigenic draft within the subtype has been revealed to support the hypothesis of limited antigenic evolution within a single subtype.     

Effects of single-point amino acid substitutions on the structure and function neuraminidase proteins in influenza A virus

In order to clarify the effect of amino acid substitutions on the structure and function of the neuraminidase (NA) protein of influenza A virus, we introduced single-point amino acid substitutions into the NA protein of the A/Tokyo/3/67 (H2N2) strain using PCR-based random mutation. The rate of tolerant random one amino acid substitutions in the NA protein was 47%. Rates of tolerant substitutions for the stalk and for the surface and inner portion of the head region of the NA protein were 79, 54, and 19%, respectively. Deleterious changes, such as those causing the NA protein to stop at the Golgi/endoplasmic reticulum, were scattered throughout the protein. On the other hand, the ratio of mutations with which the NA protein lost neuraminidase activity, but was transported to the cell surface, decreased in proportion to the distance from the structural center of enzyme active site. In order to investigate the effect of accumulated amino acid substitutions on the structural character of the N2NA protein during evolution, the same amino acid substitutions were introduced by site-directed mutagenesis at 23 homologous positions on N2 proteins of A/Tokyo/3/67, A/Bangkok/15/85 (H3N2), and A/Mie/1/2004 (H3N2). The results showed a shift, or discordance, in tolerance at some of the positions. An increase in discordance was correlated with the interval in years between virus strains, and the discordance rate was estimated to be 0.6-0.7% per year.     

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 anantigenic surface glycoprotein, which is essential to the influenza infection process, and is the target of anti-flu drugs oseltamivirand zanamivir. Currently, NA inhibitors are the pillar pharmacological strategy against seasonal and global influenza. Althoughmutations observed after NA-inhibitor treatment are characterized by changes in conserved amino acids of the enzyme catalyticsite, it is possible that specific amino acid substitutions (AASs) distant from the active site such as H274Y, could confer oseltamiviror zanamivir resistance. To better understand the molecular distribution pattern of NA AASs, we analyzed NA AASs from allavailable reported pandemic AH1N1 NA sequences, including those reported from America, Africa, Asia, Europe, Oceania, andspecifically from Mexico. The molecular distributions of the AASs were obtained at the secondary structure domain level for boththe 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 highlyconserved amino acids from the NA active site and catalytic site are indeed being affected by mutations. The reported NA AASsfollow a similar molecular distribution pattern worldwide. Although most AASs are distributed distantly from the active site, thisstudy shows the emergence of mutations affecting the previously conserved active and catalytic site. A significant number ofunique AASs were reported simultaneously on different continents.     

Antigenic drift in influenza virus H3 hemagglutinin from 1968 to 1980: multiple evolutionary pathways and sequential amino acid changes at key antigenic sites

Surveys of the antigenic properties of a wide range of variants of the H3N2 (Hong Kong) influenza virus subtype have revealed complex patterns of variants cocirculating during each of the main epidemic eras of the subtype. We determined hemagglutinin (HA) gene sequences for 14 isolates chosen to give the wildest possible spread of variant types. The addition of these data to existing HA gene sequence information for other variants provides a comprehensive picture of HA gene evolution during antigenic drift among H3N2 subtype viruses. The data reveal the existence of multiple evolutionary pathways during at least one period of development of the subtype and strikingly demonstrate that amino acid changes are limited to a small number of locations on the HA molecule during antigenic drift. The occurrence of sequential amino acid changes at key positions within these variable regions suggests that the HA structure has remained constant during subtype evolution so that only limited possibilities remain for further antigenic drift among H3N2 viruses.