A universal evolutionary index for amino acid changes

Different nonsynonymous changes may be under different selective pressure during evolution. Of the 190 possible interchanges among the 20 amino acids, only 75 can be attained by a single-base substitution. An evolutionary index (EI) can be empirically computed for each of the 75 elementary changes as the likelihood of substitutions, relative to that of synonymous changes. We used 280, 1,306, 2,488, and 309 orthologous genes from primates (human versus Old World monkey), rodents (mouse versus rat), yeast (S. cerevisiae versus S. paradoxus), and Drosophila (D. melanogaster versus D. simulans), respectively, to estimate the EIs. In each data set, EI varies more than 10-fold, and the correlation coefficients of EIs from the pairwise comparisons are high (e.g., r = 0.91 between rodent and yeast). The high correlations suggest that the amino acid properties are strong determinants of protein evolution, irrespective of the identities of the proteins or the taxa of interest. However, these properties are not well captured in conventional measures of amino acid exchangeability. We, therefore, propose a universal index of exchange (U): for any large data set, its EI can be expressed as U*R, where R is the average Ka/Ks for that data set. The codon-based, empirically determined EI (i.e., U*R) makes much better predictions on protein evolution than do previous methods.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.     

Comparison of CD spectra in the aromatic region on a series of variant proteins substituted at a unique position of tryptophan synthase alpha-subunit

CD spectra in the aromatic region of a series of the mutant alpha-subunits of tryptophan synthase from Escherichia coli, substituted at position 49 buried in the interior of the molecule, were measured at pH 7.0 and 25 degrees C. These measurements were taken to gain information on conformational change produced by single amino acid substitutions. The CD spectra of the mutant proteins, substituted by Tyr or Trp residue in place of Glu residue at position 49, showed more intense positive bands due to one additional Tyr or Trp residue at position 49. The CD spectra of other mutant proteins also differed from that of the wild-type protein, despite the fact that the substituted residues at position 49 were not aromatic. Using the spectrum of the wild-type protein (Glu49) as a standard, the spectra of the other mutants were classified into three major groups. For 10 mutant proteins substituted by Ile, Ala, Leu, Met, Val, Cys, Pro, Ser, His, or Gly, their CD values of bands (due to Tyr residues) decreased in comparison with those of the wild-type protein. The mutant protein substituted by Phe also belonged to this group. These substituted amino acid residues are more hydrophobic than the original residue, Glu. In the second group, three mutant proteins were substituted by Lys, Gln, or Asn, and the CD values of tyrosyl bands increased compared to those of the wild-type proteins. These residues are polar. In the third group, the CD values of tyrosyl bands of two mutant proteins substituted by Asp or Thr were similar to those of the wild-type protein, except for one band at 276.5 nm. These results suggested that the changes in the CD spectra for the mutant proteins were affected by the hydrophobicity of the residues at position 49.     

Models of protein sequence evolution and their applications

Homologous sequences are correlated due to their common ancestry. Probabilistic models of sequence evolution are employed routinely to properly account for these phylogenetic correlations. These increasingly realistic models provide a basis for studying evolution and for exploiting it to better understand protein structure and function. Notable recent advances have been made in the treatment of insertion and deletion events, the estimation of amino-acid replacement rates, and the detection of positive selection.     

MtArt: a new model of amino acid replacement for Arthropoda

A statistical approach was applied to select those models that best fit each individual mitochondrial (mt) protein at different taxonomic levels of metazoans. The existing mitochondrial replacement matrices, MtREV and MtMam, were found to be the best-fit models for the mt-proteins of vertebrates, with the exception of Nd6, at different taxonomic levels. Remarkably, existing mitochondrial matrices generally failed to best-fit invertebrate mt-proteins. In an attempt to better model the evolution of invertebrate mt-proteins, a new replacement matrix, named MtArt, was constructed based on arthropod mt-proteomes. The new model was found to best fit almost all analyzed invertebrate mt-protein data sets. The observed pattern of model fit across the different data sets indicates that no single replacement matrix is able to describe the general evolutionary properties of mt-proteins but rather that taxonomical biases and/or the existence of different mt-genetic codes have great influence on which model is selected.     

