Challenges of the genetic code for exploring sequence space in directed protein evolution

Directed protein evolution is the most versatile method for studying protein structure–function relationships, and for tailoring a protein's properties to the needs of industrial applications. In this review, we performed a statistical analysis on the genetic code to study the extent and consequence of the organization of the genetic code on amino acid substitution patterns generated in directed evolution experiments. In detail, we analyzed amino acid substitution patterns caused by (a) a single nucleotide (nt) exchange at each position of all 64 codons, and (b) two subsequent nt exchanges (first and second nt, first and third nt, second and third nt). Additionally, transitions and transversions mutations were compared at the level of amino acid substitution patterns. The latter analysis showed that single nucleotide substitution in a codon generates only 39.5% of the natural diversity on the protein level with 5.2–7 amino acid substitutions per codon. Transversions generate more complex amino acid substitution patterns (increased number and chemically more diverse amino acid substitutions) than transitions. Simultaneous nt exchanges at both first and second nt of a codon generates very diverse amino acid substitution patterns, achieving 83.2% of the natural diversity. The statistical analysis described in this review sets the objectives for novel random mutagenesis methods that address the consequences of the organization of the genetic code. Random mutagenesis methods that favor transversions or introduce consecutive nt exchanges can contribute in this regard.     

A statistical analysis of random mutagenesis methods used for directed protein evolution

We have developed a statistical method named MAP (mutagenesis assistant program) to equip protein engineers with a tool to develop promising directed evolution strategies by comparing 19 mutagenesis methods. Instead of conventional transition/transversion bias indicators as benchmarks for comparison, we propose to use three indicators based on the subset of amino acid substitutions generated on the protein level: (1) protein structure indicator; (2) amino acid diversity indicator with a codon diversity coefficient; and (3) chemical diversity indicator. A MAP analysis for a single nucleotide substitution was performed for four genes: (1) heme domain of cytochrome P450 BM-3 from Bacillus megaterium (EC 1.14.14.1); (2) glucose oxidase from Aspergillus niger (EC 1.1.3.4); (3) arylesterase from Pseudomonas fluorescens (EC 3.1.1.2); and (4) alcohol dehydrogenase from Saccharomyces cerevisiae (EC 1.1.1.1). Based on the MAP analysis of these four genes, 19 mutagenesis methods have been evaluated and criteria for an ideal mutagenesis method have been proposed. The statistical analysis showed that existing gene mutagenesis methods are limited and highly biased. An average amino acid substitution per residue of only 3.15-7.4 can be achieved with current random mutagenesis methods. For the four investigated gene sequences, an average fraction of amino acid substitutions of 0.5-7% results in stop codons and 4.5-23.9% in glycine or proline residues. An average fraction of 16.2-44.2% of the amino acid substitutions are preserved, and 45.6% (epPCR method) are chemically different. The diversity remains low even when applying a non-biased method: an average of seven amino acid substitutions per residue, 2.9-4.7% stop codons, 11.1-16% glycine/proline residues, 21-25.8% preserved amino acids, and 55.5% are amino acids with chemically different side-chains. Statistical information for each mutagenesis method can further be used to investigate the mutational spectra in protein regions regarded as important for the property of interest.     

Challenges of the genetic code for exploring sequence space in directed protein evolution

Directed protein evolution is the most versatile method for studying protein structure-function relationships, and for tailoring a protein's properties to the needs of industrial applications. In this review, we performed a statistical analysis on the genetic code to study the extent and consequence of the organization of the genetic code on amino acid substitution patterns generated in directed evolution experiments. In detail, we analyzed amino acid substitution patterns caused by (a) a single nucleotide (nt) exchange at each position of all 64 codons, and (b) two subsequent nt exchanges (first and second nt, first and third nt, second and third nt). Additionally, transitions and transversions mutations were compared at the level of amino acid substitution patterns. The latter analysis showed that single nucleotide substitution in a codon generates only 39.5% of the natural diversity on the protein level with 5.2-7 amino acid substitutions per codon. Transversions generate more complex amino acid substitution patterns (increased number and chemically more diverse amino acid substitutions) than transitions. Simultaneous nt exchanges at both first and second nt of a codon generates very diverse amino acid substitution patterns, achieving 83.2% of the natural diversity. The statistical analysis described in this review sets the objectives for novel random mutagenesis methods that address the consequences of the organization of the genetic code. Random mutagenesis methods that favor transversions or introduce consecutive nt exchanges can contribute in this regard.     

