Ras mutation cooperates with ��-catenin activation to drive bladder tumourigenesis

Mutations in the Ras family of proteins (predominantly in H-Ras) occur in approximately 40% of urothelial cell carcinoma (UCC). However, relatively little is known about subsequent mutations/pathway alterations that allow tumour progression. Indeed, expressing mutant H-Ras within the mouse bladder does not lead to tumour formation, unless this is expressed at high levels. The Wnt signalling pathway is deregulated in approximately 25% of UCC, so we examined if this correlated with the activation of MAPK signalling in human UCC and found a significant correlation. To test the functional significance of this association we examined the impact of combining Ras mutation (H-Ras^Q61L or K-Ras^G12D) with an activating β-catenin mutation within the mouse bladder using Cre-LoxP technology. Although alone, neither Ras mutation nor β-catenin activation led to UCC (within 12 months), mice carrying both mutations rapidly developed UCC. Mechanistically this was associated with reduced levels of p21 with dependence on the MAPK signalling pathway. Moreover, tumours from these mice were sensitive to MEK inhibition. Importantly, in human UCC there was a negative correlation between levels of p-ERK and p21 suggesting that p21 accumulation may block tumour progression following Ras mutation. Taken together these data definitively show Ras pathway activation strongly cooperates with Wnt signalling to drive UCC in vivo.     

High throughput interrogation of somatic mutations in high grade serous cancer of the ovary

Background

Epithelial ovarian cancer is the most lethal of all gynecologic malignancies, and high grade serous ovarian cancer (HGSC) is the most common subtype of ovarian cancer. The objective of this study was to determine the frequency and types of point somatic mutations in HGSC using a mutation detection protocol called OncoMap that employs mass spectrometric-based genotyping technology.

Methodology/Principal Findings

The Center for Cancer Genome Discovery (CCGD) Program at the Dana-Farber Cancer Institute (DFCI) has adapted a high-throughput genotyping platform to determine the mutation status of a large panel of known cancer genes. The mutation detection protocol, termed OncoMap has been expanded to detect more than 1000 mutations in 112 oncogenes in formalin-fixed paraffin-embedded (FFPE) tissue samples. We performed OncoMap on a set of 203 FFPE advanced staged HGSC specimens. We isolated genomic DNA from these samples, and after a battery of quality assurance tests, ran each of these samples on the OncoMap v3 platform. 56% (113/203) tumor samples harbored candidate mutations. Sixty-five samples had single mutations (32%) while the remaining samples had ≥2 mutations (24%). 196 candidate mutation calls were made in 50 genes. The most common somatic oncogene mutations were found in EGFR, KRAS, PDGRFα, KIT, and PIK3CA. Other mutations found in additional genes were found at lower frequencies (<3%).

Conclusions/Significance

Sequenom analysis using OncoMap on DNA extracted from FFPE ovarian cancer samples is feasible and leads to the detection of potentially druggable mutations. Screening HGSC for somatic mutations in oncogenes may lead to additional therapies for this patient population.     

The dynamical response of hen egg white lysozyme to the binding of a carbohydrate ligand

It has become clear that the binding of small and large ligands to proteins can invoke significant changes in side chain and main chain motion in the fast picosecond to nanosecond timescale. Recently, the use of a "dynamical proxy" has indicated that changes in these motions often reflect significant changes in conformational entropy. These entropic contributions are sometimes of the same order as the total entropy of binding. Thus, it is important to understand the connections amongst motion between the manifold of states accessible to the native state of proteins, the corresponding entropy, and how this impacts the overall energetics of protein function. The interaction of proteins with carbohydrate ligands is central to a range of biological functions. Here, we examine a classic carbohydrate interaction with an enzyme: the binding of wild-type hen egg white lysozyme (HEWL) to the natural, competitive inhibitor chitotriose. Using NMR relaxation experiments, backbone amide and side chain methyl axial order parameters were obtained across apo and chitotriose-bound HEWL. Upon binding, changes in the apparent amplitude of picosecond to nanosecond main chain and side chain motions are seen across the protein. Indeed, binding of chitotriose renders a large contiguous fraction of HEWL effectively completely rigid. Changes in methyl flexibility are most pronounced closest to the binding site, but average to only a small overall change in the dynamics across the protein. The corresponding change in conformational entropy is unfavorable and estimated to be a significant fraction of the total binding entropy.     

