Gnarled-trunk evolutionary model of influenza A virus hemagglutinin

Human influenza A viruses undergo antigenic changes with gradual accumulation of amino acid substitutions on the hemagglutinin (HA) molecule. A strong antigenic mismatch between vaccine and epidemic strains often requires the replacement of influenza vaccines worldwide. To establish a practical model enabling us to predict the future direction of the influenza virus evolution, relative distances of amino acid sequences among past epidemic strains were analyzed by multidimensional scaling (MDS). We found that human influenza viruses have evolved along a gnarled evolutionary pathway with an approximately constant curvature in the MDS-constructed 3D space. The gnarled pathway indicated that evolution on the trunk favored multiple substitutions at the same amino acid positions on HA. The constant curvature was reasonably explained by assuming that the rate of amino acid substitutions varied from one position to another according to a gamma distribution. Furthermore, we utilized the estimated parameters of the gamma distribution to predict the amino acid substitutions on HA in subsequent years. Retrospective prediction tests for 12 years from 1997 to 2009 showed that 70% of actual amino acid substitutions were correctly predicted, and that 45% of predicted amino acid substitutions have been actually observed. Although it remains unsolved how to predict the exact timing of antigenic changes, the present results suggest that our model may have the potential to recognize emerging epidemic strains.     

Accumulation of amino acid substitutions promotes irreversible structural changes in the hemagglutinin of human influenza AH3 virus during evolution

During protein evolution the amino acid substitutions accumulate with time. However, the effect of accumulation of the amino acid substitutions to structural changes has not been estimated well. We will propose that the discordance of amino acid substitution on the HA protein of influenza A virus is useful for the assessment of structural changes during evolution. Discordance value can be obtained from the experimental data of tolerance or intolerance by introducing site directed mutagenesis at the homologous positions of two HA proteins holding the same amino acid residues. The value of discordance correlated to the number of amino acid differences among proteins. In the H3HA discordance rate was calculated to be 0.45% per one amino acid change. Furthermore, discordance of amino acid substitutions suggests that tolerable amino acid substitutions in different order have a probability of promoting irreversible divergence of the HA protein to different subtypes.     

Selection of genetic variants of simian immunodeficiency virus in persistently infected rhesus monkeys.

Genetic and antigenic variation may be one means by which lentiviruses that cause AIDS avoid elimination by host immune responses. Genetic variation in the envelope gene (env) was studied by comparing the nucleotide sequences of 27 clones obtained from two rhesus monkeys infected with molecularly cloned simian immunodeficiency virus. All 27 clones differed from each other and differed from the input clone in the gp120 (SU) portion of the envelope gene. Nucleotide substitutions were shown to accumulate with time at an average rate of 8.5 per 1,000 per year in SU. Surprisingly, the majority of nucleotide substitutions (81%) resulted in amino acid changes. Variation in SU was not random but occurred predominantly in five discrete regions. Within these variable regions, a remarkable 98% of the nucleotide substitutions changed the amino acid. These results demonstrate that extensive sequence variability accumulates in vivo after infection with molecularly cloned virus and that selection occurs in vivo for changes in distinct variable regions in env.     

Reactivity of human convalescent sera with influenza virus hemagglutinin protein mutants at antigenic site A

How the antibodies of individual convalescent human sera bind to each amino acid residue at the antigenic sites of hemagglutinin (HA) of influenza viruses, and how the antigenic drift strains of influenza viruses are selected by human sera, is not well understood. In our previous study, it was found by a binding assay with a chimeric HA between A/Kamata/14/91 (Ka/91) and A/Aichi/2/68 that convalescent human sera, following Ka/91 like (H3N2) virus infection, bind to antigenic site A of Ka/91 HA. Here using chimeric HAs possessing single amino acid substitutions at site A, it was determined how those human sera recognize each amino acid residue at antigenic site A. It was found that the capacity of human sera to recognize amino acid substitutions at site A differs from one person to another and that some amino acid substitutions result in all convalescent human sera losing their binding capacity. Among these amino acid substitutions, certain ones might be selected by chance, thus creating successive antigenic drift. Phylogenetic analysis of the drift strains of Ka/91 showed amino acid substitutions at positions 133, 135 and 145 were on the main stream of the phylogenetic tree. Indeed, all of the investigated convalescent sera failed to recognize one of them.     

Insertions and deletions trigger adaptive walks in Drosophila proteins