Probing beta-lactamase structure and function using random replacement mutagenesis.

A new analytical mutagenesis technique is described that involves randomizing the DNA sequence of a short stretch of a gene (3-6 codons) and determining the percentage of all possible random sequences that produce a functional protein. A low percentage of functional random sequences in a complete library of random substitutions indicates that the region mutagenized is important for the structure and/or function of the protein. Repeating the mutagenesis over many regions throughout a protein gives a global perspective of which amino acid sequences in a protein are critical. We applied this method to 66 codons of the gene encoding TEM-1 beta-lactamase in 19 separate experiments. We found that TEM-1 beta-lactamase is extremely tolerant of amino acid substitutions: on average, 44% of all mutants with random substitutions function and 20% of the substitutions are expressed, secreted, and fold well enough to function at levels similar to those for the wild-type enzyme. We also found a few exceptional regions where only a few random sequences function. Examination of the X-ray structures of homologous beta-lactamases indicates that the regions most sensitive to substitution are in the vicinity of the active site pocket or buried in the hydrophobic core of the protein. DNA sequence analysis of functional random sequences has been used to obtain more detailed information about the amino acid sequence requirements for several regions and this information has been compared to sequence conservation among several related beta-lactamases.     

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.     

Mutational analysis of potential pore-lining amino acids in TM IV of the Na(+)/H(+) exchanger

The Na(+)/H(+) exchanger isoform 1 (NHE1) is an integral membrane protein that regulates intracellular pH by extruding an intracellular H(+) in exchange for one extracellular Na(+). In this study we examined the effect of site-specific mutagenesis on the pore-lining amino acid Phe161 and effects of mutagenesis on the charged amino acids Asp159 and Asp172. There was no absolute requirement for a carboxyl side chain at amino acid Asp159 or Asp172. Mutation of Asp159 to Asn or Gln maintained or increased the activity of the protein. Similarly, for Asp172, substitution with a Gln residue maintained activity of the protein, even though substitution with an Asn residue was inhibitory. The Asp172Glu mutant possessed normal activity after correction for its aberrant expression and surface targeting. Replacement of Phe161 with a Leu demonstrated that it was not irreplaceable in NHE1 function. However, the mutation Phe161lys inhibited NHE1 function, while the Phe161Ala mutation caused altered NHE1 targeting and expression levels. Our results show that these three amino acids, while being important in NHE1 function, are not irreplaceable. This study demonstrates that multiple substitutions at a single amino acid residue may be necessary to get a clearer picture membrane protein function.     

Comparison of epitope structures of H3HAs through protein modeling of influenza A virus hemagglutinin: mechanism for selection of antigenic variants in the presence of a monoclonal antibody

Starting with nine plaques of influenza A/Kamata/14/91(H3N2) virus, we selected mutants in the presence of monoclonal antibody 203 (mAb203). In total, amino acid substitutions were found at nine positions (77, 80, 131, 135, 141, 142, 143, 144 and 146), which localized in the antigenic site A of the hemagglutinin (HA). The escape mutants differed in the extent to which they had lost binding to mAb203. HA protein with substitutions of some amino acid residues created by site-directed mutagenesis in the escape mutants retained the ability to bind to mAb203. Changes in the amino acid character affecting charge or hydrophobicity accounted for the binding capacity to the antibody of the HA with most of the substitutions in the escape mutants and binding-positive mutants. However, the effect of some amino acid substitutions remained unexplained. A three-dimensional model of the 1991 HA was constructed and used to analyze substituted amino acids in these mutants for the accessible surface hydrophobic and hydrophilic characters. One amino acid substitution in an escape mutant and another amino acid substitution in a binding-positive mutant seemed to be explained by the changes noted on this model.     

