Evolutionary molecular engineering by random elongation mutagenesis

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

A novel framework for engineering protein loops exploring length and compositional variation

Insertions and deletions (indels) are known to affect function, biophysical properties and substrate specificity of enzymes, and they play a central role in evolution. Despite such clear significance, this class of mutation remains an underexploited tool in protein engineering with few available platforms capable of systematically generating and analysing libraries of varying sequence composition and length. We present a novel DNA assembly platform (InDel assembly), based on cycles of endonuclease restriction digestion and ligation of standardised dsDNA building blocks, that can generate libraries exploring both composition and sequence length variation. In addition, we developed a framework to analyse the output of selection from InDel-generated libraries, combining next generation sequencing and alignment-free strategies for sequence analysis. We demonstrate the approach by engineering the well-characterized TEM-1 β-lactamase Ω-loop, involved in substrate specificity, identifying multiple novel extended spectrum β-lactamases with loops of modified length and composition-areas of the sequence space not previously explored. Together, the InDel assembly and analysis platforms provide an efficient route to engineer protein loops or linkers where sequence length and composition are both essential functional parameters.

Enhancing Catalytic Efficiency of an Actinoplanes utahensis Echinocandin B Deacylase through Random Mutagenesis and Site-Directed Mutagenesis

Echinocandin B deacylase (EBDA), from Actinoplanes utahensis ZJB-08196, is capable of cleaving the linoleoyl group from echinocandin B (ECB), forming the echinocandin B nucleus (ECBN), which is a key precursor of semisynthetic antifungal antibiotics. In the present study, molecular evolution of AuEBDA by random mutagenesis combined with site-directed mutagenesis (SDM) and screening was performed. Random mutagenesis on the wild-type (WT) AuEBDA generated two beneficial substitutions of G287Q, R527V. The "best" variant AuEBDA-G287Q/R527V was obtained by combining G287Q with R527V through SDM, which was most active at 35 °C, pH 7.5, with K_m and v_max values of 0.68 mM and 395.26 U/mg, respectively. Mutation of G287Q/R527V markedly increased the catalytic efficiency k_cat/K_m by 290% compared with the WT-AuEBDA.

Phylogenetic divergence of cell biological features

Most cellular features have a range of states, but understanding the mechanisms responsible for interspecific divergence is a challenge for evolutionary cell biology. Models are developed for the distribution of mean phenotypes likely to evolve under the joint forces of mutation and genetic drift in the face of constant selection pressures. Mean phenotypes will deviate from optimal states to a degree depending on the effective population size, potentially leading to substantial divergence in the absence of diversifying selection. The steady-state distribution for the mean can even be bimodal, with one domain being largely driven by selection and the other by mutation pressure, leading to the illusion of phenotypic shifts being induced by movement among alternative adaptive domains. These results raise questions as to whether lineage-specific selective pressures are necessary to account for interspecific divergence, providing a possible platform for the establishment of null models for the evolution of cell-biological traits.

Introduction to the Field of Enzyme Immobilization and Stabilization

Enzyme stabilization is important for many biomedical or industrial application of enzymes (i.e., cell-free biotransformations and biosensors). In many applications, the goal is to provide extended active lifetime at normal environmental conditions with traditional substrates at low concentrations in buffered solutions. However, as enzymes are used for more and more applications, there is a desire to use them in extreme environmental conditions (i.e., high temperatures), in high substrate concentration or high ionic strength, and in nontraditional solvent systems. This chapter introduces the topic enzyme stabilization and the methods used for enzyme stabilization including enzyme immobilization.

Coiled-coil length: Size does matter

Protein evolution is governed by processes that alter primary sequence but also the length of proteins. Protein length may change in different ways, but insertions, deletions and duplications are the most common. An optimal protein size is a trade-off between sequence extension, which may change protein stability or lead to acquisition of a new function, and shrinkage that decreases metabolic cost of protein synthesis. Despite the general tendency for length conservation across orthologous proteins, the propensity to accept insertions and deletions is heterogeneous along the sequence. For example, protein regions rich in repetitive peptide motifs are well known to extensively vary their length across species. Here, we analyze length conservation of coiled-coils, domains formed by an ubiquitous, repetitive peptide motif present in all domains of life, that frequently plays a structural role in the cell. We observed that, despite the repetitive nature, the length of coiled-coil domains is generally highly conserved throughout the tree of life, even when the remaining parts of the protein change, including globular domains. Length conservation is independent of primary amino acid sequence variation, and represents a conservation of domain physical size. This suggests that the conservation of domain size is due to functional constraints.

Creation of thermostable and alkaline stable xylanase variants by DNA shuffling