Molecular dynamics simulation of thionated hen egg white lysozyme

Understanding of the driving forces of protein folding is a complex challenge because different types of interactions play a varying role. To investigate the role of hydrogen bonding involving the backbone, the effect of thio substitutions in a protein, hen egg white lysozyme (HEWL), was investigated through molecular dynamics simulations of native as well as partly (only residues in loops) and fully thionated HEWL using the GROMOS 54A7 force field. The results of the three simulations show that the structural properties of fully thionated HEWL clearly differ from those of the native protein, while for partly thionated HEWL they only changed slightly compared with native HEWL. The analysis of the torsional-angle distributions and hydrogen bonds in the backbone suggests that the α-helical segments of native HEWL tend to show a propensity to convert to 3(10)-helical geometry in fully thionated HEWL. A comparison of the simulated quantities with experimental NMR data such as nuclear overhauser effect (NOE) atom-atom distance bounds and (3)J((H)(N)(H)(α))-couplings measured for native HEWL illustrates that the information content of these quantities with respect to the structural changes induced by thionation of the protein backbone is rather limited.     

NMR-based localization of ions involved in salting out of hen egg white lysozyme

NaCl-induced aggregation of hen egg white lysozyme (HEWL) was monitored by NMR spectroscopy. Small, but significant, changes induced by salt addition in TOCSY spectra were attributed to the effect of local reorganization of protein backbone upon ion binding. Salt-induced variations in HN and H alpha chemical shifts were mapped on the HEWL 3D structure which allowed the construction of a scheme of the spatial localization of potential ion binding sites. It was found that in a 0.5 M NaCl solution six chloride anions and at least one sodium cation are bound to preferred sites on the HEWL surface.     

Is High Resolution Melting Analysis (HRMA) Accurate for Detection of Human Disease-Associated Mutations? A Meta Analysis

Background

High Resolution Melting Analysis (HRMA) is becoming the preferred method for mutation detection. However, its accuracy in the individual clinical diagnostic setting is variable. To assess the diagnostic accuracy of HRMA for human mutations in comparison to DNA sequencing in different routine clinical settings, we have conducted a meta-analysis of published reports.

Methodology/Principal Findings

Out of 195 publications obtained from the initial search criteria, thirty-four studies assessing the accuracy of HRMA were included in the meta-analysis. We found that HRMA was a highly sensitive test for detecting disease-associated mutations in humans. Overall, the summary sensitivity was 97.5% (95% confidence interval (CI): 96.8–98.5; I² = 27.0%). Subgroup analysis showed even higher sensitivity for non-HR-1 instruments (sensitivity 98.7% (95%CI: 97.7–99.3; I² = 0.0%)) and an eligible sample size subgroup (sensitivity 99.3% (95%CI: 98.1–99.8; I² = 0.0%)). HRMA specificity showed considerable heterogeneity between studies. Sensitivity of the techniques was influenced by sample size and instrument type but by not sample source or dye type.

Conclusions/Significance

These findings show that HRMA is a highly sensitive, simple and low-cost test to detect human disease-associated mutations, especially for samples with mutations of low incidence. The burden on DNA sequencing could be significantly reduced by the implementation of HRMA, but it should be recognized that its sensitivity varies according to the number of samples with/without mutations, and positive results require DNA sequencing for confirmation.     

Changes in temperature have opposing effects on current amplitude in α7 and α4β2 nicotinic acetylcholine receptors

We have examined the effect of temperature on the electrophysiological properties of three neuronal nicotinic acetylcholine receptor (nAChR) subtypes: the rapidly desensitizing homomeric α7 nAChR, the more slowly desensitizing heteromeric α4β2 nAChR and on α7 nAChRs containing a transmembrane mutation (L247T) that results in dramatically reduced desensitization. In all cases, the functional properties of receptors expressed in Xenopus oocytes at room temperature (RT; 21°C) were compared to those recorded at either physiological temperature (37°C) or at lower temperature (4°C). Alterations in temperature had dramatically differing effects on the amplitude of whole-cell responses detected with these three nAChR subtypes. Compared to responses at RT, the amplitude of agonist-evoked responses with α4β2 nAChRs was increased at high temperature (125±9%, n = 6, P<0.01) and reduced at low temperature (47±5%, n = 6, P<0.01), whereas the amplitude of α7 responses was reduced at high temperature (27±7%, n = 11, P<0.001) and increased at low temperatures (224±16%, n = 10, P<0.001). In contrast to the effects of temperature on α4β2 and wild type α7 nAChRs, the amplitude of α7 nAChRs containing the L247T mutation was unaffected by changes in temperature. In addition, changes in temperature had little or no effect on current amplitude when α7 nAChRs were activated by the largely non-desensitizing allosteric agonist 4BP-TQS. Despite these differing effects of temperature on the amplitude of agonist-evoked responses in different nAChRs, changes in temperature had a consistent effect on the rate of receptor desensitization on all subtypes examined. In all cases, higher temperature resulted in increased rates of desensitization. Thus, it appears that the differing effects of temperature on the amplitudes of whole-cell responses cannot be explained by temperature-induced changes in receptor desensitization rates.     