Mutational analysis of potential pore-lining amino acids in TM IV of the Na/H exchanger

The Na⁺/H⁺ exchanger isoform 1 (NHE1) is an integral membrane protein that regulates intracellular pH by extruding an intracellular H⁺ in exchange for one extracellular Na⁺. In this study we examined the effect of site-specific mutagenesis on the pore-lining amino acid Phe161 and effects of mutagenesis on the charged amino acids Asp159 and Asp172. There was no absolute requirement for a carboxyl side chain at amino acid Asp159 or Asp172. Mutation of Asp159 to Asn or Gln maintained or increased the activity of the protein. Similarly, for Asp172, substitution with a Gln residue maintained activity of the protein, even though substitution with an Asn residue was inhibitory. The Asp172Glu mutant possessed normal activity after correction for its aberrant expression and surface targeting. Replacement of Phe161 with a Leu demonstrated that it was not irreplaceable in NHE1 function. However, the mutation Phe161lys inhibited NHE1 function, while the Phe161Ala mutation caused altered NHE1 targeting and expression levels. Our results show that these three amino acids, while being important in NHE1 function, are not irreplaceable. This study demonstrates that multiple substitutions at a single amino acid residue may be necessary to get a clearer picture membrane protein function.     

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.     

Identification of amino acid residues critical for infection with ecotropic murine leukemia retrovirus.

The murine cationic amino acid transporter is also the receptor for murine ecotropic leukemia retrovirus (MuLV-E). Recently, we have cloned a human gene (H13) homologous to the murine ecotropic retroviral receptor (ERR). Although the human homolog is very similar to murine ERR in sequence (87.6% amino acid identity) and structure (14 transmembrane-spanning domains), the human protein fails to function as a receptor for MuLV-E. To identify amino acid residues critical for MuLV-E infection, we took advantage of this species difference and substituted human H13 and murine ERR amino acid residues. Mouse-human chimeric receptor molecules were generated by taking advantage of using common restriction sites. These studies demonstrated that extracellular domains 3 and/or 4 contain the critical amino acid residues. Oligonucleotide-directed mutagenesis was then used to create 13 individual ERR mutants containing one or two amino acids substitutions or insertions within these two extracellular domains. Substitution of as few as one amino acid residue (Tyr) at position 235 in ERR with the corresponding H13 amino acid residue Pro abrogates the ability to function as a receptor for MuLV-E infection. Conversely, substitution of just two amino acid residues at positions 240 and 242 or 242 and 244 in H13 with the corresponding amino acid residues in ERR endows H13 with the ability to function as the receptor. This observation can be utilized to significantly improve the safety of retrovirus-mediated gene therapy in humans.     

ENVIRON: a software package to compare protein three-dimensional structures with homologous sequences using local structural motifs

This work presents a method to compare local clusters of interactingresidues as observed in a known three-dimensional protein structurewith corresponding clusters inferred from homologous proteinsequences, assuming conserved protein folding. For this purposethe local environment of a selected residue in a known proteinstructure is defined as the ensemble of amino acids in contactwith it in the folded state. Using a multiple sequence alignmentto identify corresponding residues in homologous proteins, adetailed comparison can be performed between the local environmentof a selected amino acid in the template protein structure andthe expected local environments at the sets of equivalent residues,derived from the aligned protein sequences. The comparison makesit possible to detect conserved local features such as hydrogenbonding or complementarity in residue substitution. A globalmeasure of environmental similarity is also defined, to searchfor conserved amino acid clusters subject to functional or structural constraints. The proposed approach is useful for investigatingprotein function as well as for site-directed mutagenesis experiments,where appropriate amino acid substitutions can be suggestedby observing naturally occurring protein variants.     