Mutant xylanases, G41 and G53, were generated by random mutagenesis of Thermomyces lanuginosus xylanase DSM 5826 (xynA) in a previous study. Incubation at 90 min showed that G41 had 75% activity at 80 °C and G53 had 93% activity at pH 10. In order to create xylanase variants possessing both thermal and alkaline stability in a single enzyme, G41 and G53 served as templates for DNA shuffling using the StEP recombination method. One of the resulting StEP recombinants, S340, retained 54% stability at 80 °C and 60% stability at pH 10 with three resulting amino acid mutations. Another StEP recombinant, S325, displayed 85% stability at 80 °C and 60% stability at pH 10 and DNA sequencing showed that it inherited mutations from both parents. All thermostable variants displayed an increase in arginine content with poor enzyme activity. Thus, the StEP recombination method successfully recombined mutations into two xylanases that were more robust than their parent counterparts. Additionally, the 3D-models of the wild type T. lanuginosus xynA (xyl_ext) and its variants, G41 and S325, were predicted using I-TASSER and then subjected to molecular dynamics (MD) simulations at 300 K for a deeper understanding of their structural features. The results from the predicted 3D models show clearly the presence of α-helical regions in the N-terminal residues of the xyl_ext, G41 and S325. Moreover, the MD analysis suggests that the presence of additional residues (1-31) and point mutation induces slight structural changes with the stability of the protein being evenly distributed over the whole structure.

C-terminal flanking peptide stabilized the catalytic domain of a recombinant Bacillus subtilis endo-β-1, 4-glucanase

Three proteins, Egl330, Egl326 and Egl325, which covered the catalytic domain of a Bacillus subtilis endo-β-l, 4-glucanase were expressed in Escherichia coli and purified. Egl325 was a mutant of Egl330 with the peptide sequence Arg-Glu-Asn-Ile-Arg deleted in the C-terminus and Egl326 was another mutant of Egl330 with the peptide sequence Glu-Asn-Ile-Arg deleted in the C-terminus. These three proteins displayed same optimal reaction pH and temperature. However, the thermal stability and pH stability of Egl326 and Egl325 were diminished compared to Egl330. Results of ultra violet scanning, circular dichroism and Trp fluorescence spectrometry showed that the absence of the short peptide at the C-terminus of Egl330 resulted in the destabilization of the catalytic domain through affecting the folding of the protein.

Hydrolase stabilization via entanglement in poly(propylene sulfide) nanoparticles: stability towards reactive oxygen species

In the advancement of green syntheses and sustainable reactions, enzymatic biocatalysis offers extremely high reaction rates and selectivity that goes far beyond the reach of chemical catalysts; however, these enzymes suffer from typical environmental constraints, e.g. operational temperature, pH and tolerance to oxidative environments. A common hydrolase enzyme, diisopropylfluorophosphatase (DFPase, EC 3.1.8.2), has demonstrated a pronounced efficacy for the hydrolysis of a variety of substrates for potential toxin remediation, but suffers from the aforementioned limitations. As a means to enhance DFPase's stability in oxidative environments, enzymatic covalent immobilization within the polymeric matrix of poly(propylene sulfide) (PPS) nanoparticles was performed. By modifying the enzyme's exposed lysine residues via thiolation, DFPase is utilized as a comonomer/crosslinker in a mild emulsion polymerization. The resultant polymeric polysulfide shell acts as a 'sacrificial barrier' by first oxidizing to polysulfoxides and polysulfones, rendering DFPase in an active state. DFPase-PPS nanoparticles thus retain activity upon exposure to as high as 50 parts per million (ppm) of hypochlorous acid (HOCl), while native DFPase is observed as inactive at 500 parts per billion (ppb). This trend is also confirmed by enzyme-generated (chloroperoxidase (CPO), EC 1.11.1.10) reactive oxygen species (ROS) including both HOCl (3 ppm) and ClO(2) (100 ppm).

Introduction to the field of enzyme immobilization and stabilization

Enzyme stabilization is important for any biomedical or industrial application of enzymes. In many applications, the goal is to provide extended active lifetime at normal environmental conditions with traditional substrates at low concentrations in buffered solutions. However, as enzymes are used for more and more applications, there is a desire to use them in extreme environmental conditions (i.e., high temperatures), in high substrate concentration, and in nontraditional solvent systems. This chapter introduces the topic of enzyme stabilization and the methods used for enzyme stabilization including enzyme immobilization.

Stabilization by fusion to the C-terminus of hyperthermophile Sulfolobus tokodaii RNase HI: a possibility of protein stabilization tag

RNase HI from the hyperthermophile Sulfolobus tokodaii (Sto-RNase HI) is stabilized by its C-terminal residues. In this work, the stabilization effect of the Sto-RNase HI C-terminal residues was investigated in detail by thermodynamic measurements of the stability of variants lacking the disulfide bond (C58/145A), or the six C-terminal residues (ΔC6) and by structural analysis of ΔC6. The results showed that the C-terminal does not affect overall structure and stabilization is caused by local interactions of the C-terminal, suggesting that the C-terminal residues could be used as a “stabilization tag.” The Sto-RNase HI C-terminal residues (-IGCIILT) were introduced as a tag on three proteins. Each chimeric protein was more stable than its wild-type protein. These results suggested the possibility of a simple stabilization technique using a stabilization tag such as Sto-RNase HI C-terminal residues.

