Activation of insulin receptor signaling by a single amino acid substitution in the transmembrane domain.

The insulin receptor is a ligand-activated tyrosine kinase composed of two alpha and two beta subunits. A single transmembrane domain composed of 23 hydrophobic residues is contained in each beta subunit. We examined the role of the transmembrane domain in regulating insulin receptor signaling by inserting a negatively charged amino acid (Asp) for Val938 (V938D). Chinese hamster ovary (CHO) cells were stably transfected with a plasmid containing both the neomycin-resistance gene and either the wild-type or the mutant (V938D) insulin receptor cDNA. Insulin binding increased similarly in CHO cells stably transfected with the wild-type and the V938D-mutant insulin receptor cDNA. Insulin stimulated glucose transport and cell growth in cells expressing the normal insulin receptor. By contrast, in the absence of insulin, glucose transport and cell growth in CHO-V938D cells were as high as in insulin-stimulated control cells and no longer responsive to insulin stimulation. Phosphorylation of the beta subunit of the insulin receptor was also increased in CHO-V938D cells not exposed to insulin. These results support an essential role of the transmembrane domain of the insulin receptor in the transduction of insulin signaling.     

Transferrin DChi: The amino acid substitution

A peptide difference has been found in the neutral-band (pH 6.4) regions of tryptic digests of human transferrins C and D_Chi. The peptide has the composition Asp-Ser-Ala-Arg. Therefore, this peptide is proposed as the 2TD_Chi ^b peptide, the result of the replacement of histidine by arginine in the Tf C peptide.Supported in part by U.S. Public Health Service grants GM 09326, 5-K3 GM 18,381, and GM 00337 from the National Institutes of Health.     

A metric model of amino acid substitution

MOTIVATION: We address the question of whether there exists an effective evolutionary model of amino-acid substitution that forms a metric-distance function. There is always a trade-off between speed and sensitivity among competing computational methods of determining sequence homology. A metric model of evolution is a prerequisite for the development of an entire class of fast sequence analysis algorithms that are both scalable, O(log n) and sensitive. RESULTS: We have reworked the mathematics of the point accepted mutation model (PAM) by calculating the expected time between accepted mutations in lieu of calculating log-odds probabilities. The resulting substitution matrix (mPAM) forms a metric. We validate the application of the mPAM evolutionary model for sequence homology by executing sequence queries from a controlled yeast protein homology search benchmark. We compare the accuracy of the results of mPAM and PAM similarity matrices as well as three prior metric models. The experiment shows that mPAM significantly outperforms the other three metrics and sufficiently approaches the sensitivity of PAM250 to make it applicable to the management of protein sequence databases.     

Performance evaluation of amino acid substitution matrices

Several choices of amino acid substitution matrices are currently available for searching and alignment applications. These choices were evaluated using the BLAST searching program, which is extremely sensitive to differences among matrices, and the Prosite catalog, which lists members of hundreds of protein families. Matrices derived directly from either sequence-based or structurebased alignments of distantly related proteins performed much better overall than extrapolated matrices based on the Dayhoff evolutionary model. Similar results were obtained with the FASTA searching program. Improved performance appears to be general rather than family-specific, reflecting improved accuracy in scoring alignments. An implementation of a multiple matrix strategy was also tested. While no combination of three matrices performed as well as the single best matrix, BLOSUM 62, good results were obtained using a combination of sequence-based and structure-based matrices. This hybrid set of matrices is likely to be useful in certain situations. Our results illustrate the importance of matrix selection and value of a comprehensive approach to evaluation of protein comparison tools. © 1993 Wiley-Liss, Inc.     

Bayesian analysis of amino acid substitution models

Models of amino acid substitution present challenges beyond those often faced with the analysis of DNA sequences. The alignments of amino acid sequences are often small, whereas the number of parameters to be estimated is potentially large when compared with the number of free parameters for nucleotide substitution models. Most approaches to the analysis of amino acid alignments have focused on the use of fixed amino acid models in which all of the potentially free parameters are fixed to values estimated from a large number of sequences. Often, these fixed amino acid models are specific to a gene or taxonomic group (e.g. the Mtmam model, which has parameters that are specific to mammalian mitochondrial gene sequences). Although the fixed amino acid models succeed in reducing the number of free parameters to be estimated--indeed, they reduce the number of free parameters from approximately 200 to 0--it is possible that none of the currently available fixed amino acid models is appropriate for a specific alignment. Here, we present four approaches to the analysis of amino acid sequences. First, we explore the use of a general time reversible model of amino acid substitution using a Dirichlet prior probability distribution on the 190 exchangeability parameters. Second, we then explore the behaviour of prior probability distributions that are'centred' on the rates specified by the fixed amino acid model. Third, we consider a mixture of fixed amino acid models. Finally, we consider constraints on the exchangeability parameters as partitions,similar to how nucleotide substitution models are specified, and place a Dirichlet process prior model on all the possible partitioning schemes.     