Maps that relate all possible genotypes or phenotypes to fitness--fitness landscapes--are central to the evolution of life, but remain poorly known. An insertion or a deletion (indel) of one or several amino acids constitutes a substantial leap of a protein within the space of amino acid sequences, and it is unlikely that after such a leap the new sequence corresponds precisely to a fitness peak. Thus, one can expect an indel in the protein-coding sequence that gets fixed in a population to be followed by some number of adaptive amino acid substitutions, which move the new sequence towards a nearby fitness peak. Here, we study substitutions that occur after a frame-preserving indel in evolving proteins of Drosophila. An insertion triggers 1.03 ± 0.75 amino acid substitutions within the protein region centred at the site of insertion, and a deletion triggers 4.77 ± 1.03 substitutions within such a region. The difference between these values is probably owing to a higher fraction of effectively neutral insertions. Almost all of the triggered amino acid substitutions can be attributed to positive selection, and most of them occur relatively soon after the triggering indel and take place upstream of its site. A high fraction of substitutions that follow an indel occur at previously conserved sites, suggesting that an indel substantially changes selection that shapes the protein region around it. Thus, an indel is often followed by an adaptive walk of length that is in agreement with the theory of molecular adaptation.     

Mutagenic mapping of helical structures in the transmembrane segments of the yeast alpha-factor receptor

The alpha-mating pheromone receptor encoded by the yeast STE2 gene is a G protein coupled receptor that initiates signaling via a MAP kinase pathway that prepares haploid cells for mating. To establish the range of allowed amino acid substitutions within transmembrane segments of this receptor, we conducted extensive random mutagenesis of receptors followed by screening for receptor function. A total of 157 amino acid positions in seven different mutagenic libraries corresponding to the seven predicted transmembrane segments were analyzed, yielding 390 alleles that retain at least 60 % of normal signaling function. These alleles contained a total of 576 unique amino acid substitutions, including 61 % of all the possible amino acid changes that can arise from single base substitutions. The receptor exhibits a surprising tolerance for amino acid substitutions. Every amino acid in the mutagenized regions of the transmembrane regions could be substituted by at least one other residue. Polar amino acids were tolerated in functional receptors at 115 different positions (73 % of the total). Hydrophobic amino acids were tolerated in functional receptors at all mutagenized positions. Substitutions introducing proline residues were recovered at 53 % of all positions where they could be brought about by single base changes. Residues with charged side-chains could also be tolerated at 53 % of all positions where they were accessible through single base changes. The spectrum of allowed amino acid substitutions was characterized in terms of the hydrophobicity, radius of gyration, and charge of the allowed substitutions and mapped onto alpha-helical structures. By comparing the patterns of allowed substitutions with the recently determined structure of rhodopsin, structural features indicative of helix-helix interactions can be discerned in spite of the extreme sequence divergence between these two proteins.     

Cold-adapted tubulins in the glacier ice worm, Mesenchytraeus solifugus

Glacier ice worms, Mesenchytraeus solifugus and related species, are the only known annelids that survive obligately in glacier ice and snow. One fundamental component of cold temperature adaptation is the ability to polymerize tubulin, which typically depolymerizes at low physiological temperatures (e.g., <10 degrees C) in most temperate species. In this study, we isolated two alpha-tubulin (Msalpha1, Msalpha2) and two beta-tubulin (Msbeta1, Msbeta2) subunits from an ice worm cDNA library, and compared their predicted amino acid sequences with homologues from other cold-adapted organisms (e.g., Antarctic fish, ciliate) in an effort to identify species-specific amino acid substitutions that contribute to cold temperature-dependent tubulin polymerization. Our comparisons and predicted protein structures suggest that ice worm-specific amino acid substitutions stabilize lateral contact associations, particularly between beta-tubulin protofilaments, but these substitutions occur at different positions in comparison with other cold-adapted tubulins. The ice worm tubulin gene family appears relatively small, comprising one primary alpha- and one primary beta-tubulin monomers, though minor isoforms and psuedogenes were identified. Our analyses suggest that variation occurs in the strategies (i.e., species-specific amino acid substitutions, gene number) by which cold-adapted taxa have evolved the ability to polymerize tubulin at low physiological temperatures.     

Chain length-dependent binding of fatty acid anions to human serum albumin studied by site-directed mutagenesis