A frameshift mutation in golden retriever dogs with progressive retinal atrophy endorses SLC4A3 as a candidate gene for human retinal degenerations

Progressive retinal atrophy (PRA) in dogs, the canine equivalent of retinitis pigmentosa (RP) in humans, is characterised by vision loss due to degeneration of the photoreceptor cells in the retina, eventually leading to complete blindness. It affects more than 100 dog breeds, and is caused by numerous mutations. RP affects 1 in 4000 people in the Western world and 70% of causal mutations remain unknown. Canine diseases are natural models for the study of human diseases and are becoming increasingly useful for the development of therapies in humans. One variant, prcd-PRA, only accounts for a small proportion of PRA cases in the Golden Retriever (GR) breed. Using genome-wide association with 27 cases and 19 controls we identified a novel PRA locus on CFA37 (p_raw = 1.94×10⁻¹⁰, p_genome = 1.0×10⁻⁵), where a 644 kb region was homozygous within cases. A frameshift mutation was identified in a solute carrier anion exchanger gene (SLC4A3) located within this region. This variant was present in 56% of PRA cases and 87% of obligate carriers, and displayed a recessive mode of inheritance with full penetrance within those lineages in which it segregated. Allele frequencies are approximately 4% in the UK, 6% in Sweden and 2% in France, but the variant has not been found in GRs from the US. A large proportion of cases (approximately 44%) remain unexplained, indicating that PRA in this breed is genetically heterogeneous and caused by at least three mutations. SLC4A3 is important for retinal function and has not previously been associated with spontaneously occurring retinal degenerations in any other species, including humans.     

The naturally occurring YMDD mutation among patients chronically infected HBV and untreated with lamivudine: a systematic review and meta-analysis

Background

Several recent reports have demonstrated that tyrosine (Y)-methionine (M)-aspartic acid (D)-aspartic acid (D) (YMDD) motif mutations can naturally occur in chronic HBV patients without antiviral treatment such as lamivudine therapy. This paper aims to assess the overall spontaneous incidence and related risk factors of YMDD-motif mutations among lamivudine-naïve chronic HBV carriers, so as to provide some clue for clinical treatment of hepatitis B.

Methodology/Principal Findings

Chinese and English literatures were searched for studies reporting natural YMDD mutations among untreated chronic HBV patients from 2001 to 2010. The incidence estimates were summarized and analyzed by meta-analyses. Forty-seven eligible articles from eight countries were selected in this review (13 in English and 34 in Chinese). The pooled incidence of YMDD-motif mutation among untreated chronic HBV patients from eight countries was 12.21% (95% CI: 9.69%–14.95%). China had an incidence of 13.38% (95% CI: 10.90%–16.07%) and seven other countries had an incidence of 9.90% (95% CI: 3.28%–19.55%), respectively. Lamivudine therapy would increase the risk of mutations 5.23 times higher than the untreated patients. A higher HBV DNA copy number was associated with increased incidence of natural YMDD mutation. No significant difference was found in YMDD mutation incidence between groups of different gender, age, HBeAg status, patients'' ALT (alanine aminotransferase) level, and between the groups of HBV genotype B and C.

Conclusions

The YMDD-motif mutations can occur spontaneously with a relatively high incidence in CHB patients untreated with lamivudine. These mutations might be the consequence of accumulated base mismatch due to the nature of viral polymerase. More fundamental and clinical studies are needed to clarify the influence of YMDD mutations in hepatitis B progression and antiviral treatment.     

Mutation of His465 Alters the pH-dependent Spectroscopic Properties of Escherichia coli Glutamate Decarboxylase and Broadens the Range of Its Activity toward More Alkaline pH