Error-prone PCR of a fungal xylanase for improvement of its alkaline and thermal stability

Random mutagenesis was used to improve the alkaline and thermal stability of the xylanase (XynA) from Thermomyces lanuginosus. Error-prone PCR reactions were carried out; the PCR products were cloned into Escherichia coli and a library of 960 clones was selected on xylan-containing agar plates. The crude filtrates of positive xylanase producers were screened at 80 degrees C and tested separately at pH 10 for alkaline tolerance. The native XynA lost 80% activity after 90 min at 80 degrees C and lost 70% activity at pH 10. Conversely, the most thermostable variant, G41, retained 75% activity after 90 min at 80 degrees C and the best alkali-stable variant, G53, retained 93% activity at pH 10. Sequence analysis revealed four amino acid substitutions in G41 and a single substitution in G53. These variants, therefore, have improved thermal and alkaline stability and are ideal candidates for DNA shuffling experiments to produce a robust xylanase for industrial application.

Potential implications of availability of short amino acid sequences in proteins: an old and new approach to protein decoding and design

Three-dimensional structure of a protein molecule is primarily determined by its amino acid sequence, and thus the elucidation of general rules embedded in amino acid sequences is of great importance in protein science and engineering. To extract valuable information from sequences, we propose an analytical method in which a protein sequence is considered to be constructed by serial superimpositions of short amino acid sequences of n amino acid sets, especially triplets (3-aa sets). Using the comprehensive nonredundant protein database, we first examined "availability" of all possible combinatorial sets of 8,000 triplet species. Availability score was mathematically defined as an indicator for the relative "preference" or "avoidance" for a given short constituent sequence to be used in protein chain. Availability scores of real proteins were clearly biased against those of randomly generated proteins. We found many triplet species that occurred in the database more than expected or less than expected. Such bias was extended to longer sets, and we found that some species of pentats (5-aa sets) that occurred reasonably frequently in the randomly generated protein population did not occur at all in any real proteins known today. Availability score was dependent on species, potentially serving as a phylogenetic indicator. Furthermore, we suggest possibilities of various biotechnological applications of characteristic short sequences such as human-specific and pathogen-specific short sequences obtained from availability analysis. Availability score was also dependent on secondary structures, potentially serving as a structural indicator. Availability analysis on triplets may be combined with a comprehensive data collection on the varphi and psi peptide-bond angles of the amino acid at the center of each triplet, i.e., a collection of Ramachandran plots for each triplet. These triplet characters, together with other physicochemical data, will provide us with basic information between protein sequence and structure, by which structure prediction and engineering may be greatly facilitated. Availability analysis may also be useful in identifying word processing units in amino acid sequences based on an analogy to natural languages. Together with other approaches, availability analysis will elucidate general rules hidden in the primary sequences and eventually contributes to rebuilding the paradigm of protein science.

A comparative molecular dynamics study of thermophilic and mesophilic ribonuclease HI enzymes

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.

The conversion of 3' UTRs into coding regions

A possible origin of novel coding sequences is the removal of stop codons, leading to the inclusion of 3' untranslated regions (3' UTRs) within genes. We classified changes in the position of stop codons in closely related Saccharomyces species and in a mouse/rat comparison as either additions to or subtractions from coding regions. In both cases, the position of stop codons is highly labile, with more subtractions than additions found. The subtraction bias may be balanced by the input of new coding regions through gene duplication. Saccharomyces shows less stop codon lability than rodents, probably due to greater selective constraint. A higher proportion of 3' UTR incorporation events preserve frame in Saccharomyces. This higher proportion is consistent with the action of the [PSI(+)] prion as an evolutionary capacitor to facilitate 3' UTR incorporation in yeast.

Directed evolution of the thermostable xylanase from Thermomyces lanuginosus

The thermostability of the endo-beta-1,4-xylanase from Thermomyces lanuginosus (xynA) was improved by directed evolution using error-prone PCR. Transformants expressing the variant xylanases were first selected on 0.4% Remazol Brilliant Blue-xylan and then exposed to 80 degrees C. Whereas the wild type XynA lost 90% activity after 10 min at 80 degrees C, five mutants displayed both higher stabilities and activities than XynA. Four mutants were subjected to further mutagenesis to improve the stability and activity of the xylanase. Subsequent screening revealed three mutants with enhanced thermostability. Mutant 2B7-10 retained 71% of its activity after treatment at 80 degrees C for 60 min and had a half-life of 215 min at 70 degrees C, which is higher than that attained by XynA. Sequence analysis of second generation mutants revealed that mutations were not concentrated in any particular region of the protein and exhibited much variation. The best mutant obtained from this study was variant 2B7-10, which had a single substitution (Y58F) in beta-sheet A of the protein, which is the hydrophilic, solvent-accessible outer surface of the enzyme. Most of the mutants obtained in this study displayed a compromise between stability and activity, the only exception being mutant 2B7-10. This variant showed increased activity and thermostability.

