Advantages of a mechanistic codon substitution model for evolutionary analysis of protein-coding sequences

BACKGROUND: A mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated. RESULTS: The present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.     

ATP and other nucleotide interaction with model compounds of amino acids.

With the use of a new liquid scintillation technique, the partitioning and interfacial interaction of ³H-ATP and other adenine nucleotides to various amphiphiles has been studied. The amphiphiles consisted of lipophilic groups attached to amino acid residues or certain functional groups and were as follows: trimethylindole, octadecylamine, octadecyl alcohol, dodecylguanidine, stearylamide, acetododecyl amide, p-nonyl phenol, stearyl mercaputan, stearic acid, and imidazole. A toluene solution containing the lipophilic substance and fluor molecules was layered onto aqueous solutions containing a ³H-nucleotide, and the interaction followed by the liquid scintillation technique. A strong interaction with the nucleotides occurred with dodecylguanidine, octadecylamine, and trimethylindole, while the interaction with the other amphiphiles was slight or absent. The time course of the interaction followed first order kinetics irrespective of the number of reactive species; the rate of diffusion of the nucleotide being the rate-determining step. The technique permitted a simple analysis of formation constants; which for ATP were 1–2 orders of magnitude greater than for ADP and AMP. In the presence of Ca²⁺, but not Mg²⁺, an interaction occurred between ATP and either stearic acid, stearyl mercaputan, and p-nonyl phenol. On the other hand, in the presence of Mg²⁺, but not Ca²⁺, an interaction occurred between stearic acid and polyadenylic acid. The results are discussed in relation to ATP interactions with proteins and lipids.     

A nucleotide substitution model with nearest-neighbour interactions

MOTIVATION: It is well known that neighbouring nucleotides in DNA sequences do not mutate independently of each other. In this paper, we introduce a context-dependent substitution model and derive an algorithm to calculate the likelihood of sequences evolving under this model. We use this algorithm to estimate neighbour-dependent substitution rates, as well as rates for dinucleotide substitutions, using a Bayesian sampling procedure. The model is irreversible, giving an arrow to time, and allowing the position of the root between a pair of sequences to be inferred without using out-groups. RESULTS: We applied the model upon aligned human-mouse non-coding data. Clear neighbour dependencies were observed, including 17-18-fold increased CpG to TpG/CpA rates compared with other substitutions. Root inference positioned the root halfway the mouse and human tips, suggesting an approximately clock-like behaviour of the irreversible part of the substitution process.     

Selective separation of amino acids by reactive extraction

The method of reactive extraction with di-(2-ethylhexyl)phosphoric acid (D2EHPA) for the separation of a range of amino acids is studied. The results obtained on the individual reactive extraction indicated the possibility of the amino acids selective separation as a function of the pH value of aqueous solution and the acidic or basic character of each amino acid. Thus, using multistage extraction, the total separation of the following amino acids groups has been performed: neutral amino acids (l-glycine, l-alanine, l-tryptophan) at pH 5–5.5 (nine extraction stages), basic amino acids (l-lysine, l-arginine) and l-cysteine at pH 4–4.5 (ten extraction stages), l-histidine at pH 3–3.5 (five extraction stages), and acidic amino acids (l-aspartic acid, l-glutamic acid) at pH 2–2.5 (three extraction stages).     

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

Position-specific incorporation of dansylated non-natural amino acids into streptavidin by using a four-base codon

Novel non-natural amino acids carrying a dansyl fluorescent group were designed, synthesized, and incorporated into various positions of streptavidin by using a CGGG four-base codon in an Escherichia coli in vitro translation system. 2,6-Dansyl-aminophenylalanine (2,6-dnsAF) was found to be incorporated into the protein more efficiently than 1,5-dansyl-lysine, 2,6-dansyl-lysine, and 1,5-dansyl-aminophenylalanine. Fluorescence measurements indicate that the position-specific incorporation of the 2,6-dnsAF is a useful technique to probe protein structures. These results also indicate that well-designed non-natural amino acids carrying relatively large side chains can be accepted as substrates of the translation system.     

Four-base codon mediated mRNA display to construct peptide libraries that contain multiple nonnatural amino acids

In vitro selection and directed evolution of peptides from mRNA display are powerful strategies to find novel peptide ligands that bind to target biomolecules. In this study, we expanded the mRNA display method to include multiple nonnatural amino acids by introducing three different four-base codons at a randomly selected single position on the mRNA. Another nonnatural amino acid may be introduced by suppressing an amber codon that may appear from a (NNK)_n nucleotide sequence on the mRNA. The mRNA display was expressed in an Escherichia coli in vitro translation system in the presence of three types of tRNAs carrying different four-base anticodons and a tRNA carrying an amber anticodon, the tRNAs being chemically aminoacylated with different nonnatural amino acids. The complexity of the starting mRNA-displayed peptide library was estimated to be 1.1 × 10¹² molecules. The effectiveness of the four-base codon mediated mRNA display method was demonstrated in the selection of biocytin-containing peptides on streptavidin-coated beads. Moreover, a novel streptavidin-binding nonnatural peptide containing benzoylphenylalanine was obtained from the nonnatural peptide library. The nonnatural peptide library from the four-base codon mediated mRNA display provides much wider functional and structural diversity than conventional peptide libraries that are constituted from 20 naturally occurring amino acids.     

