Directed evolution of (betaalpha)(8)-barrel enzymes

Natural molecular evolution supplies us with manifold examples of protein engineering. The imitation of these natural processes in the design of new enzymes has led to surprising and insightful results. Well-suited for design by evolutionary methods are enzymes with the common and versatile (betaalpha)(8)-barrel fold. Studies of enzyme stability, folding and design as well as the evolution of (betaalpha)(8)-barrel enzymes are discussed.     

Stabilization of a (betaalpha)8-barrel protein by an engineered disulfide bridge.

The aim of this study was to increase the stability of the thermolabile (betaalpha)8-barrel enzyme indoleglycerol phosphate synthase from Escherichia coli by the introduction of disulfide bridges. For the design of such variants, we selected two out of 12 candidates, in which newly introduced cysteines potentially form optimal disulfide bonds. These variants avoid short-range connections, substitutions near catalytic residues, and crosslinks between the new and the three parental cysteines. The variant linking residues 3 and 189 fastens the N-terminus to the (betaalpha)8-barrel. The rate of thermal inactivation at 50 degrees C of this variant with a closed disulfide bridge is 65-fold slower than that of the reference dithiol form, but only 13-fold slower than that of the parental protein. The near-ultraviolet CD spectrum, the reactivity of parental buried cysteines with Ellman's reagent as well as the decreased turnover number indicate that the protein structure is rigidified. To confirm these data, we have solved the X-ray structure to 2.1-A resolution. The second variant was designed to crosslink the terminal modules betaalpha1 and betaalpha8. However, not even the dithiol form acquired the native fold, possibly because one of the targeted residues is solvent-inaccessible in the parental protein.     

Dissection of a (betaalpha)8-barrel enzyme into two folded halves

The (betaalpha)8-barrel, which is the most frequently encountered protein fold, is generally considered to consist of a single structural domain. However, the X-ray structure of the imidazoleglycerol phosphate synthase (HisF) from Thermotoga maritima has identified it as a (betaalpha) 8-barrel made up of two superimposable subdomains (HisF-N and HisF-C). HisF-N consists of the four N-terminal (betaalpha) units and HisF-C of the four C-terminal (betaalpha) units. It has been postulated, therefore, that HisF evolved by tandem duplication and fusion from an ancestral half-barrel. To test this hypothesis, HisF-N and HisF-C were produced in Escherichia coli, purified and characterized. Separately, HisF-N and HisF-C are folded proteins, but are catalytically inactive. Upon co-expression in vivo or joint refolding in vitro, HisF-N and HisF-C assemble to the stoichiometric and catalytically fully active HisF-NC complex. These findings support the hypothesis that the (betaalpha)8-barrel of HisF evolved from an ancestral half-barrel and have implications for the folding mechanism of the members of this large protein family.     

Stabilisation of a (betaalpha)8-barrel protein designed from identical half barrels

It has been suggested that the common (betaalpha)(8)-barrel enzyme fold has evolved by the duplication and fusion of identical (betaalpha)(4)-half barrels, followed by the optimisation of their interface. In our attempts to reconstruct these events in vitro we have previously linked in tandem two copies of the C-terminal half barrel HisF-C of imidazole glycerol phosphate synthase from Thermotoga maritima and subsequently reconstituted in the fusion construct HisF-CC a salt bridge cluster present in wild-type HisF. The resulting recombinant protein HisF-C*C, which was produced in an insoluble form and unfolded with low cooperativity at moderate urea concentrations has now been stabilised and solubilised by a combination of random mutagenesis and selection in vivo. For this purpose, Escherichia coli cells were transformed with a plasmid-based gene library encoding HisF-C*C variants fused to chloramphenicol acetyltransferase (CAT). Stable and soluble variants were identified by the survival of host cells on solid medium containing high concentrations of the antibiotic. The selected HisF-C*C proteins, which were characterised in vitro in the absence of CAT, contained eight different amino acid substitutions. One of the exchanges (Y143C) stabilised HisF-C*C by the formation of an intermolecular disulfide bond. Three of the substitutions (G245R, V248M, L250Q) were located in the long loop connecting the two HisF-C copies, whose subsequent truncation from 13 to 5 residues yielded the stabilised variant HisF-C*C Delta. From the remaining substitutions, Y143H and V234M were most beneficial, and molecular dynamics simulations suggest that they strengthen the interactions between the half barrels by establishing a hydrogen-bonding network and an extensive hydrophobic cluster, respectively. By combining the loop deletion of HisF-C*C Delta with the Y143H and V234M substitutions, the variant HisF-C**C was generated. Recombinant HisF-C**C is produced in soluble form, forms a pure monomer with its tryptophan residues shielded from solvent and unfolds with similar cooperativity as HisF. Our results show that, starting from two identical and fused half barrels, few amino acid exchanges are sufficient to generate a highly stable and compact (betaalpha)(8)-barrel protein with wild-type like structural properties.     