44-amino-acid E5 transforming protein of bovine papillomavirus requires a hydrophobic core and specific carboxyl-terminal amino acids.

The 44-amino-acid E5 protein of bovine papillomavirus type 1 is the shortest known protein with transforming activity. To identify the specific amino acids required for in vitro focus formation in mouse C127 cells, we used oligonucleotide-directed saturation mutagenesis to construct an extensive collection of mutants with missense mutations in the E5 gene. Characterization of mutants with amino acid substitutions in the hydrophobic middle third of the E5 protein indicated that efficient transformation requires a stretch of hydrophobic amino acids but not a specific amino acid sequence in this portion of the protein. Many amino acids in the carboxyl-terminal third of the protein can also undergo substitution without impairment of focus-forming activity, but the amino acids at seven positions, including two cysteine residues that mediate dimer formation, appear essential for efficient transforming activity. These essential amino acids are the most well conserved among related fibropapillomaviruses. The small size of the E5 protein, its lack of similarity to other transforming proteins, and its ability to tolerate many amino acid substitutions implies that it transforms cells via a novel mechanism.     

Many amino acid substitution variants identified in DNA repair genes during human population screenings are predicted to impact protein function

Over 520 different amino acid substitution variants have been previously identified in the systematic screening of 91 human DNA repair genes for sequence variation. Two algorithms were employed to predict the impact of these amino acid substitutions on protein activity. Sorting Intolerant from Tolerant (SIFT) classified 226 of 508 variants (44%) as "Intolerant." Polymorphism Phenotyping (PolyPhen) classed 165 of 489 amino acid substitutions (34%) as "Probably or possibly damaging." Another 9-15% of the variants were classed as "Potentially intolerant or damaging." The results from the two algorithms are highly associated, with concordance in predicted impact observed for approximately 62% of the variants. Twenty-one to thirty-one percent of the variant proteins are predicted to exhibit reduced activity by both algorithms. These variants occur at slightly lower individual allele frequency than do the variants classified as "Tolerant" or "Benign." Both algorithms correctly predicted the impact of 26 functionally characterized amino acid substitutions in the APE1 protein on biochemical activity, with one exception. It is concluded that a substantial fraction of the missense variants observed in the general human population are functionally relevant. These variants are expected to be the molecular genetic and biochemical basis for the associations of reduced DNA repair capacity phenotypes with elevated cancer risk.     

SIFT Indel: Predictions for the Functional Effects of Amino Acid Insertions/Deletions in Proteins

Indels in the coding regions of a gene can either cause frameshifts or amino acid insertions/deletions. Frameshifting indels are indels that have a length that is not divisible by 3 and subsequently cause frameshifts. Indels that have a length divisible by 3 cause amino acid insertions/deletions or block substitutions; we call these 3n indels. The new amino acid changes resulting from 3n indels could potentially affect protein function. Therefore, we construct a SIFT Indel prediction algorithm for 3n indels which achieves 82% accuracy, 81% sensitivity, 82% specificity, 82% precision, 0.63 MCC, and 0.87 AUC by 10-fold cross-validation. We have previously published a prediction algorithm for frameshifting indels. The rules for the prediction of 3n indels are different from the rules for the prediction of frameshifting indels and reflect the biological differences of these two different types of variations. SIFT Indel was applied to human 3n indels from the 1000 Genomes Project and the Exome Sequencing Project. We found that common variants are less likely to be deleterious than rare variants. The SIFT indel prediction algorithm for 3n indels is available at http://sift-dna.org/     

Expanded molecular diversity generation during directed evolution by trinucleotide exchange (TriNEx)