In vitro evolution of proteins

Consecutive rounds of diversification and selection of the fittest is believed to be the main driving force for the evolution of life. For the evolution of life to proceed, all living cells are surrounded by a lipid bilayer that separates their own genes from the external environment and from those of other organisms. In this way, the genetic information of an individual is replicated on the basis of their phenotype; thus the enrichment of the fittest will occur. Hence, evolution is based on linkage between genotype and phenotype owing to the surrounding of the genetic material with a barrier. The linkage between genotype and phenotype is also known to be essential for the directed evolution of proteins. Indeed, systems for molecular evolution, including phage display, ribosome display, and in vitro compartmentalization, all satisfy this requirement in different ways. These systems have been shown to be powerful tools for high-throughput screening for the functions of proteins, screening as many as <10(12) molecules in 1 d. These selection systems in combination with various gene libraries yield proteins with improved or altered biophysical properties, and may even allow the generation of proteins with novel functions.

HIV protease-activated molecular switches based on beta-glucuronidase and alkaline phosphatase

Our long-term goal is to direct the evolution of novel protease variants. To this end we have engineered a new type of protease-activated reporter enzyme. Many protease-activated enzymes evolved in nature, but the introduction of novel regulatory mechanisms into normally unregulated enzymes poses a difficult design challenge. Random Elongation Mutagenesis [1] was used to fuse the p6 peptide, which is recognized and cleaved by HIV protease, and twelve random sequence amino acids to the C-termini of beta-glucuronidase (GUS) and alkaline phosphatase (AP). The resulting GUS-p6-(NNN)12 and AP-p6-(NNN)12 libraries were expressed in E. coli and screened for clones that were inactivated by the C-terminal extension (tail). The inactivated clones were co-expressed with HIV protease, and those that were re-activated were isolated. The AP and GUS activities of the most responsive clones were each >3.5-fold higher when co-expressed with HIV protease, and this activation is correlated with in vivo proteolysis. It should be possible to generalize this strategy to different reporter enzymes, different target proteases, and perhaps to other types of protein-modifying enzymes.

The structure of a eukaryotic nicotinic acid phosphoribosyltransferase reveals structural heterogeneity among type II PRTases

Nicotinamide adenine dinucleotide (NAD) is an essential cofactor for cellular redox reactions and can act as an important substrate in numerous biological processes. As a result, nature has evolved multiple biosynthetic pathways to meet this high chemical demand. In Saccharomyces cerevisiae, the NAD salvage pathway relies on the activity of nicotinic acid phosphoribosyltransferase (NAPRTase), a member of the phosphoribosyltransferase (PRTase) superfamily. Here, we report the structure of a eukaryotic (yeast) NAPRTase at 1.75 A resolution (locus name: YOR209C, gene name: NPT1). The structure reveals a two-domain fold that resembles the architecture of quinolinic acid phosphoribosyltransferases (QAPRTases), but with completely different dispositions that provide evidence for structural heterogeneity among the Type II PRTases. The identification of a third domain in NAPRTases provides a structural basis and possible mechanism for the functional modulation of this family of enzymes by ATP.

Engineering novel binding proteins from nonimmunoglobulin domains

Not all adaptive immune systems use the immunoglobulin fold as the basis for specific recognition molecules: sea lampreys, for example, have evolved an adaptive immune system that is based on leucine-rich repeat proteins. Additionally, many other proteins, not necessarily involved in adaptive immunity, mediate specific high-affinity interactions. Such alternatives to immunoglobulins represent attractive starting points for the design of novel binding molecules for research and clinical applications. Indeed, through progress and increased experience in library design and selection technologies, gained not least from working with synthetic antibody libraries, researchers have now exploited many of these novel scaffolds as tailor-made affinity reagents. Significant progress has been made not only in the basic science of generating specific binding molecules, but also in applications of the selected binders in laboratory procedures, proteomics, diagnostics and therapy. Challenges ahead include identifying applications where these novel proteins can not only be an alternative, but can enable approaches so far deemed technically impossible, and delineate those therapeutic applications commensurate with the molecular properties of the respective proteins.

Programmable assembly of nanoarchitectures using genetically engineered viruses

