Identification of amino acid substitutions that alter the substrate specificity of TEM-1 beta-lactamase.

TEM-1 beta-lactamase is the most prevalent plasmid-mediated beta-lactamase in gram-negative bacteria. Recently, TEM beta-lactamase variants with amino acid substitutions in the active-site pocket of the enzyme have been identified in natural isolates with increased resistance to extended-spectrum cephalosporins. To identify other amino acid substitutions that alter the activity of TEM-1 towards extended-spectrum cephalosporins, we probed regions around the active-site pocket by random-replacement mutagenesis. This mutagenesis technique involves randomizing the DNA sequence of three to six codons in the blaTEM-1 gene to form a library containing all or nearly all of the possible substitutions for the region randomized. In total, 20 different residue positions that had been randomized were screened for amino acid substitutions that increased enzyme activity towards the extended-spectrum cephalosporin cefotaxime. Substitutions at positions 104, 168, and 238 in the TEM-1 beta-lactamase that resulted in increased enzyme activity towards extended-spectrum cephalosporins were found. In addition, small deletions in the loop containing residues 166 to 170 drastically altered the substrate specificity of the enzyme by increasing activity towards extended-spectrum cephalosporins while virtually eliminating activity towards ampicillin.     

The hydrophobic core of Escherichia coli thioredoxin shows a high tolerance to nonconservative single amino acid substitutions.

A set of single amino acid substitutions has been constructed at positions Leu42 and Leu78 in the hydrophobic core of Escherichia coli thioredoxin. This protein is required for the in vivo assembly of filamentous bacteriophages such as M13. Almost all the mutants retain this activity regardless of the change in size, hydrophobic nature, or charge of the substitution. Determination of the free energies of unfolding of the mutants containing charged residues shows that these are significantly destabilized as would be expected from simple considerations of the hydrophobic effect. Thioredoxin therefore represents a class of proteins where the often observed correlation between a particular biological activity and thermodynamic stability is not evident for single mutants in the all-or-none assay used. Native thioredoxin is very stable. Thus, structurally single mutants may not perturb the folding equilibrium or the dynamic behavior sufficiently for the effects to be sensed in vivo.     

Tracing hybrid incompatibilities to single amino acid substitutions

Deleterious interactions among genes cause reductions in fitness of interpopulation hybrids (hybrid breakdown). Identifying genes involved in hybrid breakdown has proven difficult, and few studies have addressed the molecular basis of this widespread phenomenon. Because proper function of the mitochondrial electron transport system (ETS) requires a coadapted set of nuclear and mitochondrial gene products, ETS genes present an attractive system for studying the evolution of coadapted gene complexes within isolated populations and the loss of fitness in interpopulation hybrids. Here we show the effects of single amino acid substitutions in cytochrome c (CYC) on its functional interaction with another ETS protein, cytochrome c oxidase (COX) in the intertidal copepod Tigriopus californicus. The individual and pairwise consequences of three naturally occurring amino acid substitutions in CYC are examined by site-directed mutagenesis and found to differentially effect the rates of CYC oxidation by COX variants from different source populations. In one case, we show that interpopulation hybrid breakdown in COX activity can be attributed to a single naturally occurring amino acid substitution in CYC.     

Combinatorial codon-based amino acid substitutions

Twenty Fmoc-protected trinucleotide phosphoramidites representing a complete set of codons for the natural amino acids were chemically synthesized for the first time. A pool of these reagents was incorporated into oligonucleotides at substoichiometric levels to generate two libraries of variants that randomly carry either few or many codon replacements on a region encoding nine amino acids of the bacterial enzyme TEM-1 β-lactamase. Assembly of the libraries was performed in a completely automated mode through a simple modification of ordinary protocols. This technology eliminates codon redundancy, stop codons and enables complete exploration of sequence space for single, double and triple mutations throughout a protein region spanning several residues. Sequence analysis of many non-selected clones revealed a good incorporation of the trinucleotides, producing combinations of mutations quite different from those obtained using conventional degenerate oligonucleotides. Ceftazidime-selection experiments yielded several never before reported variants containing novel amino acid combinations in the β-lactamase omega loop region.     

Multiple amino acid substitutions allow DNA polymerases to synthesize RNA