Switch from an aquaporin to a glycerol channel by two amino acids substitution

The MIP (major intrinsic protein) proteins constitute a channel family of currently 150 members that have been identified in cell membranes of organisms ranging from bacteria to man. Among these proteins, two functionally distinct subgroups are characterized: aquaporins that allow specific water transfer and glycerol channels that are involved in glycerol and small neutral solutes transport. Since the flow of small molecules across cell membranes is vital for every living organism, the study of such proteins is of particular interest. For instance, aquaporins located in kidney cell membranes are responsible for reabsorption of 150 liters of water/day in adult human. To understand the molecular mechanisms of solute transport specificity, we analyzed mutant aquaporins in which highly conserved residues have been substituted by amino acids located at the same positions in glycerol channels. Here, we show that substitution of a tyrosine and a tryptophan by a proline and a leucine, respectively, in the sixth transmembrane helix of an aquaporin leads to a switch in the selectivity of the channel, from water to glycerol.     

Removal of kinetic traps and enhanced protein folding by strategic substitution of amino acids in a model alpha-helical hairpin peptide

The presence of non-native kinetic traps in the free energy landscape of a protein may significantly lengthen the overall folding time so that the folding process becomes unreliable. We use a computational model alpha-helical hairpin peptide to calculate structural free energy landscapes and relate them to the kinetics of folding. We show how protein engineering through strategic changes in only a few amino acid residues along the primary sequence can greatly increase the speed and reliability of the folding process, as seen experimentally. These strategic substitutions also prevent the formation of long-lived misfolded configurations that can cause unwanted aggregations of peptides. These results support arguments that removal of kinetic traps, obligatory or nonobligatory, is crucial for fast folding.     

Free energy changes associated with amino acid substitution in proteins

The estimation of free energy differences from computer simulation of macromolecular systems is important for rational strategies for drug design and for protein engineering. As an example of one mutation, we have studied the free energy change resulting from the conversion of a polar group (OH) to an apolar group (CH3) in aqueous solution. We have estimated the effect of various local environments on the magnitude of the free energy difference and find that significant environmental effects are found. We have also studied the reliability of the results in detail.     

Amino acid and codon usage profiles: adaptive changes in the frequency of amino acids and codons

In the work presented, the changes in codon and amino acid contents have been studied as a function of environmental conditions by comparing pairs of homologs in a group of extremophilic/non-extremophilic genomes. Our results obtained based on such analysis highlights a number of notable observations: (i) the overall preference of amino acid usages in the proteins of a given organism is significantly affected by major environmental factors. The changes in amino acid preferences (amino acid usage profiles) in an extremophile compared to its non-extremophile relative recurs in the organisms of similar extreme habitats. (ii) On the other hand, changes in codon usage preferences in these extremophilic/non-extremophilic pairs, lack such persistency not only in different genome-pairs but also in the individual genes of a specific pair. (iii) We have noted a correlation between cellular function and codon usage profiles of the genes in the studied pairs. (iv) Based on this correlation, we could obtain a decent prediction of cellular functions solely based on codon usage profile data. (v) Comparisons made between two sets of randomly generated genomes suggest that different patterns of codon usage changes in genes of different functional categories result in a partial resistance towards the changes in the concentration of a given amino acid. This buffering capacity might explain the observed differences in codon usage trends in genes of different functions. In the end, we suggest codon usage and amino acid profiles as powerful tools that can be utilized to improve function predictions and genome-environment mappings.     

Selective constraints in yeast genes with differential expressivity: codon pair usage and mRNA stability perspectives

Protein translation has been elucidated to be dictated by evolutionary constraints, namely, variations in tRNA availabilities and/or variations in codon-anticodon binding that is manifested in biased codon usage. Taking advantage of publicly available mRNA expression and protein abundance data for Saccharomyces cerevisiae, we have performed a comprehensive analysis of the diverse factors guiding translation leading to desired protein levels irrespective of the corresponding high or low mRNA levels. It has been elucidated in this study that different combinations of most abundant/non abundant tRNA isoacceptors are selected for in S. cerevisiae that helps in achieving the optimum speed and accuracy in the protein translation process. This is also accompanied by the strategic location of codon pairs in coherence to mRNA secondary structure folding stability for the above mentioned combinations of tRNA isoacceptors. We thus find that codon pair contextual effects; in addition to tRNA abundance and mRNA folding stability during translation elongation process play plausible roles in maintaining translation accuracy and speed that can achieve desired protein levels.     

The ''polysemous'' codon--a codon with multiple amino acid assignment caused by dual specificity of tRNA identity.

