Previous Article | Next Article ![]()
Journal of Virology, September 2007, p. 9950-9955, Vol. 81, No. 18
0022-538X/07/$08.00+0 doi:10.1128/JVI.00468-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Faculty of Medicine, Siriraj Hospital,1 Faculty of Science,2 Faculty of Tropical Medicine,3 Mahidol University, Bangkok, Thailand; Bureau of Epidemiology,4 Department of Disease Control, Bangkok, Thailand,5 Faculty of Science, Kasetsart University, Bangkok, Thailand,6 College of Life and Health Sciences, Chubu University, Kasugai, Japan7
Received 5 March 2007/ Accepted 18 June 2007
|
|
|---|
2,3-galactose (SA
2,3Gal), whereas human influenza viruses preferentially recognize SA
2,6Gal. A conversion to SA
2,6Gal specificity is believed to be one of the changes required for the introduction of new hemagglutinin (HA) subtypes to the human population, which can lead to pandemics. Avian influenza H5N1 virus is a major threat for the emergence of a pandemic virus. As of 12 June 2007, the virus has been reported in 45 countries, and 312 human cases with 190 deaths have been confirmed. We describe here substitutions at position 129 and 134 identified in a virus isolated from a fatal human case that could change the receptor-binding preference of HA of H5N1 virus from SA
2,3Gal to both SA
2,3Gal and SA
2,6Gal. Molecular modeling demonstrated that the mutation may stabilize SA
2,6Gal in its optimal cis conformation in the binding pocket. The mutation was found in approximately half of the viral sequences directly amplified from a respiratory specimen of the patient. Our data confirm the presence of H5N1 virus with the ability to bind to a human-type receptor in this patient and suggest the selection and expansion of the mutant with human-type receptor specificity in the human host environment. |
|
|---|
2,6-galactose (SA
2,6Gal) receptor with high affinity (17). It is believed that this receptor binding property is the major factor preventing the H5N1 virus from efficiently transmitting from person to person and causing a pandemic (17). The receptor binding preference of H5N1 viruses can be altered by only a few amino acid substitutions in the hemagglutinin (HA) protein. Mutations that change the receptor binding preference from the avian to the human type could potentially enable the virus to transmit efficiently in the human population and cause a catastrophic pandemic. Monitoring of the viral changes is therefore extremely important in the current situation where H5N1 viruses are spreading progressively.
A previous study showed that mutations at positions 226 and 228 (H3 numbering) (Q226L and G228S), which are the adaptive mutations for H2 and H3 (5, 19), could reduce the binding affinity to sialic acid-
2,3-galactose (SA
2,3Gal) of an H5N1 virus isolated in 1997 (5, 19). Human H5N1 isolates from Hong Kong that were isolated in 2003, which contain a mutation at position 227 (S227N) (H3 numbering), were shown to have a reduced binding affinity toward SA
2,3Gal and an increased affinity toward SA
2,6Gal (1). Another report demonstrated previously that neither these mutations nor the mutations that could adapt H1 viruses to the human receptor (E190D and G225D) (H3 numbering) could completely convert a Vietnam H5N1 virus isolated in 2004 to the 2,6-type receptor specificity (16). Except for S227N, these mutations have not been found in H5N1 viruses isolated from humans or animals. Recently, N182K and Q192R mutations were shown to enhance the binding of a Vietnam H5N1 HA isolated in 2004 to the SA
2,6Gal receptor (20). Although it was described previously that the N182K mutation was found in Kan-1 HA (20), the original Kan-1 virus as well as the HA sequence of Kan-1 in the GenBank Database do not contain this mutation. It is not clear from where this mutation was derived. The Q192R mutation was found in a clone that was present as a minor population in a viral isolate and was identified after plaque purification. In contrast, here, we show a naturally occurring mutant with an ability to bind SA
2,6Gal that was directly identified in a human nasopharyngeal specimen.
|
|
|---|
Hemagglutination assay.
A 10% suspension of goose (Anser cygnoides) red blood cells was prepared in phosphate-buffered saline (PBS), and a 50-µl aliquot was treated with 1.25 U of
2,3-sialidase cloned from Salmonella enterica serovar Typhimurium LT2 (7) (Takara, Japan) for 1 h at 37°C. Untreated and treated red blood cells were used in the hemagglutination assay.
HA-receptor binding assay.
The receptor binding preference was analyzed by a solid-phase direct binding assay as previously described (20) by using the sialylglycopolymer, which contains N-acetylneuraminic acid linked to galactose through either an
2,3 or an
2,6 bond (Neu5Ac
2, 3LacNAcb-pAP, and Neu5Ac
2,6LacNAcb-pAP) (18). Serial dilutions of each sialylglycopolymer were prepared in PBS, and 100 µl was added to the wells of 96-well microtiter plates (Polystyrene Universal-Bind Microplate; Corning) and allowed to attach overnight at 4°C. The plate containing sialylglycopolymer was irradiated under UV light at 254 nm for 10 min and then washed five times with 200 µl PBS. The plate was blocked for 8 h at 4°C with PBS containing 2% skim milk powder. After washing five times with 200 µl PBS containing 0.1% Tween 20, virus culture supernatant containing 128 HA units was allowed to attach onto the plate on ice overnight. After the incubation, the plate was washed five times with ice-cold PBS-0.1% Tween 20, and 50 µl of anti-HA goat hyperimmune serum at a dilution of 1:2,000 was added to each well and was allowed to incubate on ice for 2 h. The plate was then washed again with ice-cold PBS-0.1% Tween 20, and 50 µl of the polyclonal rabbit anti-goat immunoglobulin-horseradish peroxidase conjugate (DakoCytomation) at a dilution of 1:2,000 was added to each well. After incubation on ice for 2 h, the wells were extensively washed with PBS-0.1% Tween 20, and 100 µl per well of premixed tetramethylbenzidine-H2O2 substrate was added. After incubation at room temperature for 15 min, the reaction was stopped with 50 µl of 1 M H2SO4, and the absorbance at 450/630 nm was read.
Molecular dynamic simulation.
The crystal structures of HA from A/Duck/Singapore/3/97 (H5N1) (4) were used as templates for the simulations of HA binding to SA
2,3Gal (Protein Data Bank [PDB] accession number 1JSN) and SA
2,6Gal (PDB accession number 1JSO). In the structural template under PDB accession number 1JSO, the sialic residue has no galactose unit connected; therefore, we added it with the torsion angle of 55°. Both glycosides in the two structures were terminated with a methoxy group and used as the input for molecular dynamic simulations. In homology modeling, the wild-type and mutant HA (L129V/A134V) were three-dimensionally aligned via the SWISS-MODEL server (12) using PDB accession numbers 1JSN and 1JSO and used as the initial input structures for molecular dynamic simulations. In order to provide a control of classical SA
2,6Gal-tropic HA, we ran a human influenza virus H1 HA structure (PDB accession number 1RVZ) (2) in a similar simulation.
All structures were solvated using the TIP5P water model (10), and we performed energy minimization to relieve bad contacts caused by unreasonable distances in the initial structures and then equilibrated for 100 ps before a 3-ns productive run at 300K using the SANDER module in the AMBER9 program (University of California, San Francisco) with a Glycam04 parameter (http://glycam.ccrc.uga.edu). Xmgrace (http://plasma-gate.weizmann.ac.il/Grace/), VMD (8), and AMBER tools running on UNIX were exploited to visualize and manipulate all figures.
|
|
|---|
2,3Gal-specific sialidase (7). Theoretically, the sialidase digestion should abolish hemagglutination by SA
2,3Gal-specific viruses, whereas viruses that can bind to SA
2,6Gal should maintain hemagglutination activity with the treated red blood cells. The sialidase treatment did not affect the hemagglutination titer of human influenza viruses. An example of H1N1 virus (A/Thailand/Siriraj-12/06), the hemagglutination titer of which was not affected by the treatment, is shown in Fig. 1a. We have tested four human H5N1 isolates [A/Thailand/1(KAN-1)/2004, A/Thailand/3(SP-83)/2004, A/Thailand/5(KK-494)/2004, and A/Thailand/676/2005]. The sialidase treatment completely abolished the hemagglutination activity of all the viruses with more than a 256-fold reduction in the hemagglutination titer except for A/Thailand/676/2005 (Th676), which partially maintained its hemagglutination activity in sialidase-treated red blood cells with some reduction in the hemagglutination titer (Fig. 1a). This suggested an enhanced binding of this virus to SA
2,6Gal. The virus was isolated from a 5-year-old boy in Thailand. The patient had a progressive viral pneumonia that led to respiratory failure and death by 12 days after onset of illness.
![]() View larger version (31K): [in a new window] |
FIG. 1. (a) Abolition of hemagglutination titers in SA 2,3Gal-specific sialidase-treated red blood cells compared to those in untreated cells indicated SA 2,3Gal monospecificity for H5N1 viruses except for Th676 and the L129V/A134V reverse genetic virus, the titers of which were only modestly affected, suggesting an enhanced binding to SA 2,6Gal. (b) Electropherogram of sequencing reaction of the Th676 virus shows a T-to-G mutation changing the amino acid from leucine to valine at position 129 and a C-to-T mutation changing the amino acid from alanine to valine at position 134 (arrowheads).
|
2,6Gal (20), our Kan-1 HA clone (GenBank accession number EF107522) as well as the Kan-1 viral isolate (GenBank accession number AY555150) do not contain the N182K mutation, which was described to enhance binding to SA
2,6Gal (20). The hemagglutination pattern of the L129V/A134V reverse genetic virus indicated an enhanced SA
2,6Gal binding, whereas that of the L129V virus did not (Fig. 1a). The lesser degree of reduction of the hemagglutination titer of the L129V/A134V reverse genetic virus (fourfold) compared to that of Th676 (eightfold) might reflect the fact that Th676 was a mixture of wild-type and mutant viruses. We also used a recently described direct binding assay using sialylglycopolymers (3) and confirmed the effect of the L129V/A134V mutation on the receptor binding preference. The reverse genetic virus carrying wild-type Kan-1 HA showed SA
2,3Gal specificity and did not bind significantly to SA
2,6Gal, whereas a human H1N1 virus, A/Thailand/Siriraj-12/06, bound only to SA
2,6Gal (Fig. 2a and b). Also, while the L129V mutation alone did not change the receptor binding preference, the HA containing both L129V and A134V mutations bound equally well to both SA
2,3Gal and SA
2,6Gal (Fig. 2c and d).
![]() View larger version (26K): [in a new window] |
FIG. 2. Direct binding assay using sialylglycopolymers. Human H1N1 virus (A/Thailand/Siriraj-12/06) showed receptor preference for SA 2,6Gal (a), whereas the reverse genetic virus with Kan-1 wild-type HA showed receptor preference for SA 2,3Gal (b). The reverse genetic virus carrying HA with the L129V mutation showed a receptor preference similar to that of the wild-type virus (c), whereas the virus carrying HA with L129V/A134V mutations showed receptor preference for both SA 2,3Gal and SA 2,6Gal (d). Diamonds represent absorbencies of SA 2,3Gal binding, and squares represent absorbencies of SA 2,6Gal binding.
|
2,3Gal and SA
2,6Gal receptors were shown to have a specific conformation, either cis or trans (4). Torsion angle (
) is the angle between two planes containing O6, C1 of the SA unit, and O3 (or O6), C3 of Gal unit (Fig. 3a). The angle indicates whether the glycoside is in the cis (
= 56°) or trans (
= –55°) conformation (Fig. 3a). Within its bound state to H5 in an X-ray cocrystal structure (4), SA
2,3Gal was found in the trans conformation. The trans conformation of SA
2,3Gal allows 4-OH (O4) and the glycosidic oxygen (O3) of the Gal unit to interact optimally with Gln222 of H5 (4). In contrast, if bound in the trans conformation, the 4-OH in SA
2,6Gal would be too far away from Gln222, and if bound in the cis conformation, the glycosidic oxygen (O6) in SA
2,6Gal would be too far away from Gln222. The lack of interactions between 4-OH and Gln222 in the trans conformation and between O6 and Gln222 in the cis conformation makes the binding of H5 to SA
2,6Gal unstable in both conformations. The conformation of SA
2,6Gal in solution was shown to have a cis-to-trans ratio of 9:1 (11), which suggests that the optimal binding conformation for SA
2,6Gal is cis because no additional energy would be required for the conformational change. Furthermore, recent molecular modeling has predicted SA
2,6Gal in H5 binding to be in the cis conformation (9). In our simulation, human H1 (PDB accession number 1RVZ) showed an average SA
2,6Gal
angle of 63°, which indicated a cis conformation (Fig. 3b). This provides a positive control for SA
2,6Gal binding in our simulation and indicates that the
angle can be used as an indicator of the receptor preference. In order to have structural insight into the altered receptor specificity of the mutant HA, we performed molecular dynamic simulations using two sialic acid-H5 cocrystals as reference structures for the two types of glycosidic linkages, SA
2,3Gal (PDB accession number 1JSN) and SA
2,6Gal (PDB accession number 1JSO) (4). In our modeling, both the wild-type and the mutant L129V/A134V HAs shared similar binding patterns within the sialic binding pocket as previously reported for both SA
2,3Gal and SA
2,6Gal binding (14, 19). However, it was observed that over the period of simulations, hydrophobicity and spatial constraint changes at residues 129 and 134 due to the different alkyl side chains, especially the A134V mutation, which happened near the glycosidic linkage, caused some changes in glycoside binding patterns. The
angle indicated that SA
2,6Gal in the L129V/A134V HA spent most of the time in the low-energy cis conformation, whereas those in the wild-type Kan-1 HA and the template H5 were forced to change from their cis conformation in the initial input structures to the trans conformation (Fig. 3c). Although SA
2,6Gal in both wild-type Kan-1 HA and the template HA was in the trans conformation at the end of the simulation period, SA
2,6Gal in Kan-1 HA stayed longer in the cis conformation than that in the template HA. This suggested that wild-type Kan-1 HA might have a slightly increased affinity for SA
2,6Gal compared to the template HA of A/Duck/Singapore/3/97. The increased stability in the cis conformation of the binding to L129V/A134V HA was due to an alternative interaction between Gly221 and 4-OH or 3-OH (O3) of Gal caused by a displacement of Gly221 by the A134V mutation (Fig. 3d). In the simulations of SA
2,3Gal binding to L129V/A134V HA, the
angle significantly moved away from the preferred angle (from –55° to –30°) (Fig. 3c). Although SA
2,3Gal was still in the trans conformation, the A134V mutation pushed Gln222 slightly away from the optimal binding condition. The widening of the gap between Gln222 and the glycosidic oxygen (O3) caused by larger hydrophobic side chain of Val134 might reduce the binding affinity of the mutant HA for SA
2,3Gal.
![]() View larger version (44K): [in a new window] |
FIG. 3. The two glycoside backbone conformations, SA 2,6Gal and SA 2,3Gal, cut from the X-ray cocrystals under PDB accession numbers 1JSI and 1JSIN (4) with different values for torsion ( ) angles that were used as the cis or trans conformation characteristics (a). For a human H1 virus (PDB accession number 1RVZ), SA 2,6Gal showed an average angle of 63°, indicating a cis conformation in the molecular dynamic simulation (b). The amount of time that each receptor analog spent having a particular angle for SA 2,3Gal (top) and SA 2,6Gal (bottom) in the binding sites of the template H5 (A/Duck/Singapore/3/97) (black), wild-type Kan-1 (green), and the L129V/A134V mutant (red) as the results of molecular dynamics simulations indicated an altered binding conformation in the L129V/A134V mutant (c). The modeling showed that SA 2,6Gal in a cis conformation in the wild-type Kan-1 HA pockets had long distances between O6 and Gln222 (3.47 Å) and between 4-OH and Gly221 (7.72 Å) (left). However, in the L129V/A134V mutant, the glycoside was stabilized with the Gly221 and Gal interactions as observed by the shorter distance between 4-OH and Gln221 (2.02 Å) (right) (d).
|
|
|
|---|
2,6Gal receptor is expressed mainly in the upper airway, while the SA
2,3Gal receptor is expressed in alveoli and the terminal bronchiole (13). A virus with good affinity to both SA
2,3Gal and SA
2,6Gal receptors may be a very dangerous one, which could both infect efficiently via its binding to SA
2,6Gal in the upper airway and cause severe infection in the lung via its binding to SA
2,3Gal. This hypothesis is supported by the fact that one of the two well-characterized HA genes from the H1N1 1918 pandemic virus binds efficiently to both SA
2,3Gal and SA
2,6Gal (3, 15). Although receptor binding preference is a major factor determining host species tropism, we do not know whether the alteration in receptor binding properties is sufficient to enable the virus to transmit from person to person efficiently. Nevertheless, our finding indicated that an adaptation of H5N1 virus to the human host by receptor binding site modification could and did indeed happen. Since the identification of this patient, there have been three confirmed human cases in the country, but none of the viral isolates from these cases contained the L129V and A134V substitutions. It is likely that this particular mutant has been eliminated by infection control measures.
Our report demonstrates that avian influenza H5N1 virus could be naturally adapted to the human-type receptor. We need to intensify our effort to detect such viruses as early as possible. Our data also provide a genetic marker that can be used to screen for an H5N1 virus with pandemic potential. H5N1 viruses may have several ways to adapt their receptor binding properties, but mutations that were found directly in patients, such as L129V/A134V, are more likely than those artificially generated and tested mutations to be the mutation that could cause a pandemic.
We thank Atthapon Eiamudomkan, Prasongsak Nakhornkwang, Manee Phonpasee, Sanya Kittisoontaropas, Vivat Pvongpasert, Wannasatre Rattanalum, Chansuda Sukbumrung, and Sompong Boonsuepchat of the NakhonNayok Provincial Public Health Office. The reverse genetic plasmid pHw2000 with A/PR/8/34 genomic segments was kindly provided by R. G. Webster and E. Hoffmann. The activity was a part of the newly established Thailand Avian Influenza Monitoring Network (TAIM Net).
Published ahead of print on 11 July 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»