DNA and RNA polymerase exhibit similarities in structures and catalytic mechanisms, suggesting that both classes of enzymes are evolutionarily related. To probe the biochemical and structure-function relationship between the two classes of polymerases, a large library (200,000 members) of mutant Thermus aquaticus DNA polymerase I (Taq pol I) was created containing random substitutions within a portion of the dNTP binding site (motif A; amino acids 605-617), and a fraction of all selected active Taq pol I (291 of 8000) was tested for the ability to incorporate successive ribonucleotides; 23 unique mutants that added rNTPs into a growing polynucleotide chain were identified and sequenced. These mutants, each containing one to four substitutions, incorporate ribonucleotides at a efficiency approaching 10(3)-fold greater than that of wild type Taq pol I. Several mutants added successive ribonucleotides and thus can catalyze the synthesis of RNA. Sequence analysis of these mutants demonstrates that at least two amino acid residues are involved in excluding ribonucleotides from the active site. Interestingly, wild type DNA polymerases from several distinct families selectively discriminate against rUTP. This study suggests that current DNA and RNA polymerases could have evolved by divergent evolution from an ancestor that shared a common mechanism for polynucleotide synthesis.     

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.     

Prediction of protein disorder on amino acid substitutions

Intrinsically disordered regions of proteins are known to have many functional roles in cell signaling and regulatory pathways. The altered expression of these proteins due to mutations is associated with various diseases. Currently, most of the available methods focus on predicting the disordered proteins or the disordered regions in a protein. On the other hand, methods developed for predicting protein disorder on mutation showed a poor performance with a maximum accuracy of 70%. Hence, in this work, we have developed a novel method to classify the disorder-related amino acid substitutions using amino acid properties, substitution matrices, and the effect of neighboring residues that showed an accuracy of 90.0% with a sensitivity and specificity of 94.9 and 80.6%, respectively, in 10-fold cross-validation. The method was evaluated with a test set of 20% data using 10 iterations, which showed an average accuracy of 88.9%. Furthermore, we systematically analyzed the features responsible for the better performance of our method and observed that neighboring residues play an important role in defining the disorder of a given residue in a protein sequence. We have developed a prediction server to identify disorder-related mutations, and it is available at http://www.iitm.ac.in/bioinfo/DIM_Pred/.     

Selective amino acid substitutions convert the creatine transporter to a gamma-aminobutyric acid transporter

The creatine transporter (CRT) is a member of a large family of sodium-dependent neurotransmitter and amino acid transporters. The CRT is closely related to the gamma-aminobutyric acid (GABA) transporter, GAT-1, yet GABA is not an effective substrate for the CRT. The high resolution structure of a prokaryotic homologue, LeuT has revealed precise details of the substrate binding site for leucine (Yamashita, A., Singh, S. K., Kawate, T., Jin, Y., and Gouaux, E. (2005) Nature 437, 215-223). We have now designed mutations based on sequence comparisons of the CRT with GABA transporters and the LeuT structural template in an attempt to alter the substrate specificity of the CRT. Combinations of two or three amino acid substitutions at four selected positions resulted in the loss of creatine transport activity and gain of a specific GABA transport function. GABA transport by the "gain of function" mutants was sensitive to nipecotic acid, a competitive inhibitor of GABA transporters. Our results show LeuT to be a good structural model to identify amino acid residues involved in the substrate and inhibitor selectivity of eukaryotic sodium-dependent neurotransmitter and amino acid transporters. However, modification of the binding site alone appears to be insufficient for efficient substrate translocation. Additional residues must mediate the conformational changes required for the diffusion of substrate from the binding site to the cytoplasm.     

A transition probability model for amino acid substitutions from blocks.

Substitution matrices have been useful for sequence alignment and protein sequence comparisons. The BLOSUM series of matrices, which had been derived from a database of alignments of protein blocks, improved the accuracy of alignments previously obtained from the PAM-type matrices estimated from only closely related sequences. Although BLOSUM matrices are scoring matrices now widely used for protein sequence alignments, they do not describe an evolutionary model. BLOSUM matrices do not permit the estimation of the actual number of amino acid substitutions between sequences by correcting for multiple hits. The method presented here uses the Blocks database of protein alignments, along with the additivity of evolutionary distances, to approximate the amino acid substitution probabilities as a function of actual evolutionary distance. The PMB (Probability Matrix from Blocks) defines a new evolutionary model for protein evolution that can be used for evolutionary analyses of protein sequences. Our model is directly derived from, and thus compatible with, the BLOSUM matrices. The model has the additional advantage of being easily implemented.     

Enhancing antimicrobial activity of mastoparan-B by amino acid substitutions