Interconverting the catalytic activities of (betaalpha)(8)-barrel enzymes from different metabolic pathways: sequence requirements and molecular analysis

The (betaalpha)(8)-barrel enzymes N'-[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide isomerase (tHisA) and imidazole glycerol phosphate synthase (tHisF) from Thermotoga maritima catalyze two successive reactions in the biosynthesis of histidine. In both enzymes, aspartate residues at the C-terminal end of beta-strand 1 (Asp8 in tHisA and Asp11 in tHisF) and beta-strand 5 (Asp127 in tHisA and Asp130 in tHisF) are essential for catalytic activity. It was demonstrated earlier that in tHisA the substitution of Asp127 by valine (tHisA-D127V) generates phosphoribosylanthranilate isomerase (TrpF) activity, a related (betaalpha)(8)-barrel enzyme participating in tryptophan biosynthesis. It is shown here that in tHisF the corresponding substitution of Asp130 by valine (tHisF-D130V) also generates TrpF activity. To determine the effectiveness of individual amino acid exchanges in these conversions, each of the 20 standard amino acid residues was introduced at position 127 of tHisA and 130 of tHisF by saturation random mutagenesis. The tHisA-D127X and tHisF-D130X variants with TrpF activity were identified by selection in vivo, and the proteins purified and characterized. The results obtained show that removal of the negatively charged carboxylate side-chain at the C-terminal end of beta-strand 5 is sufficient to establish TrpF activity in tHisA and tHisF, presumably because it allows the binding of the negatively charged TrpF substrate, phosphoribosylanthranilate. In contrast, the double mutants tHisA-D8N+D127V and tHisF-D11N+D130V did not show detectable activity, demonstrating that the aspartate residues at the C-terminal end of beta-strand 1 are essential for catalysis of the TrpF reaction. The ease with which TrpF activity can be established on both the tHisA and tHisF scaffolds supports the evolutionary relationship of these three enzymes and highlights the functional plasticity of the (betaalpha)(8)-barrel enzyme fold.     

Evolution of function in (beta/alpha)8-barrel enzymes

The (beta/alpha)(8)-barrel is the most common fold in structurally characterized enzymes. Whether the functionally diverse enzymes that share this fold are the products of either divergent or convergent evolution (or both) is an unresolved question that will probably be answered as the sequence databases continue to expand. Recent work has examined natural, designed, and directed evolution of function in several superfamilies of (beta/alpha)(8)-barrel containing enzymes.     

Two (betaalpha)(8)-barrel enzymes of histidine and tryptophan biosynthesis have similar reaction mechanisms and common strategies for protecting their labile substrates

The enzymes N'-[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide isomerase (HisA) and phosphoribosylanthranilate isomerase (TrpF) are sugar isomerases that are involved in histidine and tryptophan biosynthesis, respectively. Both enzymes have the (betaalpha)(8)-barrel fold and catalyze Amadori rearrangements of a thermolabile aminoaldose into the corresponding aminoketose. To identify those amino acids that are essential for catalysis, conserved residues at the active sites of both HisA and TrpF from the hyperthermophile Thermotoga maritima were replaced by site-directed mutagenesis, and the purified variants were investigated by steady-state enzyme kinetics. Aspartate 8, aspartate 127, and threonine 164 appeared to be important for the HisA reaction, whereas cysteine 7 and aspartate 126 appeared to be important for the TrpF reaction. On the basis of these results and the X-ray structure of a complex between TrpF and a bound product analogue, a reaction mechanism involving general acid-base catalysis and a Schiff base intermediate is proposed for both enzymes. A comparison of the HisA and TrpF enzymes from T. maritima and Escherichia coli showed that, at the physiological temperatures of 80 and 37 degrees C, respectively, the enzymes from the hyperthermophile have significantly higher catalytic efficiencies than the corresponding enzymes from mesophiles. These results suggest that HisA and TrpF have similar chemical reaction mechanisms and use the same strategy to prevent the loss of their thermolabile substrates.     

Stability, catalytic versatility and evolution of the (beta alpha)(8)-barrel fold