Different versions of the Dayhoff rate matrix

Many phylogenetic inference methods are based on Markov models of sequence evolution. These are usually expressed in terms of a matrix (Q) of instantaneous rates of change but some models of amino acid replacement, most notably the PAM model of Dayhoff and colleagues, were originally published only in terms of time-dependent probability matrices (P(t)). Previously published methods for deriving Q have used eigen-decomposition of an approximation to P(t). We show that the commonly used value of t is too large to ensure convergence of the estimates of elements of Q. We describe two simpler alternative methods for deriving Q from information such as that published by Dayhoff and colleagues. Neither of these methods requires approximation or eigen-decomposition. We identify the methods used to derive various different versions of the Dayhoff model in current software, perform a comparison of existing and new implementations, and, to facilitate agreement among scientists using supposedly identical models, recommend that one of the new methods be used as a standard.     

Constructing amino acid residue substitution classes maximally indicative of local protein structure

Using an information theoretic formalism, we optimize classes of amino acid substitution to be maximally indicative of local protein structure. Our statistically-derived classes are loosely identifiable with the heuristic constructions found in previously published work. However, while these other methods provide a more rigid idealization of physicochemically constrained residue substitution, our classes provide substantially more structural information with many fewer parameters. Moreover, these substitution classes are consistent with the paradigmatic view of the sequence-to-structure relationship in globular proteins which holds that the three-dimensional architecture is predominantly determined by the arrangement of hydrophobic and polar side chains with weak constraints on the actual amino acid identities. More specific constraints are imposed on the placement of prolines, glycines, and the charged residues. These substitution classes have been used in highly accurate predictions of residue solvent accessibility. They could also be used in the identification of homologous proteins, the construction and refinement of multiple sequence alignments, and as a means of condensing and codifying the information in multiple sequence alignments for secondary structure prediction and tertiary fold recognition. © 1996 Wiley-Liss, Inc.     

Parameters of in vitroEvolution of Proteins and Peptides

In vitroevolution is used to study protein sequences, structures, and interactions and to obtain proteins with new properties. To analyze the specific features of this process in phage display experiments, we studied the amino acid composition of selected sequences, constructed a matrix of amino acid substitutions, and identified pairs of coadaptive substitutions. The amino acid frequency proved to be tightly associated with the number of corresponding codons; numerous correlated substitutions were found.     

Markovian and non-Markovian protein sequence evolution: aggregated Markov process models

Over the years, there have been claims that evolution proceeds according to systematically different processes over different timescales and that protein evolution behaves in a non-Markovian manner. On the other hand, Markov models are fundamental to many applications in evolutionary studies. Apparent non-Markovian or time-dependent behavior has been attributed to influence of the genetic code at short timescales and dominance of physicochemical properties of the amino acids at long timescales. However, any long time period is simply the accumulation of many short time periods, and it remains unclear why evolution should appear to act systematically differently across the range of timescales studied. We show that the observed time-dependent behavior can be explained qualitatively by modeling protein sequence evolution as an aggregated Markov process (AMP): a time-homogeneous Markovian substitution model observed only at the level of the amino acids encoded by the protein-coding DNA sequence. The study of AMPs sheds new light on the relationship between amino acid-level and codon-level models of sequence evolution, and our results suggest that protein evolution should be modeled at the codon level rather than using amino acid substitution models.     

An expectation maximization algorithm for training hidden substitution models

We derive an expectation maximization algorithm for maximum-likelihood training of substitution rate matrices from multiple sequence alignments. The algorithm can be used to train hidden substitution models, where the structural context of a residue is treated as a hidden variable that can evolve over time. We used the algorithm to train hidden substitution matrices on protein alignments in the Pfam database. Measuring the accuracy of multiple alignment algorithms with reference to BAliBASE (a database of structural reference alignments) our substitution matrices consistently outperform the PAM series, with the improvement steadily increasing as up to four hidden site classes are added. We discuss several applications of this algorithm in bioinformatics.     

氨基酸的置换数与生物大分子进化改变量的测度