Biological systems possess inherent molecular recognition and self-assembly capabilities and are attractive templates for constructing complex material structures with molecular precision. Here we report the assembly of various nanoachitectures including nanoparticle arrays, hetero-nanoparticle architectures, and nanowires utilizing highly engineered M13 bacteriophage as templates. The genome of M13 phage can be rationally engineered to produce viral particles with distinct substrate-specific peptides expressed on the filamentous capsid and the ends, providing a generic template for programmable assembly of complex nanostructures. Phage clones with gold-binding motifs on the capsid and streptavidin-binding motifs at one end are created and used to assemble Au and CdSe nanocrytals into ordered one-dimensional arrays and more complex geometries. Initial studies show such nanoparticle arrays can further function as templates to nucleate highly conductive nanowires that are important for addressing/interconnecting individual nanostructures.

Testing the relationship between foldability and the early folding events of dihydrofolate reductase from Escherichia coli

A "folding element" is a contiguous peptide segment crucial for a protein to be foldable and is a new concept that could assist in our understanding of the protein-folding problem. It is known that the presence of the complete set of folding elements of dihydrofolate reductase (DHFR) from Escherichia coli is essential for the protein to be foldable. Since almost all of the amino acid residues known to be involved in the early folding events of DHFR are located within the folding elements, a close relationship between the folding elements and early folding events is hypothesized. In order to test this hypothesis, we have investigated whether or not the early folding events are preserved in circular permutants and topological mutants of DHFR, in which the order of the folding elements is changed but the complete set of folding elements is present. The stopped-flow circular dichroism (CD) measurements show that the CD spectra at the early stages of folding are similar among the mutants and the wild-type DHFR, indicating that the presence of the complete set of folding elements is sufficient to preserve the early folding events. We have further examined whether or not sequence perturbation on the folding elements by a single amino acid substitution affects the early folding events of DHFR. The results show that the amino acid substitutions inside of the folding elements can affect the burst-phase CD spectra, whereas the substitutions outside do not. Taken together, these results indicate that the above hypothesis is true, suggesting a close relationship between the foldability of a protein and the early folding events. We propose that the folding elements interact with each other and coalesce to form a productive intermediate(s) early in the folding, and these early folding events are important for a protein to be foldable.

The insertion of palindromic repeats in the evolution of proteins

The current theory of protein evolution is that all contemporary proteins are derived from an ancestral subset. However, each new sequenced genome exhibits many genes with no detectable homologues in other species, leading to the paradoxical picture of a universal ancestor with more genes than any of its progeny. Standard explanations indicate that fast evolving genes might disappear into the 'twilight zone' of sequence similarity. Regardless of the size of the original ancestral subset, its origin and the potential mechanisms of its subsequent enlargement are rarely addressed. Sequencing of Rickettsia conorii genome recently led to the discovery of three families of repeat-mobile elements frequently inserted into the middle of protein coding genes. Although not yet identified in other species of bacteria, this discovery has provided the first clear evidence for the de novo creation of long protein segments (up to 50 amino acid residues) by repeat insertion. Based on previous results and theories on the coding potential of palindromic elements, we speculate that their insertion and mobility might have played a significant role in the early stages of protein evolution.

Directed evolution of barnase stability using proteolytic selection

We report the construction of a phage-displayed repertoire of mutants of the ribonuclease barnase from Bacillus amyloliquefaciens. The construction was guided by the natural variability between two closely related ribonucleases, barnase and binase from Bacillus intermedius. This repertoire was selected using a proteolytic selection method, allowing sorting of the library according to the resistance of the mutants toward proteolysis. Susceptibility toward proteolysis has been correlated with flexibility and unfolding, and is thus expected to yield mutants with increased thermal stability. Enrichment of barnase mutants with specific combinations of amino acid residues at four of the randomised positions was observed. Three of these enriched amino acid residues are present in neither barnase nor binase. For some of the mutations, the improvement in proteolytic stability does not lead to a pronounced improvement in thermodynamic stability, indicating that the factors governing the proteolytic stability in some cases may be different from those governing the thermodynamic stability, e.g. propensity to local unfolding.The results obtained add important knowledge to a novel use of phage display technology for selection of thermodynamically stable proteins. Only by carefully establishing the parameters that can be adjusted, and recognising the influence this will have on the outcome of selection, will it be possible to realise the powerful technique of proteolytic selection.

Molecular evolution in static and dynamical landscapes

In an attempt to understand protein evolution, we address the issues ofhow much variety in the sequences is needed to prompt the evolution ofan enzyme from random polypeptides and how does cellular interactionaffect the dynamics of molecular evolution to allow genetic diversity inpopulation. The experimental evolution of phage-displayed randompolypeptides of about 140 amino acid residues panned with transition stateanalogue for an esterase reaction showed that even with a population sizeas small as ten, not only could significant varieties be found but also therandom polypeptides in each of the generation had great promise towardsdeveloping into functional proteins. Hence, it is evident that the enzymeevolution is prompted even within a small local area of the static landscapeof the sequence space. Considering that interaction among living cells is aninevitable event in natural evolution, its role was investigated through threeconsecutive rounds of random mutagenesis on the glutamine synthetasegene and chemostat culture of the transformed Escherichia colicellscontaining the mutated genes. The molecular phylogeny and populationdynamics show the coexistence of some mutants having different level ofglutamine synthetase at each generation. In addition, it was confirmed thatcellular interaction via the medium influences the stability of the coexistenceand bring forth fitness change to the coexisting members of the population,thereby, leading to a dynamical landscape. Based on experimental resultsreflecting the extent of interaction among members in population, here, Iproposed that protein evolution could change its mode from theoptimization on static landscape to diversification on dynamicallandscape.