Trinucleotide exchange (TriNEx) is a method for generating novel molecular diversity during directed evolution by random substitution of one contiguous trinucleotide sequence for another. Single trinucleotide sequences were deleted at random positions in a target gene using the engineered transposon MuDel that were subsequently replaced with a randomized trinucleotide sequence donated by the DNA cassette termed SubSeq(NNN). The bla gene encoding TEM-1 beta-lactamase was used as a model to demonstrate the effectiveness of TriNEx. Sequence analysis revealed that the mutations were distributed throughout bla, with variants containing single, double and triple nucleotide changes. Many of the resulting amino acid substitutions had significant effects on the in vivo activity of TEM-1, including up to a 64-fold increased activity toward ceftazidime and up to an 8-fold increased resistance to the inhibitor clavulanate. Many of the observed amino acid substitutions were only accessible by exchanging at least two nucleotides per codon, including charge-switch (R164D) and aromatic substitution (W165Y) mutations. TriNEx can therefore generate a diverse range of protein variants with altered properties by combining the power of site-directed saturation mutagenesis with the capacity of whole-gene mutagenesis to randomly introduce mutations throughout a gene.     

Heterotachy and functional shift in protein evolution

Study of structure/function relationships constitutes an important field of research, especially for modification of protein function and drug design. However, the fact that rational design (i.e. the modification of amino acid sequences by means of directed mutagenesis, based on knowledge of the three-dimensional structure) appears to be much less efficient than irrational design (i.e. random mutagenesis followed by in vitro selection) clearly indicates that we understand little about the relationships between primary sequence, three-dimensional structure and function. The use of evolutionary approaches and concepts will bring insights to this difficult question. The increasing availability of multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico evolutionary methods to predict details of protein function in duplicated (paralogous) proteins. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogs. It has been proposed that the positions that show switches in substitution rate over time--i.e., 'heterotachous sites'--are good indicators of functional divergence. However, it appears that heterotachy is a much more general process, since most variable sites of homologous proteins with no evidence of functional shift are heterotachous. Similarly, it appears that switches in substitution rate are as frequent when paralogous sequences are compared as when orthologous sequences are compared. Heterotachy, instead of being indicative of functional shift, may more generally reflect a less specific process related to the many intra- and inter-molecular interactions compatible with a range of more or less equally viable protein conformations. These interactions will lead to different constraints on the nature of the primary sequences, consistently with theories suggesting the non-independence of substitutions in proteins. However, a specific type of amino acid variation might constitute a good indicator of functional divergence: substitutions occurring at positions that are generally slowly evolving. Such substitutions at constrained sites are indeed much more frequent soon after gene duplication. The identification and analysis of these sites by complementing structural information with evolutionary data may represent a promising direction to future studies dealing with the functional characterization of an ever increasing number of multi-gene families identified by complete genome analysis.     

Effect of point mutations in the native simian virus 40 tumor antigen, and in synthetic peptides corresponding to the H-2Db-restricted epitopes, on antigen presentation and recognition by cytotoxic T lymphocyte clones.

The effects on CTL recognition of individual amino acid substitutions within epitopes I, II, and III of SV40 tumor Ag (T Ag) were examined. Epitope I spans amino acids 207 to 215, and epitope II/III is within residues 223 to 231 of SV40 T Ag. An amino acid substitution at position 207 (Ala----Val) or 214 (Lys----Glu) of SV40 T Ag expressed in transformed cells resulted in loss of epitope I, recognized by CTL clone Y-1. The amino acid substitution at residue 214 in the corresponding synthetic peptide, LT207-215(214-Lys----Glu), also led to loss of recognition by CTL clone Y-1. The recognition, by CTL clone Y-1, of peptides LT207-215 and LT207-217 with an Ala----Val substitution at position 207 was severely affected. Peptides LT205-215 and LT205-219 with the Ala----Val substitution at residue 207 were, however, recognized by CTL clone Y-1, suggesting that residues 205 and 206 may be involved in presentation of site I. Alteration of residue 224 (Lys----Glu) in the native T Ag resulted in loss of recognition by both CTL clones Y-2 and Y-3. However, a peptide corresponding to epitope II/III with an identical amino acid substitution at residue 224 provided a target for CTL clone Y-3 but not clone Y-2. A change of Lys----Gln at residue 224 in both the native protein and a synthetic peptide caused loss of recognition by CTL clone Y-2 but not CTL clone Y-3. Further, an amino acid substitution of Lys----Arg at position 224 of the native T Ag decreased recognition of epitope II/III by CTL clones Y-2 and Y-3 but had no effect on recognition of a synthetic peptide bearing the same substitution. These results indicate that the mutagenesis approach, resulting in identical amino acid substitutions in the native protein and in the synthetic peptides, may provide insight into the role of individual residues in the processing, presentation, and recognition of CTL recognition epitopes.     

