BACKGROUND
avian influenza h5n <dig> virus uses its hemagglutinin  protein to bind with a host receptor before entering the cell.
this protein binds avidly to the avian-type receptor.
however, a major health concern is that ha mutation could alter the binding preference to that of human receptor, which could occur before the virus is completely adapted to its new host  <cit> .
the incidences of human infection by h5n <dig> virus and the spectrum of h5n <dig> mutations are increasing  <cit> .
some of the mutated viruses could potentially infect humans and be spread person-to-person causing an outbreak  <cit> .
the host cell selectivity of influenza a viruses is mediated by the interaction of particular viral ha variants to different host cell receptor types.
the cell receptor that is bound by ha is a penta-saccharide chain.
the first sugar unit is sialic acid , followed by galactose , n-acetylglucosamine , gal, and glucose  units.
however, the available x-ray crystal structures of the host cell receptor are in the tri-saccharide form as shown in figure  <dig>  precluding accurate simulation of the full-length receptor.
two types of receptor are bound by influenza virus: the first type has the α-linkage between the first two units to form siaαgal glycosides.
the other receptor type contains an α-linkage to form siaαgal.
in avian viruses, the preferred ha receptors are of the siaαgal type, while most human viruses interact with siaαgal glycoside receptors.
normally, the avian influenza virus h5n <dig> infects birds rather than humans or other mammalian hosts because their ha binds better to the avian-type siaαgal glycoside receptor  <cit> .
despite the number of ha variants that have been reported in the protein database  <cit> , we do not know enough about the binding mechanism to predict which ha variants can bind efficiently to the human receptor.
the most reliable source of information to study the binding mechanism is from x-ray crystal structure.
however, only a small number of h <dig> ha-receptor co-crystals are available in public databases  <cit> .
another approach to study which ha variants can bind preferentially to human receptor is by binding assay experiments which systematically screen interactions between ha variants and receptor analogs  <cit> .
nevertheless, to produce and screen many proteins in order to search for human-receptor binding has is impractical.
alternatively, insight into the binding mechanism can be observed by ha protein-receptor simulation to predict the binding potential between different ha variants and human receptor.
by this means, one can effectively screen and assign priority to a small number of ha variants for further in vitro experiments.
surveillance of h <dig> ha mutations over a number of years has led to the discovery of viruses with changes in receptor binding selectivity caused by mutation of the receptor binding site.
the h <dig> ha x-ray structure from a/duck/singapore/3/ <dig>  shows preferential binding to siaαgal receptor, owing to q222l and g224s mutations  <cit> .
a sing- <dig> descendent that infected a human in  <dig>  a/thailand/1/ <dig> , has a single mutation in the ha binding pocket .
this variant, however, still binds preferentially with the avian receptor  <cit> .
a kan- <dig> derived strain with further ha mutations l129v and a134v is considered a quasi-species which exhibited higher selectivity toward human cell type receptor  <cit> .
it should be noted that the mutated amino acids considered in the study are not the main binding residues and the mutations do not greatly alter the resulting complex structure.
however, some of these mutations do indeed change the receptor type selectivity by contributing to changes in binding mechanism.
this modification in host type preference could be caused by differential binding of mutually exclusive conformations of the cell receptor in different ha binding environments.
it is hypothesized, that by comparing the conformations of different receptor types bound to different ha binding sites with known binding selectivity, we could distinguish different bound conformations.
a recent study by xu and coworkers  <cit>  has shown that receptor binding preference of different influenza viral types can be modelled by simulation measuring receptor torsion angles.
in this study, it is shown how the information from the available structures of ha complexed with avian and human receptors can be used to predict the ha binding selectivity from different ha variants in receptor-based conformational analysis by measuring a single torsion angle during molecular dynamic simulation.
furthermore, the predictions of receptor preference agree with the available in vitro binding data.
materials and methods
crystallographic datasets and ha variants used for simulation
x-ray crystallographic datasets of ha variants complexed with tri-saccharide receptor analogs were obtained from the pdb .
in the 1jso dataset, the electron density is ill-defined for the galβglcnac sugar residues.
therefore, to allow accurate comparisons between the different crystallographic templates in md simulations, the structures of the equivalent sugar residues from the receptor in the 1rvz structural template  <cit> , the closest receptor structure available, were inserted into the 1jso template.
all glycosides were terminated with a methoxy group and were used as the input for molecular dynamics simulations.
to prepare ha variants for md simulation for which no crystallographic data are available, homology modelling was performed by three-dimensional alignment with the most similar structural template  variants) using the swiss-model server  <cit> .
a number of ha strains of known host selectivity were used to test whether siaαgal would bind to ha with avian receptor preference in trans, and whether siaαgal would bind to ha with human receptor preference in cis.
the md simulations on seven known ha variants with avian and human receptor types were performed: sing- <dig>  puerto- <dig>  kan- <dig>  mutated kan- <dig>  and three mutated sing- <dig> systems.
among these systems, according to the available binding assay data, two of them, sing- <dig> and kan- <dig>  bind siaαgal stronger than siaαgal  <cit> .
one of the mutated sing- <dig> strains  is known to bind siaαgal stronger, suggesting that these mutations could be important for changing the receptor preference from avian to human.
the other strains except puerto- <dig> can equally bind to both types of the receptor.
for puerto- <dig> there are no experimental data from binding assay, but we assumed strong human receptor binding for this is the strain since it caused a human pandemic  <cit> .
our prediction scheme was tested on two sets of data.
the first set comprised four kan- <dig> ha variants each with single mutations m226t, k189g, k218s, and l190p obtained as quasi-species from the source where kan- <dig> was found  <cit> .
the reason for testing these variants was to determine whether these strains are potentially harmful, i.e.
stronger binding to human receptor than kan- <dig>  the second set contained three mutated puerto- <dig> ha variants, q222l, g224s, and q222l/g224s .
these variants were chosen to determine what effect these mutations had on a different ha type  to sing- <dig> .
molecular dynamic simulations
to permit comparison with previous experimental results  <cit> , all the simulations were done using the same glycam <dig>  <cit>  parameters for receptor and amber  <dig> force field for protein.
initial structure was first solvated using the tip5p water model  <cit>  in the truncated octahedron box.
energy minimization was then performed to relieve bad contacts caused by unreasonable distances in the structure by keeping the protein and receptor restrained.
the whole system was relaxed at  <dig> k with  <dig> Å non-bonded cutoff.
the temperature of the system was then set to  <dig> k and equilibrated for  <dig> picoseconds with weak restraint on both receptor and protein, where bonds involving hydrogen are constrained using the shake algorithm  <cit> .
torsion angles  were monitored throughout the simulation to determine the conformation of the glycosides.
the Φ angle is defined as the angle between the o6-c <dig> bond of sia and the glycosidic bond of gal units .
note that the Φ angle defined by xu and colleagues  <cit>  refers to a different plane of rotation of the receptor.
the Φ torsion angle described in this study was not considered in xu et al.
to determine the binding preference, the Φ angle was monitored in order to reflect the receptor type selectivity.
the  <dig> nanosecond production run was performed at a constant temperature and pressure with  <dig>  picoseconds time step  using the sander module in the amber <dig> program  <cit> .
the structures stabilized after  <dig>  nanoseconds as shown in figure  <dig>  the highest degree of fluctuation was observed for residues located at the terminal chains of the structures.
these residues were inserted in ha chain b and were left unrestrained during the simulations; however, they did not disturb the binding site  between residues  <dig> and  <dig> was less than  <dig>  Å).
the utility programs, xmgrace  <cit> , and vmd  <cit>  were used to visualize and render all the figures presented in this paper.
RESULTS
data from available co-crystal structures of cell receptor analog and ha were used to establish relationships between bound receptor conformation and host type preferences.
comparison between these structures revealed that the Φ torsion angle, has different values for siaαgal and siaαgal binding .
the observed values show siaαgal has a Φ angle approximately - <dig> degrees in the h <dig> binding pocket, i.e.
the receptor is in trans conformation.
meanwhile, siaαgal exhibits a Φ angle of approximately + <dig> degrees, i.e.
the receptor is in cis conformation  <cit> .
in other words, siaαgal seems to have an optimal binding geometry when the gal and sia units are bound in the trans conformation, while for the siaαgal they are instead bound in the cis conformation.
according to this, we hypothesized that in solution, both conformations are in equilibrium, as suggested by the available crystallographic data  <cit> .
upon binding of ha to the receptor, it is hypothesized that one conformation is favoured; thus, binding drives the equilibrium to this conformation without any molecular readjustment or thermodynamic cost.
in order to use the relationships between torsion angles and binding preferences to explain host selectivity of the unknown influenza virus structures, homology modelling and molecular dynamics  simulations were employed.
during each md simulation, Φ was monitored and interpreted in terms of the binding preference.
md simulations were performed with different combinations of avian and human-type receptor analogs with ha variants.
the duration of time in which the receptor analog spent in the cis or trans conformation when bound to each ha type varied according to the type of glycosidic linkage and the number of sugar units.
figure  <dig> presents plots showing different cis/trans conformations.
for the binding studies of avian receptor analogs with α-linkage, the glycosidic linkage in trans conformation was observed during most of the simulation in all ha systems tested , albeit with some transient fluctuations to the cis conformation for the q222l, g224s, and q222l/g224s sing- <dig> ha variants.
in the simulation of human receptor analogs with α-linkage, the trans configuration  was observed in the majority during the simulation with some transient fluctuations to the cis conformation when bound to sing- <dig>  sing- <dig> mutants q222l, g224s and kan- <dig> systems,.
in contrast, the cis conformation was observed  for the majority of the simulation time when the analogs bound with puerto- <dig> human influenza, l129v/a134v kan- <dig>  and q222l/g224s sing- <dig> has.
discussion
from the experimental results, it can be concluded that each of the h <dig> ha variants exhibited different binding behaviours to different receptor analogs.
the available x-ray structures contain only the tri-saccharide part of receptor analogs bound to the binding site, and the last two sugar units are missing or unsolved  <cit> .
nonetheless, the full-length receptor structures can be constructed using molecular modelling software  <cit> .
in this study, modelling was restricted to the tri-saccharide receptor system to minimize the error in the simulations.
in our previous work  <cit> , kan- <dig> was predicted to bind fairly well to human receptor for a significantly long period; however, this prediction is at variance with in vitro experimental data .
from the above experimental results, the tri-saccharide simulation model can better estimate the host type binding preference of different h <dig> ha variants than di-saccharide.
in vitro
puerto-34
sing-97
kan-1
na - no data available, +++ most preferred, ++ moderately preferred, + least preferred
* - siaαgalβglcnac not siaαgalβglcnac
the crystallographic data in ha et al  <cit>  showed that the typical h <dig> ha binds preferentially to avian receptor with an α-linkage in the trans conformation, whereas the typical h <dig> ha binds preferentially to human receptor with an α-linkage in the cis conformation.
the simulations presented here predicted that:
sing- <dig> and kan- <dig> bound better to α than to α, since the observed predominant conformation of receptor was trans for both receptor types
l129v/a134v kan- <dig> could bind to both receptor types, since the trans conformation was predominant for α binding and the cis conformation for α binding
q222l, g224s, q222l/g224s sing- <dig> mutants appear to have a weaker preference for α than non-mutated sing- <dig> because fluctuations from the trans to cis conformation were observed.
the α simulation of mutated sing- <dig> implied that the q222l/g224l variant had markedly greater binding affinity toward human cell receptor as it bound in cis with α all the time.
human virus puerto- <dig> ha bound preferentially with α, since the observed conformation of receptor was in cis configuration.
prediction of the relative binding selectivity
based on the ha binding conformation preferences from md simulation , predictions of relative binding selectivity  can be made as follows.
the selectivity toward siaαgal binding was similar among the three ha variants puerto- <dig>  sing <dig> and kan- <dig>  the order of selectivity toward siaαgal binding was puerto- <dig> > l129v/a134v kan- <dig> ≅ q222l/g224s sing- <dig> > sing- <dig> ha ≅ kan- <dig>  these tendencies were in good agreement with the in vitro binding assays  <cit> , in terms of order of preference.
therefore, the duration of the cis conformation during the simulation may be correlated with the selectivity of the docked ha.
according to our prediction scheme, the l129v/a134v kan- <dig> variants, i.e.
mutations m226t, k189g, k218s, and l190p, may have increased the selectivity slightly toward human receptor, since the receptor was present in cis conformation for some of the simulation .
the results from the three mutated puerto- <dig> systems also showed some changes in their binding behaviors compared to non-mutated puerto- <dig> .
the two single mutations, q222l and g224s cause a loss in human receptor affinity as shown by fluctuations to the trans-conformation, while the double mutation q222l/g224s maintained its preference for human receptor as the bound α glycosides were in the cis conformation.
for the α receptor, all the puerto- <dig> systems except for g224s appear to interact weakly since the receptor is in the cis-conformation.
the trans conformation is observed for g224s, although the average torsion angle is increased from - <dig> to - <dig> degrees, suggesting that it may not be optimal for binding.
the results show that mutations at residues  <dig> and  <dig> have minor impact on host preference for puerto- <dig> ha
CONCLUSIONS
we have shown that our cis-trans conformational analysis scheme could predict the host type selectivity of ha variants.
our cis-trans conformation hypothesis also worked well under another ha system where it predicted the cis conformation and revealed the similar mechanism in puerto- <dig> simulation.
the binding patterns and mechanisms of the adopted receptor model, siaαgal and siaαgal, to wild-type and mutated kan- <dig> ha, and sing- <dig> ha were proposed.
the results could be used to explain why the l129v/a134v kan- <dig> and q222l/g224s sing- <dig> could bind better to human receptor analog in in vitro assays.
the underlying proposed mechanism that made h <dig> bind to human host without mutation at residue  <dig> or  <dig> involved the interaction between residue  <dig> side-chain and gln <dig>  it is proposed that mutations change the ha binding preference from siaαgal to siaαgal.
our study also suggested that even mutations outside of key binding residues  <cit> , e.g.
residue  <dig> or  <dig>  have consequences on altering receptor type and should not be ignored.
furthermore, our procedure is useful for predicting host type, which can be tested by in vitro binding assays.
competing interests
the authors declare that they have no competing interests.
authors' contributions
nj, dc and wm performed modeling, md simulations and generated figures.
st, pjs and cs helped with data analysis and drafting of the manuscript.
st and cs obtained funding.
all authors read and approved the final manuscript.
note
other papers from the meeting have been published as part of bmc bioinformatics volume  <dig> supplement  <dig>  2009: eighth international conference on bioinformatics : bioinformatics, available online at http://www.biomedcentral.com/1471-2105/10?issue=s <dig>

 avian influenza h5n virus uses its hemagglutinin protein to bind with a host receptor before entering the cell. this protein binds avidly to the avian-type receptor. however, a major health concern is that ha mutation could alter the binding preference to that of human receptor, which could occur before the virus is completely adapted to its new host. the incidences of human infection by h5n virus and the spectrum of h5n mutations are increasing. some of the mutated viruses could potentially infect humans and be spread person-to-person causing an outbreak. the host cell selectivity of influenza a viruses is mediated by the interaction of particular viral ha variants to different host cell receptor types. the cell receptor that is bound by ha is a penta-saccharide chain. the first sugar unit is sialic acid, followed by galactose, n-acetylglucosamine, gal, and glucose units. however, the available x-ray crystal structures of the host cell receptor are in the tri-saccharide form as shown in figure precluding accurate simulation of the full-length receptor. two types of receptor are bound by influenza virus: the first type has the α-linkage between the first two units to form siaαgal glycosides. the other receptor type contains an α-linkage to form siaαgal. in avian viruses, the preferred ha receptors are of the siaαgal type, while most human viruses interact with siaαgal glycoside receptors. normally, the avian influenza virus h5n infects birds rather than humans or other mammalian hosts because their ha binds better to the avian-type siaαgal glycoside receptor. despite the number of ha variants that have been reported in the protein database, we do not know enough about the binding mechanism to predict which ha variants can bind efficiently to the human receptor. the most reliable source of information to study the binding mechanism is from x-ray crystal structure. however, only a small number of h ha-receptor co-crystals are available in public databases. another approach to study which ha variants can bind preferentially to human receptor is by binding assay experiments which systematically screen interactions between ha variants and receptor analogs. nevertheless, to produce and screen many proteins in order to search for human-receptor binding has is impractical. alternatively, insight into the binding mechanism can be observed by ha protein-receptor simulation to predict the binding potential between different ha variants and human receptor. by this means, one can effectively screen and assign priority to a small number of ha variants for further in vitro experiments. surveillance of h ha mutations over a number of years has led to the discovery of viruses with changes in receptor binding selectivity caused by mutation of the receptor binding site. the h ha x-ray structure from a/duck/singapore/3/ shows preferential binding to siaαgal receptor, owing to q222l and g224s mutations. a sing- descendent that infected a human in a/thailand/1/, has a single mutation in the ha binding pocket. this variant, however, still binds preferentially with the avian receptor. a kan- derived strain with further ha mutations l129v and a134v is considered a quasi-species which exhibited higher selectivity toward human cell type receptor. it should be noted that the mutated amino acids considered in the study are not the main binding residues and the mutations do not greatly alter the resulting complex structure. however, some of these mutations do indeed change the receptor type selectivity by contributing to changes in binding mechanism. this modification in host type preference could be caused by differential binding of mutually exclusive conformations of the cell receptor in different ha binding environments. it is hypothesized, that by comparing the conformations of different receptor types bound to different ha binding sites with known binding selectivity, we could distinguish different bound conformations. a recent study by xu and coworkers has shown that receptor binding preference of different influenza viral types can be modelled by simulation measuring receptor torsion angles. in this study, it is shown how the information from the available structures of ha complexed with avian and human receptors can be used to predict the ha binding selectivity from different ha variants in receptor-based conformational analysis by measuring a single torsion angle during molecular dynamic simulation. furthermore, the predictions of receptor preference agree with the available in vitro binding data. materials and methods crystallographic datasets and ha variants used for simulation x-ray crystallographic datasets of ha variants complexed with tri-saccharide receptor analogs were obtained from the pdb. in the 1jso dataset, the electron density is ill-defined for the galβglcnac sugar residues. therefore, to allow accurate comparisons between the different crystallographic templates in md simulations, the structures of the equivalent sugar residues from the receptor in the 1rvz structural template, the closest receptor structure available, were inserted into the 1jso template. all glycosides were terminated with a methoxy group and were used as the input for molecular dynamics simulations. to prepare ha variants for md simulation for which no crystallographic data are available, homology modelling was performed by three-dimensional alignment with the most similar structural template variants) using the swiss-model server. a number of ha strains of known host selectivity were used to test whether siaαgal would bind to ha with avian receptor preference in trans, and whether siaαgal would bind to ha with human receptor preference in cis. the md simulations on seven known ha variants with avian and human receptor types were performed: sing- puerto- kan- mutated kan- and three mutated sing- systems. among these systems, according to the available binding assay data, two of them, sing- and kan- bind siaαgal stronger than siaαgal. one of the mutated sing- strains is known to bind siaαgal stronger, suggesting that these mutations could be important for changing the receptor preference from avian to human. the other strains except puerto- can equally bind to both types of the receptor. for puerto- there are no experimental data from binding assay, but we assumed strong human receptor binding for this is the strain since it caused a human pandemic. our prediction scheme was tested on two sets of data. the first set comprised four kan- ha variants each with single mutations m226t, k189g, k218s, and l190p obtained as quasi-species from the source where kan- was found. the reason for testing these variants was to determine whether these strains are potentially harmful, i.e. stronger binding to human receptor than kan- the second set contained three mutated puerto- ha variants, q222l, g224s, and q222l/g224s. these variants were chosen to determine what effect these mutations had on a different ha type to sing-. molecular dynamic simulations to permit comparison with previous experimental results, all the simulations were done using the same glycam parameters for receptor and amber force field for protein. initial structure was first solvated using the tip5p water model in the truncated octahedron box. energy minimization was then performed to relieve bad contacts caused by unreasonable distances in the structure by keeping the protein and receptor restrained. the whole system was relaxed at k with Å non-bonded cutoff. the temperature of the system was then set to k and equilibrated for picoseconds with weak restraint on both receptor and protein, where bonds involving hydrogen are constrained using the shake algorithm. torsion angles were monitored throughout the simulation to determine the conformation of the glycosides. the Φ angle is defined as the angle between the o6-c bond of sia and the glycosidic bond of gal units. note that the Φ angle defined by xu and colleagues refers to a different plane of rotation of the receptor. the Φ torsion angle described in this study was not considered in xu et al. to determine the binding preference, the Φ angle was monitored in order to reflect the receptor type selectivity. the nanosecond production run was performed at a constant temperature and pressure with picoseconds time step using the sander module in the amber program. the structures stabilized after nanoseconds as shown in figure the highest degree of fluctuation was observed for residues located at the terminal chains of the structures. these residues were inserted in ha chain b and were left unrestrained during the simulations; however, they did not disturb the binding site between residues and was less than Å). the utility programs, xmgrace, and vmd were used to visualize and render all the figures presented in this paper.  data from available co-crystal structures of cell receptor analog and ha were used to establish relationships between bound receptor conformation and host type preferences. comparison between these structures revealed that the Φ torsion angle, has different values for siaαgal and siaαgal binding. the observed values show siaαgal has a Φ angle approximately - degrees in the h binding pocket, i.e. the receptor is in trans conformation. meanwhile, siaαgal exhibits a Φ angle of approximately + degrees, i.e. the receptor is in cis conformation. in other words, siaαgal seems to have an optimal binding geometry when the gal and sia units are bound in the trans conformation, while for the siaαgal they are instead bound in the cis conformation. according to this, we hypothesized that in solution, both conformations are in equilibrium, as suggested by the available crystallographic data. upon binding of ha to the receptor, it is hypothesized that one conformation is favoured; thus, binding drives the equilibrium to this conformation without any molecular readjustment or thermodynamic cost. in order to use the relationships between torsion angles and binding preferences to explain host selectivity of the unknown influenza virus structures, homology modelling and molecular dynamics simulations were employed. during each md simulation, Φ was monitored and interpreted in terms of the binding preference. md simulations were performed with different combinations of avian and human-type receptor analogs with ha variants. the duration of time in which the receptor analog spent in the cis or trans conformation when bound to each ha type varied according to the type of glycosidic linkage and the number of sugar units. figure presents plots showing different cis/trans conformations. for the binding studies of avian receptor analogs with α-linkage, the glycosidic linkage in trans conformation was observed during most of the simulation in all ha systems tested, albeit with some transient fluctuations to the cis conformation for the q222l, g224s, and q222l/g224s sing- ha variants. in the simulation of human receptor analogs with α-linkage, the trans configuration was observed in the majority during the simulation with some transient fluctuations to the cis conformation when bound to sing- sing- mutants q222l, g224s and kan- systems,. in contrast, the cis conformation was observed for the majority of the simulation time when the analogs bound with puerto- human influenza, l129v/a134v kan- and q222l/g224s sing- has. discussion from the experimental results, it can be concluded that each of the h ha variants exhibited different binding behaviours to different receptor analogs. the available x-ray structures contain only the tri-saccharide part of receptor analogs bound to the binding site, and the last two sugar units are missing or unsolved. nonetheless, the full-length receptor structures can be constructed using molecular modelling software. in this study, modelling was restricted to the tri-saccharide receptor system to minimize the error in the simulations. in our previous work, kan- was predicted to bind fairly well to human receptor for a significantly long period; however, this prediction is at variance with in vitro experimental data. from the above experimental results, the tri-saccharide simulation model can better estimate the host type binding preference of different h ha variants than di-saccharide. in vitro puerto-34 sing-97 kan-1 na - no data available, +++ most preferred, ++ moderately preferred, + least preferred * - siaαgalβglcnac not siaαgalβglcnac the crystallographic data in ha et al showed that the typical h ha binds preferentially to avian receptor with an α-linkage in the trans conformation, whereas the typical h ha binds preferentially to human receptor with an α-linkage in the cis conformation. the simulations presented here predicted that: sing- and kan- bound better to α than to α, since the observed predominant conformation of receptor was trans for both receptor types l129v/a134v kan- could bind to both receptor types, since the trans conformation was predominant for α binding and the cis conformation for α binding q222l, g224s, q222l/g224s sing- mutants appear to have a weaker preference for α than non-mutated sing- because fluctuations from the trans to cis conformation were observed. the α simulation of mutated sing- implied that the q222l/g224l variant had markedly greater binding affinity toward human cell receptor as it bound in cis with α all the time. human virus puerto- ha bound preferentially with α, since the observed conformation of receptor was in cis configuration. prediction of the relative binding selectivity based on the ha binding conformation preferences from md simulation, predictions of relative binding selectivity can be made as follows. the selectivity toward siaαgal binding was similar among the three ha variants puerto- sing and kan- the order of selectivity toward siaαgal binding was puerto- > l129v/a134v kan- ≅ q222l/g224s sing- > sing- ha ≅ kan- these tendencies were in good agreement with the in vitro binding assays, in terms of order of preference. therefore, the duration of the cis conformation during the simulation may be correlated with the selectivity of the docked ha. according to our prediction scheme, the l129v/a134v kan- variants, i.e. mutations m226t, k189g, k218s, and l190p, may have increased the selectivity slightly toward human receptor, since the receptor was present in cis conformation for some of the simulation. the results from the three mutated puerto- systems also showed some changes in their binding behaviors compared to non-mutated puerto-. the two single mutations, q222l and g224s cause a loss in human receptor affinity as shown by fluctuations to the trans-conformation, while the double mutation q222l/g224s maintained its preference for human receptor as the bound α glycosides were in the cis conformation. for the α receptor, all the puerto- systems except for g224s appear to interact weakly since the receptor is in the cis-conformation. the trans conformation is observed for g224s, although the average torsion angle is increased from - to - degrees, suggesting that it may not be optimal for binding. the results show that mutations at residues and have minor impact on host preference for puerto- ha  we have shown that our cis-trans conformational analysis scheme could predict the host type selectivity of ha variants. our cis-trans conformation hypothesis also worked well under another ha system where it predicted the cis conformation and revealed the similar mechanism in puerto- simulation. the binding patterns and mechanisms of the adopted receptor model, siaαgal and siaαgal, to wild-type and mutated kan- ha, and sing- ha were proposed. the results could be used to explain why the l129v/a134v kan- and q222l/g224s sing- could bind better to human receptor analog in in vitro assays. the underlying proposed mechanism that made h bind to human host without mutation at residue or involved the interaction between residue side-chain and gln it is proposed that mutations change the ha binding preference from siaαgal to siaαgal. our study also suggested that even mutations outside of key binding residues, e.g. residue or have consequences on altering receptor type and should not be ignored. furthermore, our procedure is useful for predicting host type, which can be tested by in vitro binding assays. competing interests the authors declare that they have no competing interests. authors' contributions nj, dc and wm performed modeling, md simulations and generated figures. st, pjs and cs helped with data analysis and drafting of the manuscript. st and cs obtained funding. all authors read and approved the final manuscript. note other papers from the meeting have been published as part of bmc bioinformatics volume supplement 2009: eighth international conference on bioinformatics: bioinformatics, available online at http://www.biomedcentral.com/1471-2105/10?issue=s 