Improved partitioning in aqueous two-phase system of tyrosine-tagged recombinant lactate dehydrogenase

The partitioning of Bacillus stearothermophilus lactate dehydrogenase (LDH) in an aqueous two-phase system was studied. Particularly, the influence of tyrosine tags on the partitioning was evaluated. The hydrophobic effect, caused by the addition of tyrosine residues, was determined in a system based on dextran and the thermoseparating ethylene oxide-propylene oxide random copolymer (EO30PO70). Five different LDH variants were constructed with N-terminal tags containing tyrosines (Y3 and Y6), tyrosines and prolines (Y3P2 and Y6P2), and only prolines (P2). LDH fused with tags containing tyrosines increased the partitioning coefficient, and the more tyrosines added to the protein, the larger improvement in partitioning. When prolines were added between the tyrosine-rich tag and the protein, a further increased partitioning was obtained. The enhanced partitioning was attributed to the rigid structure of the proline, which in turn led to an increase in the exposure of the tag to the surroundings. The best tyrosine tag, Y6P2, increased the partition coefficient four times and additionally, a higher thermostability was observed.

Evolvability of random polypeptides through functional selection within a small library

A directed evolution with phage-displayed random polypeptides of about 140 amino acid residues was followed until the sixth generation under a selection based on affinity to a transition state analog for an esterase reaction. The experimental design deliberately limits the observation to only 10 clones per generation. The first generation consists of three soluble random polypeptides and seven arbitrarily chosen clones from a previously constructed library. The clone showing the highest affinity in a generation was selected and subjected to random mutagenesis to generate variants for the next generation. Even within only 10 arbitrarily chosen polypeptides in each of the generations, there are enough variants in accord to capacity of binding affinity. In addition, the binding capacity of the selected polypeptides showed a gradual continuous increase over the generation. Furthermore, the purified selected random polypeptides exhibited a gradual but significant increase in esterase activity. The ease of the functional development within a small sequence variety implies that enzyme evolution is prompted even within a small population of random polypeptides.

Oxidative burst by acellular haemoglobin and neurotransmitters

Acellular haemoglobin (Hb) has intrinsic toxicity to the tissues since harmful reactive oxygen species are readily produced during auto-oxidation of Hb. On the other hand, Hb is known to have peroxidase-like activity monovalently oxidizing various peroxidase substrates. Thus, monovalently oxidized organic free radical species may be produced. This may relay the radical reactions leading to the production of reactive oxygen species such as superoxide. Such substrates possibly generating superoxide, include aromatic monoamines such as neurotransmitters and their precursors rich in neural a tissues. Based on our knowledge on the reactivity of haemoproteins against phenolics and aromatic monoamines, we proposed a hindered danger in use of Hb as a reperfusion agent. Clinical use of recently developing Hb-based blood substitutes must be reconsidered.

Thermoadaptation of alpha-galactosidase AgaB1 in Thermus thermophilus

The evolutionary potential of a thermostable alpha-galactosidase, with regard to improved catalytic activity at high temperatures, was investigated by employing an in vivo selection system based on thermophilic bacteria. For this purpose, hybrid alpha-galactosidase genes of agaA and agaB from Bacillus stearothermophilus KVE39, designated agaA1 and agaB1, were cloned into an autonomously replicating Thermus vector and introduced into Thermus thermophilus OF1053GD (DeltaagaT) by transformation. This selector strain is unable to metabolize melibiose (alpha-galactoside) without recombinant alpha-galactosidases, because the native alpha-galactosidase gene, agaT, has been deleted. Growth conditions were established under which the strain was able to utilize melibiose as a single carbohydrate source when harboring a plasmid-encoded agaA1 gene but unable when harboring a plasmid-encoded agaB1 gene. With incubation of the agaB1 plasmid-harboring strain under selective pressure at a restrictive temperature (67 degrees C) in a minimal melibiose medium, spontaneous mutants as well as N-methyl-N'-nitro-N-nitrosoguanidine-induced mutants able to grow on the selective medium were isolated. The mutant alpha-galactosidase genes were amplified by PCR, cloned in Escherichia coli, and sequenced. A single-base substitution that replaces glutamic acid residue 355 with glycine or valine was found in the mutant agaB1 genes. The mutant enzymes displayed the optimum hydrolyzing activity at higher temperatures together with improved catalytic capacity compared to the wild-type enzyme and furthermore showed an enhanced thermal stability. To our knowledge, this is the first report of an in vivo evolution of glycoside-hydrolyzing enzyme and selection within a thermophilic host cell.