Effects of single-point amino acid substitutions on the structure and function neuraminidase proteins in influenza A virus

In order to clarify the effect of amino acid substitutions on the structure and function of the neuraminidase (NA) protein of influenza A virus, we introduced single-point amino acid substitutions into the NA protein of the A/Tokyo/3/67 (H2N2) strain using PCR-based random mutation. The rate of tolerant random one amino acid substitutions in the NA protein was 47%. Rates of tolerant substitutions for the stalk and for the surface and inner portion of the head region of the NA protein were 79, 54, and 19%, respectively. Deleterious changes, such as those causing the NA protein to stop at the Golgi/endoplasmic reticulum, were scattered throughout the protein. On the other hand, the ratio of mutations with which the NA protein lost neuraminidase activity, but was transported to the cell surface, decreased in proportion to the distance from the structural center of enzyme active site. In order to investigate the effect of accumulated amino acid substitutions on the structural character of the N2NA protein during evolution, the same amino acid substitutions were introduced by site-directed mutagenesis at 23 homologous positions on N2 proteins of A/Tokyo/3/67, A/Bangkok/15/85 (H3N2), and A/Mie/1/2004 (H3N2). The results showed a shift, or discordance, in tolerance at some of the positions. An increase in discordance was correlated with the interval in years between virus strains, and the discordance rate was estimated to be 0.6-0.7% per year.     

Eigenvalue analysis of amino acid substitution matrices reveals a sharp transition of the mode of sequence conservation in proteins

The pattern of amino acid substitutions and sequence conservation over many structure-based alignments of protein sequences was analyzed as a function of percentage sequence identity. The statistics of the amino acid substitutions were converted into the form of log-odds amino acid substitution matrices to which eigenvalue decomposition was applied. It was found that the most important component of the substitution matrices exhibited a sharp transition at the sequence identity of 30-35%, which coincides with the twilight zone. Above the transition point, the most dominant component is related to the mutability of amino acids and it acts to disfavor any substitutions, whereas below the transition point, the most dominant component is related to the hydrophobicity of amino acids and substitutions between residues of similar hydrophobic character are positively favored. Implications for protein evolution and sequence analysis are discussed.     

Thermostabilization of Escherichia coli ribonuclease HI by replacing left-handed helical Lys95 with Gly or Asn.

From the systematic replacements of amino acid residues of Escherichia coli ribonuclease HI with those of its thermophilic counterpart, the basic protrusion domain including region 6 (R6) from residues 91 to 95 was found to increase the structural stability of the mutant protein (Kimura, S., Nakamura, H., Hashimoto, T., Oobatake, M., and Kanaya, S. (1992) J. Biol. Chem. 267, 21535-21542). Further mutagenesis concentrating in the R6 region has revealed that replacements of Lys95 at the left-handed structure with Gly or Asn essentially enhances the protein stability. Gly and Asn substitutions stabilize the protein up to 1.9 kcal/mol and 0.9 kcal/mol in the free energy changes of unfolding, respectively. We propose that the amino acid substitution of left-handed non-Gly residue with Gly or Asn residue can be used as one of the general strategies to enhance protein stability, when such a non-Gly residue itself does not seriously contribute to protein stability.