The (beta alpha)(8)-barrel is a versatile single-domain protein fold that is adopted by a large number of enzymes. The (beta alpha)(8)-barrel fold has been used as a model to elucidate the structural basis of protein thermostability and in studies to interconvert catalytic activities or substrate specificities by rational design or directed evolution. Recently, the (beta alpha)(4)-half-barrel was identified as a possible structural subdomain.     

Crystal structure of an archaeal class I aldolase and the evolution of (betaalpha)8 barrel proteins

Fructose-1,6-bisphosphate aldolase (FBPA) catalyzes the reversible cleavage of fructose 1,6-bisphosphate to glyceraldehyde 3-phosphate and dihydroxyacetone phosphate in the glycolytic pathway. FBPAs from archaeal organisms have recently been identified and characterized as a divergent family of proteins. Here, we report the first crystal structure of an archaeal FBPA at 1.9-A resolution. The structure of this 280-kDa protein complex was determined using single wavelength anomalous dispersion followed by 10-fold non-crystallographic symmetry averaging and refined to an R-factor of 14.9% (Rfree 17.9%). The protein forms a dimer of pentamers, consisting of subunits adopting the ubiquitous (betaalpha)8 barrel fold. Additionally, a crystal structure of the archaeal FBPA covalently bound to dihydroxyacetone phosphate was solved at 2.1-A resolution. Comparison of the active site residues with those of classical FBPAs, which share no significant sequence identity but display the same overall fold, reveals a common ancestry between these two families of FBPAs. Structural comparisons, furthermore, establish an evolutionary link to the triosephosphate isomerases, a superfamily hitherto considered independent from the superfamily of aldolases.     

Homology among (betaalpha)(8) barrels: implications for the evolution of metabolic pathways

We provide statistically reliable sequence evidence indicating that at least 12 of 23 SCOP (betaalpha)(8) (TIM) barrel superfamilies share a common origin. This includes all but one of the known and predicted TIM barrels found in central metabolism. The statistical evidence is complemented by an examination of the details of protein structure, with certain structural locations favouring catalytic residues even though the nature of their molecular function may change. The combined analysis of sequence, structure and function also enables us to propose a phylogeny of TIM barrels. Based on these data, we are able to examine differing theories of pathway and enzyme evolution, by mapping known TIM barrel folds to the pathways of central metabolism. The results favour widespread recruitment of enzymes between pathways, rather than a "backwards evolution" model, and support the idea that modern proteins may have arisen from common ancestors that bound key metabolites.     

Role of the N-terminal extension of the (betaalpha)8-barrel enzyme indole-3-glycerol phosphate synthase for its fold, stability, and catalytic activity

Indole-3-glycerol phosphate synthase (IGPS) catalyzes the fifth step in the biosynthesis of tryptophan. It belongs to the large and versatile family of (betaalpha)(8)-barrel enzymes but has an unusual N-terminal extension of about 40 residues. Limited proteolysis with trypsin of IGPS from both Sulfolobus solfataricus (sIGPS) and Thermotoga maritima (tIGPS) removes about 25 N-terminal residues and one of the two extra helices contained therein. To assess the role of the extension, the N-terminally truncated variants sIGPSDelta(1-26) and tIGPSDelta(1-25) were produced recombinantly in Escherichia coli, purified, and characterized in comparison to the wild-type enzymes. Both sIGPSDelta(1-26) and tIGPSDelta(1-25) have unchanged oligomerization states and turnover numbers. In contrast, their Michaelis constants for the substrate 1-(o-carboxyphenylamino)-1-deoxyribulose 5-phosphate are increased, and their resistance toward unfolding induced by heat and guanidinium chloride is decreased. sIGPSDelta(1-26) was crystallized, and its X-ray structure was solved at 2.8 A resolution. The comparison with the known structure of sIGPS reveals small differences that account for its reduced substrate affinity and protein stability. The structure of the core of sIGPSDelta(1-26) is, however, unchanged compared to sIGPS, explaining its retained catalytic activity and consistent with the idea that it evolved from the same ancestor as the phosphoribosyl anthranilate isomerase and the alpha-subunit of tryptophan synthase. These (betaalpha)(8)-barrel enzymes catalyze the reactions preceding and following IGPS in tryptophan biosynthesis but lack an N-terminal extension.     

Evolutionary markers in the (beta/alpha)8-barrel fold

Enzymes with the (beta/alpha)(8)-barrel fold are involved in the catalysis of a wide variety of biochemical reactions. The active sites of these enzymes are located on the C-terminal face of the central beta-barrel. Conserved amino acid sequence, as well as secondary, tertiary and quaternary structure patterns are providing a rich body of data to support the premise of a common ancestry of many members of the (beta/alpha)(8)-barrel fold family of enzymes. Recent data indicate that there is at least one example of a bienzyme that functions as an ammonia channel, adding a new level of functional diversity to the (beta/alpha)(8)-barrel fold. These proteins have become ideal tools that can be used in conjunction with directed evolution techniques to engineer novel catalytic activities.     