the avian receptor was represented by siaαgal substructure and human receptor by siaαgal. the aim of this study is to further explore whether computational simulation could be used as a prediction tool for host type selectivity in emerging variants. 7– september asia pacific bioinformatics network eighth international conference on bioinformatics singapore the simulation results showed all hemagglutinin proteins used the same set of amino acid residues to bind with the glycoside; however, some mutations alter linkage preferences. according to the conformation analysis of the bound receptors, we could predict the relative selectivity in accordance with in vitro experimental data when disaccharides receptor analogs were used. this hypothesis has been proven correct, since viral isolates from vietnam and thailand have been found which have increased selectivity toward the human cell receptor. it is known that the highly pathogenic avian influenza a virus h5n binds strongly and with high specificity to the avian-type receptor by its hemagglutinin surface protein. preference toward human-type receptors is associated with a positive torsion angle, while avian-type receptor preference is associated with a negative torsion angle.
 it is known that the highly pathogenic avian influenza a virus h5n binds strongly and with high specificity to the avian-type receptor by its hemagglutinin surface protein. this specificity is normally a barrier to viral transmission from birds to humans. however, strains may emerge with mutated hemagglutinin, potentially changing the receptor binding preference from avian to human-type. this hypothesis has been proven correct, since viral isolates from vietnam and thailand have been found which have increased selectivity toward the human cell receptor. the change in binding preference is due to mutation, which can be computationally modelled. the aim of this study is to further explore whether computational simulation could be used as a prediction tool for host type selectivity in emerging variants.  molecular dynamics simulation was employed to study the interactions between receptor models and hemagglutinin proteins from h5n strains a/duck/singapore/3/ mutated a/duck/singapore/3/, a/thailand/1/ and mutated a/thailand/1/. the avian receptor was represented by siaαgal substructure and human receptor by siaαgal. the glycoside binding conformation was monitored throughout the simulations since high selectivity toward a particular host occurs when the sialoside bound with the near-optimized conformation.  the simulation results showed all hemagglutinin proteins used the same set of amino acid residues to bind with the glycoside; however, some mutations alter linkage preferences. preference toward human-type receptors is associated with a positive torsion angle, while avian-type receptor preference is associated with a negative torsion angle. according to the conformation analysis of the bound receptors, we could predict the relative selectivity in accordance with in vitro experimental data when disaccharides receptor analogs were used. 7– september asia pacific bioinformatics network eighth international conference on bioinformatics singapore 