Human serum albumin is the most abundant protein in the circulatory system, and one of its principal functions is to transport fatty acids. Binding of octanoate, decanoate, laurate and myristate was studied by a rate-of-dialysis technique. The primary association constants increased, but not linearly, with chain length. The number of high-affinity sites also increased with chain length; octanoate and decanoate bind to one such site, whereas laurate and myristate most probably bind to two sites. Albumin is composed of three homologous helical domains (I-III), which can be subdivided into two subdomains (A and B). For getting information about the positions of the high-affinity sites we produced 13 recombinant isoforms mutated in four different subdomains. Results obtained with these albumins are in accordance with the following model: octanoate and decanoate bind to a single site in subdomain IIIA, laurate binds to sites in subdomains IIIA and IIIB, whereas myristate binds in subdomains IB and IIIB. The results also showed that primary fatty acid binding is sensitive to amino acid substitutions in other parts of the protein. This is in contrast to the effect of amino acid substitutions of genetic albumin variants (alloalbumins). Usually these substitutions, which are situated at the surface of the protein, have no effect on fatty acid binding. Binding of fatty acid anions to different high-affinity sites and the sensitivity of these sites to amino acid substitutions elsewhere in the protein (and perhaps also to other types of modifications) are important factors that could effect simultaneous binding of other ligands, e.g. in patients treated with albumin-binding drugs.     

Amino acid substitutions in structurally related proteins. A pattern recognition approach. Determination of a new and efficient scoring matrix

Amino acid substitutions in evolutionarily related proteins have been studied from a structural point of view. We consider here that an amino acid al in a protein p1 has been replaced by the amino acid a2 in the structurally similar protein p2 if, after superposition of the p1 and p2 structures, the a1 and a2 C alpha atoms are no more than 1.2 A apart. Thirty-two proteins, grouped in 11 classes, have been analysed by this method. This produced 2860 amino acid pairs (substitutions), which were analysed by multi-dimensional statistical methods. The main results are as follows: (1) according to the observed exchangeability of amino acid side-chains, only four groups (strong clusters) could be delineated; (i) Ile and Val, (ii) Leu and Met, (iii) Lys, Arg and Gln, and (iv) Tyr and Phe. The other residues could not be classified. (2) The matrix of distances between amino acids, or scoring matrix, determined from this study, is different from any other published matrix. (3) Except for the distance matrices based on the chemical properties of amino acid side-chains, which can be grouped together, all other published matrices are different from one another. (4) The distance matrix determined in this study seems to be very efficient for aligning distantly related protein sequences.     

A loss-of-function mutation in the rice KNOX type homeobox gene, OSH3

KNOX homeodomain (HD) proteins encoded by KNOTTED1-like homeobox genes (KNOX genes) are thought to work as switches for cells to change from an indeterminate to a determinate state, although their direct functions are not clear. In the process of isolating KNOX genes from rice, we found that one gene, named OSH3, has two amino acid substitutions in three of the invariant amino acid residues in the HD of KNOX proteins. These amino acid substitutions are not universal in rice: two of the cultivars from the Indica variety of rice do not carry those substitutions but two of the cultivars from Japonica variety do. We tested the effect of these amino acid substitutions on their ability to form dimers and to induce abnormal morophologies when overexpressed in transgenic plants. We found that OSH3 without those substitutions can form dimers and can induce an abnormal phenotype in overexpression studies, and that OSH3 with those amino acid substitutions is defective in both. Based on these observations, we concluded that OSH3 from two of the cultivars from the Japonica variety could have lost its original function, or could have acquired a novel function by modifying the action of HD, or both.     

The influence of amino acid substitutions on the conformational energy of cytochrome c.

Conformational energies have been evaluated for each of the staggered side-chain conformations associated with the 261 amino acid substitutions known to occur among 60 eucaryotic species. At least 86% of these substitutions can be sterically accommodated (one at a time) within the structure of horse-heart cytochrome c resulting from conformational energy refinement. Simultaneous incorporation of all pertinent amino acid substitutions found in eight representative species into the refined horse-heart structure is also shown to be sterically possible, with few exceptions. In two cases (Pekin duck cytochrome with 10 substitutions and Samia cynthia cytochrome with 24 substitutions), all substitutions could be readily incorporated, and the total energies associated with their computed structures differed by less than 10 kcal/mol from that of horse-heart cytochrome c. In the cytochromes from rattlesnake (22 substitutions), tuna (18 substitutions), and Neurospora crassa (36 substitutions), tyrosine could not be substituted for phenylalanine at position 46, within the constraints of the calculations. However, when all of the remaining substitutions were incorporated into these three cytochromes, their computed conformational energies differed by less than 30 kcal/mol from that of horse-heart cytochrome c. Between two and four amino acid substitutions cause high energies in the cytochromes from human, baker's yeast, and cotton seed, but all of the remaining substitutions are consistent with a low energy conformation. These results suggest that the structures of homologous proteins may be even more similar than has previously been recognized. Substitutions of all possible amino acid types at the invariant positions (where all eucaryotic cytochromes c bear the same amino acid) have revealed some cases where different amino acids can be accommodated, thus demonstrating that the biological constraints on amino acid substitutions are often different from the purely steric constraints investigated in this work.     