In some Candida species, the universal CUG leucine codon is translated as serine. However, in most cases, the serine tRNAs responsible for this non-universal decoding (tRNA(Ser)CAG) accept in vitro not only serine, but also, to some extent, leucine. Nucleotide replacement experiments indicated that m1G37 is critical for leucylation activity. This finding was supported by the fact that the tRNA(Ser)CAGs possessing the leucylation activity always have m1G37, whereas that of Candida cylindracea, which possesses no leucylation activity, has A37. Quantification of defined aminoacetylated tRNAs in cells demonstrated that 3% of the tRNA(Ser)CAGs possessing m1G37 were, in fact, charged with leucine in vivo. A genetic approach using an auxotroph mutant of C.maltosa possessing this type of tRNA(Ser)CAG also suggested that the URA3 gene inactivated due to the translation of CUG as serine was rescued by a slight incorporation of leucine into the polypeptide, which demonstrated that the tRNA charged with multiple amino acids could participate in the translation. These findings provide the first evidence that two distinct amino acids are assigned by a single codon, which occurs naturally in the translation process of certain Candida species. We term this novel type of codon a 'polysemous codon'.     

Genome-wide patterns of nucleotide substitution reveal stringent functional constraints on the protein sequences of thermophiles

To test the hypothesis that the proteins of thermophilic prokaryotes are subject to unusually stringent functional constraints, we estimated the numbers of synonymous and nonsynonymous nucleotide substitutions per site between 17,957 pairs of orthologous genes from 22 pairs of closely related species of Archaea and Bacteria. The average ratio of nonsynonymous to synonymous substitutions was significantly lower in thermophiles than in nonthermophiles, and this effect was observed in both Archaea and Bacteria. There was no evidence that this difference could be explained by factors such as nucleotide content bias. Rather, the results support the hypothesis that proteins of thermophiles are subject to unusually strong purifying selection, leading to a reduced overall level of amino acid evolution per mutational event. The results show that genome-wide patterns of sequence evolution can be influenced by natural selection exerted by a species' environment and shed light on a previous observation that relatively few of the mutations arising in a thermophilic archaeon were nucleotide substitutions in contrast to indels.     

A similarity ring for amino acids based on their evolutionary substitution rates

A similarity ring for amino acids is obtained by reordering the rows and columns of their substitution probability matrix so as to minimize the variance of the matric elements about the main diagonal.     

Conversion of monocyte chemoattractant protein-1 into a neutrophil attractant by substitution of two amino acids.

The small cytokine monocyte chemoattractant protein-1 has structural similarity to the neutrophil chemoattractant interleukin-8, but each protein is specific in attracting its own target cell. To investigate the structural basis of this cell type specificity, we have developed an Escherichia coli expression system for the monocyte chemoattractant and mutagenized selected amino acid residues to ones found at the corresponding positions of interleukin-8. We find that a double mutation of tyrosine 28 and arginine 30 to leucine and valine, respectively, causes a drastic decrease in chemotactic activity toward monocytes with the appearance of a novel (interleukin-8-like) neutrophil chemotactic activity. Computer graphic analysis predicts that, with the double substitution, a putative receptor binding groove of the monocyte chemoattractant protein would become topographically similar to that of interleukin-8. We therefore postulate that one or both of these amino acid residues are part of the binding contact of these small cytokines and their receptors.     

Studies on the active sites ofBacillus cereus sphingomyelinase substitution of some amino acids by site-directed mutagenesis

Summary Chemical modifications suggested that acidic amino acids such as aspartic and glutamic acids are involved in the active sites ofBacillus cereus sphingomyelinase. Among aspartic acid residues in the conserved regions of this enzyme, Asp-126, Asp-156, Asp-233 and Asp-295 were converted to glycine by site-directed mutagenesis. According to prediction on structural similarity to pancreatic DNase I, His-151 and His-296 were also converted to alanine. The Asp and His mutants, D126G, D156G, D233G, D295G, H151A and H296A, were produced inBacillus brevis 47, a protein-hyperproducing strain. The catalytic activities of D295G, H151A and H296A were completely abolished, and sphingomyelin-hydrolyzing activity of D126G or D156G was reduced by more than 50%. The activity of D126G towardp-NPPC was comparable to that of the wild-type, while D156G catalyzed the hydrolysis of HNP andp-NPPC more efficiently than the wild-type. Hemolytic activities of the mutants were parallel to their sphingomyelin-hydrolyzing activities.     

Four-base codon/anticodon strategy and non-enzymatic aminoacylation for protein engineering with non-natural amino acids

Techniques for position-specific incorporation of non-natural amino acids in an in vitro protein synthesizing system are described. First, a PNA-assisted non-enzymatic tRNA aminoacylation with a variety of natural and non-natural amino acids is described. With this technique, one can aminoacylate a specific tRNA simply by adding a preformed amino acid activated ester-PNA conjugate into an in vitro protein biosynthesizing system. Second, the genetic code is expanded by introducing 4-base codons that can be exclusively translated to non-natural amino acids. The most advantageous point of the 4-base codon strategy is to introduce multiple amino acids into specific positions in single proteins by using mutually orthogonal 4-base codons and orthogonal tRNAs. An easy and quick method for preparation of tRNAs possessing 4-base anticodons is also described. Combination of the non-enzymatic aminoacylation and the 4-base codon/anticodon strategy gives an easy and widely applicable technique for incorporating a variety of non-natural amino acids into proteins in vitro.