This study evaluated antimicrobial and hemolytic activities of mastoparan-B (MP-B) isolated from the venom of the hornet, Vespa basalis, and its analogs after substituting certain amino acid (aa) residues. MP-B exhibited significantly different antimicrobial activities against bacteria species/strains tested, especially two Escherichia coli strains, Staphylococcus xylosus and Citrobacter koseri at low dosages, and was non-specific against certain Gram-positive and -negative bacteria. Our results indicated that hydrophobicity modification by single aa substitution may enhance their antimicrobial activities. An aa substituted MP-B, viz., MP-B-1, in which Trp substituted for Leu³, became more effective than the original peptide at inhibiting or killing the bacterial species tested, especially Klebsiella pneumonia, Salmonella typhimurium, and Salmonella Cholerasuis, and even up to 8 times more effective in some cases. However, MP-B-2 was virtually similar to MP-B against each bacterial species assayed, while MP-B-3 reduced its effectiveness greatly compared to others. On the other hand, MP-B and its analogs were not effective against the beneficial probiotics and they were not hemolytic to erythrocytes at the dosages tested. Our results suggested that MP-B becomes more potent against specific pathogenic bacteria and safe to the probiotics after undergoing appropriate amino acid substitutions.     

Two types of amino acid substitutions in protein evolution

Summary The frequency of amino acid substitutions, relative to the frequency expected by chance, decreases linearly with the increase in physico-chemical differences between amino acid pairs involved in a substitution. This correlation does not apply to abnormal human hemoglobins. Since abnormal hemoglobins mostly reflect the process of mutation rather than selection, the correlation manifest during protein evolution between substitution frequency and physico-chemical difference in amino acids can be attributed to natural selection. Outside of abnormal proteins, the correlation also does not apply to certain regions of proteins characterized by rapid rates of substitution. In these cases again, except for the largest physico-chemical differences between amino acid pairs, the substitution frequencies seem to be independent of the physico-chemical parameters. The limination of the substituents involving the largest physicochemical differences can once more be attributed to natural selection. For smaller physico-chemical differences, natural selection, if it is operating in the polypeptide regions, must be based on parameters other than those examined.     

Uneven distribution of amino acid substitutions throughout the amino acid sequence of homologous proteins]

A set of aligned homologous protein sequences is divided into two groups consisting of m and n most related sequences. The value of position variability for homologous protein sequences is defined as a number of failures to coincide in the intergroup comparison of all possible m*n pairs of amino acid residues in that position divided by m*n. The position variability value plotted versus the sequence position number with a window of 10 positions gives the intergroup local variability profile. Area S of the figure included between the local variability profile and the straight line corresponding to the mean local variability value is compared with the average area Sr for 1000 random homologous protein families. If S is greater than Sr by more than 2 standard deviation units sigma r, the local variability profile is assumed to contain peaks and hollows corresponding to significant variable and conservative regions of the sequences. The profile extrema containing the area surplus delta S = S-(Sr+ 2 sigma r) are cut off by two straight lines to locate significant regions. The difference (S-Sr) given in standard deviation units sigma r is believed to be the amino acid substitution overall irregularity along the homologous protein sequences OI = (S-Sr)/sigma r. The significant conservative and variable regions of six homologous sequence families (phospholipase A2, cytochromes b, alpha-subunits of Na,K-ATPase, L- and M-subunits of photosynthetic bacteria photoreaction centre and human rhodopsins) were identified. It was shown that for artificial homologous protein sequences derived by k-fold lengthening of natural protein sequences, the OI value rises as square root of k. To compare the degree of substitution irregularity in homologous protein sequence families of different lengths L the value of standard substitution overall irregularity for L = 250 is proposed.     

Accumulation pattern of amino acid substitutions in protein evolution

Summary A simple method for the evolutionary analysis of amino acid sequence data is presented and used to examine whether the number of variable sites (NVS) of a protein is constant during its evolution. The NVSs for hemoglobin and for mitochondrial cytochrome c are each found to be almost constant, and the ratio between the NVSs is close to the ratio between the unit evolutionary periods. This indicates that the substitution rate per variable site is almost uniform for these proteins, as the neutral theory claims. An advantage of the present analysis is that it can be done without knowledge of paleontological divergence times and can be extended to bacterial proteins such as bacterial c-type cytochromes. It is suggested that the NVS of cytochrome c has been almost constant even over the long period (ca. 3.0 billion years) of bacterial evolution but that at least two different substitution rates are necessary to describe the accumulated changes in the sequence. This two clock interpretation is consistent with fossil evidence for the appearance times of photosynthetic bacteria and eukaryotes.     

Evaluating antioxidative activities of amino acid substitutions on mastoparan-B