Molecular conversion of NAD kinase to NADH kinase through single amino acid residue substitution

NAD kinase phosphorylates NAD+ to form NADP+ and is strictly specific to NAD+, whereas NADH kinase phosphorylates both NAD+ and NADH, thereby showing relaxed substrate specificity. Based on their primary and tertiary structures, the difference in the substrate specificities between NAD and NADH kinases was proposed to be caused by one aligned residue: Gly or polar amino acid (Gln or Thr) in five NADH kinases and a charged amino acid (Arg) in two NAD kinases. The substitution of Arg with Gly in the two NAD kinases relaxed the substrate specificity (i.e. converted the NAD kinases to NADH kinases). The substitution of Arg in one NAD kinase with polar amino acids also relaxed the substrate specificity, whereas substitution with charged and hydrophobic amino acids did not show a similar result. In contrast, the substitution of Gly with Arg in one NADH kinase failed to convert it to NAD kinase. These results suggest that a charged or hydrophobic amino acid residue in the position of interest is crucial for strict specificity of NAD kinases to NAD+, whereas Gly or polar amino acid residue is not the sole determinant for the relaxed substrate specificity of NADH kinases. The significance of the conservation of the residue at the position in 207 NAD kinase homologues is also discussed.     

Improved amino acid flexibility parameters

Protein molecules exhibit varying degrees of flexibility throughout their three-dimensional structures, with some segments showing little mobility while others may be so disordered as to be unresolvable by techniques such as X-ray crystallography. Atomic displacement parameters, or B-factors, from X-ray crystallographic studies give an experimentally determined indication of the degree of mobility in a protein structure. To provide better estimators of amino acid flexibility, we have examined B-factors from a large set of high-resolution crystal structures. Because of the differences among structures, it is necessary to normalize the B-factors. However, many proteins have segments of unusually high mobility, which must be accounted for before normalization can be performed. Accordingly, a median-based method from quality control studies was used to identify outliers. After removal of outliers from, and normalization of, each protein chain, the B-factors were collected for each amino acid in the set. It was found that the distribution of normalized B-factors followed a Gumbel, or extreme value distribution, and the location parameter, or mode, of this distribution was used as an estimator of flexibility for the amino acid. These new parameters have a higher correlation with experimentally determined B-factors than parameters from earlier methods.     

Switching DNA-binding specificity by unnatural amino acid substitution

The specificity of protein–nucleic acid recognition is believed to originate largely from hydrogen bonding between protein polar atoms, primarily side-chain and polar atoms of nucleic acid bases. One way to design new nucleic acid binding proteins of novel specificity is by structure-guided alterations of the hydrogen bonding patterns of a nucleic acid–protein complex. We have used cI repressor of bacteriophage λ as a model system. In the λ-repressor–DNA complex, the -NH₂ group (hydrogen bond donor) of lysine-4 of λ-repressor forms hydrogen bonds with the amide carbonyl atom of asparagine-55 (acceptor) and the O6 (acceptor) of CG6 of operator site O_L1. Substitution of lysine-4 (two donors) by iso-steric S-(2-hydroxyethyl)-cysteine (one donor and one acceptor), by site-directed mutagenesis and chemical modification, leads to switch of binding specificity of λ-repressor from C:G to T:A at position 6 of O_L1. This suggests that unnatural amino acid substitutions could be a simple way of generating nucleic acid binding proteins of altered specificity.     

Higher rates of amino acid substitution in rodents than in humans.

An analysis of 54 protein sequences from humans and rodents (mice or rats), with the chicken as an outgroup, indicates that, from the common ancestor of primates and rodents, 35 of the proteins have evolved faster in the lineage to mouse or rat (rodent lineage) whereas only 12 proteins have evolved faster in the lineage to humans (human lineage). The average rate of amino acid substitution is significantly faster in the rodent lineage than in the human lineage. In addition, the average rate of insertion/deletion is also faster in rodents than in humans and there is a positive correlation between the rate of amino acid substitution and the rate of insertion/deletion in a protein sequence.     

Estimating amino acid substitution models for metazoan evolutionary studies

Amino acid substitution models represent the substitution rates among amino acids during the evolution of protein sequences. The models are a prerequisite for maximum likelihood or Bayesian methods to analyse the phylogenetic relationships among species based on their protein sequences. Estimating amino acid substitution models requires large protein datasets and intensive computation. In this paper, we presented the estimation of both time-reversible model (Q.met) and time non-reversible model (NQ.met) for multicellular animals (Metazoa). Analyses showed that the Q.met and NQ.met models were significantly better than existing models in analysing metazoan protein sequences. Moreover, the time non-reversible model NQ.met enables us to reconstruct the rooted phylogenetic tree for Metazoa. We recommend researchers to employ the Q.met and NQ.met models in analysing metazoan protein sequences.     