Neutralization epitopes on rotavirus SA11 4fM outer capsid proteins.

The VP7 and VP4 genes of seven antigenic mutants of simian rotavirus SA11 4fM (serotype 3) selected after 39 passages in the presence of SA11 4fM hyperimmune antiserum, were sequenced. Nucleotide sequence analysis indicated the following. (i) Twice as many amino acid substitutions occurred in the VP7 protein than in VP4, which has a molecular weight twice that of VP7. (ii) Most amino acid changes that occurred clustered in six variable regions of VP7 and in two variable regions of VP4; these variable regions may represent immunodominant epitopes. (iii) Most amino acid substitutions that occurred in VP7 and VP4 of these mutants were also observed in antigenic mutants selected with neutralizing monoclonal antibodies (NMAbs); however, some amino acid substitutions occurred that were not selected for NMAbs. (iv) On VP7, some of the neutralization epitopes appeared to be interrelated because amino acid substitution in one site affected binding of specific NMAbs to other sites, while other neutralization epitopes on VP7 appeared to be independent, in that amino acid substitution in one site did not affect the binding of NMAbs to another distant site.     

Guanidine hydrochloride denaturation studies of mutant forms of staphylococcal nuclease

Several mutant forms of staphylococcal nuclease with one or two defined amino acid substitutions have been purified, and the effects of the altered amino acid sequence on the stability of the folded conformation have been analyzed by guanidine hydrochloride denaturation. Two nuc- mutations, which greatly reduced the level of enzyme activity accumulated in E coli colonies carrying a recombinant plasmid with the mutant nuc gene (ie, a NUC- phenotype), both result in protein unfolding at significantly lower guanidine hydrochloride concentrations than the wild-type protein, whereas three sup mutations isolated on the basis of their ability to suppress partially the NUC- phenotype of the above two mutations result in unfolding at significantly higher guanidine hydrochloride concentrations. Characterization of nuclease molecules with two different amino acid substitutions, either nuc- + sup pairs or sup + sup pairs, suggests that the effect of an amino acid substitution on the stability of the native conformation, as measured by the value of delta delta GD, may not be a constant, but rather a variable that is sensitive to the presence of other substitutions at distant sites in the same molecule. Surprisingly, the slopes of the log Kapp vs guanidine hydrochloride concentration plots vary by as much as 35% among the different proteins.     

Nitrate reductase of Neurospora crassa: the functional role of individual amino acids in the heme domain as examined by site-directed mutagenesis

The enzyme nitrate reductase, which catalyzes the reduction of nitrate to nitrite, is a multi-redox center homodimeric protein. Each polypeptide subunit is approximately 100 kDa in size and contains three separate domains, one each for a flavin, a heme-iron, and a molybdopterin cofactor. The heme-iron domain of nitrate reductase has homology with the simple redox protein, cytochrome b5, whose crystal structure was used to predict a three-dimensional structure for the heme domain. Two histidine residues have been identified that appear to coordinate the iron of the heme moiety, while other residues may be important in the folding or the function of the heme pocket. Site-directed mutagenesis was employed to obtain mutants that encode nitrate reductase derivatives with eight different single amino acid substitutions within the heme domain, including the two central histidine residues. Replacement of one of these histidines by alanine resulted in a completely nonfunctional enzyme whereas replacement of the other histidine resulted in a stable and functional enzyme with a lower affinity for heme. Certain amino acid substitutions appeared to cause a rapid turnover of the heme domain, whereas other substitutions were tolerated and yielded a stable and fully active enzyme. Three different single amino acid replacements within the heme domain led to a dramatic change in regulation of nitrate reductase synthesis, with significant expression of the enzyme even in the absence of nitrate induction.     

Nitrate reductase of Neurospora crassa: the functional role of individual amino acids in the heme domain as examined by site-directed mutagenesis