Glutamate decarboxylase (GadB) from Escherichia coli is a hexameric, pyridoxal 5′-phosphate-dependent enzyme catalyzing CO₂ release from the α-carboxyl group of l-glutamate to yield γ-aminobutyrate. GadB exhibits an acidic pH optimum and undergoes a spectroscopically detectable and strongly cooperative pH-dependent conformational change involving at least six protons. Crystallographic studies showed that at mildly alkaline pH GadB is inactive because all active sites are locked by the C termini and that the 340 nm absorbance is an aldamine formed by the pyridoxal 5′-phosphate-Lys²⁷⁶ Schiff base with the distal nitrogen of His⁴⁶⁵, the penultimate residue in the GadB sequence. Herein we show that His⁴⁶⁵ has a massive influence on the equilibrium between active and inactive forms, the former being favored when this residue is absent. His⁴⁶⁵ contributes with n ≈ 2.5 to the overall cooperativity of the system. The residual cooperativity (n ≈ 3) is associated with the conformational changes still occurring at the N-terminal ends regardless of the mutation. His⁴⁶⁵, dispensable for the cooperativity that affects enzyme activity, is essential to include the conformational change of the N termini into the cooperativity of the whole system. In the absence of His⁴⁶⁵, a 330-nm absorbing species appears, with fluorescence emission spectra more complex than model compounds and consisting of two maxima at 390 and 510 nm. Because His⁴⁶⁵ mutants are active at pH well above 5.7, they appear to be suitable for biotechnological applications.     

Inclusion of Many-Body Effects in the Additive CHARMM Protein CMAP Potential Results in Enhanced Cooperativity of α-Helix and β-Hairpin Formation

Folding simulations on peptides and proteins using empirical force fields have demonstrated the sensitivity of the results to details of the backbone potential. A recently revised version of the additive CHARMM protein force field, which includes optimization of the backbone CMAP potential to achieve good balance between different types of secondary structure, correcting the α-helical bias present in the former CHARMM22/CMAP energy function, is shown to result in improved cooperativity for the helix-coil transition. This is due to retention of the empirical corrections introduced in the original CMAP to reproduce folded protein structures—corrections that capture many-body effects missing from an energy surface fitted to gas phase calculations on dipeptides. The experimental temperature dependence of helix formation in (AAQAA)₃ and parameters for helix nucleation and elongation are in much better agreement with experiment than those obtained with other recent force fields. In contrast, CMAP parameters derived by fitting to a vacuum quantum mechanical surface for the alanine dipeptide do not reproduce the enhanced cooperativity, showing that the empirical backbone corrections, and not some other feature of the force field, are responsible. We also find that the cooperativity of β-hairpin formation is much improved relative to other force fields we have studied. Comparison with (ϕ,ψ) distributions from the Protein Data Bank further justifies the inclusion of many-body effects in the CMAP. These results suggest that the revised energy function will be suitable for both simulations of unfolded or intrinsically disordered proteins and for investigating protein-folding mechanisms.     

Effects of dithiothreitol on the amyloid fibrillogenesis of hen egg-white lysozyme

At least 25 human proteins can fold abnormally to form pathological deposits that are associated with several degenerative diseases. Despite extensive investigation on amyloid fibrillation, the detailed molecular mechanisms remain rather elusive and there are currently no effective cures for treating these amyloid diseases. The present study examined the effects of dithiothreitol on the fibrillation of hen egg-white lysozyme (HEWL). Our results revealed that the fibrillation of hen lysozyme was significantly inhibited by reduced dithiothreitol (DTT^red) while oxidized dithiothreitol (DTT^ox) had no anti-aggregating activity. Effective inhibitory activity against hen lysozyme fibrillation was observed only when DTT^red was added within 8 days of incubation. Our results showed that the initial addition of DTT^red interacted with HEWL, leading to a loss in conformational stability. It was concluded from our findings that DTT^red-induced attenuation of HEWL fibrillation may be associated with disulfide disruption and extensive structural unfolding of HEWL. Our data may contribute to rational design of effective therapeutic strategies for amyloid diseases.     

Effective charge measurements reveal selective and preferential accumulation of anions, but not cations, at the protein surface in dilute salt solutions