(Beta alpha)8-barrel proteins of tryptophan biosynthesis in the hyperthermophile Thermotoga maritima.

To better understand the evolution of a key metabolic pathway, we have sequenced the trpCFBA gene cluster of the hyperthermophilic bacterium Thermotoga maritima. The genes were cloned by complementation in vivo of trp deletion strains of Escherichia coli. The new sequences, together with earlier findings, establish that the trp operon of T.maritima has the order trpE(G.D)CFBA, which might represent the ancestral organization of the tryptophan operon. Heterologous expression of the trp(G.D) and trpC genes in E.coli and N-terminal sequencing of their polypeptide products showed that their translation is initiated at the rate start codons TTG and ATC, respectively. Consequently, the N-terminus of the trp(G.D) fusion protein is 43 residues shorter than previously postulated. Amino acid composition and sequence analyses of the protein products of T.maritima trpC (indoleglycerol phosphate synthase), trpF (phosphoribosyl anthranilate isomerase) and trpA (alpha-subunit of tryptophan synthase) suggest that these thermostable (beta alpha)8-barrel proteins may be stabilized by additional salt bridges, compared with the mesostable forms. Another notable feature is the predicted lack of the N-terminal helix alpha 0 in the alpha-subunit of tryptophan synthase.     

The (betaalpha)(8) glycosidases: sequence and structure analyses suggest distant evolutionary relationships.

There are currently at least nine distinct glycosidase sequence families which are all known to adopt a TIM barrel fold [Henrissat,B. and Davies,G. (1997) CURR: Opin. Struct. Biol., 7, 637-644]. To explore the relationships between these enzymes and their evolution, comprehensive sequence and structure comparisons were performed, generating four distinct clusters. The first cluster, S1, comprises the alpha-amylase related enzymes, all with the retention mechanism (axial-->axial). The second cluster, S2, included two functional subgroups, one composed of various kinds of glucosidases all with the retention mechanism (equatorial-->equatorial) (the so-called 4/7 superfamily), and the other subgroup including the beta-amylases with the inversion mechanism (axial--> equatorial). The third cluster, S3, with the retention mechanism (equatorial-->equatorial), could be subdivided, based on the catalytic residues and mechanisms, into two functional subgroups: the chitinase group, catalysed by two acidic residues on the C-termini of beta-4 and beta-6, and the hevamine group, using two acidic residues on the C-termini of beta-4 for catalysis. The fourth cluster, S4, is composed of chitobiase with the retention mechanism (equatorial--> equatorial). These clusters are compared with the sequence families derived by Henrissat and coworkers. PSI-BLAST profiles and multiple-alignments of tertiary structures suggest that S1 and S2 are distantly related, as are S3 and S4, which have N-acetylated substrates. This work highlights the difficulties of untangling distant evolutionary relationships in ubiquitous folds such as the TIM barrel.     

The structure/function relationship of a dual-substrate (betaalpha)8-isomerase

Two structures of phosphoribosyl isomerase A (PriA) from Streptomyces coelicolor, involved in both histidine and tryptophan biosynthesis, were solved at 1.8 Å resolution. A closed conformer was obtained, which represents the first complete structure of PriA, revealing hitherto unnoticed molecular interactions and the occurrence of conformational changes. Inspection of these conformers, including ligand-docking simulations, allowed identification of residues involved in substrate recognition, chemical catalysis and conformational changes. These predictions were validated by mutagenesis and functional analysis. Arg¹⁹ and Ser⁸¹ were shown to play critical roles within the carboxyl and amino phosphate-binding sites, respectively; the catalytic residues Asp¹¹ and Asp¹³⁰ are responsible for both activities; and Thr¹⁶⁶ and Asp¹⁷¹, which make an unusual contact, are likely to elicit the conformational changes needed for adopting the active site architectures. This represents the first report of the structure/function relationship of this (βα)₈-isomerase.     

Structural (betaalpha)8 TIM barrel model of 3-hydroxy-3-methylglutaryl-coenzyme A lyase

