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.

[1]  F. Arnold,et al.  Patterns of adaptation in a laboratory evolved thermophilic enzyme. , 2001, Biochimica et biophysica acta.

[2]  B. Connolly,et al.  Low-fidelity Pyrococcus furiosus DNA polymerase mutants useful in error-prone PCR. , 2004, Nucleic acids research.

[3]  K. R. Marshall,et al.  P450 BM3: the very model of a modern flavocytochrome. , 2002, Trends in biochemical sciences.

[4]  R D Schmid,et al.  Directed evolution of the fatty-acid hydroxylase P450 BM-3 into an indole-hydroxylating catalyst. , 2000, Chemistry.

[5]  T. S. Wong,et al.  The diversity challenge in directed protein evolution. , 2006, Combinatorial chemistry & high throughput screening.

[6]  S. Freeland,et al.  The Case for an Error Minimizing Standard Genetic Code , 2003, Origins of life and evolution of the biosphere.

[7]  Massimo Di Giulio,et al.  The origin of the genetic code: theories and their relationships, a review. , 2005 .

[8]  R. Myers,et al.  A general method for saturation mutagenesis of cloned DNA fragments. , 1985, Science.

[9]  A new method for random mutagenesis of complete genes: enzymatic generation of mutant libraries in vitro. , 1988, Protein engineering.

[10]  K. Choi,et al.  Cloning and nucleotide sequence of an esterase gene from Pseudomonas fluorescens and expression of the gene in Escherichia coli. , 1990, Agricultural and Biological Chemistry.

[11]  J. Punnonen,et al.  Development of novel vaccines using DNA shuffling and screening strategies. , 2004, Current opinion in molecular therapeutics.

[12]  J. Bennetzen,et al.  The primary structure of the Saccharomyces cerevisiae gene for alcohol dehydrogenase. , 1982, The Journal of biological chemistry.

[13]  F. Arnold,et al.  Designed evolution of enzymatic properties. , 2000, Current opinion in biotechnology.

[14]  C R Woese,et al.  Order in the genetic code. , 1965, Proceedings of the National Academy of Sciences of the United States of America.

[15]  George Georgiou,et al.  Preparative expression of secreted proteins in bacteria: status report and future prospects. , 2005, Current opinion in biotechnology.

[16]  P. Stougaard,et al.  A new glucose oxidase from Aspergillus niger: characterization and regulation studies of enzyme and gene , 1996, Applied Microbiology and Biotechnology.

[17]  Hiroshi Murakami,et al.  Random insertion and deletion of arbitrary number of bases for codon-based random mutation of DNAs , 2002, Nature Biotechnology.

[18]  Jean-Louis Reymond,et al.  Recent advances in enzyme assays. , 2004, Trends in biotechnology.

[19]  R D Schmid,et al.  Rational evolution of a medium chain-specific cytochrome P-450 BM-3 variant. , 2001, Biochimica et biophysica acta.

[20]  F. Arnold,et al.  How enzymes adapt: lessons from directed evolution , 2001, Trends in biochemical sciences.

[21]  C. Schmidt-Dannert,et al.  A continuous spectrophotometric assay for P450 BM-3, a fatty acid hydroxylating enzyme, and its mutant F87A. , 1999, Analytical biochemistry.

[22]  M. Zacharias,et al.  Sensitive Assay for Laboratory Evolution of Hydroxylases toward Aromatic and Heterocyclic Compounds , 2005, Journal of biomolecular screening.

[23]  L. Loeb,et al.  A Single Highly Mutable Catalytic Site Amino Acid Is Critical for DNA Polymerase Fidelity* , 2001, The Journal of Biological Chemistry.

[24]  Jingding Huang,et al.  A new approach to random mutagenesis in vitro , 2004, Biotechnology and bioengineering.

[25]  Jürgen Pleiss,et al.  A Molecular Mechanism of Enantiorecognition of Tertiary Alcohols by Carboxylesterases , 2003, Chembiochem : a European journal of chemical biology.

[26]  Manel Camps,et al.  Targeted gene evolution in Escherichia coli using a highly error-prone DNA polymerase I , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[27]  J. Reymond Spectrophotometric Enzyme Assays for High-Throughput Screening , 2004 .

[28]  Guy Sella,et al.  On the Evolution of Redundancy in Genetic Codes , 2001, Journal of Molecular Evolution.

[29]  R. Siegel,et al.  Generation of large libraries of random mutants in Bacillus subtilis by PCR-based plasmid multimerization. , 1997, BioTechniques.

[30]  D. M. Brown,et al.  An approach to random mutagenesis of DNA using mixtures of triphosphate derivatives of nucleoside analogues. , 1996, Journal of molecular biology.

[31]  Sergey V. Balashov,et al.  Specificity of spontaneous mutations induced in mutA mutator cells. , 2004, Mutation research.

[32]  U. Schwaneberg,et al.  A P450 BM-3 mutant hydroxylates alkanes, cycloalkanes, arenes and heteroarenes. , 2001, Journal of biotechnology.

[33]  Jian-Kang Chen,et al.  Transfection of an Active Cytochrome P450 Arachidonic Acid Epoxygenase Indicates That 14,15-Epoxyeicosatrienoic Acid Functions as an Intracellular Second Messenger in Response to Epidermal Growth Factor* , 1999, The Journal of Biological Chemistry.

[34]  Jean-Louis Reymond,et al.  Enzyme assays for high-throughput screening. , 2004, Current opinion in biotechnology.

[35]  A. Fulco,et al.  Coding nucleotide, 5' regulatory, and deduced amino acid sequences of P-450BM-3, a single peptide cytochrome P-450:NADPH-P-450 reductase from Bacillus megaterium. , 1989, The Journal of biological chemistry.

[36]  Exon Inser,et al.  PCR-Based Random Mutagenesis Using Manganese and Reduced dNTP Concentration , 1997 .

[37]  J. Falck,et al.  An Active Site Substitution, F87V, Converts Cytochrome P450 BM-3 into a Regio- and Stereoselective (14S,15R)-Arachidonic Acid Epoxygenase* , 1997, The Journal of Biological Chemistry.

[38]  S. Wain-Hobson,et al.  Hypermutagenic PCR involving all four transitions and a sizeable proportion of transversions. , 1996, Nucleic acids research.

[39]  F. Arnold,et al.  Directed evolution of enzyme catalysts. , 1997, Trends in biotechnology.

[40]  Ryota Fujii,et al.  One-step random mutagenesis by error-prone rolling circle amplification. , 2004, Nucleic acids research.