Ideal amino acid exchange forms for approximating substitution matrices

We have analyzed 29 published substitution matrices (SMs) and five statistical protein contact potentials (CPs) for comparison. We find that popular, 'classical' SMs obtained mainly from sequence alignments of globular proteins are mostly correlated by at least a value of 0.9. The BLOSUM62 is the central element of this group. A second group includes SMs derived from alignments of remote homologs or transmembrane proteins. These matrices correlate better with classical SMs (0.8) than among themselves (0.7). A third group consists of intermediate links between SMs and CPs - matrices and potentials that exhibit mutual correlations of at least 0.8. Next, we show that SMs can be approximated with a correlation of 0.9 by expressions c(0) + x(i)x(j) + y(i)y(j) + z(i)z(j), 1相似文献   

Genome bias influences amino acid choices: analysis of amino acid substitution and re-compilation of substitution matrices exclusive to an AT-biased genome

The genomic era has seen a remarkable increase in the number of genomes being sequenced and annotated. Nonetheless, annotation remains a serious challenge for compositionally biased genomes. For the preliminary annotation, popular nucleotide and protein comparison methods such as BLAST are widely employed. These methods make use of matrices to score alignments such as the amino acid substitution matrices. Since a nucleotide bias leads to an overall bias in the amino acid composition of proteins, it is possible that a genome with nucleotide bias may have introduced atypical amino acid substitutions in its proteome. Consequently, standard matrices fail to perform well in sequence analysis of these genomes. To address this issue, we examined the amino acid substitution in the AT-rich genome of Plasmodium falciparum, chosen as a reference and reconstituted a substitution matrix in the genome's context. The matrix was used to generate protein sequence alignments for the parasite proteins that improved across the functional regions. We attribute this to the consistency that may have been achieved amid the target and background frequencies calculated exclusively in our study. This study has important implications on annotation of proteins that are of experimental interest but give poor sequence alignments with standard conventional matrices.     

A single amino acid substitution of Leu130Ile in snake DNases I contributes to the acquisition of thermal stability.

We purified pancreatic deoxyribonucleases I (DNases I) from three snakes, Elaphe quadrivirgata, Elaphe climacophora and Agkistrodon blomhoffii, and cloned their cDNAs. Each mature snake DNase I protein comprised 262 amino acids. Wild-type snake DNases I with Leu130 were more thermally unstable than wild-type mammalian and avian DNases I with Ile130. After substitution of Leu130Ile, the thermal stabilities of the snake enzymes were higher than those of their wild-type counterparts and similar to mammalian wild-type enzyme levels. Conversely, substituting Ile130Leu of mammalian DNases I made them more thermally unstable than their wild-type counterparts. Therefore, a single amino acid substitution, Leu130Ile, might be involved in an evolutionally critical change in the thermal stabilities of vertebrate DNases I. Amphibian DNases I have a Ser205 insertion in a Ca2+-binding site of mammalian and avian enzymes that reduces their thermal stabilities [Takeshita, H., Yasuda, T., Iida, R., Nakajima, T., Mori, S., Mogi, K., Kaneko, Y. & Kishi, K. (2001) Biochem. J.357, 473-480]. Thus, it is plausible that the thermally stable wild-type DNases I of the higher vertebrates, such as mammals and birds, have been generated by a single Leu130Ile substitution of reptilian enzymes through molecular evolution following Ser205 deletion from amphibian enzymes. This mechanism may reflect one of the evolutionary changes from cold-blooded to warm-blooded vertebrates.     

An amino acid substitution matrix for protein conformation identification

Amino acid substitution matrices play an essential role in protein sequence alignment, a fundamental task in bioinformatics. Most widely used matrices, such as PAM matrices derived from homologous sequences and BLOSUM matrices derived from aligned segments of PROSITE, did not integrate conformation information in their construction. There are a few structure-based matrices, which are derived from limited data of structure alignment. Using databases PDB_SELECT and DSSP, we create a database of sequence-conformation blocks which explicitly represent sequence-structure relationship. Members in a block are identical in conformation and are highly similar in sequence. From this block database, we derive a conformation-specific amino acid substitution matrix CBSM60. The matrix shows an improved performance in conformational segment search and homolog detection.     

Effect of amino acid substitution by sited-directed mutagenesis on the carbohydrate recognition and stability of human 14-kDa beta-galactoside-binding lectin.