The enzyme nitrate reductase, which catalyzes the reduction of nitrate to nitrite, is a multi-redox center homodimeric protein. Each polypeptide subunit is approximately 100 kDa in size and contains three separate domains, one each for a flavin, a heme-iron, and a molybdopterin cofactor. The heme-iron domain of nitrate reductase has homology with the simple redox protein, cytochrome b5, whose crystal structure was used to predict a three-dimensional structure for the heme domain. Two histidine residues have been identified that appear to coordinate the iron of the heme moiety, while other residues may be important in the folding or the function of the heme pocket. Site-directed mutagenesis was employed to obtain mutants that encode nitrate reductase derivatives with eight different single amino acid substitutions within the heme domain, including the two central histidine residues. Replacement of one of these histidines by alanine resulted in a completely nonfunctional enzyme whereas replacement of the other histidine resulted in a stable and functional enzyme with a lower affinity for heme. Certain amino acid substitutions appeared to cause a rapid turnover of the heme domain, whereas other substitutions were tolerated and yielded a stable and fully active enzyme. Three different single amino acid replacements within the heme domain led to a dramatic change in regulation of nitrate reductase synthesis, with significant expression of the enzyme even in the absence of nitrate induction.     

Folding defects caused by single amino acid substitutions in a subunit are not alleviated by assembly

Significant stabilization of a protein often occurs when it is assembled into an oligomer. Bacteriophage P22 contains 420 monomers of coat protein that are stabilized by the assembly and maturation processes. The effects of eight single amino acid substitutions in coat protein that each cause a temperature-sensitive-folding defect were investigated to determine if the conformational differences previously observed in the monomers could be alleviated by assembly or maturation. Several techniques including differential scanning calorimetry, heat-induced expansion, urea denaturation, and sensitivity to protease digestion were used to explore the effects of the amino acid substitutions on the conformation of coat protein, once assembled. Each of the amino acid substitutions caused a change in the conformation as compared to wild-type coat protein, observed by at least one of the probes used. Thus, neither assembly nor expansion entirely corrected the conformational defects in the monomeric subunits of the folding mutants.     

Six amino acid substitutions in the carboxyl-transferase domain of the plastidic acetyl-CoA carboxylase gene are linked with resistance to herbicides in a Lolium rigidum population

The molecular basis of an acetyl-CoA carboxylase (ACCase) target-based resistant Lolium rigidum population (WLR 96) was studied here. The carboxyl-transferase domain of the plastidic ACCase gene from resistant individuals was amplified by PCR and sequenced. The DNA sequences were aligned and compared with a susceptible population. Six amino acid substitutions were identified in the resistant population. The substitution Ile-2041-Asn, known to confer resistance to ACCase-inhibiting herbicides aryloxyphenoxypropionate (APP) in Alopecurus myosuroides, was identified in most resistant plants but it is always linked with other amino acid substitutions. This was confirmed by a cleaved amplified polymorphism (CAP) marker and an allele-specific PCR. The sole amino acid substitution Ile-2041-Asn was not found in this population. It is likely this mutation evolved later among individuals already possessing the other substitutions. Three haplotypes were identified from the resistant population based on the six amino acid combinations, and two are linked with herbicide resistance in this population. The multiple amino acid substitutions including the Ile-2041-Asn form the molecular basis endowing a high degree of resistance to ACCase-inhibiting herbicides in this L. rigidum population.     

V3 loop of the human immunodeficiency virus type 1 Env protein: interpreting sequence variability.

Two different states of human immunodeficiency virus type 1 are apparent in the asymptomatic and late stages of infection. Important determinants associated with these two states have been found within the V3 loop of the viral Env protein. In this study, two large data sets of published V3 sequences were analyzed to identify patterns of sequence variability that would correspond to these two states of the virus. We were especially interested in the pattern of basic amino acid substitutions, since the presence of basic amino acids in V3 has been shown to change virus tropism in cell culture. Four features of the sequence heterogeneity in V3 were observed: (i) approximately 70% of all nonconservative basic substitutions occur at four positions in V3, and V3 sequences with a basic substitution in at least one of these four positions contain approximately 95% of all nonconservative basic substitutions; (ii) substitution patterns within V3 are influenced by the identity of the amino acid at position 25; (iii) sequence polymorphisms account for a significant fraction of uncharged amino acid substitutions at several positions in V3, and sequence heterogeneity other than these polymorphisms is most significant at two positions near the tip of V3; and (iv) sequence heterogeneity in V3 (in addition to the basic amino acid substitutions) is approximately twofold greater in V3 sequences that contain basic amino acid substitutions. By using this sequence analysis, we were able to identify distinct groups of V3 sequences in infected patients that appear to correspond to these two virus states. The identification of these discrete sequence patterns in vivo demonstrates how the V3 sequence can be used as a genetic marker for studying the two states of human immunodeficiency virus type 1.