Specific-ion effects are ubiquitous in nature; however, their underlying mechanisms remain elusive. Although Hofmeister-ion effects on proteins are observed at higher (>0.3M) salt concentrations, in dilute (<0.1M) salt solutions nonspecific electrostatic screening is considered to be dominant. Here, using effective charge (Q*) measurements of hen-egg white lysozyme (HEWL) as a direct and differential measure of ion-association, we experimentally show that anions selectively and preferentially accumulate at the protein surface even at low (<100 mM) salt concentrations. At a given ion normality (50 mN), the HEWL Q* was dependent on anion, but not cation (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, GdnH⁺, and Ca²⁺), identity. The Q* decreased in the order F⁻ > Cl⁻ > Br⁻ > NO₃⁻ ∼ I⁻ > SCN⁻ > ClO₄⁻ ≫ SO₄²⁻, demonstrating progressively greater binding of the monovalent anions to HEWL and also show that the SO₄²⁻ anion, despite being strongly hydrated, interacts directly with the HEWL surface. Under our experimental conditions, we observe a remarkable asymmetry between anions and cations in their interactions with the HEWL surface.     

Comparison of the effects of 2,2,2-trifluoroethanol on peptide and protein structure and function

The co-solvent 2,2,2-trifluoroethanol (TFE) has been often used to aid formation of secondary structure in solution peptides or alternately as a denaturant within protein folding studies. Hen egg white lysozyme (HEWL) and a synthetic model peptide defining HEWL helix-4 were used as comparative model systems to systematically investigate the effect of increasing TFE concentrations on the structure of proteins and peptides. HEWL was analyzed using NMR, far-UV CD and fluorescence spectroscopy; with correlation of these results towards changes in enzymatic activity and the helix-4 peptide was analysed using NMR. Data illustrates two conflicting modes of interaction: Low TFE concentrations stabilize tertiary structure, observed from an increase in the number of NMR NOE contacts. Higher TFE concentrations denatured HEWL with the loss of lysozyme tertiary structure. The effects of TFE upon secondary structural elements within HEWL are distinct from those observed for the helix-4 peptide. This illustrates a dissimilar interaction of TFE towards both protein and peptide at equivalent TFE concentrations. The concentration that TFE promotes stabilization over denaturation is likely to be protein dependent although the structural action can be extrapolated to other protein systems with implications for the use of TFE in structural stability studies.     

TGFBR2 and BAX mononucleotide tract mutations, microsatellite instability, and prognosis in 1072 colorectal cancers

Background

Mononucleotide tracts in the coding regions of the TGFBR2 and BAX genes are commonly mutated in microsatellite instability-high (MSI-high) colon cancers. The receptor TGFBR2 plays an important role in the TGFB1 (transforming growth factor-β, TGF-β) signaling pathway, and BAX plays a key role in apoptosis. However, a role of TGFBR2 or BAX mononucleotide mutation in colorectal cancer as a prognostic biomarker remains uncertain.

Methodology/Principal Findings

We utilized a database of 1072 rectal and colon cancers in two prospective cohort studies (the Nurses'' Health Study and the Health Professionals Follow-up Study). Cox proportional hazards model was used to compute mortality hazard ratio (HR), adjusted for clinical, pathological and molecular features including the CpG island methylator phenotype (CIMP), LINE-1 methylation, and KRAS, BRAF and PIK3CA mutations. MSI-high was observed in 15% (162/1072) of all colorectal cancers. TGFBR2 and BAX mononucleotide mutations were detected in 74% (117/159) and 30% (48/158) of MSI-high tumors, respectively. In Kaplan-Meier analysis as well as univariate and multivariate Cox regression analyses, compared to microsatellite stable (MSS)/MSI-low cases, MSI-high cases were associated with superior colorectal cancer-specific survival [adjusted HR, 0.34; 95% confidence interval (CI), 0.20–0.57] regardless of TGFBR2 or BAX mutation status. Among MSI-high tumors, TGFBR2 mononucleotide mutation was associated with CIMP-high independent of other variables [multivariate odds ratio, 3.57; 95% CI, 1.66–7.66; p = 0.0011].

Conclusions

TGFBR2 or BAX mononucleotide mutations are not associated with the patient survival outcome in MSI-high colorectal cancer. Our data do not support those mutations as prognostic biomarkers (beyond MSI) in colorectal carcinoma.     

Kinetic epitope mapping of the chicken lysozyme.HyHEL-10 Fab complex: delineation of docking trajectories.