Increasing the thermal stability of an oligomeric protein, beta-glucuronidase

The reporter enzyme beta-glucuronidase was mutagenized and evolved for thermostability. After four cycles of screening the best variant was more active than the wild-type enzyme, and retained function at 70 degrees C, whereas the wild-type enzyme lost function at 65 degrees C. Variants derived from sequential mutagenesis were shuffled together, and re-screened for thermostability. The best variants retained activities at even higher temperatures (80 degrees C), but had specific activities that were now less than that of the wild-type enzyme. The mutations clustered near the tetramer interface of the enzyme, and many of the evolved variants showed much greater resistance to quaternary structure disruption at high temperatures, which is also a characteristic of naturally thermostable enzymes. Together, these results suggest a pathway for the evolution of thermostability in which enzymes initially become stable at high temperatures without loss of activity at low temperatures, while further evolution leads to enzymes that have kinetic parameters that are optimized for high temperatures.

In vitro adaptation of a ligase ribozyme for activity under a low-pH condition

A ligase ribozyme accelerating a ligation reaction with oligonucleotide under a low-pH condition was selected by in vitro adaptation. A ribozyme active at pH 7 was randomly mutated, and the resultant RNA library was subjected to in vitro adaptation under a low-pH reaction condition. At pH 4, the adapted RNAs reacted with the oligonucleotide substrates about 200 times faster than the original ribozyme. When the ribozyme was cloned and sequenced, 10 of the 30 clones sequenced had identical sequences. The differences in sequence from the original ribozyme were found at four positions in the middle region and at the 3' end. A few sequential differences dominated the activity of the ribozyme under the extreme condition. The adapted ribozyme had one repeating sequence that was critical for the activity.

Molecular cloning and sequence analysis of the manganese catalase gene from Thermoleophilum album NM

The manganese catalase gene (mnct) from Thermoleophilum album NM, a thermophilic bacterium, was cloned and its nucleotide sequence was analyzed. The gene consists of 885 bp (65.4% GC content) encoding 294 amino acids with a molecular mass of 32,500 Da. The deduced amino acid sequence shows similarities to those of Thermus species strain YS 8-13 (a thermophilic bacterium) and Bacillus halodurans (an alkaliphilic bacterium) with 61 and 54% identities, respectively.

Generation of protein lineages with new sequence spaces by functional salvage screen

A variety of different methods to generate diverse proteins, including random mutagenesis and recombination, are currently available and most of them accumulate the mutations on the target gene of a protein, whose sequence space remains unchanged. On the other hand, a pool of diverse genes, which is generated by random insertions, deletions and exchange of the homologous domains with different lengths in the target gene, would present the protein lineages resulting in new fitness landscapes. Here we report a method to generate a pool of protein variants with different sequence spaces by employing green fluorescent protein (GFP) as a model protein. This process, designated functional salvage screen (FSS), comprises the following procedures: a defective GFP template expressing no fluorescence is first constructed by genetically disrupting a predetermined region(s) of the protein and a library of GFP variants is generated from the defective template by incorporating the randomly fragmented genomic DNA from Escherichia coli into the defined region(s) of the target gene, followed by screening of the functionally salvaged, fluorescence-emitting GFPs. Two approaches, sequence-directed and PCR-coupled methods, were attempted to generate the library of GFP variants with new sequences derived from the genomic segments of E.coli. The functionally salvaged GFPs were selected and analyzed in terms of the sequence space and functional properties. The results demonstrate that the functional salvage process not only can be a simple and effective method to create protein lineages with new sequence spaces, but also can be useful in elucidating the involvement of a specific region(s) or domain(s) in the structure and function of protein.

The catalytic domain of a Piromyces rhizinflata cellulase expressed in Escherichia coli was stabilized by the linker peptide of the enzyme