This study describes three novel homozygous missense mutations (S75R, S201Y, and D204N) in the 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) lyase gene, which caused 3-hydroxy-3-methylglutaric aciduria in patients from Germany, England, and Argentina. Expression studies in Escherichia coli show that S75R and S201Y substitutions completely abolished the HMG-CoA lyase activity, whereas D204N reduced catalytic efficiency to 6.6% of the wild type. We also propose a three-dimensional model for human HMG-CoA lyase containing a (betaalpha)8 (TIM) barrel structure. The model is supported by the similarity with analogous TIM barrel structures of functionally related proteins, by the localization of catalytic amino acids at the active site, and by the coincidence between the shape of the substrate (HMG-CoA) and the predicted inner cavity. The three novel mutations explain the lack of HMG-CoA lyase activity on the basis of the proposed structure: in S75R and S201Y because the new amino acid residues occlude the substrate cavity, and in D204N because the mutation alters the electrochemical environment of the active site. We also report the localization of all missense mutations reported to date and show that these mutations are located in the beta-sheets around the substrate cavity.     

Divergent evolution of chloroplast-type ferredoxins

The TOL plasmid pWW0 of Pseudomonas putida encodes a set of enzymes required for the oxidation of toluene to Krebs cycle intermediates. The structural genes for these enzymes are encoded in two operons which comprise the xylCMABN and xylXYZLTEGFJQKIH genes, respectively. The function of the xylT gene has not yet been identified. The nucleotide sequence of xylT was determined in this study and putative gene product was shown to contain a sequence characteristic for chloroplast-type ferredoxins. The nahT gene, the homologue of xylT, present on NAH plasmid NAH7 encoding naphthalene-degrading enzymes, was also sequenced. The sequence conservation between xylT and nahT strongly suggests that both gene products have some physiological function. Chloroplast-type ferredoxins have been discovered in photosynthetic organisms (plants, algae, cyanobacteria and Rhodobacter) and Halobacterium species. Furthermore, chloroplast-type ferredoxin-like sequences have been found in the electron-transfer components of some oxygenases. The sequences of XylT and NahT were compared with those of the previously identified chloroplast-type ferredoxins, in order to examine their evolutionary relationships.     

Evaluation of the sequence template method for protein structure prediction. Discrimination of the (beta/alpha)8-barrel fold.

A multiple alignment of five (beta/alpha)8-barrel enzymes has been derived from their structure. The eight beta-strands and eight alpha-helices of the (beta/alpha)8-barrel are correctly aligned and the equivalenced residues in these regions fulfil similar structural roles. Each beta-strand has a central core of usually four residues, two residues contribute side-chains to the barrel core and the other two residues are involved in beta-strand/alpha-helix contacts. However, the fold imposes no constraints on the volumes of the residues at either a local or global level: the volume of the beta-barrel core varies between 1088 A3 in glycolate oxidase and 1571 A3 in taka-amylase. Sequence motifs derived from the multiple alignment were scanned against a database of 124 protein sequences, including 17 (beta/alpha)8-barrel enzymes. The results were evaluated in terms of the discrimination of (beta/alpha)8-barrel sequences and the quality of the alignments obtained. One motif was able to identify the top 12% of high scoring sequences as forming (beta/alpha)8-barrels with 50% accuracy and the bottom 50% of sequences as not being (beta/alpha)8-barrel proteins with 100% accuracy. However, in most instances the alignments were poor. The reasons for this are discussed with reference to the (beta/alpha)8-barrel proteins and the sequence motif method in general.     

Divergent evolution of oxidosqualene cyclases in plants

Triterpenes are one of the largest classes of plant metabolites and have important functions. A diverse array of triterpenoid skeletons are synthesized via the isoprenoid pathway by enzymatic cyclization of 2,3-oxidosqualene. The genomes of the lower plants Chlamydomonas reinhardtii and moss (Physcomitrella patens) contain just one oxidosqualene cyclase (OSC) gene (for sterol biosynthesis), whereas the genomes of higher plants contain nine to 16 OSC genes. Here we carry out functional analysis of rice OSCs and rigorous phylogenetic analysis of 96 OSCs from higher plants, including Arabidopsis thaliana, Oryza sativa, Sorghum bicolor and Brachypodium distachyon. The functional analysis identified an amino acid sequence for isoarborinol synthase (OsIAS) (encoded by Os11g35710/OsOSC11) in rice. Our phylogenetic analysis suggests that expansion of OSC members in higher plants has occurred mainly through tandem duplication followed by positive selection and diversifying evolution, and consolidated the previous suggestion that dicot triterpene synthases have been derived from an ancestral lanosterol synthase instead of directly from their cycloartenol synthases. The phylogenetic trees are consistent with the reaction mechanisms of the protosteryl and dammarenyl cations which parent a wide variety of triterpene skeletal types, allowing us to predict the functions of the uncharacterized OSCs.