The rate constants, k(on), for the formation of hen (chicken) lysozyme (HEWL). Fab-10 complexes have been determined for wild-type (WT) and epitope-mutated lysozymes by a homogeneous solution method based on the 95% reduced enzymatic activity of the complex. The values fall within a narrow 10-fold range [(0.18 to 1.92) x 10(6) M(-1)s(-l)]. The affinity constants, K(D), cover a broader, 440-fold, range from 0.075 to 33 nM. Values of K(D) as high as 7 microM were obtained for the complexes prepared from some mutations at HEWL positions 96 and 97, but the associated kinetic constants could not be determined. The values of k(on) are negatively correlated with side-chain volume at position 101HEWL, but are essentially independent of this parameter for position 21HEWL substitutions. The multiple mutations made at positions 21HEWL and 101HEWL provide sufficient experimental data on complex formation to evaluate phi values [phi = (deltadeltaGon)/(deltadeltaG(D))] at these two positions to begin to define trajectories for protein-protein association. The data, when interpreted within the concept of a two-step association sequence embracing a metastable encounter complex intermediate, argue that the rate determining step at position 21HEWL (phiavg = 0.2) is encounter complex formation, but the larger phi(avg) value of 0.36 experienced for most position 101HEWL mutations indicates a larger contribution from the post-encounter annealing process at this site for these replacements.     

DNA sequence profiles of the colorectal cancer critical gene set KRAS-BRAF-PIK3CA-PTEN-TP53 related to age at disease onset

The incidence of colorectal cancer (CRC) increases with age and early onset indicates an increased likelihood for genetic predisposition for this disease. The somatic genetics of tumor development in relation to patient age remains mostly unknown. We have examined the mutation status of five known cancer critical genes in relation to age at diagnosis, and compared the genomic complexity of tumors from young patients without known CRC syndromes with those from elderly patients. Among 181 CRC patients, stratified by microsatellite instability status, DNA sequence changes were identified in KRAS (32%), BRAF (16%), PIK3CA (4%), PTEN (14%) and TP53 (51%). In patients younger than 50 years (n = 45), PIK3CA mutations were not observed and TP53 mutations were more frequent than in the older age groups. The total gene mutation index was lowest in tumors from the youngest patients. In contrast, the genome complexity, assessed as copy number aberrations, was highest in tumors from the youngest patients. A comparable number of tumors from young (<50 years) and old patients (>70 years) was quadruple negative for the four predictive gene markers (KRAS-BRAF-PIK3CA-PTEN); however, 16% of young versus only 1% of the old patients had tumor mutations in PTEN/PIK3CA exclusively. This implies that mutation testing for prediction of EGFR treatment response may be restricted to KRAS and BRAF in elderly (>70 years) patients. Distinct genetic differences found in tumors from young and elderly patients, whom are comparable for known clinical and pathological variables, indicate that young patients have a different genetic risk profile for CRC development than older patients.     

Evolution of H3N2 influenza virus in a guinea pig model

Studies of influenza virus evolution under controlled experimental conditions can provide a better understanding of the consequences of evolutionary processes with and without immunological pressure. Characterization of evolved strains assists in the development of predictive algorithms for both the selection of subtypes represented in the seasonal influenza vaccine and the design of novel immune refocused vaccines. To obtain data on the evolution of influenza in a controlled setting, naïve and immunized Guinea pigs were infected with influenza A/Wyoming/2003 (H3N2). Virus progeny from nasal wash samples were assessed for variation in the dominant and other epitopes by sequencing the hemagglutinin (HA) gene to quantify evolutionary changes. Viral RNA from the nasal washes from infection of naïve and immune animals contained 6% and 24.5% HA variant sequences, respectively. Analysis of mutations relative to antigenic epitopes indicated that adaptive immunity played a key role in virus evolution. HA mutations in immunized animals were associated with loss of glycosylation and changes in charge and hydrophobicity in and near residues within known epitopes. Four regions of HA-1 (75–85, 125–135, 165–170, 225–230) contained residues of highest variability. These sites are adjacent to or within known epitopes and appear to play an important role in antigenic variation. Recognition of the role of these sites during evolution will lead to a better understanding of the nature of evolution which help in the prediction of future strains for selection of seasonal vaccines and the design of novel vaccines intended to stimulated broadened cross-reactive protection to conserved sites outside of dominant epitopes.