The roles of selected amino acid residues of human 14-kDa beta-galactoside-binding lectin were studied by site-directed mutagenesis. Ten mutant lectin proteins were produced, in each of which one of the residues regarded as possibly related to the stability of the lectin (6 cysteine residues) or one of those highly conserved in the vertebrate beta-galactoside-binding lectin family (Asn46, Trp68, Glu71, and Arg73), was substituted. All the mutant lectins in which one of the cysteine residues had been substituted with serine (C2S, C16S, C42S, C60S, C88S, and C130S) proved to have sugar binding ability comparable with that of the wild-type lectin. In addition, one of the mutants in which Cys2 was substituted (C2S) was found to have become considerably more stable under non-reducing conditions. It retained asialofetuin binding activity for over a week in the absence of beta-mercaptoethanol, while the wild-type lectin lost it within a day. This suggests that oxidation of Cys2 could be a key process in the inactivation of human 14-kDa lectin. Substitution of highly conservative Trp68 to tyrosine (W68Y) slightly reduced lactose binding ability, but the mutant was still adsorbed strongly on asialofetuin-agarose. Other mutant lectins in which conservative hydrophilic amino acids were substituted (N46D, E71Q, and R73H) failed to bind to the asialofetuin agarose, with no sign of retardation. Thus, conservative hydrophilic residues proved to be more important in carbohydrate recognition than the cysteine and tryptophan residues, contrary to the widely accepted concept that these latter residues are essential.     

Nature of amino acid side chain and alpha-helix stability.

In order to investigate the ability of neutral amino acids to support the α-helix conformation, the coil–helix transition of poly(L -lysine) and of lysine copolymers with these amino acids was studied in water/methanol using circular dichroism. The transtions were recorded at constant pH adding buffer to the methanol/water mixtures. With poly(L -lysine), experiments were performed at several constant pH's; the transition midpoint on the water (methanol) concentration scale was found to depend strongly upon pH; the helix stability region is shifted towards higher water concentrations, when the pH is increased. Copolymers of lysine and several neutral amino acids revealed the same effect in that increasing amounts of, for example, norleucine also shifted the transition midpoint to higher water concentrations. A series of copolymers containing L -lysine as the host and different hydrophobic amino acids were synthesized and the helix–coil transition in water/methanol was observed at constant pH. Different copolymers of equal composition showed significant differences with respect to the nature of the amino acid incorporated into polylysine. From these studies an α-helix-philic scale (in decreasing order): Leu, Nle, Ile, Ala, Phe, Val, Gly is deduced and discussed; the results obtained were compared with those of different procedures.     

A single amino acid substitution in SecY stabilizes the interaction with SecA.

The SecYEG complex constitutes a protein conducting channel across the bacterial cytoplasmic membrane. It binds the peripheral ATPase SecA to form the translocase. When isoleucine 278 in transmembrane segment 7 of the SecY subunit was replaced by a unique cysteine, SecYEG supported an increased preprotein translocation and SecA translocation ATPase activity, and allowed translocation of a preprotein with a defective signal sequence. SecY(I278C)EG binds SecA with a higher affinity than normal SecYEG, in particular in the presence of ATP. The increased translocation activity of SecY(I278C)EG was confirmed in a purified system consisting of SecYEG proteoliposomes, while immunoprecipitation in detergent solution reveal that translocase-preprotein complexes are more stable with SecY(I278C) than with normal SecY. These data imply an important role for SecY transmembrane segment 7 in SecA binding. As improved SecA binding to SecY was also observed with the prlA4 suppressor mutation, it may be a general mechanism underlying signal sequence suppression.     

The genomic rate of adaptive amino acid substitution in Drosophila

The proportion of amino acid substitutions driven by adaptive evolution can potentially be estimated from polymorphism and divergence data by an extension of the McDonald-Kreitman test. We have developed a maximum-likelihood method to do this and have applied our method to several data sets from three Drosophila species: D. melanogaster, D. simulans, and D. yakuba. The estimated number of adaptive substitutions per codon is not uniformly distributed among genes, but follows a leptokurtic distribution. However, the proportion of amino acid substitutions fixed by adaptive evolution seems to be remarkably constant across the genome (i.e., the proportion of amino acid substitutions that are adaptive appears to be the same in fast-evolving and slow-evolving genes; fast-evolving genes have higher numbers of both adaptive and neutral substitutions). Our estimates do not seem to be significantly biased by selection on synonymous codon use or by the assumption of independence among sites. Nevertheless, an accurate estimate is hampered by the existence of slightly deleterious mutations and variations in effective population size. The analysis of several Drosophila data sets suggests that approximately 25% +/- 20% of amino acid substitutions were driven by positive selection in the divergence between D. simulans and D. yakuba.