Comparing the structures of H3, H5 and H9 subtype haemagglutinins, we deduced a structural basis for including all 15 influenza subtypes in four clades. H3, H5 and H9 represent three of these clades; we now report the structure of an H7 HA as a representative of the fourth clade. We confirm the structure of the turn at the N-terminus of the conserved central alpha-helix of HA2, and the combination of ionisable residues near the "fusion peptide" as clade-specific features. We compare the structures of three H1 HAs with H5 HA in the same clade, to refine our previous classification and we confirm the division of the clades into two groups of two. We also show the roles of carbohydrate side chains in the esterase-fusion domain boundaries in the formation of clade-specific structural markers.

There are 15 subtypes of influenza A virus (H1-H15), all of which are found in avian species. Three caused pandemics in the last century: H1 in 1918 (and 1977), H2 in 1957 and H3 in 1968. In 1997, an H5 avian virus and in 1999 an H9 virus caused outbreaks of respiratory disease in Hong Kong. We have determined the three-dimensional structures of the haemagglutinins (HAs) from H5 avian and H9 swine viruses closely related to the viruses isolated from humans in Hong Kong. We have compared them with known structures of the H3 HA from the virus that caused the 1968 H3 pandemic and of the HA?esterase?fusion (HEF) glycoprotein from an influenza C virus. Structure and sequence comparisons suggest that HA subtypes may have originated by diversification of properties that affected the metastability of HAs required for their membrane fusion activities in viral infection.

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.

A single amino acid substitution, Asp-63 to Asn-63, was detected in the hemagglutinin of an antigenic variant of the 1968 Hong Kong (H3) influenza virus that was selected by growth of the wild-type virus in the presence of a monoclonal antibody. The mutation generates an oligosaccharide attachment site, Asn-Cys-Thr at residues 63-65, that is glycosylated. Immunoprecipitation experiments with extracts from variant virus-infected cells prepared in the presence or absence of tunicamycin, which inhibits glycosylation, demonstrate that addition of the new oligosaccharide side chain is required to prevent reaction with the monoclonal antibody. Similar experiments with the virus of the 1969 Hong Kong influenza epidemic, A/England/878/69, which also contains a hemagglutinin glycosylated at residue 63, support this conclusion and provide evidence for the epidemiological significance of carbohydrate-mediated modifications of hemagglutinin antigenicity.

We have determined the structure of the HA of an avian influenza virus, A/duck/Ukraine/63, a member of the same antigenic subtype, H3, as the virus that caused the 1968 Hong Kong influenza pandemic, and a possible progenitor of the pandemic virus. We find that structurally significant differences between the avian and the human HAs are restricted to the receptor-binding site particularly the substitutions Q226L and G228S that cause the site to open and residues within it to rearrange, including the conserved residues Y98, W153, and H183. We have also analyzed complexes formed by the HA with sialopentasaccharides in which the terminal sialic acid is in either alpha2,3- or alpha2,6-linkage to galactose. Comparing the structures of complexes in which an alpha2,3-linked receptor analog is bound to the H3 avian HA or to an H5 avian HA leads to the suggestion that all avian influenza HAs bind to their preferred alpha2,3-linked receptors similarly, with the analog in a trans conformation about the glycosidic linkage. We find that alpha2,6-linked analogs are bound by both human and avian HAs in a cis conformation, and that the incompatibility of an alpha2,6-linked receptor with the alpha2,3-linkage-specific H3 avian HA-binding site is partially resolved by a small change in the position and orientation of the sialic acid. We discuss our results in relation to the mechanism of transfer of influenza viruses between species.

SUMMARY: StrBioLib is a library of Java classes useful for developing software for computational structural biology research. StrBioLib contains classes to represent and manipulate protein structures, biopolymer sequences, sets of biopolymer sequences, and alignments between biopolymers based on either sequence or structure. Interfaces are provided to interact with commonly used bioinformatics applications, including (psi)-blast, modeller, muscle and Primer3, and tools are provided to read and write many file formats used to represent bioinformatic data. The library includes a general-purpose neural network object with multiple training algorithms, the Hooke and Jeeves non-linear optimization algorithm, and tools for efficient C-style string parsing and formatting. StrBioLib is the basis for the Pred2ary secondary structure prediction program, is used to build the astral compendium for sequence and structure analysis, and has been extensively tested through use in many smaller projects. Examples and documentation are available at the site below. AVAILABILITY: StrBioLib may be obtained under the terms of the GNU LGPL license from http://strbio.sourceforge.net/

The observed distribution of protein structures can give us important clues about the underlying evolutionary process, imposing important constraints on possible models. The availability of results from an increasing number of genome projects has made the development of these models an active area of research. Models explaining the observed distribution of structures have focused on the inherent functional capabilities and structural properties of different folds and on the evolutionary dynamics. Increasingly, these elements are being combined.

Evolutionary traces of thermophilic adaptation are manifest, on the whole-genome level, in compositional biases toward certain types of amino acids. However, it is sometimes difficult to discern their causes without a clear understanding of underlying physical mechanisms of thermal stabilization of proteins. For example, it is well-known that hyperthermophiles feature a greater proportion of charged residues, but, surprisingly, the excess of positively charged residues is almost entirely due to lysines but not arginines in the majority of hyperthermophilic genomes. All-atom simulations show that lysines have a much greater number of accessible rotamers than arginines of similar degree of burial in folded states of proteins. This finding suggests that lysines would preferentially entropically stabilize the native state. Indeed, we show in computational experiments that arginine-to-lysine amino acid substitutions result in noticeable stabilization of proteins. We then hypothesize that if evolution uses this physical mechanism as a complement to electrostatic stabilization in its strategies of thermophilic adaptation, then hyperthermostable organisms would have much greater content of lysines in their proteomes than comparably sized and similarly charged arginines. Consistent with that, high-throughput comparative analysis of complete proteomes shows extremely strong bias toward arginine-to-lysine replacement in hyperthermophilic organisms and overall much greater content of lysines than arginines in hyperthermophiles. This finding cannot be explained by genomic GC compos