Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis
Identifieur interne : 000B42 ( Pmc/Corpus ); précédent : 000B41; suivant : 000B43Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis
Auteurs : Wentian Chen ; Shisheng Sun ; Zheng LiSource :
- PLoS ONE [ 1932-6203 ] ; 2012.
Abstract
Increasing numbers of H5N1 influenza viruses (IVs) are responsible for human deaths, especially in North Africa and Southeast Asian. The binding of hemagglutinin (HA) on the viral surface to host sialic acid (SA) receptors is a requisite step in the infection process. Phylogenetic analysis reveals that H5N1 viruses can be divided into 10 clades based on their HA sequences, with most human IVs centered from clade 1 and clade 2.1 to clade 2.3. Protein sequence alignment in various clades indicates the high conservation in the receptor-binding domains (RBDs) is essential for binding with the SA receptor. Two glycosylation sites, 158N and 169N, also participate in receptor recognition. In the present work, we attempted to construct a serial H5N1 HA models including diverse glycosylated HAs to simulate the binding process with various SA receptors
Url:
DOI: 10.1371/journal.pone.0038794
PubMed: 22719948
PubMed Central: 3375263
Links to Exploration step
PMC:3375263Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis</title><author><name sortKey="Chen, Wentian" sort="Chen, Wentian" uniqKey="Chen W" first="Wentian" last="Chen">Wentian Chen</name><affiliation><nlm:aff id="aff1"><addr-line>Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an, People’s Republic of China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Sun, Shisheng" sort="Sun, Shisheng" uniqKey="Sun S" first="Shisheng" last="Sun">Shisheng Sun</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Pathology, Clinical Chemistry Division, Johns Hopkins University, Baltimore, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Li, Zheng" sort="Li, Zheng" uniqKey="Li Z" first="Zheng" last="Li">Zheng Li</name><affiliation><nlm:aff id="aff1"><addr-line>Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an, People’s Republic of China</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">22719948</idno><idno type="pmc">3375263</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3375263</idno><idno type="RBID">PMC:3375263</idno><idno type="doi">10.1371/journal.pone.0038794</idno><date when="2012">2012</date><idno type="wicri:Area/Pmc/Corpus">000B42</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000B42</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis</title><author><name sortKey="Chen, Wentian" sort="Chen, Wentian" uniqKey="Chen W" first="Wentian" last="Chen">Wentian Chen</name><affiliation><nlm:aff id="aff1"><addr-line>Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an, People’s Republic of China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Sun, Shisheng" sort="Sun, Shisheng" uniqKey="Sun S" first="Shisheng" last="Sun">Shisheng Sun</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Pathology, Clinical Chemistry Division, Johns Hopkins University, Baltimore, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Li, Zheng" sort="Li, Zheng" uniqKey="Li Z" first="Zheng" last="Li">Zheng Li</name><affiliation><nlm:aff id="aff1"><addr-line>Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an, People’s Republic of China</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2012">2012</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Increasing numbers of H5N1 influenza viruses (IVs) are responsible for human deaths, especially in North Africa and Southeast Asian. The binding of hemagglutinin (HA) on the viral surface to host sialic acid (SA) receptors is a requisite step in the infection process. Phylogenetic analysis reveals that H5N1 viruses can be divided into 10 clades based on their HA sequences, with most human IVs centered from clade 1 and clade 2.1 to clade 2.3. Protein sequence alignment in various clades indicates the high conservation in the receptor-binding domains (RBDs) is essential for binding with the SA receptor. Two glycosylation sites, 158N and 169N, also participate in receptor recognition. In the present work, we attempted to construct a serial H5N1 HA models including diverse glycosylated HAs to simulate the binding process with various SA receptors <italic>in silico</italic>. As the SA-α-2,3-Gal and SA-α-2,6-Gal receptor adopted two distinctive topologies, straight and fishhook-like, respectively, the presence of N-glycans at 158N would decrease the affinity of HA for all of the receptors, particularly SA-α-2,6-Gal analogs. The steric clashes of the huge glycans shown at another glycosylation site, 169N, located on an adjacent HA monomer, would be more effective in preventing the binding of SA-α-2,3-Gal analogs.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Zachary, S" uniqKey="Zachary S">S Zachary</name></author><author><name sortKey="Rahul, R" uniqKey="Rahul R">R Rahul</name></author><author><name sortKey="Karthik, V" uniqKey="Karthik V">V Karthik</name></author><author><name sortKey="Ram, S" uniqKey="Ram S">S Ram</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author><author><name sortKey="Wiley, Dc" uniqKey="Wiley D">DC Wiley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wiley, Dc" uniqKey="Wiley D">DC Wiley</name></author><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Senne, Da" uniqKey="Senne D">DA Senne</name></author><author><name sortKey="Panigrahy, B" uniqKey="Panigrahy B">B Panigrahy</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name></author><author><name sortKey="Pearson, Je" uniqKey="Pearson J">JE Pearson</name></author><author><name sortKey="Suss, J" uniqKey="Suss J">J Suss</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sammalkorpi, M" uniqKey="Sammalkorpi M">M Sammalkorpi</name></author><author><name sortKey="Lazaridis, T" uniqKey="Lazaridis T">T Lazaridis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gambaryan, As" uniqKey="Gambaryan A">AS Gambaryan</name></author><author><name sortKey="Tuzikov, Ab" uniqKey="Tuzikov A">AB Tuzikov</name></author><author><name sortKey="Pazynina, Gv" uniqKey="Pazynina G">GV Pazynina</name></author><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name></author><author><name sortKey="Matrosovich, Mn" uniqKey="Matrosovich M">MN Matrosovich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shinya, K" uniqKey="Shinya K">K Shinya</name></author><author><name sortKey="Ebina, M" uniqKey="Ebina M">M Ebina</name></author><author><name sortKey="Yamada, S" uniqKey="Yamada S">S Yamada</name></author><author><name sortKey="Ono, E" uniqKey="Ono E">E Ono</name></author><author><name sortKey="Kasai, N" uniqKey="Kasai N">N Kasai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Karthik, V" uniqKey="Karthik V">V Karthik</name></author><author><name sortKey="Xiao, Yk" uniqKey="Xiao Y">YK Xiao</name></author><author><name sortKey="Aarthi, C" uniqKey="Aarthi C">C Aarthi</name></author><author><name sortKey="Claudia, P" uniqKey="Claudia P">P Claudia</name></author><author><name sortKey="Rahul, R" uniqKey="Rahul R">R Rahul</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, D" uniqKey="Xu D">D Xu</name></author><author><name sortKey="Newhouse, Ei" uniqKey="Newhouse E">EI Newhouse</name></author><author><name sortKey="Amaro, Re" uniqKey="Amaro R">RE Amaro</name></author><author><name sortKey="Pao, Hc" uniqKey="Pao H">HC Pao</name></author><author><name sortKey="Cheng, Ls" uniqKey="Cheng L">LS Cheng</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Iwata, T" uniqKey="Iwata T">T Iwata</name></author><author><name sortKey="Fukuzawa, K" uniqKey="Fukuzawa K">K Fukuzawa</name></author><author><name sortKey="Nakajima, K" uniqKey="Nakajima K">K Nakajima</name></author><author><name sortKey="Aida Hyugaji, S" uniqKey="Aida Hyugaji S">S Aida-Hyugaji</name></author><author><name sortKey="Mochizuki, Y" uniqKey="Mochizuki Y">Y Mochizuki</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vines, A" uniqKey="Vines A">A Vines</name></author><author><name sortKey="Wells, K" uniqKey="Wells K">K Wells</name></author><author><name sortKey="Matrosovich, M" uniqKey="Matrosovich M">M Matrosovich</name></author><author><name sortKey="Castrucci, Mr" uniqKey="Castrucci M">MR Castrucci</name></author><author><name sortKey="Ito, T" uniqKey="Ito T">T Ito</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Wj" uniqKey="Wang W">WJ Wang</name></author><author><name sortKey="Lu, B" uniqKey="Lu B">B Lu</name></author><author><name sortKey="Zhou, H" uniqKey="Zhou H">H Zhou</name></author><author><name sortKey="Suguitan, Al" uniqKey="Suguitan A">AL Suguitan</name></author><author><name sortKey="Cheng, X" uniqKey="Cheng X">X Cheng</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stevens, J" uniqKey="Stevens J">J Stevens</name></author><author><name sortKey="Blixt, O" uniqKey="Blixt O">O Blixt</name></author><author><name sortKey="Chen, Lm" uniqKey="Chen L">LM Chen</name></author><author><name sortKey="Donis, Ro" uniqKey="Donis R">RO Donis</name></author><author><name sortKey="Paulson, Jc" uniqKey="Paulson J">JC Paulson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ohuchi, R" uniqKey="Ohuchi R">R Ohuchi</name></author><author><name sortKey="Ohuchi, M" uniqKey="Ohuchi M">M Ohuchi</name></author><author><name sortKey="Garten, W" uniqKey="Garten W">W Garten</name></author><author><name sortKey="Klenk, Hd" uniqKey="Klenk H">HD Klenk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author><author><name sortKey="Stevens, Dj" uniqKey="Stevens D">DJ Stevens</name></author><author><name sortKey="Daniels, Rs" uniqKey="Daniels R">RS Daniels</name></author><author><name sortKey="Douglas, Ar" uniqKey="Douglas A">AR Douglas</name></author><author><name sortKey="Knossow, M" uniqKey="Knossow M">M Knossow</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Deshpande, Kl" uniqKey="Deshpande K">KL Deshpande</name></author><author><name sortKey="Fried, Va" uniqKey="Fried V">VA Fried</name></author><author><name sortKey="Ando, M" uniqKey="Ando M">M Ando</name></author><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Das, Sr" uniqKey="Das S">SR Das</name></author><author><name sortKey="Puigbo, P" uniqKey="Puigbo P">P Puigbo</name></author><author><name sortKey="Hensley, Se" uniqKey="Hensley S">SE Hensley</name></author><author><name sortKey="Hurt, De" uniqKey="Hurt D">DE Hurt</name></author><author><name sortKey="Bennink, Jr" uniqKey="Bennink J">JR Bennink</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name></author><author><name sortKey="Ma, C" uniqKey="Ma C">C Ma</name></author><author><name sortKey="Wong, C" uniqKey="Wong C">C Wong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sun, S" uniqKey="Sun S">S Sun</name></author><author><name sortKey="Wang, Q" uniqKey="Wang Q">Q Wang</name></author><author><name sortKey="Zhao, F" uniqKey="Zhao F">F Zhao</name></author><author><name sortKey="Chen, W" uniqKey="Chen W">W Chen</name></author><author><name sortKey="Li, Z" uniqKey="Li Z">Z Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Khurana, S" uniqKey="Khurana S">S Khurana</name></author><author><name sortKey="Verma, S" uniqKey="Verma S">S Verma</name></author><author><name sortKey="Verma, N" uniqKey="Verma N">N Verma</name></author><author><name sortKey="Crevar, Cc" uniqKey="Crevar C">CC Crevar</name></author><author><name sortKey="Carter, Dm" uniqKey="Carter D">DM Carter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, K" uniqKey="Wang K">K Wang</name></author><author><name sortKey="Holtz, Km" uniqKey="Holtz K">KM Holtz</name></author><author><name sortKey="Anderson, K" uniqKey="Anderson K">K Anderson</name></author><author><name sortKey="Chubet, R" uniqKey="Chubet R">R Chubet</name></author><author><name sortKey="Mahmoud, W" uniqKey="Mahmoud W">W Mahmoud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rebecca, Md" uniqKey="Rebecca M">MD Rebecca</name></author><author><name sortKey="Jose, M" uniqKey="Jose M">M Josè</name></author><author><name sortKey="Gonzalo, Im" uniqKey="Gonzalo I">IM Gonzalo</name></author><author><name sortKey="Yuriana, Oa" uniqKey="Yuriana O">OA Yuriana</name></author><author><name sortKey="Stacey, Sc" uniqKey="Stacey S">SC Stacey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, C" uniqKey="Wang C">C Wang</name></author><author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name></author><author><name sortKey="Tseng, Y" uniqKey="Tseng Y">Y Tseng</name></author><author><name sortKey="Hsu, C" uniqKey="Hsu C">C Hsu</name></author><author><name sortKey="Hung, Y" uniqKey="Hung Y">Y Hung</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, M" uniqKey="Li M">M Li</name></author><author><name sortKey="Wang, B" uniqKey="Wang B">B Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bao, Y" uniqKey="Bao Y">Y Bao</name></author><author><name sortKey="Bolotov, P" uniqKey="Bolotov P">P Bolotov</name></author><author><name sortKey="Dernovoy, D" uniqKey="Dernovoy D">D Dernovoy</name></author><author><name sortKey="Kiryutin, B" uniqKey="Kiryutin B">B Kiryutin</name></author><author><name sortKey="Zaslavsky, L" uniqKey="Zaslavsky L">L Zaslavsky</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Larkin, Ma" uniqKey="Larkin M">MA Larkin</name></author><author><name sortKey="Blackshields, G" uniqKey="Blackshields G">G Blackshields</name></author><author><name sortKey="Brown, Np" uniqKey="Brown N">NP Brown</name></author><author><name sortKey="Chenna, R" uniqKey="Chenna R">R Chenna</name></author><author><name sortKey="Mcgettigan, Pa" uniqKey="Mcgettigan P">PA McGettigan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tamura, K" uniqKey="Tamura K">K Tamura</name></author><author><name sortKey="Peterson, D" uniqKey="Peterson D">D Peterson</name></author><author><name sortKey="Peterson, N" uniqKey="Peterson N">N Peterson</name></author><author><name sortKey="Stecher, G" uniqKey="Stecher G">G Stecher</name></author><author><name sortKey="Nei, M" uniqKey="Nei M">M Nei</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Delano, Wl" uniqKey="Delano W">WL DeLano</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ekiert, Dc" uniqKey="Ekiert D">DC Ekiert</name></author><author><name sortKey="Bhabha, G" uniqKey="Bhabha G">G Bhabha</name></author><author><name sortKey="Elsliger, Ma" uniqKey="Elsliger M">MA Elsliger</name></author><author><name sortKey="Friesen, Rh" uniqKey="Friesen R">RH Friesen</name></author><author><name sortKey="Jongeneelen, M" uniqKey="Jongeneelen M">M Jongeneelen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Humphrey, W" uniqKey="Humphrey W">W Humphrey</name></author><author><name sortKey="Dalke, A" uniqKey="Dalke A">A Dalke</name></author><author><name sortKey="Schulten, K" uniqKey="Schulten K">K Schulten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Andreas, Bl" uniqKey="Andreas B">BL Andreas</name></author><author><name sortKey="Von Der Lieth, Cw" uniqKey="Von Der Lieth C">CW von der Lieth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Phillips, Jc" uniqKey="Phillips J">JC Phillips</name></author><author><name sortKey="Braun, R" uniqKey="Braun R">R Braun</name></author><author><name sortKey="Wang, W" uniqKey="Wang W">W Wang</name></author><author><name sortKey="Gumbart, J" uniqKey="Gumbart J">J Gumbart</name></author><author><name sortKey="Tajkhorshid, E" uniqKey="Tajkhorshid E">E Tajkhorshid</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mackerell, Ad" uniqKey="Mackerell A">AD MacKerell</name></author><author><name sortKey="Feig, M" uniqKey="Feig M">M Feig</name></author><author><name sortKey="Brooks, Cl" uniqKey="Brooks C">CL Brooks</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Raman, Ep" uniqKey="Raman E">EP Raman</name></author><author><name sortKey="Guvench, O" uniqKey="Guvench O">O Guvench</name></author><author><name sortKey="Mackerell, Ad" uniqKey="Mackerell A">AD Mackerell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kirschner, Kn" uniqKey="Kirschner K">KN Kirschner</name></author><author><name sortKey="Yongye, Ab" uniqKey="Yongye A">AB Yongye</name></author><author><name sortKey="Tschampel, Sm" uniqKey="Tschampel S">SM Tschampel</name></author><author><name sortKey="Gonzalez Outeiri O, J" uniqKey="Gonzalez Outeiri O J">J González-Outeiriño</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vendome, J" uniqKey="Vendome J">J Vendome</name></author><author><name sortKey="Posy, S" uniqKey="Posy S">S Posy</name></author><author><name sortKey="Jin, X" uniqKey="Jin X">X Jin</name></author><author><name sortKey="Bahna, F" uniqKey="Bahna F">F Bahna</name></author><author><name sortKey="Goran, Ahlsen" uniqKey="Goran A">Ahlsen Goran</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bohne, A" uniqKey="Bohne A">A Bohne</name></author><author><name sortKey="Lang, E" uniqKey="Lang E">E Lang</name></author><author><name sortKey="Von Der Lieth, Cw" uniqKey="Von Der Lieth C">CW von der Lieth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Morris, Gm" uniqKey="Morris G">GM Morris</name></author><author><name sortKey="Huey, R" uniqKey="Huey R">R Huey</name></author><author><name sortKey="Lindstrom, W" uniqKey="Lindstrom W">W Lindstrom</name></author><author><name sortKey="Sanner, Mf" uniqKey="Sanner M">MF Sanner</name></author><author><name sortKey="Belew, Rk" uniqKey="Belew R">RK Belew</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chandrasekaran, A" uniqKey="Chandrasekaran A">A Chandrasekaran</name></author><author><name sortKey="Srinivasan, A" uniqKey="Srinivasan A">A Srinivasan</name></author><author><name sortKey="Raman, R" uniqKey="Raman R">R Raman</name></author><author><name sortKey="Viswanathan, K" uniqKey="Viswanathan K">K Viswanathan</name></author><author><name sortKey="Raguram, S" uniqKey="Raguram S">S Raguram</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ducatez, Mf" uniqKey="Ducatez M">MF Ducatez</name></author><author><name sortKey="Olinger, Cm" uniqKey="Olinger C">CM Olinger</name></author><author><name sortKey="Owoade, Aa" uniqKey="Owoade A">AA Owoade</name></author><author><name sortKey="Tarnagda, Z" uniqKey="Tarnagda Z">Z Tarnagda</name></author><author><name sortKey="Tahita, Mc" uniqKey="Tahita M">MC Tahita</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, R" uniqKey="Chen R">R Chen</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Ha, Y" uniqKey="Ha Y">Y Ha</name></author><author><name sortKey="Stevens, Dj" uniqKey="Stevens D">DJ Stevens</name></author><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author><author><name sortKey="Wiley, Dc" uniqKey="Wiley D">DC Wiley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stevens, J" uniqKey="Stevens J">J Stevens</name></author><author><name sortKey="Blixt, O" uniqKey="Blixt O">O Blixt</name></author><author><name sortKey="Tumpey, Tm" uniqKey="Tumpey T">TM Tumpey</name></author><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name></author><author><name sortKey="Paulson, Jc" uniqKey="Paulson J">JC Paulson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yamada, S" uniqKey="Yamada S">S Yamada</name></author><author><name sortKey="Suzuki, Y" uniqKey="Suzuki Y">Y Suzuki</name></author><author><name sortKey="Suzuki, T" uniqKey="Suzuki T">T Suzuki</name></author><author><name sortKey="Le, Mq" uniqKey="Le M">MQ Le</name></author><author><name sortKey="Nidom, Ca" uniqKey="Nidom C">CA Nidom</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sui, J" uniqKey="Sui J">J Sui</name></author><author><name sortKey="Hwang, Wc" uniqKey="Hwang W">WC Hwang</name></author><author><name sortKey="Perez, S" uniqKey="Perez S">S Perez</name></author><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Aird, D" uniqKey="Aird D">D Aird</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bohne, A" uniqKey="Bohne A">A Bohne</name></author><author><name sortKey="Lieth, Cw" uniqKey="Lieth C">CW Lieth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Moustakas, Dt" uniqKey="Moustakas D">DT Moustakas</name></author><author><name sortKey="Lang, Pt" uniqKey="Lang P">PT Lang</name></author><author><name sortKey="Pegg, S" uniqKey="Pegg S">S Pegg</name></author><author><name sortKey="Pettersen, Et" uniqKey="Pettersen E">ET Pettersen</name></author><author><name sortKey="Kuntz, Id" uniqKey="Kuntz I">ID Kuntz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nicholls, Jm" uniqKey="Nicholls J">JM Nicholls</name></author><author><name sortKey="Chan, R" uniqKey="Chan R">R Chan</name></author><author><name sortKey="Russell, Rj" uniqKey="Russell R">RJ Russell</name></author><author><name sortKey="Air, Gm" uniqKey="Air G">GM Air</name></author><author><name sortKey="Peiris, M" uniqKey="Peiris M">M Peiris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Leo, F" uniqKey="Leo F">F Leo</name></author><author><name sortKey="Asakuma, S" uniqKey="Asakuma S">S Asakuma</name></author><author><name sortKey="Fukuda, K" uniqKey="Fukuda K">K Fukuda</name></author><author><name sortKey="Senda, A" uniqKey="Senda A">A Senda</name></author><author><name sortKey="Urashima, T" uniqKey="Urashima T">T Urashima</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kaverin, Nv" uniqKey="Kaverin N">NV Kaverin</name></author><author><name sortKey="Rudneva, Ia" uniqKey="Rudneva I">IA Rudneva</name></author><author><name sortKey="Llyushina, Na" uniqKey="Llyushina N">NA Llyushina</name></author><author><name sortKey="Varich, Nl" uniqKey="Varich N">NL Varich</name></author><author><name sortKey="Lipatov, As" uniqKey="Lipatov A">AS Lipatov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Salahuddin, P" uniqKey="Salahuddin P">P Salahuddin</name></author><author><name sortKey="Khan, Au" uniqKey="Khan A">AU Khan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turner, Bg" uniqKey="Turner B">BG Turner</name></author><author><name sortKey="Summers, Mf" uniqKey="Summers M">MF Summers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rota, Pa" uniqKey="Rota P">PA Rota</name></author><author><name sortKey="Oberste, Ms" uniqKey="Oberste M">MS Oberste</name></author><author><name sortKey="Monroe, Ss" uniqKey="Monroe S">SS Monroe</name></author><author><name sortKey="Nix, Wa" uniqKey="Nix W">WA Nix</name></author><author><name sortKey="Campagnoli, R" uniqKey="Campagnoli R">R Campagnoli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Neu, U" uniqKey="Neu U">U Neu</name></author><author><name sortKey="Woellner, K" uniqKey="Woellner K">K Woellner</name></author><author><name sortKey="Gauglitz, G" uniqKey="Gauglitz G">G Gauglitz</name></author><author><name sortKey="Stehle, T" uniqKey="Stehle T">T Stehle</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tsai, B" uniqKey="Tsai B">B Tsai</name></author><author><name sortKey="Gilbert, Jm" uniqKey="Gilbert J">JM Gilbert</name></author><author><name sortKey="Stehle, T" uniqKey="Stehle T">T Stehle</name></author><author><name sortKey="Lencer, W" uniqKey="Lencer W">W Lencer</name></author><author><name sortKey="Benjamin, Tl" uniqKey="Benjamin T">TL Benjamin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, T" uniqKey="Lin T">T Lin</name></author><author><name sortKey="Wang, G" uniqKey="Wang G">G Wang</name></author><author><name sortKey="Li, A" uniqKey="Li A">A Li</name></author><author><name sortKey="Zhang, Q" uniqKey="Zhang Q">Q Zhang</name></author><author><name sortKey="Wu, C" uniqKey="Wu C">C Wu</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">22719948</article-id><article-id pub-id-type="pmc">3375263</article-id><article-id pub-id-type="publisher-id">PONE-D-12-04973</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0038794</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Computational Biology</subject><subj-group><subject>Macromolecular Structure Analysis</subject><subj-group><subject>Macromolecular Complex Analysis</subject></subj-group></subj-group><subj-group><subject>Population Modeling</subject><subj-group><subject>Infectious Disease Modeling</subject></subj-group></subj-group><subj-group><subject>Evolutionary Modeling</subject></subj-group></subj-group><subj-group><subject>Evolutionary Biology</subject><subj-group><subject>Evolutionary Processes</subject><subj-group><subject>Genetic Drift</subject></subj-group></subj-group></subj-group><subj-group><subject>Microbiology</subject><subj-group><subject>Virology</subject><subj-group><subject>Viral Transmission and Infection</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Medicine</subject><subj-group><subject>Infectious Diseases</subject><subj-group><subject>Viral Diseases</subject><subj-group><subject>Influenza</subject></subj-group></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis</article-title><alt-title alt-title-type="running-head">Glycosylation Affect the Preference of H5N1 Virus</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Chen</surname><given-names>Wentian</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Sun</surname><given-names>Shisheng</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Li</surname><given-names>Zheng</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an, People’s Republic of China</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Department of Pathology, Clinical Chemistry Division, Johns Hopkins University, Baltimore, United States of America</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Brown</surname><given-names>Earl G.</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1">University of Ottawa, Canada</aff><author-notes><corresp id="cor1">* E-mail: <email>zhengli@nwu.edu.cn</email></corresp><fn fn-type="con"><p>Conceived and designed the experiments: WC SS ZL. Performed the experiments: WC SS. Analyzed the data: WC SS ZL. Contributed reagents/materials/analysis tools: WC. Wrote the paper: WC ZL.</p></fn></author-notes><pub-date pub-type="collection"><year>2012</year></pub-date><pub-date pub-type="epub"><day>14</day><month>6</month><year>2012</year></pub-date><volume>7</volume><issue>6</issue><elocation-id>e38794</elocation-id><history><date date-type="received"><day>20</day><month>2</month><year>2012</year></date><date date-type="accepted"><day>10</day><month>5</month><year>2012</year></date></history><permissions><copyright-statement>Chen et al.</copyright-statement><copyright-year>2012</copyright-year><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p></license></permissions><abstract><p>Increasing numbers of H5N1 influenza viruses (IVs) are responsible for human deaths, especially in North Africa and Southeast Asian. The binding of hemagglutinin (HA) on the viral surface to host sialic acid (SA) receptors is a requisite step in the infection process. Phylogenetic analysis reveals that H5N1 viruses can be divided into 10 clades based on their HA sequences, with most human IVs centered from clade 1 and clade 2.1 to clade 2.3. Protein sequence alignment in various clades indicates the high conservation in the receptor-binding domains (RBDs) is essential for binding with the SA receptor. Two glycosylation sites, 158N and 169N, also participate in receptor recognition. In the present work, we attempted to construct a serial H5N1 HA models including diverse glycosylated HAs to simulate the binding process with various SA receptors <italic>in silico</italic>. As the SA-α-2,3-Gal and SA-α-2,6-Gal receptor adopted two distinctive topologies, straight and fishhook-like, respectively, the presence of N-glycans at 158N would decrease the affinity of HA for all of the receptors, particularly SA-α-2,6-Gal analogs. The steric clashes of the huge glycans shown at another glycosylation site, 169N, located on an adjacent HA monomer, would be more effective in preventing the binding of SA-α-2,3-Gal analogs.</p></abstract><counts><page-count count="14"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>The influenza viruses (IVs) have had an important effect on human history. Three influenza A pandemics have taken place in 1918–1919 (H1N1), 1957 (H2N2) and 1968 (H3N2), which killed up to 50 million people all over the world. According to phylogenetic analysis, IVs originated in waterfowl; due to antigenic mutations, including antigenic drift and antigenic shift, the mutation in the receptor-binding domains (RBDs) could result in a human pandemic similar to H1N1 <xref rid="pone.0038794-Zachary1" ref-type="bibr">[1]</xref>. Since the last century, there are increasing numbers of reports about avian IVs, such as the H5N1, H7N7 and H9N2 virus, that have infected humans <xref rid="pone.0038794-Obenauer1" ref-type="bibr">[2]</xref>, <xref rid="pone.0038794-Fouchier1" ref-type="bibr">[3]</xref>. As of 02 May 2012, 603 confirmed human cases with avian influenza H5N1 virus were reported to the WHO since 2003; of these patients, 356 died (fatality rate: 59.03%) <xref rid="pone.0038794-WHOOIEFAO1" ref-type="bibr">[4]</xref>.</p><p>The recognition, attachment and infection of IVs require the participation of a kind of enveloped glycoprotein, named hemagglutinin (HA), which is a homotrimeric transmembrane protein. Three RBDs at the tip of HA binding to the host sialoglycans are essential for endocytosis, each monomer can be cleaved by proteases in phagosome into a globular head and a stem region, known as HA1 and HA2, respectively <xref rid="pone.0038794-Skehel1" ref-type="bibr">[5]</xref>, <xref rid="pone.0038794-Wiley1" ref-type="bibr">[6]</xref>. The whole HA unfold and expose their interior HA2 in the acid environment, then the fusion peptides in HA2 insert themselves into the host membrane and lock the membrane bilayers. Finally, uncharacterized changes to its conformation pull the two membranes together <xref rid="pone.0038794-Senne1" ref-type="bibr">[7]</xref>, <xref rid="pone.0038794-Sammalkorpi1" ref-type="bibr">[8]</xref>.</p><p>As the key role in the infection process, HA is thought to be crucial in host switching. HA interacts with the distal sialic acid (SA) and the penultimate β-D-galactose (Gal) on the host glycans. Commonly, the human IVs prefer the SA-α-2,6-Gal terminal glycan, whereas the avian IVs prefer SA-α-2,3-Gal <xref rid="pone.0038794-Gambaryan1" ref-type="bibr">[9]</xref>. SA-α-2,6-Gal is the characteristic glycan pattern in the upper respiratory airway of humans, but a small amount of SA-α-2,3-Gal also appears in the lower respiratory airway and breast milk in humans <xref rid="pone.0038794-Shinya1" ref-type="bibr">[10]</xref>, <xref rid="pone.0038794-Karthik1" ref-type="bibr">[11]</xref>. To cross the species barrier and result in the human-to-human infections, HA must change its receptor-binding preference from SA-α-2,3-Gal to SA-α-2,6-Gal <xref rid="pone.0038794-Xu1" ref-type="bibr">[12]</xref>.</p><p>Studies of H5N1 virus demonstrate that receptor recognition occurs in the RBD. The RBD consists of three secondary structure elements and four conserved residues: Helix190, Loop130 and Loop220 combined with 98Tyr, 153Trp, 183His, and 195Tyr <xref rid="pone.0038794-Iwata1" ref-type="bibr">[13]</xref>. The Q226L and G228S mutations in HA of H3N2 and H2N2 virus, correlated with the host switching from SA-α-2,3-Gal to SA-α-2,6-Gal receptor specificity <xref rid="pone.0038794-Vines1" ref-type="bibr">[14]</xref>. Actually, as the pocket is shallow, the mutation of one or two residues in the RBD would be responsible for host switching. The S227N mutation resulted in binding to SA-α-2,3-Gal receptors with lower affinity, which caused several human infections in Hong Kong and Egypt <xref rid="pone.0038794-Wang1" ref-type="bibr">[15]</xref>. The combination of the Q226L and G228S mutations would be more efficient in binding SA-α-2,6-Gal, though there is no such mutation in H5N1 <xref rid="pone.0038794-Stevens1" ref-type="bibr">[16]</xref>.</p><p>Glycosylation plays an important role in the viral life cycle. Among the different subtypes of IV HA, there are extensive variations in the glycosylation sites of the head region, whereas the stem region are more conserved <xref rid="pone.0038794-Ohuchi1" ref-type="bibr">[17]</xref>. Glycans near the antigenic epitopes interfere with antibody recognition to escape immune pressure from hosts, and those near HA0 cleavage site protect their infectivity <xref rid="pone.0038794-Skehel2" ref-type="bibr">[18]</xref>, <xref rid="pone.0038794-Deshpande1" ref-type="bibr">[19]</xref>. The gain or loss of glycosylation sites near the RBD can also affect viral receptor recognition <xref rid="pone.0038794-Das1" ref-type="bibr">[20]</xref>, <xref rid="pone.0038794-Chen1" ref-type="bibr">[21]</xref>. Sun <italic>et al.</italic> describe several alternative models of glycosylation sites that affect the influenza virus, such as increasing the number of glycosylation sites and changing their positional conversion. One conspicuous feature is that, when a new human epidemic virus emerges, the virus would reveal a high affinity to the host target. As co-evolution occurs with the host, the virulence of the virus would decrease, accompanied by a lower affinity and increasing numbers of glycosylation site <xref rid="pone.0038794-Sun1" ref-type="bibr">[22]</xref>.</p><p>The N-glycan synthesis of viral glycoprotein depended on the host cell. Non-glycosylated HA can be obtained by expression in <italic>Escherichia coli</italic><xref rid="pone.0038794-Khurana1" ref-type="bibr">[23]</xref>. HA expressed in insect cells, such as S2 cells, differs from the HA in mammalian cells, with heterogeneous mannose-type glycans in the N-glycan terminus <xref rid="pone.0038794-Wang2" ref-type="bibr">[24]</xref>. Additionally, high-mannose-type glycans can be produced in human cells deficient in a glycosyltransferase named N-acetylglucosaminyltransferase I (GnTI-HEK293 cells) <xref rid="pone.0038794-Rebecca1" ref-type="bibr">[25]</xref>. Wang <italic>et al.</italic> detect the HAs with gradually increasing glycosylation by glycan microarrays and find that their binding ability is not only related to different SA candidates but also to the glycosylation around the RBD. The truncation of the N-glycans on HA increased its binding affinity while decreasing its specificity toward disparate SA receptors <xref rid="pone.0038794-Wang3" ref-type="bibr">[26]</xref>.</p><p>To further explore the binding affinity between HA and SA receptors, computational methods were adopted <xref rid="pone.0038794-Xu1" ref-type="bibr">[12]</xref>, <xref rid="pone.0038794-Iwata1" ref-type="bibr">[13]</xref>, <xref rid="pone.0038794-Li1" ref-type="bibr">[27]</xref>. Commonly, researchers focused on the interactions between HA and different SA receptors, with less consideration of the N-glycosylation on HA. In the present study, the contribution of the glycosylation of HA to host recognition was analyzed by available H5N1 HA sequences and crystal structures. A phylogenetic analysis was used to investigate the key amino-acid residues and glycosylation sites in the RBD. Further construction of models, molecular dynamic (MD) simulations and docking analysis were utilized to investigate the structure and composition of the N-glycans and their effects on HA activity.</p></sec><sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title><sec id="s2a"><title>Protein Sequence Data and Phylogenetic Analysis</title><p>A dataset of the protein sequences for H5N1 HA was retrieved from the National Center for Biotechnology Information (NCBI) flu database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/genomes/FLU">http://www.ncbi.nlm.nih.gov/genomes/FLU</ext-link>) accessed in Sep. 2011 <xref rid="pone.0038794-Bao1" ref-type="bibr">[28]</xref>, containing approximately 330 human reports in 3602 sequences. An alignment of whole sequences was performed by ClustalW 2.0 (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>) <xref rid="pone.0038794-Larkin1" ref-type="bibr">[29]</xref>. To investigate the distribution of N-glycosylation sites in different clades, a smaller dataset was used for further analysis with representative virus, including human-infecting viruses such as A/Hong Kong/486/1997, A/Hong Kong/213/2003, A/Indonesia/560H/2006 and A/Cambodia/V0219301/2011, as well as candidate vaccine viruses, such as A/Viet Nam/1203/2004, A/bar-headed goose/Qinghai/1A/2005 and A/Egypt/N03072/2010. Unrooted phylogenetic tree was constructed using the Neighbor-Joining method and the Poisson correction model in MEGA 5.05 <xref rid="pone.0038794-Tamura1" ref-type="bibr">[30]</xref>. The internal branching probabilities were determined by bootstrap analysis with 1,000 replicates.</p></sec><sec id="s2b"><title>MD Simulations for HA Models</title><p>The crystal structure of the H5N1 HA (Protein Data Bank/PDB code: 3GBM) was used as the original coordinate file. Only the chain A and B in 3GBM, which represent HA1 and HA2, respectively, were extracted by PyMOL 0.99 and used as an HA monomer <xref rid="pone.0038794-DeLano1" ref-type="bibr">[31]</xref>, <xref rid="pone.0038794-Ekiert1" ref-type="bibr">[32]</xref>. This new HA monomer was denoted by the A/Viet Nam/1203/2004 strain (04VN). Based on the alignment analysis, a 03HK-RBD-like HA model was modified by S159N, T160A and N227S, replayed with the mutator 1.3 plugin of VMD <xref rid="pone.0038794-Humphrey1" ref-type="bibr">[33]</xref>, and denoted as 03HK.</p><p>To explore the influence of N-glycans near the RBD, four glycosylated HA models were constructed by the “GlycoProt” program (<ext-link ext-link-type="uri" xlink:href="http://www.glycoscience.de">http://www.glycoscience.de</ext-link>) <xref rid="pone.0038794-Andreas1" ref-type="bibr">[34]</xref>. In consideration of existing studies, four glycans were added only at the site 158N in 04VN (<xref ref-type="fig" rid="pone-0038794-g001">Figure 1A and 1B</xref>). The names of four models described their glycan feature: mono N-acetyl-D-glucosamine HA (MG) had only one N-acetyl-D-glucosamine (GlcNAc) residue on 158N, high-mannosylation HA (HM) represented a mannose-rich glycan and the native full-glycosylation HA (FG) was commonly expressed from the mammalian cell system. Comparatively, the de-sialic acid HA (DS) model was a typical production treated with neuramidinase <xref rid="pone.0038794-Wang3" ref-type="bibr">[26]</xref>.</p><fig id="pone-0038794-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g001</object-id><label>Figure 1</label><caption><title>Schematic diagrams for HA and SA receptors models.</title><p>(A) Five initial models of HA are represented by 04VN, MG, HM, DS and FG respectively. Only one glycan is added on 158N near the RBD. (B) Four glycans on glycosylated HAs are represented by MG, HM, DS and FG respectively. (C) Eight sialoglycans and one pentaglucose are used for docking assays. The abbreviations for each glycan are as follows: lactosialyltetraoses (LSTa/LSTc), disialyllacto-N-tetraose (3DSLNT/6DSLNT) and bisialyantennas based on two monosaccharides (bisialyantennary mannose, BM3/BM6 and bisialyantennary GlcNAc, BG3/BG6). The sequence of monosaccharides and glycosidic bonds are illustrated using Consortium for Functional Glycomics nomenclature.</p></caption><graphic xlink:href="pone.0038794.g001"></graphic></fig><p>To perform the MD simulation, six HA models were solvated in rectangular water boxes, and the largest system had 70,000 atoms (including protein, glycan and waters) in a 140 Å×70 Å×70 Å cuboid of TIP3P water molecules. Minimization and equilibration were performed using the NAMD 2.8 program <xref rid="pone.0038794-Phillips1" ref-type="bibr">[35]</xref> with the CHARMM22 all-atom force field for the protein <xref rid="pone.0038794-MacKerell1" ref-type="bibr">[36]</xref>; the missing parameters for the glycans were adopted from the GLYCAM06 and top_all36_carb V1.83 force field <xref rid="pone.0038794-Raman1" ref-type="bibr">[37]</xref>, <xref rid="pone.0038794-Kirschner1" ref-type="bibr">[38]</xref>. A periodic boundary condition was used along with the particle mesh Ewald method for electrostatic interaction and a 12.5 Å cutoff for Van der Waals interaction. All of the systems were first energy-minimized for 0.5 ns with the protein and glycan atoms fixed and then for another 0.5 ns with all of the atoms relaxed. After the minimization, the system was heated from 0 to 300 K in 1 ns and then equilibrated for 4 ns with the pressure and temperature controlled (NVT). The temperature was held at 300 K using Langevin dynamics, and the pressure was held at 1 atm by the Langevin piston method. Production runs were subsequently performed for another 10 ns for each HA model.</p><p>The root-mean-square deviation (RMSD) was used to evaluate their stability, which was defined as the active range of the non-hydrogen atoms in the protein backbone and the pyranose ring in the glycans during the production runs. The root-mean-square fluctuation (RMSF) was used to compare the flexibility of the glycan residues and amino-acid residues. The RMSF reflects the fluctuation of the α-C atoms in amino-acid residues or the C2 atoms in glycans. The RMSF of the protein backbone was calculated by sampling every 2 ps, whereas the RMSD was extracted at 2 ps interval during the production runs <xref rid="pone.0038794-Vendome1" ref-type="bibr">[39]</xref>.</p></sec><sec id="s2c"><title>Docking Analysis of HAs and SA Receptor Complexes</title><p>To compare the interaction between HAs and SA receptors, four SA-α-2,3-Gal analogs and four SA-α-2,6-Gal analogs were designed and constructed to mimic the sialoglycans on human and avian cells. These SA receptors could be divided into four types: lactosialyltetraoses (LSTa/LSTc), disialyllacto-N-tetraose (3DSLNT/6DSLNT; 6DSLNT does not actually exist in humans) and bisialyantennary based on two monosaccharides (bisialyantennary mannose, BM3/BM6 and bisialyantennary GlcNAc, BG3/BG6). Additionally, a 5-glucose-pentose (5GLC) was used as a negative quality control (<xref ref-type="fig" rid="pone-0038794-g001">Figure 1C</xref>). All of the receptors were constructed by the SWEET2 program and optimized by the MM3 force field <xref rid="pone.0038794-Bohne1" ref-type="bibr">[40]</xref>.</p><p>Six reasonable HA models gained after the 16 ns MD simulation were used for the flexible docking. Each model was docked with 9 receptors by turns. AutoDock 4.2 was used to perform the flexible docking in which the grid spacing was set to 0.375 Å and each grid map consisted of 80×80×80 grid points in three dimensions <xref rid="pone.0038794-Morris1" ref-type="bibr">[41]</xref>. For each docking analysis, the grid center coordinates were set as the mean coordinates of the RBD. All of the glycosidic linkages from the receptors were treated as flexible linkages and allowed to rotate in 20-degree increments, together with key amino-acid residues in the RBD. With the other options set to the default values, 100 docking runs were performed during each docking experiment, and the results were analyzed by cluster analysis in their largest cluster. For each flexible docking, the three binding models with the lowest binding energies were selected for further statistical analysis. The flexible docking calculation and complementary analysis were performed on a Windows XP workstation.</p></sec><sec id="s2d"><title>MD Simulations for Free Sialoglycans and Sialoglycan-HA Complexes</title><p>While the knowledge about topological structures of free or docked sialoglycans is limited, MD simulations make it possible to explore the adpoted characteristics. The explicitly solvated MD simulation of sialoglycans for 50 ns, together with 5 ns docked MD simulation in trimer HA have provided essential dynamic information about free sialoglycans conformations. The θ angle parameter, which define the angle between the C2 atom of SA1 and the C1 atoms of Gal2 and GlcNAc3 (or Man3), was used to describe the topologies adopted by various sialoglycans <xref rid="pone.0038794-Chandrasekaran1" ref-type="bibr">[42]</xref>. As the θ angle was maked up of the last three residues, it could reflect the dynamic changes during simulation complementally.</p><p>The initial conformations of 5 ns docked MD simulation were attained from optimal trimer sialoglycan-FG docking complexes by symmetry transformation script, then three complex N-glycans were added on 169N in the trimer docked complexes <xref rid="pone.0038794-Humphrey1" ref-type="bibr">[33]</xref>, <xref rid="pone.0038794-Andreas1" ref-type="bibr">[34]</xref>. For each trimer complexes, the MD simulation were performed through the same procedure as described previously.</p><p>The volumetric analysing method, which presented the overall spatial volume sampled by the glycans, used for monitoring the N-glycans and sialoglycans. As the powerful visual method to judge the interaction between the sialoglycans and N-glycans, the volumetric methods were constructed by averaging atomic density of the extracted MD frames using the VMD Volmap Tool plugin with 1 Å resolution grids <xref rid="pone.0038794-Humphrey1" ref-type="bibr">[33]</xref>.</p><p>All the MD simulations (NAMD_2.8b1_Linux-x86_64-CUDA) were performed on an ubuntu 10.04-based workstation, and the benchmark maintained 1 ns/day with the maximal system.</p></sec></sec><sec id="s3"><title>Results</title><sec id="s3a"><title>Phylogenetic Analysis of Glycosylation Sites of HA in H5N1 Virus</title><p>There have been several works about the phylogenetic analysis in H5N1. The evolutionary dynamics of eleven genes in H5N1 virus are complicated, but the substitution rates of HA and neuramidinase (NA) are significantly higher than the rates of the internal protein segments (such as PA or M1 ) <xref rid="pone.0038794-Ducatez1" ref-type="bibr">[43]</xref>, <xref rid="pone.0038794-Chen2" ref-type="bibr">[44]</xref>. According to a nomenclature system devised for H5N1 in 2008, the HAs can be divided into ten phylogenetic clades (0–9), of which three clades, namely clade 0, clade 1 and clade 2, are responsible for most human infections <xref rid="pone.0038794-WHOOIEFAO2" ref-type="bibr">[45]</xref>. There was a marked divergence in the last four years, as several new subclades emerged and predominated in the H5N1 epidemics. When some groups in one monophyletic clade meet the divergence criteria (bootstrap values >60), they can split into second, third or even fourth-order clades <xref rid="pone.0038794-WHOOIEFAO3" ref-type="bibr">[46]</xref>. The clade 2 began to spread westward during 2005; reached the Middle East, Europe and Africa in 2006; and caused hundreds of deaths within the past five years. Within clade 2, clades 2.1, 2.2 and 2.3 were responsible for most human cases in North Africa and Southeastern Asian. In addition, clades 2.4 and 2.5 had no relationship with human infections. Another major group, clade 1, which was isolated mainly in China, Cambodia, Thailand and Vietnam, was responsible for most documented human infections in Asia during 2004–2005, Thailand during 2006 and Cambodia during 2010–2012. Remarkably, after a 5-year quiet period, clade 1.1 reemerged in Cambodia with a 100% human fatality rate <xref rid="pone.0038794-WHOOIEFAO1" ref-type="bibr">[4]</xref>. However, several previously circulating clades were not detected for the past 4 years, and it was likely that these clades were supplanted by the new clades. Based on an updated unified nomenclature system, the current circulating H5N1 virus in the past two years are centralized primarily in clades 1.1, 2.2.1, 2.2.1.1, 2.2.2, 2.3.2.1, 2.3.4, 2.3.4.2 and 7.1 <xref rid="pone.0038794-WHOOIEFAO3" ref-type="bibr">[46]</xref>.</p><p>The homology of the 3602 HA dataset was rather low, 25 out 325 amino-acid residues in HA1 revealed an identity rate of 100% (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>). Those conserved residues were closely associated with the essential functions of HA; four residues located at the bottom of the RBD, 95Tyr, 153Trp, 183His and 195Trp exhibited an identity rate of 100%, and other residues in the RBD were also conservative relatively (<xref ref-type="fig" rid="pone-0038794-g002">Figure 2</xref>). One serious human isolate of H5N1 HA (A/Hong Kong/213/2003) had a single amino-acid residue difference in the RBD compared with most avian and human IVs. The S227N would result in an enhancement of its ability to bind to a human host <xref rid="pone.0038794-Stevens1" ref-type="bibr">[16]</xref>. Actually, since 2003, nine human records appeared in ten S227N mutations, such as the A/Cambodia/V0606311/2011 (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>).</p><fig id="pone-0038794-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g002</object-id><label>Figure 2</label><caption><title>The amino-acid residues in the RBD of H5N1 HA.</title><p>(A) The RBD in H5N1 HA consists of three secondary structure elements, Loop130, Loop220 and Helix190, together with four 100% conserved residues at the bottom (red). The residues with a conservation rate higher than 99% are labeled in yellow. (B) The statistics of the conserved amino-acid residues in the RBD. The predominant amino-acid residues are underlined, and the number of different types of mutations are shown as various blocks. Those with conservation rates lower than 99% are labeled. All of the residue numbers were adopted from the H3 HA numbering system.</p></caption><graphic xlink:href="pone.0038794.g002"></graphic></fig><p>As is known, N-X-T/S (X cannot be a proline) is the motif of an N-glycosylation site. The statistical analysis of 3,602 records elucidated that seven conservative glycosylation sites appeared in H5N1 HA. The conservation rates of glycosylation at 15N, 27N and 290N in HA1 and 154N and 213N in HA2 were 99.86%, 99.75%, 99.69%, 99.88%, and 99.73%, respectively (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>). Another two glycosylation sites appear to be more complex. The conservation rate of glycosylation on 169N was 95.73%. All of the 169N glycosylation site deficiencies presented an avian specificity, and the earliest deficient record was A/Turkey/England/1991, the root of the tree in <xref ref-type="fig" rid="pone-0038794-g003">Figure 3</xref>. This glycosylation site was conserved until 2006; henceforth, sporadic deficiency cases have been reported in China, Vietnam and Egypt. Most of these cases were in clade 2.2.1 or clade 7.1. The conservation rate of the glycosylation on 158N was 47.61%. Strains with glycosylated and non-glycosylated 158N have co-existed extensively even before 1997. According to our phylogenetic analysis, viruses with 158N glycosylation site deficiency existed in nine clades extensively (all except clade 7, <xref ref-type="fig" rid="pone-0038794-g003">Figure 3</xref>). Most recent human records belonging to clade 2.2.1 in North Africa stayed a stable 158N glycosylation site deficiency. In contrast, the human virus belonging to clade 1.1 or clade 2.3.4 in Southeast Asia were conservative for the 158N glycosylation site. This discrepancy confirms that no close relationship exists between the loss of the 158N glycosylation site and host switching. Moreover, approximately 55 dual-glycosylation site (158N and 169N) deficiency virus were responsible for avian infections, most of which clustered in clade 2.2 during 2007 and 2008.</p><fig id="pone-0038794-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g003</object-id><label>Figure 3</label><caption><title>Phylogenetic tree of 103 HAs from H5N1 strains.</title><p>Phylogenetic trees are inferred from protein sequences by the Neighbor-Joining method and rooted using A/turkey/England/1991(red text). Estimates of the statistical significance of the phylogenies are calculated by performing 1,000 bootstrap replicates. The lengths of the horizontal line are proportional to the numbers of protein sequence differences, as indicated by the scale bars. Different clades classified by the WHO are shown as grayish-white bars. Virus with glycosylation site 158N deficiency are labeled in red, whereas glycosylation site 169N deficiency are labeled in green, and the dual deficiencies are labeled in purple. Several important mutation strains mentioned in the article are labeled in blue.</p></caption><graphic xlink:href="pone.0038794.g003"></graphic></fig><p>To confirm the evolutionary position of the HAs in the PDB and remodel the comparable HA models, further investigation was undertaken in H5N1 HAs from the PDB. Seven H5 HA structure files were available, of which three incomplete structures, 1JSO, 1JSM and 1JSN, represent the A/Duck/Singapore/3/97 monomer <xref rid="pone.0038794-Ha1" ref-type="bibr">[47]</xref>. The files of 2FKO and 3GBM correspond to A/Viet Nam/1203/2004, whereas 2IBX and 3FKU represented A/Viet Nam/1194/2004 <xref rid="pone.0038794-Stevens2" ref-type="bibr">[48]</xref>, <xref rid="pone.0038794-Yamada1" ref-type="bibr">[49]</xref>, <xref rid="pone.0038794-Sui1" ref-type="bibr">[50]</xref>. The similarity of the last four HA sequence was 98.27%, with the identical RBD (<xref ref-type="fig" rid="pone-0038794-g004">Figure 4</xref>). Accordingly, 3GBM would represent four HAs in subsequent experiments due to its minimal structural resolution. Although two strains mentioned above are responsible for human infections, Yamada <italic>et al</italic>. showed that A/Viet Nam/1194/2004 still had a high affinity for SA-α-2,3-Gal receptors <xref rid="pone.0038794-Yamada1" ref-type="bibr">[49]</xref>. Coincidentally, A/Viet Nam/1194/2004 and A/Viet Nam/1203/2004 belong to clade 1 and both are chosen as candidate vaccine viruses. A/Hong Kong/213/2003 is another definite human virus that also belongs to clade 1. As <xref ref-type="fig" rid="pone-0038794-g004">Figure 4</xref> showed, the RBDs of various H5N1 clades were highly conserved, the only differences among the three strains were residue 227 and the glycosylation site at 158N. A T160A mutation in “158NST” would result in the deficiency of this glycosylated site.</p><fig id="pone-0038794-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g004</object-id><label>Figure 4</label><caption><title>Sequence alignment of A/Hong Kong/486/97, A/Viet Nam/1203/2004(3GBM), A/Hong Kong/213/2003, A/Cambodia/S1211394/2008, A/Egypt/2321-NAMRU3/2007, A/Anhui/1/2005 and A/chicken/Shanxi/10/2006.</title><p>The RBDs are shown in blue lines with the identical amino acids marked with asterisks. The glycosylation site 158N are boxed in red, whereas the four most conserved residues at the bottom of the RBD are boxed in green.</p></caption><graphic xlink:href="pone.0038794.g004"></graphic></fig></sec><sec id="s3b"><title>MD Simulation for HA Models Reveals Stable RBD and Relative Flexible Glycans</title><p>Glycosylation exerts a significant influence on envelope glycoproteins. To perform the MD simulation for glycoprotein, the complete coordinate file is necessary. Unfortunately, complete larger glycans are too flexible to yield sufficient electron density; consequently, their 3D structures cannot be derived using X-ray crystallography. Many oligosaccharide NOEs (Nuclear Overhouser Effects) also cannot be resolved and are difficult to assign by NMR. Additionally, there are often too few inter-residue NOEs to render an unambiguous 3D structure determination possible <xref rid="pone.0038794-Bohne2" ref-type="bibr">[51]</xref>. Given this situation, glycosylated HAs must be reconstructed manually.</p><p>The individual RMSFs in the RBD showed significant difference. Those relatively active residues, which have a higher RMSF and less impacted by steric clashes, may change their conformation to accommodate the receptors. During the 16 ns MD simulations, the energetic and structural properties were monitored for the 6 systems. The RMSDs and the energies converged in all of the systems, which indicated well-behaved simulations (<xref ref-type="fig" rid="pone-0038794-g005">Figure 5A</xref>). We found that the helical structure and the overall protein backbones were stable during the whole equilibrium, whereas the loop experienced larger fluctuations relatively. Loop220 may fluctuate more than Loop130 and Helix190, based on their respective RMSFs. Considering the results from the phylogenetic analysis, four most conserved residues, 95Tyr, 153Trp, 183His and 195Tyr, located at the bottom of the RBD, had lower RMSFs compare with the non-conserved residues (<xref ref-type="fig" rid="pone-0038794-g005">Figure 5B</xref>). By analyzing the MD trajectories, the replaced residues in 03HK would affect their acroteric amino-acid residues. The RMSF of 222Lys in 03HK decreased, which was responsible for the interaction of the basic group in 222Lys with the newly introduced, weakly acidic amide group of the S227N residue, a motion that made the two sides of Loop220 become closer. Moreover, by changing S159N and T160A in 03HK, the RMSFs of the adjacent residues were also changed, which indicated the steric alteration.</p><fig id="pone-0038794-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g005</object-id><label>Figure 5</label><caption><title>The dynamical properties of the HA models during MD simulation.</title><p>(A) Six convergent RMSDs of the HA models indicate the stable states, while highly glycosylated HAs experienced larger fluctuation against the initial conformations. (B) The RMSF values reflect the fluctuation of those residues in the RBD. The S227N mutation in 03HK would result in a relatively stable 222Lys (with a lower RMSF value in 222). The glycans at 158N of affect their acroteric amino-acid residues. (C) The RMSF of the three glycans on HM, DS, and FG present the complicated praxiology that the glycan residues in the core are more stable compared with the terminal glycan residues. Columns 2–5 represent the RMSF in the core of the N-glycan compared with the glycan residues at the branches and terminus.</p></caption><graphic xlink:href="pone.0038794.g005"></graphic></fig><p>The behavior of N-glycans on HAs was also investigated. The RMSDs of the longer glycans in HM, DS and FG had similar features: at the minimization stage, the RMSDs rose vigorously; however, the whole glycan became stable and convergent immediately. Each monosaccharide in the N-glycan had a corresponding RMSF that reflects its activity, such as the mannoses in HM, GalNAc in DS and SA in FG; as was shown in <xref ref-type="fig" rid="pone-0038794-g005">Figure 5C</xref>, the higher RMSFs in the distal monosaccharide residues revealed their active fluctuation.</p></sec><sec id="s3c"><title>Energy Analysis from HAs and SA Receptor Docking Complexes</title><p>The interaction between HA and SA receptors could be deemed as a kind of docking mechanism. Several software programs for molecular docking analysis have been developed in recent years, such as Dock and AutoDock <xref rid="pone.0038794-Moustakas1" ref-type="bibr">[52]</xref>, <xref rid="pone.0038794-Morris1" ref-type="bibr">[41]</xref>. AutoDock 4.2 adopts the Lamarckian genetic algorithm to estimate the binding energies and explore the most probable docking conformation, of which the docking conformation with the minimal value may be closely related to a natural state.</p><p>Based on the binding energies provided from the output files, all HAs showed a conspicuous affinity for α-2,3-sialoglycans or α-2,6-sialoglycans (<xref ref-type="fig" rid="pone-0038794-g006">Figure 6</xref>). Interestingly, as the negative quality control, the binding energy from HA-5GLC complexes presented a positive value of approximately +3- +6 kcal/mol, which corresponded to a nonbinding property between HA and glucoses. In the docking assays, all the HAs presented relatively lower energies with α-2,3-sialoglycans and higher values with α-2,6-sialoglycans, which indicated that the HA of A/Viet Nam/1203/2004 would still have a preference for avian sialoglycans, although 04VN was derived from a human infection. One explanation for this paradox could be that humans exposed to a high concentration of A/Viet Nam/1203/2004 strains would also be infected <xref rid="pone.0038794-Yamada1" ref-type="bibr">[49]</xref>.</p><fig id="pone-0038794-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g006</object-id><label>Figure 6</label><caption><title>The binding energies collected from flexible docking.</title><p>The binding energies are calculated from the three lowest energies provided in the largest clusters. A notable tendency is the increasing binding energy accompanying with increasing glycosylation.</p></caption><graphic xlink:href="pone.0038794.g006"></graphic></fig><p>A notable tendency was the increasing binding energy accompanying increasing glycosylation, which demonstrated that the binding affinities were impaired by complicated glycans. The average binding energies collected from the glycosylated HAs were higher than those from 04VN and 03HK.</p><p>Another remarkable feature in the 03HK-sialoglycans docked complexes which puzzled us, was that there is no obvious numerical difference between the two types of sialoglycan docking assays (<xref ref-type="fig" rid="pone-0038794-g006">Figure 6</xref>). Although a standard error of approximately 2–3 kcal/mol in predicted binding energies would be acceptable <xref rid="pone.0038794-Morris1" ref-type="bibr">[41]</xref>, we supposed that the existing amino-acid residue mutations in A/Hong Kong/213/2003 did not enhance the binding affinity to α-2,6-sialoglycans but decreased its ability to bind α-2,3-sialoglycans. An analogical description in H3N2 HA noted that the trans conformation of SA-α-2,3-Gal directs the Gal ring away from the receptor-binding domain surface, and the amide group of Gln226 formed an H-bond with the axial 4-hydroxyl group of the penultimate Gal <xref rid="pone.0038794-Li1" ref-type="bibr">[27]</xref>. In this case, the total binding energy would be impaired by the 227Ser with a hydroxyl group being replaced by the Asn with a larger amide group.</p></sec><sec id="s3d"><title>Conformation Analysis from HAs and SA Receptor Docking Complexes</title><p>Carboxylate, hydroxyl, N-acetyl and glycerol group are four functional groups in SA residue that participate in numerous hydrogen bonds (H-bonds) and hydrophobic interactions with the RBD. Those H-bonds formed from the hydroxyl groups in the penultimate Gal and Loop220 are responsible for host preference <xref rid="pone.0038794-Nicholls1" ref-type="bibr">[53]</xref>.</p><p>Commonly, the SA in α-2,3-sialoglycans and α-2,6-sialoglycans adopt the trans or cis conformation, respectively, which is oriented to the 158N glycosylation site <xref rid="pone.0038794-Xu1" ref-type="bibr">[12]</xref>, <xref rid="pone.0038794-Li1" ref-type="bibr">[27]</xref>. In all of the docking complexes, the carboxylate, hydroxyl, N-acetyl and glycerol groups in SA buried inside the shallow pocket (<xref ref-type="fig" rid="pone-0038794-g007">Figure 7</xref>). This reflected that HA, as a kind of glycan-binding protein (GBP), had a natural tendency to bind SA. The penultimate Gal were distinctive based on the preference of HA, as not all of the Gals would fit into the RBD due to the steric clashes of amino-acid residues or N-glycan, much less the antepenultimate GlcNAc/Man in receptors. Most antepenultimate residues would form a β1-3 or β1-4 glycosidic linkage with the penultimate β-D-Gal. As a result, the α-2,3-sialoglycan and α-2,6-sialoglycan conformations were predominantly straight and fishhook-like, respectively (<xref ref-type="supplementary-material" rid="pone.0038794.s001">Figure S1</xref>, <xref ref-type="supplementary-material" rid="pone.0038794.s002">S2</xref>). The topologies of the α-2,3-sialoglycans was extrorse, whereas the α-2,6-sialoglycans, which lean toward Helix190, were ental.</p><fig id="pone-0038794-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g007</object-id><label>Figure 7</label><caption><title>The docking complexes between six HAs and eight sialoglycans.</title><p>SA-α-2,3-Gal receptors are superimposed with distal SA residues at left, whereas the SA-α-2,6-Gal receptors are superimposed with distal SA residues at right. LSTa/LSTc are shown in blue, 3DSLNT/6DSLNT are shown in green, BM3/BM6 are shown in red and BG3/BG6 are shown in yellow. The types of receptors deploy distinctive topologies: most SA-α-2,3-Gal receptors are straight and extrorse, whereas the SA-α-2,6-Gal receptors are fishhook-like and ental. (A) 03HK-sialoglycan docking complexes, (B) 04VN-sialoglycan docking complexes, (C) MG-sialoglycan docking complexes, (D) HM-sialoglycan docking complexes, (E) DS-sialoglycan docking complexes, (F) FG-sialoglycan docking complexes.</p></caption><graphic xlink:href="pone.0038794.g007"></graphic></fig><p>As the long glycans, LSTa and LSTc would interact with the RBD in a pre-description <xref rid="pone.0038794-Xu1" ref-type="bibr">[12]</xref>. Due to the fishhook-like topology, the remainder of the residues in LSTc exhibited more contact with Helix190 compared with LSTa. The branched 3DSLNT is the most abundant sialyl human milk oligosaccharide; in contrast, 6DSLNT does not exist in human milk <xref rid="pone.0038794-Leo1" ref-type="bibr">[54]</xref>. Another α-2,6 glycosidic linkage was formed between a second SA and a GlcNAc, however, this α-2,6-linked SA could not fit into the RBD. Although there were two SA-α-2,6-Gal or SA-α-2,3-Gal linkages in biantennary sialoglycans, these did not enhance the binding affinity. Subsequent conformation analysis showed that only the SA in one antenna would fit into the RBD, whereas the remaining components interacted slightly with several residues. Actually, one biantennary sialoglycan molecule cannot dock with two RBDs in a trimer HAs due to the short distance between the two antennas. The estimated curvilinear distance between two RBDs is approximately 60 Å, whereas the maximal distances between two SAs in most bi- or tri-antennary sialoglycans are shorter (30–40 Å).</p><p>The SA receptors would interact with the RBD through two forces:hydrogen bond (H-bond) and hydrophobic interaction. H-bond was the essential interaction in stabilizing ligand-receptor complexes. There were copious polar groups in the receptors and RBD that facilitate hydrogen-bond formation. During the course of dynamic bonding, several persistent and transient hydrogen bonds could be observed: most persistent H-bonds were formed from the hydroxyl groups offered by SA and Gal with the side-chain groups. In all of the docking complexes, the four most conserved residues located at the bottom of the RBD have formed stable H-bonds with SA: the glycerol group in SA interacted with the p-hydroxyphenyl group of 195Tyr, the imidazole group of 183His interacted with the glycerol group of SA and the carboxyl group in SA formed a stable H-bond with the p-hydroxyphenyl group of 95Tyr. Obviously, Loop220 played an important part in the formation of H-bonds with the penultimate Gal. Many observable transient H-bonds were formed from the remainder of sialoglycans and non-conserved residues around the RBD. One typical example was the residue 194Leu or 190Glu, located in Helix190, which would form an H-bond with the antepenultimate GlcNAc/Man, especially in α-2,6-sialoglycans. In various hydrophobic interactions, one common stacking force was formed by the N-acetyl group in SA and 153Trp, whereas other aromatic groups are far from the docking pocket. Several residues, such as Ser226, Trp153, His183, or Gln191, would also form hydrophobic interactions with sialoglycans, as in the previous descriptions <xref rid="pone.0038794-Li1" ref-type="bibr">[27]</xref>.</p><p>The Ser has a hydroxy group at the γ-C atom that can form an H-bond with other polar groups, such as water or the receptor. In the docking conformations collected previously, however, the residue 227Ser did not participate directly in stabilizing interactions with the penultimate Gal. A persistent H-bond between 226Gln and Gal formed in several HA-α-2,3-sialoglycan docking complexes but rarely in HA-α-2,6-sialoglycans docking complexes <xref rid="pone.0038794-Xu1" ref-type="bibr">[12]</xref>. After the introduction of the S227N mutation, an amide group would interact with the adjacent basic group in 222Lys, forming a more compact Loop220, which moved 226Gln away from its previous position. This change would decrease the probability of the above H-bond formation, which caused most S227N strains exhibiting a decreased affinity for α-2,3-sialoglycans. Moreover, almost all S227N mutations are followed by the loss of the 158N glycosylation site (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>), the deficiency of this glycosylation site is conducive to bind the α-2,6-sialoglycans. Hence this slight change may be one reason responsible for more human infections <xref rid="pone.0038794-Wang3" ref-type="bibr">[26]</xref>.</p><p>We found that all of the glycosylated HAs except MG exhibited a clearly weakened binding ability. One factor was the steric clashes that would block the complex, bi- or tri-antennary sialoglycans to the target proteins. The topologies of glycans, the distribution of hydroxy groups in different saccharide residues, the amide group in GalNAc and GlcNAc or the groups in SA would also affect the binding affinity. Analysis of binding energies and conformations indicated that N-glycans on 158N have more influence on fishhook-like α-2,6-sialoglycans.</p></sec><sec id="s3e"><title>Glycosylation as a Subsidiary Contribution in Receptor Recognition</title><p>Several epitope sites in the tip of HA1 are near the RBD, and one function of glycosylation is to provide the shield to escape the host immune system <xref rid="pone.0038794-Kaverin1" ref-type="bibr">[55]</xref>, <xref rid="pone.0038794-Salahuddin1" ref-type="bibr">[56]</xref>. As a type of hydrophilic long-chain and high-mass molecule, N-glycans have only one covalent bond to the Asn, which discloses that the N-glycans are rather mobile and flexible. Undoubtedly, the presence of huge N-glycans would impact the binding efficacy of the RBD.</p><p>One factor that cannot be neglected in trimeric HAs is that some glycosylation sites appear to be far from one RBD (>25 Å), but their glycans still affect the binding affinity for the adjoining RBD in another monomer. Such as the 169N glycosylation site: the distance between the center of the RBD in chain A and 169N is 34.18 Å, whereas the distance to 158N is 18.72 Å; however, it would be 20.05 Å between the RBD and the adjoining 169N (measured by the nitrogen atom in Asn and the weight center of the RBD, from a trimeric HA coordinate file coded 2IBX). A previous survey showed that the conservation rate of the 169N glycosylation site is 95.73%, whereas an S171A change would result in the loss of this glycosylation site (<xref ref-type="supplementary-material" rid="pone.0038794.s004">File S1</xref>).</p><p>As observed from the trimer HA in <xref ref-type="fig" rid="pone-0038794-g008">Figure 8</xref>, we found that one RBD in chain A was located among three glycosylated sites: 158N in chain A and 169N in chain E, with the distance from the center of the RBD to the center of each glycan maintained at approximately 18 to 21 Å. Additionally, the distance between the center of the RBD and the third glycosylated site, 158N in chain E, is approximately 35–40 Å. As summarized from the MD simulation and docking complexes, the SA-α-2,3-Gal and SA-α-2,6-Gal receptors deployed two distinctive topologies (<xref ref-type="supplementary-material" rid="pone.0038794.s001">Figure S1</xref>, <xref ref-type="supplementary-material" rid="pone.0038794.s002">S2</xref> and <xref ref-type="supplementary-material" rid="pone.0038794.s003">S3</xref>). It could be hypothesized that the shapes and components of the glycans on HA have various influences: huge glycans around the RBD would restrict the affinity and specificity of HA.</p><fig id="pone-0038794-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0038794.g008</object-id><label>Figure 8</label><caption><title>The influence of vicinal N-glycans on the RBD of HA trimer.</title><p>(A) The volumetric topologies of the N-glycans near the RBD during the 5 ns MD simulation. The complex N-glycans on the sites 158N and 169N would swing dramatically. (B) The distances from the weight center of the RBD to the three topological centers of the N-glycans are calculated at 10 ps intervals, three colors represent the corresponding distances in (A).</p></caption><graphic xlink:href="pone.0038794.g008"></graphic></fig><p>The N-glycans and sialoglycans topologies were examined using the volumes sampled. The orientations and shapes of the volumetric maps varied dramatically for all the glycans: three N-glycans on 158N and 169N were swinging dramatically while the sialoglycans constrained in the RBD, especially for the SA residues. Straight and extrorse topologies were observed in the free and docked SA-α-2,3-Gal receptors. Bound SA-α-2,6-Gal receptors adopted tightly folded and fishhook-like conformations and bented toward the Helix190. These two distinctive conformations resulted in different preference of contaction to the N-glycans on 158N or 169N (<xref ref-type="supplementary-material" rid="pone.0038794.s001">Figure S1</xref>, <xref ref-type="supplementary-material" rid="pone.0038794.s002">S2</xref> and <xref ref-type="supplementary-material" rid="pone.0038794.s003">S3</xref>). Obviously, the loss of the 169N glycosylation site has an advantage in binding the straight and extrorse α-2,3-sialoglycans, whereas the loss of the 158N glycosylation site facilitates the recognition of the fishhook-like and ental α-2,6-sialoglycans. With more complicated glycosylation, the reasonable sialoglycans for docking would be more restricted.</p></sec></sec><sec id="s4"><title>Discussion</title><p>There is one specified group of glycoproteins appearing on viral envelope that participates in the recognition and invasion, such as HA in IVs, gp160 in HIV and the spike (S) glycoprotein in the SARS coronaviruses <xref rid="pone.0038794-Turner1" ref-type="bibr">[57]</xref>, <xref rid="pone.0038794-Rota1" ref-type="bibr">[58]</xref>. At the very first stage, these proteins recognize a specified glycan pattern on the surface of the host cell; from this perspective, these glycoproteins can be deemed as a kind of lectin or GBP. Many viruses target the distal SA, including Simian Virus 40 (SV40) and murine polyomavirus (Polyoma) <xref rid="pone.0038794-Neu1" ref-type="bibr">[59]</xref>, <xref rid="pone.0038794-Tsai1" ref-type="bibr">[60]</xref>.</p><p>HA have a high specificity in host recognition. The mechanism in H1N1 HA is relatively clear: outbreak and epidemic HAs have a high preference for SA-α-2,6-Gal, which is distinctive in humans compared with SA-α-2,3-Gal in avian hosts. The E190D and G225D mutations in the RBD of H1N1 would result in host switching <xref rid="pone.0038794-Lin1" ref-type="bibr">[61]</xref>. Similarly, the Q226L and G228S mutations in H2N2 and H3N2 would also lead to host switching. However, H5N1, as one of the highly pathogenic avian IVs, has not been reported to have such definite mutations until now. In the past five years, several subclades have predominated in the H5N1 virus, and some subclades exhibited a preference for human infection and were responsible for a large number of human deaths. In 2011, a 100% human fatality rate of clade 1.1 in Cambodia caused a great panic; however, genetic and <italic>in vitro</italic> binding analysis indicated that the receptor specificity of the H5N1 virus is still of the avian type, which suggested that the restriction by receptor specificity may not be as stringent as previously thought in certain avian viruses <xref rid="pone.0038794-Lin1" ref-type="bibr">[61]</xref>.</p><p>Though there are not sufficient cumulative data to analyze H5N1 compared with H1N1 or H3N2, prospective work can be performed by computer simulation. One common method is to explore the relationship between host switching and protein sequences. However, secondary sequences are not enough to further investigate the reaction mechanism. By conformation analysis of a three-dimensional structure, a concrete target can be found readily. The RBD in H5N1 HA mainly consisted of three conserved secondary structures; moreover, acroteric modifications, such as glycosylation, would also affect the binding preference. Conformation analysis showed that the acquisition or deficiency of an N-glycosylation site would be more effective than simple amino-acid residue mutations. Furthermore, the composition and shapes of the N-glycans on HA are dependent on the host N-glycan synthetic system in endoplasmic reticulum. Thus, the modification of heterogeneous N-glycans make the HA more intricate <xref rid="pone.0038794-Chen1" ref-type="bibr">[21]</xref>.</p><p>In this study, we presented a series of methods and models to simulate the effect of glycosylation on receptor binding <italic>in silico</italic>. The previous conclusion revealed that the RBD pocket of HA in H5N1 was shallow and specified; one terminal residue, SA, could fit into the RBD easily in different HAs. Detailed and complete studies have been conducted with several methods, and the mechanism of host switching and praxiology of HAs was clear. No covalent bonds existed between HA and its receptors, but the abundant H-bonds and hydrophobic bonds facilitated the docking conformation. A re-constructed 03HK-RBD-like HA model revealed that 03HK has a negative influence on the binding of SA-α-2,3-Gal. This change could be one reason that 03HK has a dual-host affinity. As mentioned previously, the 04VN model derived from A/Viet Nam/1203/2004, a vaccine candidate for humans, but still had a higher affinity SA-α-2,3-Gal.</p><p>Docking assays from four glycosylated HAs revealed similar outcomes to that of 04VN, but with numerical disparities. Generally, the HA models with increasing glycans have a gradually decreasing binding ability; larger glycans on 158N in the prominent antigenic loop on the edge of the RBD would obviously decrease the binding affinity for the SA-α-2,6-Gal analogs. Recent H5N1 viruses (especially in clade 2.2), which are characterized by a deficiency of the 158N glycosylation site, exhibit an increased propensity for acquiring human receptor specificity. This effect could be due to the large glycans on 158N, which would hamper the affinity for the fishhook-like α-2,6-sialoglycans, based on our results.</p><p>A typical monomer affinity in GBP is measured at the millimolar level, which is quite weak compared with enzymatic reactions, so-called “high affinity but low avidity”. To achieve the high avidity, like other GBPs, a complete HA consists with three monomers. The topologies of the three RBDs in the HA trimer face to three symmetry directions would facilitate the binding of the host SA glycans. As the primary interaction occurs in the RBD and SA and Gal, the remainder of the components in sialoglycans would also interact with the acroteric amino-acid residues. It is often overlooked that glycosylation at 169N of an adjacent monomer would also affect the binding affinity, especially for the straight SA-α-2,3-Gal analogs. Coincidentally, the deficiency of the 169N glycosylation site only appears in avian infection virus. The result showed that the N-glycans on two glycosylation sites would form the extensive steric clashes and restrict receptor specificity.</p><p>Actually, the factors for cross-species and human-to-human transmission in IVs are a composite mechanism. The preference switching of HA is only one necessary step; another envelope glycoprotein, NA, the function of which is related to the releasing of the virus, also participates in human infections. A definite human infection would be affected by virus, host and environment. To prevent a potential human epidemic caused by H5N1 virus, more investigations must be performed.</p></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0038794.s001"><label>Figure S1</label><caption><p><bold>The superimpositions of sialoglycans during 50</bold><bold>ns explicitly solvated MD simulation reflect the free</bold> S<bold>A-α-2,3-Gal and SA-α-2,6-Gal receptors adopted two distinctive topologies.</bold> (A) As the initial sialyglacans are shown as the cyan CPK models, nine representative conformations of sialoglycan from 5 ns intervals are superimposed with distal SA residues. It clearly revealed the remainder of all SA-α-2,3-Gal receptors are leftward, compare to the rightward SA-α-2,6-Gal receptors, which result in two orientations.</p><p>(TIF)</p></caption><media xlink:href="pone.0038794.s001.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0038794.s002"><label>Figure S2</label><caption><p><bold>Topological θ angle plots of the free sialoglycans during 50 ns.</bold> (A) the θ angle in four SA-α-2,3-Gal receptors swinging between 110° and 180°. (B) the θ angle in four SA-α-2,6-Gal receptors swinging between 60° and 110°.</p><p>(TIF)</p></caption><media xlink:href="pone.0038794.s002.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0038794.s003"><label>Figure S3</label><caption><p><bold>The volumetric topologies of glycans in the Sialoglycans and trimer HA complexes.</bold> As it indicated that the θ angles of sialoglycans are stable during 5 ns MD simultion in the trimer HA, all the SA-α-2,3-Gal receptors maintain a straight-like topology (110°<θ<180°) and SA-α-2,6-Gal receptors maintain a fishhook-like topology (60°<θ<110°) resepectively. The orientations and shapes of the volumetric maps vary dramatically for the N-glycans on 158N and 169N while the sialoglycans varied in the smaller spatial volume in RBD. More complicated N-glycans would sterically hinder the receptor binding, even with different preference.</p><p>(TIF)</p></caption><media xlink:href="pone.0038794.s003.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0038794.s004"><label>File S1</label><caption><p><bold>The alignment file containing the 3,602 protein sequences of H5N1 HA1.</bold> As is shown in the file, the glycosylation sites are marked in red, whereas the 100% conserved residues are marked in cyan.</p><p>(XLS)</p></caption><media xlink:href="pone.0038794.s004.xls"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was supported by the International Science & Technology Cooperation Program (2009DFA32730) and the Key Science and Technology Program (2010K12-02(5)) of the Shaanxi Province of China. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group><ref-list><title>References</title><ref id="pone.0038794-Zachary1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zachary</surname><given-names>S</given-names></name><name><surname>Rahul</surname><given-names>R</given-names></name><name><surname>Karthik</surname><given-names>V</given-names></name><name><surname>Ram</surname><given-names>S</given-names></name></person-group><year>2009</year><article-title>Context-specific target definition in influenza A virus hemagglutinin-glycan receptor interactions.</article-title><source>Chem Biol</source><volume>16</volume><fpage>803</fpage><lpage>814</lpage><pub-id pub-id-type="pmid">19716471</pub-id></element-citation></ref><ref id="pone.0038794-Obenauer1"><label>2</label><element-citation publication-type="journal"><collab>Obenauer JC, Denson J, Mehta PK, Su XP, et a1</collab><year>2006</year><article-title>Large-scale sequence analysis of avian influenza isolates.</article-title><source>Science 311</source><volume>(5767)</volume><fpage>1576</fpage><lpage>1580</lpage></element-citation></ref><ref id="pone.0038794-Fouchier1"><label>3</label><element-citation publication-type="journal"><collab>Fouchier RA, Schneeberger PM, Rozendaal FW, Broekman JM, et a1</collab><year>2004</year><article-title>Avian influenza A virus (H7N7) associated with human conjunctivitis and a fatal case of acute respiratory distress syndrome.</article-title><source>PNAS</source><volume>101(5)</volume><fpage>1356</fpage><lpage>1361</lpage><pub-id pub-id-type="pmid">14745020</pub-id></element-citation></ref><ref id="pone.0038794-WHOOIEFAO1"><label>4</label><element-citation publication-type="other"><collab>WHO/OIE/FAO H5N1 Evolution Working Group</collab><year>2012</year><article-title>Accessed 2012 May 5.</article-title><comment>Cumulative number of confirmed human cases of avian influenza A/(H5N1): <ext-link ext-link-type="uri" xlink:href="http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/index.html">http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/index.html</ext-link>.</comment></element-citation></ref><ref id="pone.0038794-Skehel1"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Skehel</surname><given-names>JJ</given-names></name><name><surname>Wiley</surname><given-names>DC</given-names></name></person-group><year>2000</year><article-title>Receptor binding and membrane fusion in virus entry: the influenza hemagglutinin.</article-title><source>Annu Rev Biochem</source><volume>69</volume><fpage>531</fpage><lpage>569</lpage><pub-id pub-id-type="pmid">10966468</pub-id></element-citation></ref><ref id="pone.0038794-Wiley1"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wiley</surname><given-names>DC</given-names></name><name><surname>Skehel</surname><given-names>JJ</given-names></name></person-group><year>1987</year><article-title>The structure and function of the hemagglutinin membrane glycoprotein of influenza virus, Annu Rev Biochem</article-title><volume>56</volume><fpage>365</fpage><lpage>395</lpage></element-citation></ref><ref id="pone.0038794-Senne1"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Senne</surname><given-names>DA</given-names></name><name><surname>Panigrahy</surname><given-names>B</given-names></name><name><surname>Kawaoka</surname><given-names>Y</given-names></name><name><surname>Pearson</surname><given-names>JE</given-names></name><name><surname>Suss</surname><given-names>J</given-names></name><etal></etal></person-group><year>1996</year><article-title>Survey of the hemagglutinin (HA) cleavage site sequence of H5 and H7 avian influenza viruses: Amino acid sequence at the HA cleavage site as a marker of pathogenicity potential.</article-title><source>Avian Dis</source><volume>40</volume><fpage>425</fpage><lpage>437</lpage><pub-id pub-id-type="pmid">8790895</pub-id></element-citation></ref><ref id="pone.0038794-Sammalkorpi1"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sammalkorpi</surname><given-names>M</given-names></name><name><surname>Lazaridis</surname><given-names>T</given-names></name></person-group><year>2007</year><article-title>Configuration of influenza hemagglutinin fusion peptide monomers and oligomers in membranes.</article-title><source>Biochem Bioph Res Co</source><volume>1768</volume><fpage>30</fpage><lpage>38</lpage></element-citation></ref><ref id="pone.0038794-Gambaryan1"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gambaryan</surname><given-names>AS</given-names></name><name><surname>Tuzikov</surname><given-names>AB</given-names></name><name><surname>Pazynina</surname><given-names>GV</given-names></name><name><surname>Webster</surname><given-names>RG</given-names></name><name><surname>Matrosovich</surname><given-names>MN</given-names></name><etal></etal></person-group><year>2004</year><article-title>H5N1 chicken influenza viruses display a high binding affinity for Neu5Acalpha2,3Galbeta1–4(6-HSO3)Glc-NAc-containing receptors.</article-title><source>Virology</source><volume>326</volume><fpage>310</fpage><lpage>316</lpage><pub-id pub-id-type="pmid">15302215</pub-id></element-citation></ref><ref id="pone.0038794-Shinya1"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shinya</surname><given-names>K</given-names></name><name><surname>Ebina</surname><given-names>M</given-names></name><name><surname>Yamada</surname><given-names>S</given-names></name><name><surname>Ono</surname><given-names>E</given-names></name><name><surname>Kasai</surname><given-names>N</given-names></name><etal></etal></person-group><year>2006</year><article-title>Avian flu: influenza virus receptors in the human airway.</article-title><source>Nature</source><volume>440</volume><fpage>435</fpage><lpage>436</lpage><pub-id pub-id-type="pmid">16554799</pub-id></element-citation></ref><ref id="pone.0038794-Karthik1"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Karthik</surname><given-names>V</given-names></name><name><surname>Xiao</surname><given-names>YK</given-names></name><name><surname>Aarthi</surname><given-names>C</given-names></name><name><surname>Claudia</surname><given-names>P</given-names></name><name><surname>Rahul</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>Determinants of glycan receptor specificity of H2N2 influenza A virus hemagglutinin.</article-title><source>PloS One</source><volume>5</volume><fpage>e13768</fpage></element-citation></ref><ref id="pone.0038794-Xu1"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname><given-names>D</given-names></name><name><surname>Newhouse</surname><given-names>EI</given-names></name><name><surname>Amaro</surname><given-names>RE</given-names></name><name><surname>Pao</surname><given-names>HC</given-names></name><name><surname>Cheng</surname><given-names>LS</given-names></name><etal></etal></person-group><year>2009</year><article-title>Distinct glycan topology for avian and human sialopentasaccharide receptor analogues upon binding different hemagglutinins: A molecular dynamics perspective.</article-title><source>J Mol Biol</source><volume>387</volume><fpage>465</fpage><lpage>491</lpage><pub-id pub-id-type="pmid">19356594</pub-id></element-citation></ref><ref id="pone.0038794-Iwata1"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iwata</surname><given-names>T</given-names></name><name><surname>Fukuzawa</surname><given-names>K</given-names></name><name><surname>Nakajima</surname><given-names>K</given-names></name><name><surname>Aida-Hyugaji</surname><given-names>S</given-names></name><name><surname>Mochizuki</surname><given-names>Y</given-names></name><etal></etal></person-group><year>2008</year><article-title>Theoretical analysis of binding specificity of influenza viral hemagglutinin to avian and human receptors based on the fragment molecular orbital method.</article-title><source>Comput Biol Chem</source><volume>32</volume><fpage>198</fpage><lpage>211</lpage><pub-id pub-id-type="pmid">18485828</pub-id></element-citation></ref><ref id="pone.0038794-Vines1"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vines</surname><given-names>A</given-names></name><name><surname>Wells</surname><given-names>K</given-names></name><name><surname>Matrosovich</surname><given-names>M</given-names></name><name><surname>Castrucci</surname><given-names>MR</given-names></name><name><surname>Ito</surname><given-names>T</given-names></name><etal></etal></person-group><year>1998</year><article-title>The role of influenza A virus hemagglutinin residues 226 and 228 in receptor specificity and host range restriction.</article-title><source>J Virol</source><volume>72</volume><fpage>7626</fpage><lpage>7631</lpage><pub-id pub-id-type="pmid">9696865</pub-id></element-citation></ref><ref id="pone.0038794-Wang1"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>WJ</given-names></name><name><surname>Lu</surname><given-names>B</given-names></name><name><surname>Zhou</surname><given-names>H</given-names></name><name><surname>Suguitan</surname><given-names>AL</given-names></name><name><surname>Cheng</surname><given-names>X</given-names></name><etal></etal></person-group><year>2010</year><article-title>Glycosylation at 158N of the hemagglutinin protein and receptor binding specificity synergistically affect the antigenicity and immunogenicity of a live attenuated H5N1 A/Vietnam/1203/2004 vaccine virus in ferrets.</article-title><source>J Virol</source><volume>84</volume><fpage>6570</fpage><lpage>6577</lpage><pub-id pub-id-type="pmid">20427525</pub-id></element-citation></ref><ref id="pone.0038794-Stevens1"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stevens</surname><given-names>J</given-names></name><name><surname>Blixt</surname><given-names>O</given-names></name><name><surname>Chen</surname><given-names>LM</given-names></name><name><surname>Donis</surname><given-names>RO</given-names></name><name><surname>Paulson</surname><given-names>JC</given-names></name><etal></etal></person-group><year>2008</year><article-title>Recent avian H5N1 viruses exhibit increased propensity for acquiring human receptor specificity.</article-title><source>J Mol Biol</source><volume>381</volume><fpage>1382</fpage><lpage>1394</lpage><pub-id pub-id-type="pmid">18672252</pub-id></element-citation></ref><ref id="pone.0038794-Ohuchi1"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ohuchi</surname><given-names>R</given-names></name><name><surname>Ohuchi</surname><given-names>M</given-names></name><name><surname>Garten</surname><given-names>W</given-names></name><name><surname>Klenk</surname><given-names>HD</given-names></name></person-group><year>1997</year><article-title>Oligosaccharides in the stem region maintain the influenza virus hemagglutinin in the metastable form required for fusion activity.</article-title><source>J Virol</source><volume>71</volume><fpage>3719</fpage><lpage>3725</lpage><pub-id pub-id-type="pmid">9094646</pub-id></element-citation></ref><ref id="pone.0038794-Skehel2"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Skehel</surname><given-names>JJ</given-names></name><name><surname>Stevens</surname><given-names>DJ</given-names></name><name><surname>Daniels</surname><given-names>RS</given-names></name><name><surname>Douglas</surname><given-names>AR</given-names></name><name><surname>Knossow</surname><given-names>M</given-names></name><etal></etal></person-group><year>1984</year><article-title>A carbohydrate side-chain on hemagglutinins of Hong Kong influenza-viruses inhibits recognition by a monoclonal-antiboby.</article-title><source>Immunology</source><volume>81</volume><fpage>1779</fpage><lpage>1783</lpage></element-citation></ref><ref id="pone.0038794-Deshpande1"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deshpande</surname><given-names>KL</given-names></name><name><surname>Fried</surname><given-names>VA</given-names></name><name><surname>Ando</surname><given-names>M</given-names></name><name><surname>Webster</surname><given-names>RG</given-names></name></person-group><year>1987</year><article-title>Glycosylation affects cleavage of an H5N2 influenza-virus hemagglutinin and regulates virulence.</article-title><source>Biochemistry</source><volume>84</volume><fpage>36</fpage><lpage>40</lpage></element-citation></ref><ref id="pone.0038794-Das1"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Das</surname><given-names>SR</given-names></name><name><surname>Puigbo</surname><given-names>P</given-names></name><name><surname>Hensley</surname><given-names>SE</given-names></name><name><surname>Hurt</surname><given-names>DE</given-names></name><name><surname>Bennink</surname><given-names>JR</given-names></name><etal></etal></person-group><year>2010</year><article-title>Glycosylation focuses sequence variation in the influenza A virus H1 hemagglutinin globular domain.</article-title><source>PloS Pathogens</source><volume>6</volume><fpage>e1001211</fpage><pub-id pub-id-type="pmid">21124818</pub-id></element-citation></ref><ref id="pone.0038794-Chen1"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>J</given-names></name><name><surname>Ma</surname><given-names>C</given-names></name><name><surname>Wong</surname><given-names>C</given-names></name></person-group><year>2011</year><article-title>Vaccine design of hemagglutinin glycoprotein against influenza.</article-title><source>Trends Biotechnol</source><volume>29</volume><fpage>426</fpage><lpage>434</lpage><pub-id pub-id-type="pmid">21640418</pub-id></element-citation></ref><ref id="pone.0038794-Sun1"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname><given-names>S</given-names></name><name><surname>Wang</surname><given-names>Q</given-names></name><name><surname>Zhao</surname><given-names>F</given-names></name><name><surname>Chen</surname><given-names>W</given-names></name><name><surname>Li</surname><given-names>Z</given-names></name></person-group><year>2011</year><article-title>Glycosylation site alteration in the evolution of influenza A (H1N1) viruses.</article-title><source>PLoS One</source><volume>6</volume><fpage>e22844</fpage><pub-id pub-id-type="pmid">21829533</pub-id></element-citation></ref><ref id="pone.0038794-Khurana1"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Khurana</surname><given-names>S</given-names></name><name><surname>Verma</surname><given-names>S</given-names></name><name><surname>Verma</surname><given-names>N</given-names></name><name><surname>Crevar</surname><given-names>CC</given-names></name><name><surname>Carter</surname><given-names>DM</given-names></name><etal></etal></person-group><year>2010</year><article-title>Properly folded bacterially expressed H1N1 hemagglutinin globular head and ectodomain vaccines protect ferrets against H1N1 pandemic influenza virus.</article-title><source>PLoS One</source><volume>5</volume><fpage>e11548</fpage><pub-id pub-id-type="pmid">20634959</pub-id></element-citation></ref><ref id="pone.0038794-Wang2"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>K</given-names></name><name><surname>Holtz</surname><given-names>KM</given-names></name><name><surname>Anderson</surname><given-names>K</given-names></name><name><surname>Chubet</surname><given-names>R</given-names></name><name><surname>Mahmoud</surname><given-names>W</given-names></name><etal></etal></person-group><year>2006</year><article-title>Expression and purification of an influenza hemagglutinin-one step closer to a recombinant protein-based influenza vaccine.</article-title><source>Vaccine</source><volume>24</volume><fpage>2176</fpage><lpage>2185</lpage><pub-id pub-id-type="pmid">16310896</pub-id></element-citation></ref><ref id="pone.0038794-Rebecca1"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rebecca</surname><given-names>MD</given-names></name><name><surname>Josè</surname><given-names>M</given-names></name><name><surname>Gonzalo</surname><given-names>IM</given-names></name><name><surname>Yuriana</surname><given-names>OA</given-names></name><name><surname>Stacey</surname><given-names>SC</given-names></name><etal></etal></person-group><year>2010</year><article-title>The influenza A virus hemagglutinin glycosylation state affects receptor-binding specificity.</article-title><source>J Virol</source><volume>85</volume><fpage>865</fpage><lpage>872</lpage><pub-id pub-id-type="pmid">21068239</pub-id></element-citation></ref><ref id="pone.0038794-Wang3"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>C</given-names></name><name><surname>Chen</surname><given-names>J</given-names></name><name><surname>Tseng</surname><given-names>Y</given-names></name><name><surname>Hsu</surname><given-names>C</given-names></name><name><surname>Hung</surname><given-names>Y</given-names></name><etal></etal></person-group><year>2009</year><article-title>Glycans on influenza hemagglutinin affect receptor binding and immune response.</article-title><source>PANS</source><volume>106</volume><fpage>18137</fpage><lpage>18142</lpage></element-citation></ref><ref id="pone.0038794-Li1"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>M</given-names></name><name><surname>Wang</surname><given-names>B</given-names></name></person-group><year>2006</year><article-title>Computational studies of H5N1 hemagglutinin binding with SA-α-2,3-Gal and SA-α-2, 6-Gal.</article-title><source>Biochem Bioph Res Co</source><volume>347</volume><fpage>662</fpage><lpage>668</lpage></element-citation></ref><ref id="pone.0038794-Bao1"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname><given-names>Y</given-names></name><name><surname>Bolotov</surname><given-names>P</given-names></name><name><surname>Dernovoy</surname><given-names>D</given-names></name><name><surname>Kiryutin</surname><given-names>B</given-names></name><name><surname>Zaslavsky</surname><given-names>L</given-names></name><etal></etal></person-group><year>2008</year><article-title>The Influenza virus resource at the national center for biotechnology information.</article-title><source>J Virol</source><volume>82</volume><fpage>596</fpage><lpage>601</lpage><pub-id pub-id-type="pmid">17942553</pub-id></element-citation></ref><ref id="pone.0038794-Larkin1"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Larkin</surname><given-names>MA</given-names></name><name><surname>Blackshields</surname><given-names>G</given-names></name><name><surname>Brown</surname><given-names>NP</given-names></name><name><surname>Chenna</surname><given-names>R</given-names></name><name><surname>McGettigan</surname><given-names>PA</given-names></name><etal></etal></person-group><year>2007</year><article-title>ClustalW and ClustalX version 2.</article-title><source>Bioinformatics</source><volume>23</volume><fpage>2947</fpage><lpage>2948</lpage><pub-id pub-id-type="pmid">17846036</pub-id></element-citation></ref><ref id="pone.0038794-Tamura1"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tamura</surname><given-names>K</given-names></name><name><surname>Peterson</surname><given-names>D</given-names></name><name><surname>Peterson</surname><given-names>N</given-names></name><name><surname>Stecher</surname><given-names>G</given-names></name><name><surname>Nei</surname><given-names>M</given-names></name><etal></etal></person-group><year>2011</year><article-title>MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.</article-title><source>Mol Biol Evol</source><volume>28</volume><fpage>2731</fpage><lpage>2739</lpage><pub-id pub-id-type="pmid">21546353</pub-id></element-citation></ref><ref id="pone.0038794-DeLano1"><label>31</label><element-citation publication-type="other"><person-group person-group-type="author"><name><surname>DeLano</surname><given-names>WL</given-names></name></person-group><year>2002</year><article-title>Pymol: An open-source molecular graphics tool.</article-title><comment><ext-link ext-link-type="uri" xlink:href="http://www.pymol.org/">http://www.pymol.org/</ext-link></comment></element-citation></ref><ref id="pone.0038794-Ekiert1"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ekiert</surname><given-names>DC</given-names></name><name><surname>Bhabha</surname><given-names>G</given-names></name><name><surname>Elsliger</surname><given-names>MA</given-names></name><name><surname>Friesen</surname><given-names>RH</given-names></name><name><surname>Jongeneelen</surname><given-names>M</given-names></name><etal></etal></person-group><year>2009</year><article-title>Antibody recognition of a highly conserved influenza virus epitope: implications for universal prevention and therapy.</article-title><source>Science</source><volume>324</volume><fpage>246</fpage><lpage>251</lpage><pub-id pub-id-type="pmid">19251591</pub-id></element-citation></ref><ref id="pone.0038794-Humphrey1"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Humphrey</surname><given-names>W</given-names></name><name><surname>Dalke</surname><given-names>A</given-names></name><name><surname>Schulten</surname><given-names>K</given-names></name></person-group><year>1996</year><article-title>VMD-visual molecular dynamics.</article-title><source>J Mol Graph Model</source><volume>1</volume><fpage>33</fpage><lpage>38</lpage></element-citation></ref><ref id="pone.0038794-Andreas1"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Andreas</surname><given-names>BL</given-names></name><name><surname>von der Lieth</surname><given-names>CW</given-names></name></person-group><year>2005</year><article-title>GlyProt: in silico glycosylation of proteins.</article-title><source>Nucleic Acids RES</source><volume>33</volume><fpage>W214</fpage><lpage>W219</lpage><pub-id pub-id-type="pmid">15980456</pub-id></element-citation></ref><ref id="pone.0038794-Phillips1"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Phillips</surname><given-names>JC</given-names></name><name><surname>Braun</surname><given-names>R</given-names></name><name><surname>Wang</surname><given-names>W</given-names></name><name><surname>Gumbart</surname><given-names>J</given-names></name><name><surname>Tajkhorshid</surname><given-names>E</given-names></name><etal></etal></person-group><year>2005</year><article-title>Scalable molecular dynamics with NAMD.</article-title><source>J Comput Chem</source><volume>26</volume><fpage>1781</fpage><lpage>1802</lpage><pub-id pub-id-type="pmid">16222654</pub-id></element-citation></ref><ref id="pone.0038794-MacKerell1"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>MacKerell</surname><given-names>AD</given-names></name><name><surname>Feig</surname><given-names>M</given-names></name><name><surname>Brooks</surname><given-names>CL</given-names></name></person-group><year>2004</year><article-title>Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations.</article-title><source>J Comput Chem</source><volume>25</volume><fpage>1400</fpage><lpage>1415</lpage><pub-id pub-id-type="pmid">15185334</pub-id></element-citation></ref><ref id="pone.0038794-Raman1"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Raman</surname><given-names>EP</given-names></name><name><surname>Guvench</surname><given-names>O</given-names></name><name><surname>Mackerell</surname><given-names>AD</given-names></name></person-group><year>2009</year><article-title>CHARMM additive all-atom force field for glycosidic linkages in carbohydrates involving furanoses.</article-title><source>J Chem Theory Comput</source><volume>5</volume><fpage>23</fpage><lpage>53</lpage><pub-id pub-id-type="pmid">26609816</pub-id></element-citation></ref><ref id="pone.0038794-Kirschner1"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kirschner</surname><given-names>KN</given-names></name><name><surname>Yongye</surname><given-names>AB</given-names></name><name><surname>Tschampel</surname><given-names>SM</given-names></name><name><surname>González-Outeiriño</surname><given-names>J</given-names></name><etal></etal></person-group><year>2008</year><article-title>GLYCAM06: a generalizable biomolecular force field. Carbohydrates.</article-title><source>J Comput Chem</source><volume>29(4)</volume><fpage>622</fpage><lpage>655</lpage><pub-id pub-id-type="pmid">17849372</pub-id></element-citation></ref><ref id="pone.0038794-Vendome1"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vendome</surname><given-names>J</given-names></name><name><surname>Posy</surname><given-names>S</given-names></name><name><surname>Jin</surname><given-names>X</given-names></name><name><surname>Bahna</surname><given-names>F</given-names></name><name><surname>Goran</surname><given-names>Ahlsen</given-names></name><etal></etal></person-group><year>2011</year><article-title>Molecular design principles underlying β-strand swapping in the adhesive dimerization of cadherins.</article-title><source>Nat Struct Mol Biol</source><volume>18</volume><fpage>693</fpage><lpage>700</lpage><pub-id pub-id-type="pmid">21572446</pub-id></element-citation></ref><ref id="pone.0038794-Bohne1"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bohne</surname><given-names>A</given-names></name><name><surname>Lang</surname><given-names>E</given-names></name><name><surname>von der Lieth</surname><given-names>CW</given-names></name></person-group><year>1999</year><article-title>SWEET - WWW-based rapid 3D construction of oligo- and polysaccharides.</article-title><source>Bioinformatics</source><volume>15</volume><fpage>767</fpage><lpage>768</lpage><pub-id pub-id-type="pmid">10498779</pub-id></element-citation></ref><ref id="pone.0038794-Morris1"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Morris</surname><given-names>GM</given-names></name><name><surname>Huey</surname><given-names>R</given-names></name><name><surname>Lindstrom</surname><given-names>W</given-names></name><name><surname>Sanner</surname><given-names>MF</given-names></name><name><surname>Belew</surname><given-names>RK</given-names></name><etal></etal></person-group><year>2009</year><article-title>AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility.</article-title><source>J Comput Chem</source><volume>0</volume><fpage>1</fpage><lpage>7</lpage></element-citation></ref><ref id="pone.0038794-Chandrasekaran1"><label>42</label><element-citation publication-type="other"><person-group person-group-type="author"><name><surname>Chandrasekaran</surname><given-names>A</given-names></name><name><surname>Srinivasan</surname><given-names>A</given-names></name><name><surname>Raman</surname><given-names>R</given-names></name><name><surname>Viswanathan</surname><given-names>K</given-names></name><name><surname>Raguram</surname><given-names>S</given-names></name><etal></etal></person-group><year>2008</year><article-title>Glycan topology determines human adaptation of avian H5N1 virus hemagglutinin.</article-title><source>Nat Biotechnol 26, 107–113</source></element-citation></ref><ref id="pone.0038794-Ducatez1"><label>43</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ducatez</surname><given-names>MF</given-names></name><name><surname>Olinger</surname><given-names>CM</given-names></name><name><surname>Owoade</surname><given-names>AA</given-names></name><name><surname>Tarnagda</surname><given-names>Z</given-names></name><name><surname>Tahita</surname><given-names>MC</given-names></name><etal></etal></person-group><year>2007</year><article-title>Molecular and antigenic evolution and geographical spread of H5N1 highly pathogenic avian influenza viruses in western Africa.</article-title><source>J Gen Virol</source><volume>88</volume><fpage>2297</fpage><lpage>2306</lpage><pub-id pub-id-type="pmid">17622635</pub-id></element-citation></ref><ref id="pone.0038794-Chen2"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>R</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name></person-group><year>2006</year><article-title>Avian influenza virus exhibits rapid evolutionary dynamics.</article-title><source>Mol Biol Evol</source><volume>23</volume><fpage>2336</fpage><lpage>2341</lpage><pub-id pub-id-type="pmid">16945980</pub-id></element-citation></ref><ref id="pone.0038794-WHOOIEFAO2"><label>45</label><element-citation publication-type="journal"><collab>WHO/OIE/FAO H5N1 Evolution Working Group</collab><year>2008</year><article-title>Toward a unified nomenclature system for highly pathogenic avian influenza virus (H5N1).</article-title><source>Emerg Infect Dis</source><volume>14(7)</volume><fpage>e1</fpage></element-citation></ref><ref id="pone.0038794-WHOOIEFAO3"><label>46</label><element-citation publication-type="journal"><collab>WHO/OIE/FAO H5N1 Evolution Working Group</collab><year>2012</year><article-title>Continued evolution of highly pathogenic avian influenza A (H5N1): updated nomenclature.</article-title><source>Influenza and Other Respiratory Viruses</source><volume>6</volume><fpage>1</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">22035148</pub-id></element-citation></ref><ref id="pone.0038794-Ha1"><label>47</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ha</surname><given-names>Y</given-names></name><name><surname>Stevens</surname><given-names>DJ</given-names></name><name><surname>Skehel</surname><given-names>JJ</given-names></name><name><surname>Wiley</surname><given-names>DC</given-names></name></person-group><year>2001</year><article-title>X-ray structures of H5 avian and H9 swine influenza virus hemagglutinins bound to avian and human receptor analogs.</article-title><source>PNAS</source><volume>98</volume><fpage>11181</fpage><lpage>11186</lpage><pub-id pub-id-type="pmid">11562490</pub-id></element-citation></ref><ref id="pone.0038794-Stevens2"><label>48</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stevens</surname><given-names>J</given-names></name><name><surname>Blixt</surname><given-names>O</given-names></name><name><surname>Tumpey</surname><given-names>TM</given-names></name><name><surname>Taubenberger</surname><given-names>JK</given-names></name><name><surname>Paulson</surname><given-names>JC</given-names></name><etal></etal></person-group><year>2007</year><article-title>Structure and receptor specificity of the hemagglutinin from an H5N1 influenza virus.</article-title><source>Science</source><volume>312</volume><fpage>404</fpage><lpage>410</lpage></element-citation></ref><ref id="pone.0038794-Yamada1"><label>49</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yamada</surname><given-names>S</given-names></name><name><surname>Suzuki</surname><given-names>Y</given-names></name><name><surname>Suzuki</surname><given-names>T</given-names></name><name><surname>Le</surname><given-names>MQ</given-names></name><name><surname>Nidom</surname><given-names>CA</given-names></name><etal></etal></person-group><year>2006</year><article-title>Haemagglutinin mutations responsible for the binding of H5N1 influenza A viruses to human-type receptors.</article-title><source>Nature</source><volume>444</volume><fpage>378</fpage><lpage>382</lpage><pub-id pub-id-type="pmid">17108965</pub-id></element-citation></ref><ref id="pone.0038794-Sui1"><label>50</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sui</surname><given-names>J</given-names></name><name><surname>Hwang</surname><given-names>WC</given-names></name><name><surname>Perez</surname><given-names>S</given-names></name><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Aird</surname><given-names>D</given-names></name><etal></etal></person-group><year>2009</year><article-title>Structural and functional bases for broad-spectrum neutralization of avian and human influenza A viruses.</article-title><source>Nat Struct Mol Biol</source><volume>16</volume><fpage>265</fpage><lpage>273</lpage><pub-id pub-id-type="pmid">19234466</pub-id></element-citation></ref><ref id="pone.0038794-Bohne2"><label>51</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bohne</surname><given-names>A</given-names></name><name><surname>Lieth</surname><given-names>CW</given-names></name></person-group><year>2002</year><article-title>Glycosylation of proteins: a computer based method for the rapid exploration of conformational space of N-glycans.</article-title><source>Pacific Symposium on Biocomputing</source><volume>7</volume><fpage>285</fpage><lpage>292</lpage></element-citation></ref><ref id="pone.0038794-Moustakas1"><label>52</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moustakas</surname><given-names>DT</given-names></name><name><surname>Lang</surname><given-names>PT</given-names></name><name><surname>Pegg</surname><given-names>S</given-names></name><name><surname>Pettersen</surname><given-names>ET</given-names></name><name><surname>Kuntz</surname><given-names>ID</given-names></name><etal></etal></person-group><year>2006</year><article-title>Development and validation of a modular, extensible docking program: DOCK 5.</article-title><source>J Comput Aided Mol Des</source><volume>20</volume><fpage>601</fpage><lpage>609</lpage><pub-id pub-id-type="pmid">17149653</pub-id></element-citation></ref><ref id="pone.0038794-Nicholls1"><label>53</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nicholls</surname><given-names>JM</given-names></name><name><surname>Chan</surname><given-names>R</given-names></name><name><surname>Russell</surname><given-names>RJ</given-names></name><name><surname>Air</surname><given-names>GM</given-names></name><name><surname>Peiris</surname><given-names>M</given-names></name></person-group><year>2008</year><article-title>Evolving complexities of influenza virus and its receptors.</article-title><source>Trends Microbiol</source><volume>16</volume><fpage>149</fpage><lpage>157</lpage><pub-id pub-id-type="pmid">18375125</pub-id></element-citation></ref><ref id="pone.0038794-Leo1"><label>54</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Leo</surname><given-names>F</given-names></name><name><surname>Asakuma</surname><given-names>S</given-names></name><name><surname>Fukuda</surname><given-names>K</given-names></name><name><surname>Senda</surname><given-names>A</given-names></name><name><surname>Urashima</surname><given-names>T</given-names></name></person-group><year>2010</year><article-title>Determination of sialyl and neutral oligosaccharides levels in transition and mature milks of Samoan women, using anthranilic derivation followed by reverse phase high performance liquid chromatography.</article-title><source>Biosci Biotech Bioch</source><volume>74</volume><fpage>298</fpage><lpage>303</lpage></element-citation></ref><ref id="pone.0038794-Kaverin1"><label>55</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kaverin</surname><given-names>NV</given-names></name><name><surname>Rudneva</surname><given-names>IA</given-names></name><name><surname>Llyushina</surname><given-names>NA</given-names></name><name><surname>Varich</surname><given-names>NL</given-names></name><name><surname>Lipatov</surname><given-names>AS</given-names></name><etal></etal></person-group><year>2002</year><article-title>Structure of antigenic sites on the haemagglutinin molecule of H5 avian influenza virus and phenotypic variation of escape mutants.</article-title><source>J Gen Virol</source><volume>83</volume><fpage>2497</fpage><lpage>2505</lpage><pub-id pub-id-type="pmid">12237433</pub-id></element-citation></ref><ref id="pone.0038794-Salahuddin1"><label>56</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salahuddin</surname><given-names>P</given-names></name><name><surname>Khan</surname><given-names>AU</given-names></name></person-group><year>2009</year><article-title>Identification of sequence mutations affecting hemagglutinin specificity to sialic acid receptor in influenza A virus subtypes.</article-title><source>Bioinformatics</source><volume>4</volume><fpage>30</fpage><lpage>35</lpage></element-citation></ref><ref id="pone.0038794-Turner1"><label>57</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turner</surname><given-names>BG</given-names></name><name><surname>Summers</surname><given-names>MF</given-names></name></person-group><year>1999</year><article-title>Structural Biology of HIV.</article-title><source>J Mol Biol</source><volume>285</volume><fpage>1</fpage><lpage>32</lpage><pub-id pub-id-type="pmid">9878383</pub-id></element-citation></ref><ref id="pone.0038794-Rota1"><label>58</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rota</surname><given-names>PA</given-names></name><name><surname>Oberste</surname><given-names>MS</given-names></name><name><surname>Monroe</surname><given-names>SS</given-names></name><name><surname>Nix</surname><given-names>WA</given-names></name><name><surname>Campagnoli</surname><given-names>R</given-names></name><etal></etal></person-group><year>2003</year><article-title>Characterization of a Novel Coronavirus Associated with Severe Acute Respiratory Syndrome.</article-title><source>Science</source><volume>300</volume><fpage>1394</fpage><lpage>1399</lpage><pub-id pub-id-type="pmid">12730500</pub-id></element-citation></ref><ref id="pone.0038794-Neu1"><label>59</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Neu</surname><given-names>U</given-names></name><name><surname>Woellner</surname><given-names>K</given-names></name><name><surname>Gauglitz</surname><given-names>G</given-names></name><name><surname>Stehle</surname><given-names>T</given-names></name></person-group><year>2008</year><article-title>Structural basis of GM1 ganglioside recognition by simian virus40.</article-title><source>PNAS</source><volume>105</volume><fpage>5219</fpage><lpage>5224</lpage><pub-id pub-id-type="pmid">18353982</pub-id></element-citation></ref><ref id="pone.0038794-Tsai1"><label>60</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tsai</surname><given-names>B</given-names></name><name><surname>Gilbert</surname><given-names>JM</given-names></name><name><surname>Stehle</surname><given-names>T</given-names></name><name><surname>Lencer</surname><given-names>W</given-names></name><name><surname>Benjamin</surname><given-names>TL</given-names></name><etal></etal></person-group><year>2003</year><article-title>Gangliosides are receptors for murine polyoma virus and SV40.</article-title><source>Embo J</source><volume>22</volume><fpage>4346</fpage><lpage>4355</lpage><pub-id pub-id-type="pmid">12941687</pub-id></element-citation></ref><ref id="pone.0038794-Lin1"><label>61</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname><given-names>T</given-names></name><name><surname>Wang</surname><given-names>G</given-names></name><name><surname>Li</surname><given-names>A</given-names></name><name><surname>Zhang</surname><given-names>Q</given-names></name><name><surname>Wu</surname><given-names>C</given-names></name><etal></etal></person-group><year>2009</year><article-title>The hemagglutinin structure of an avian H1N1 influenza A virus.</article-title><source>Virology</source><volume>392</volume><fpage>73</fpage><lpage>81</lpage><pub-id pub-id-type="pmid">19628241</pub-id></element-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/H2N2V1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000B42 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000B42 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= H2N2V1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:3375263
|texte= Two Glycosylation Sites in H5N1 Influenza Virus Hemagglutinin That Affect Binding Preference by Computer-Based Analysis
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:22719948" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a H2N2V1
| This area was generated with Dilib version V0.6.33. |  |