Mastoparan-B is a peptide toxin isolated from the venom of Vespa basalis, the most dangerous hornet found in Taiwan. This study is aimed to evaluate the antioxidative activities of several amino acid substitutions on MP-B, and examined the influences of mast cell degranulation and hemolytic activities in parallel with antioxidative activities. The correlations between the biological function and amino acid sequence were assessed. Our study shows original MP-B is a valuable antioxidant at low concentration in competing with nitric-oxide for oxygen molecules and possesses good antioxidative enzyme activities resembled to superoxidase dismutase and glutathione peroxidase. And there are no predominant rates of mast cell degranulation and hemolytic effects in such condition. With proper substitutions, the reducing power, DPPH scavenging activity and glutathione reductase-like enzyme activity of MP-B can increase clearly. The results demonstrate that MP-B analogs are very potential to be applicable antioxidants for other antioxidative usages.     

Monoclonal antibodies recognizing single amino acid substitutions in hemoglobin

Four monoclonal antibodies (mAb) to non-human primate hemoglobin referred to as Cap-4, Cap-5, Rh-2, and Rh-4, and two mAb to human hemoglobin, referred to as H-1 and H-3 were isolated and were partially characterized. Binding studies with these mAb on a panel of hemoglobins and isolated alpha and beta globin chains revealed a unique reactivity pattern for each mAb. Amino acid sequence analysis of the antigens used to generate the binding data suggests that the specific recognition of certain hemoglobin antigens by each mAb is controlled by the presence of a particular amino acid at a specific position within the epitope. The use of synthetic peptides as antigens confirmed this observation for five of the mAb. No synthetic peptides were tested with the sixth mAb, Rh-2. The amino acids required for binding of mAb Cap-4, Cap-5, Rh-4, and Rh-2 to hemoglobin are alanine at beta 5, threonine at beta 13, glutamine at beta 125, and leucine at alpha 68. The non-human primate hemoglobin antibodies require a specific amino acid that is not present in human hemoglobin. The amino acid required for binding of Cap-4, Cap-5, and Rh-4 could arise by a single base change in the beta globin gene, whereas the amino acid required for Rh-2 binding would only occur if two base changes occurred. Thus these mAb are candidate probes for a somatic cell mutation assay on the basis of the detection of peripheral blood red cells that possess single amino acid substituted hemoglobin as a result of single base substitutions in the globin genes of precursor cells.     

Reciprocal amino acid substitutions in the evolution of homologous peptides

Amino acid sequences of peptides are often inferred from their amino acid compositions by comparison with homologous peptides of known sequence. The probabilities are considered that by such an approach errors are made due to the occurrence of balanced double changes, i.e. reciprocal substitutions, between two homologous peptides of identical compositions. Formulae are derived for the calculation of these probabilities, depending on peptide length and evolutionary distance. However, such calculations requiring too much computer time, the probabilities for reciprocal substitutions are estimated by simulation of evolutionary changes in peptides. It can be concluded from the resulting data that for many purposes the possible errors in amino acid sequences partially inferred from amino acid compositions are acceptably small.     

Nature of amino acid substitutions in homologous proteins during evolution

Replacement substitutions of mitochondrial cytochrome $\text{c}$ and α- and β-chains of haemoglobin have been studied by considering the structural similarity among amino acid residues at the secondary and tertiary structural levels. Secondary structural similarity explains ～70% while tertiary structural similarity explains ～50% of observed replacements for most of the cases. These structural similarities could not account for all the replacement substitutions. The study was extended to consider the composition of codons, and the chemical nature and polarity of the replacing and replaced residues. These also could not individually account for all the affected replacements. In general, no property of amino acid residues is conserved for substitutions occurring at any single position during evolution of proteins.     

Five PDGF B amino acid substitutions convert PDGF A to a PDGF B-like transforming molecule.

We used site-directed mutagenesis to determine the minimum number of PDGF B residues needed to convert PDGF A to a potently transforming PDGF B-like molecule. Substitution of two PDGF B subdomains, 106-115 and 135-144, were found to be critical. These substitutions were sufficient to broaden the ability of PDGF A to activate beta as well as alpha platelet-derived growth factor (PDGF) receptors and increase its transforming efficiency to that of PDGF B. Within subdomain I, either PDGF B residues Arg-109 and Asn-115 or Arg-109, Leu-110, and Arg-113, in combination with subdomain II PDGF B residues Asn-136, Arg-137, and Arg-142 were identified as being essential. Those mutants with transforming ability comparable with PDGF B showed significantly lower efficiencies of beta receptor triggering. Thus, our studies identify a small number of PDGF B amino acids indispensable for beta PDGF receptor interaction and suggest that a low level of beta PDGF receptor activation is sufficient to dramatically increase PDGF transforming efficiency in NIH 3T3 cells.