以蛋白质分子的氨基酸置换数或核酸分子的核苷酸置换数为衡量尺度 ,说明生物大分子随时间的改变 (即分子进化速率 )保持相对恒定     

Divergence of the F1-ATP Synthase Complex in the Ice Worm, Mesenchytraeus solifugus (Annelida, Clitellata, Enchytraeidae)

The F₁-ATP synthase complex constitutes the catalytic component of F₁F₀-ATP synthase, the primary ATP synthetic enzyme in the cell. Previous studies indicate that the glacier ice worm, Mesenchytraeus solifugus, maintains unusually high ATP levels that continue to rise as temperatures decline, suggesting that molecular changes within ice worm F₁-ATP synthase subunits may contribute to this energetic anomaly. In this report, we compared ice worm F₁-ATP synthase subunits (α, β, γ) with homologues across metazoan phyla (arthropod, chordate, nematode) and among a group of clitellate annelids (Enchytraeus albidus, Enchytraeus buchholzi, Lumbriculus variegatus, Theromyzon tessulatum). Amino acid alignments indicated that ice worm F₁-ATP α and F₁-ATP β subunits share strong sequence homology with their mesophilic counterparts, respectively, but that ATP γ has diverged more rapidly. Moreover, F₁-ATP α and F₁-ATP β displayed amino acid compositional changes consistent with trends observed in other cold adapted proteins, while F₁-ATP γ diverged in unexpected directions (e.g., gains in size, charged residues). Several ice worm-specific amino acid substitutions map to positions near the F₁-ATP β catalytic site while others occur near subunit contact sites.     

The Rate of Cu,Zn Superoxide Dismutase Evolution

The rate of amino acid replacement in Cu, Zn SOD greatly departs from the expectations of the molecular clock. We examine 27 Cu, Zn SOD sequences available and conclude that: (I) the SOD enzymes from different mammal families differ from each other by roughly the same number of replacements, which is consistent with a simultaneous mammalian radiation; (2) over the most recent 60 million years (MY) the rate of SOD evolution is fairly high (15aa/100aa/100MYR) and may be considered constant; (3) the rate of accumulation of amino acid replacements since the divergence of fungi. plants and animals to the present is inconstant along different branches of the evolutionary tree; moreover it steadily decreases with time, to the same extent in all lineages; (4) some comparisons exhibit divergences that are in any case greater than expected from a Poisson process on the assumption of a molecular clock; (5) plant chloroplast enzymes display fewer differences from each other than cytoplasmic ones; (6) bacteriocuprein (from Photobacterium leiognathi), fluke and human extracellular SOD are all three extremely remotely related to one another and to the SOD of other eukaryotes.

The process of consistent decline of the rate of evolution of Cu. Zn SOD can be described by a number of mathematical functions. We explore simple models that assume constant rates and might be applicable to other proteins or genes that apparently evolve at disparate rates.     

Conformational effects of amino acid substitutions at positions 10, 12, and 13 in the P21 protein

Substitutions of amino acids for Gly 12 or Gly 13 in theras oncogene-encoded P21 proteins have been demonstrated to produce unique structural changes in these proteins that correlate with their ability to produce cell transformation. For example, the P21 proteins with Arg 12 or Val 13 are both known to be actively transforming. Recent site-specific mutagenesis experiments on the transforming Arg 12 protein have found that the substitution of Val for Gly 10 has no effect on transforming activity whereas the substitution of Val for Gly 13 led to a loss of transforming activity. In this study, we examine the structural effects of these substitutions on the amino terminal hydrophobic decapeptide (Leu 6-Gly 15) of P21 using conformational energy analysis. The results show that the transforming proteins with Gly 10 and Arg 12 or Val 10 and Arg 12 can both adopt the putative malignancy-causing conformation, whereas, for the nontransforming protein with Arg 12 and Val 13, this conformation is energetically disallowed. These results further support the theory that due to structural changes the transforming P21 proteins are unable to bind to some regulatory cellular element which may be the recently identified binding protein responsible for the induction of increased GTPase activity in normal P21 compared with transforming mutants.