Analysis of a carboxymethyl-cellulase clone isolated from the cDNA library of the ruminal fungus, Piromyces rhizinflata 2301, revealed that the clone encoded a polypeptide containing a cellulase catalytic domain, designated CelAcd. CelAcd was flanked by a 28-amino acid linker peptide at the N-terminus and linked to a dockerin domain by a 7-amino acid linker at the C-terminus. CelAcd showed homology with the family 5 of glycosyl hydrolases. CelAcd plus the linker peptides at both termini (designated CelcdN'C') was expressed in Escherichia coli and the purified enzyme was characterized. The feature of particular interest of CelcdN'C' was its bifunctional endo- and exo-glucanase activity, demonstrated by its ability to hydrolyse carboxymethyl cellulose (CMC) and pNP-beta-D-cellobioside. Furthermore, CelcdN'C' exhibited relatively high ability to degrade both microcrystalline Avicel and filter paper. CelcdN'C' also showed activity against barley beta-glucan, Lichenin and oat spelt xylan. The optimal activity conditions for CelcdN'C' with CMC as the substrate were pH 5.5 and 50 degrees C. Fifty percent of the enzyme activity was lost when CelcdN'C' was treated at 55 degrees C for 10 min. CelcdN'C' retained more than 10% enzyme activity after being heated under 90 degrees C for 10 min. Deletion of the N-terminal flanking linker of CelcdN'C' (the resulting protein was designated CelcdC') did not alter the enzymatic function of the catalytic domain. However, the thermal stability of CelcdC' was dramatically reduced. We conclude that the N-terminal flanking linker of CelAcd stabilizes the enzyme protein.

Structure-based chimeric enzymes as an alternative to directed enzyme evolution: phytase as a test case

Thermostability is a key feature for commercially attractive variants of the fungal enzyme phytase. In an initial set of experiments, we restored ionic interactions and hydrogen bonds on the surface of Aspergillus terreus phytase, which are present in the homologous but more thermostable enzyme from A. niger. Since these mutations turned out to be neutral, we replaced-in the same region and based on the crystal structure of A. niger phytase-entire secondary structure elements. The replacement of one alpha-helix on the surface of A. terreus phytase by the corresponding stretch of A. niger phytase resulted in an enzyme with improved thermostability and unaltered enzymatic activity. Surprisingly, the thermostability of this hybrid protein was very similar to that of A. niger phytase, although the fusion protein contained only a 31 amino acid stretch of the more stable parent enzyme. This report provides evidence that structure-based chimeric enzymes can be used to exploit the evolutionary information within a sequence alignment. We propose this method as an alternative to directed enzyme evolution if due to expression constraints the screening of large mutant populations is not feasible.

Enhancement of the enzymatic activity of ribonuclease HI from Thermus thermophilus HB8 with a suppressor mutation method

A genetic method for isolating a mutant enzyme of ribonuclease HI (RNase HI) from Thermus thermophilus HB8 with enhanced activity at moderate temperatures was developed. T. thermophilus RNase HI has an ability to complement the RNase H-dependent temperature-sensitive (ts) growth phenotype of Escherichia coli MIC3001. However, this complementation ability was greatly reduced by replacing Asp(134), which is one of the active site residues, with His, probably due to a reduction in the catalytic activity. Random mutagenesis of the gene encoding the resultant D134H enzyme, followed by screening for second-site revertants, allowed us to isolate three single mutations (Ala(12) --> Ser, Lys(75) --> Met, and Ala(77) --> Pro) that restore the normal complementation ability to the D134H enzyme. These mutations were individually or simultaneously introduced into the wild-type enzyme, and the kinetic parameters of the resultant mutant enzymes for the hydrolysis of a DNA-RNA-DNA/DNA substrate were determined at 30 degrees C. Each mutation increased the k(cat)/K(m) value of the wild-type enzyme by 2.1-4.8-fold. The effects of the mutations on the enzymatic activity were roughly cumulative, and the combination of these three mutations increased the k(cat)/K(m) value of the wild-type enzyme by 40-fold (5.5-fold in k(cat)). Measurement of thermal stability of the mutant enzymes with circular dichroism spectroscopy in the presence of 1 M guanidine hydrochloride and 1 mM dithiothreitol showed that the T(m) value of the triple mutant enzyme, in which all three mutations were combined, was comparable to that of the wild-type enzyme (75.0 vs 77.4 degrees C). These results demonstrate that the activity of a thermophilic enzyme can be improved without a cost of protein stability.

Directed molecular evolution of cytochrome c peroxidase

Cytochrome c peroxidase (CCP) from Saccharomyces cerevisiae was subjected to directed molecular evolution to generate mutants with increased activity against the classical peroxidase substrate guaiacol, thus changing the substrate specificity of CCP from the protein cytochrome c to a small organic molecule. After three rounds of DNA shuffling and screening, mutants were isolated which possessed a 300-fold increased activity against guaiacol and an up to 1000-fold increased specificity for this substrate relative to that for the natural substrate. In all of the selected mutants, the distal arginine (Arg48), which is fully conserved in the superfamily of peroxidases, was mutated to histidine, showing that this mutation plays a key role in the significant increase in activity against phenolic substrates. The results suggest that, in addition to stabilizing the reactive intermediate compound I, the distal arginine plays an important role as a gatekeeper in the active site of CCP, controlling the access to the ferryl oxygen and the distal histidine. Other isolated mutations increase the general reactivity of the peroxidase or increase the intracellular concentration of the active holo form, allowing their selection under the employed screening conditions. The results illustrate the ability of directed molecular evolution technologies to deliver solutions to biochemical problems that would not be readily predicted by rational design.