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.     

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.     

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.     

Divergent evolution of (betaalpha)8-barrel enzymes.

The (betaalpha)8-barrel is the most versatile and most frequently encountered fold among enzymes. It is an interesting question how the contemporary (betaalpha)8-barrels are evolutionarily related and by which mechanisms they evolved from more simple precursors. Comprehensive comparisons of amino acid sequences and three-dimensional structures suggest that a large fraction of the known (betaalpha)8-barrels have divergently evolved from a common ancestor. The mutational interconversion of enzymatic activities of several (betaalpha)8-barrels further supports their common evolutionary origin. Moreover, the high structural similarity between the N- and C-terminal (betaalpha)4 units of two (betaalpha)8-barrel enzymes from histidine biosynthesis indicates that the contemporary proteins evolved by tandem duplication and fusion of the gene of an ancestral 'half-barrel' precursor. In support of this hypothesis, recombinantly produced 'half-barrels' were shown to be folded, dimeric proteins.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Coevolving residues of (beta/alpha)(8)-barrel proteins play roles in stabilizing active site architecture and coordinating protein dynamics

Indole-3-glycerol phosphate synthase (IGPS) is a representative of (β/α)₈-barrel proteins—the most common enzyme fold in nature. To better understand how the constituent amino-acids work together to define the structure and to facilitate the function, we investigated the evolutionary and dynamical coupling of IGPS residues by combining statistical coupling analysis (SCA) and molecular dynamics (MD) simulations. The coevolving residues identified by the SCA were found to form a network which encloses the active site completely. The MD simulations showed that these coevolving residues are involved in the correlated and anti-correlated motions. The correlated residues are within van der Waals contact and appear to maintain the active site architecture; the anti-correlated residues are mainly distributed on opposite sides of the catalytic cavity and coordinate the motions likely required for the substrate entry and product release. Our findings might have broad implications for proteins with the highly conserved (βα)₈-barrel in assessing the roles of amino-acids that are moderately conserved and not directly involved in the active site of the (β/α)₈-barrel. The results of this study could also provide useful information for further exploring the specific residue motions for the catalysis and protein design based on the (β/α)₈-barrel scaffold.     

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.     

Specificity modulation of barley alpha-amylase through biased random mutagenesis involving a conserved tripeptide in beta --> alpha loop 7 of the catalytic (beta/alpha)(8)-barrel domain.

The relative specificity and bond cleavage pattern of barley alpha-amylase 1 (AMY1) were dramatically changed by mutation in F(286)VD that connected beta-strand 7 of the catalytic (beta/alpha)(8)-barrel to a succeeding 3(10)-helix. This conserved tripeptide of the otherwise variable beta --> alpha segment 7 lacked direct ligand contact, but the nearby residues His290 and Asp291 participated in transition-state stabilization and catalysis. On the basis of sequences of glycoside hydrolase family 13, a biased random mutagenesis protocol was designed which encoded 174 putative F(286)VD variants of C95A-AMY1, chosen as the parent enzyme to avoid inactivating glutathionylation by the yeast host. The FVG, FGG, YVD, LLD, and FLE mutants showed 12-380 and 1.8-33% catalytic efficiency (k(cat)/K(m)) toward 2-chloro-4-nitrophenyl beta-D-maltoheptaoside and amylose DP17, respectively, and 0.5-50% activity for insoluble starch compared to that of C95A-AMY1. K(m) and k(cat) were decreased 2-9- and 1.3-83-fold, respectively, for the soluble substrates. The starch:oligosaccharide and amylose:oligosaccharide specificity ratios were 13-172 and 2.4-14 for mutants and 520 and 27 for C95A-AMY1, respectively. The FVG mutant released 4-nitrophenyl alpha-D-maltotrioside (PNPG(3)) from PNPG(5), whereas C95A-AMY1 produced PNPG and PNPG(2). The mutation thus favored interaction with the substrate aglycon part, while products from PNPG(6) reflected the fact that the mutation restored binding at subsite -6 which was lost in C95A-AMY1. The outcome of this combined irrational and rational protein engineering approach was evaluated considering structural accommodation of mutant side chains. FVG and FGG, present in the most active variants, represented novel sequences. This emphasized the worth of random mutagenesis and launched flexibility as a goal for beta --> alpha loop 7 engineering in family 13.     

D-Ribulose 5-phosphate 3-epimerase: functional and structural relationships to members of the ribulose-phosphate binding (beta/alpha)8-barrel superfamily

The "ribulose phosphate binding" superfamily defined by the Structural Classification of Proteins (SCOP) database is considered the result of divergent evolution from a common (beta/alpha)(8)-barrel ancestor. The superfamily includes d-ribulose 5-phosphate 3-epimerase (RPE), orotidine 5'-monophosphate decarboxylase (OMPDC), and 3-keto-l-gulonate 6-phosphate decarboxylase (KGPDC), members of the OMPDC suprafamily, as well as enzymes involved in histidine and tryptophan biosynthesis that utilize phosphorylated metabolites as substrates. We now report studies of the functional and structural relationships of RPE to the members of the superfamily. As suggested by the results of crystallographic studies of the RPEs from rice [Jelakovic, S., Kopriva, S., Suss, K. H., and Schulz, G. E. (2003) J. Mol. Biol. 326, 127-35] and Plasmodium falciparum [Caruthers, J., Bosch, J., Bucker, F., Van Voorhis, W., Myler, P., Worthey, E., Mehlin, C., Boni, E., De Titta, G., Luft, J., Kalyuzhniy, O., Anderson, L., Zucker, F., Soltis, M., and Hol, W. G. J. (2006) Proteins 62, 338-42], the RPE from Streptococcus pyogenes is activated by Zn(2+) which binds with a stoichiometry of one ion per polypeptide. Although wild type RPE has a high affinity for Zn(2+) and inactive apoenzyme cannot be prepared, the affinity for Zn(2+) is decreased by alanine substitutions for the two histidine residues that coordinate the Zn(2+) ion (H34A and H67A); these mutant proteins can be prepared in an inactive, metal-free form and activated by exogenous Zn(2+). The crystal structure of the RPE was solved at 1.8 A resolution in the presence of d-xylitol 5-phosphate, an inert analogue of the d-xylulose 5-phosphate substrate. This structure suggests that the 2,3-enediolate intermediate in the 1,1-proton transfer reaction is stabilized by bidentate coordination to the Zn(2+) that also is liganded to His 34, Asp 36, His 67, and Asp 176; the carboxylate groups of the Asp residues are positioned also to function as the acid/base catalysts. Although the conformation of the bound analogue resembles those of ligands bound in the active sites of OMPDC and KGPDC, the identities of the active site residues that coordinate the essential Zn(2+) and participate as acid/base catalysts are not conserved. We conclude that only the phosphate binding motif located at the ends of the seventh and eighth beta-strands of the (beta/alpha)(8)-barrel is functionally conserved among RPE, OMPDC, and KGPDC, consistent with the hypothesis that the members of the "ribulose phosphate binding" (beta/alpha)(8)-barrel "superfamily" as defined by SCOP have not evolved by evolutionary processes involving the intact (beta/alpha)(8)-barrel. Instead, this "superfamily" may result from assembly from smaller modules, including the conserved phosphate binding motif associated with the C-terminal (beta/alpha)(2)-quarter barrel.     

Manipulation of the active site loops of D-hydantoinase, a (beta/alpha)8-barrel protein, for modulation of the substrate specificity

We previously proposed that the stereochemistry gate loops (SGLs) constituting the substrate binding pocket of D-hydantoinase, a (beta/alpha)(8)-barrel enzyme, might be major structural determinants of the substrate specificity [Cheon, Y. H., et al. (2002) Biochemistry 41, 9410-9417]. To construct a mutant D-hydantoinase with favorable substrate specificity for the synthesis of commercially important non-natural amino acids, the SGL loops of the enzyme were rationally manipulated on the basis of the structural analysis and sequence alignment of three hydantoinases with distinct substrate specificities. In the SGLs of D-hydantoinase from Bacillus stearothermophilus SD1, mutations of hydrophobic and bulky residues Met 63, Leu 65, Phe 152, and Phe 159, which interact with the exocyclic substituent of the substrate, induced remarkable changes in the substrate specificities. In particular, the substrate specificity of mutant F159A toward aromatic substrate hydroxyphenylhydantoin (HPH) was enhanced by approximately 200-fold compared with that of the wild-type enzyme. Saturation mutagenesis at position 159 revealed that k(cat) for aromatic substrates increased gradually as the size of the amino acid side chain decreased, and this seems to be due to reduced steric hindrance between the bulky exocyclic group of the substrate and the amino acid side chains. When site-directed random mutagenesis of residues 63 and 65 was conducted with the wild type and mutant F159A, the selected enzymes (M63F/L65V and L65F/F159A) exhibited approximately 10-fold higher k(cat) values for HPH than the wild-type counterpart, which is likely to result from reorganization of the active site for efficient turnover. These results indicate that the amino acid residues of SGLs forming the substrate binding pocket are critical for the substrate specificity of D-hydantoinase, and the results also imply that substrate specificities of cyclic amidohydrolase family enzymes can be modulated by rational design of these SGLs.     

(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 equilibrium unfolding pathway of a (beta/alpha)8 barrel

The (beta/alpha)(8) barrel is the most commonly occurring fold among enzymes. A key step towards rationally engineering (beta/alpha)(8) barrel proteins is to understand their underlying structural organization and folding energetics. Using misincorporation proton-alkyl exchange (MPAX), a new tool for solution structural studies of large proteins, we have performed a native-state exchange analysis of the prototypical (beta/alpha)(8) barrel triosephosphate isomerase. Three cooperatively unfolding subdomains within the structure are identified, as well as two partially unfolded forms of the protein. The C-terminal domain coincides with domains reported to exist in four other (beta/alpha)(8) barrels, but the two N-terminal domains have not been observed previously. These partially unfolded forms may represent sequential intermediates on the folding pathway of triosephosphate isomerase. The methods reported here should be applicable to a variety of other biological problems involving protein conformational changes.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

The evolution of alpha/beta barrel enzymes

Of the enzymes whose structure is known, roughly one out of ten has an eight-stranded alpha/beta barrel domain. Structural and chemical evidence suggests that all of these domains may have diverged from a common ancestor.     

Evolution of alpha-amylases: architectural features and key residues in the stabilization of the (beta/alpha)(8) scaffold

We provide a comprehensive analysis of the current enzymes with alpha-amylase activity (AAMYs) that belong to family 13 glycoside hydrolase (GH-13; 144 Archaea, Bacteria, and Eukaryota sequences from 87 different species). This study aims to further knowledge of the evolutionary molecular relationships among the sequences of their A and B domains with special emphasis on the correlation between what is observed in the structures and protein evolution. Multialignments for the A domain distinguish two clusters for sequences from Archaea organisms, eight for sequences from Bacteria organisms, and three for sequences from Eukaryota organisms. The clusters for Bacteria do not follow any strict taxonomic pathway; in fact, they are rather scattered. When we compared the A domains of sequences belonging to different kingdoms, we found that various pairs of clusters were significantly similar. Using either sequence similarity with crystallized structures or secondary-structure prediction methods, we identified in all AAMYs the eight putative beta-strands that constitute the beta-sheet in the TIM barrel of the A domain and studied the packing in its interior. We also discovered a "hidden homology" in the TIM barrel, an invariant Gly located upstream in the sequence before the conserved Asp in beta-strand 3. This Gly precedes an alpha-helix and is actively involved in capping its N-terminal end with a capping box. In all cases, a Schellman motif caps the C-terminal end of this helix.     

Folding of beta/alpha-unit scrambled forms of S. cerevisiae triosephosphate isomerase: Evidence for autonomy of substructure formation and plasticity of hydrophobic and hydrogen bonding interactions in core of (beta/alpha)8-barrel

We show that residues at the interfaces of protein-protein complexes have higher side-chain energy than other surface residues. Eight different sets of protein complexes were analyzed. For each protein pair, the complex structure was used to identify the interface residues in the unbound monomer structures. Side-chain energy was calculated for each surface residue in the unbound monomer using our previously developed scoring function.1 The mean energy was calculated for the interface residues and the other surface residues. In 15 of the 16 monomers, the mean energy of the interface residues was higher than that of other surface residues. By decomposing the scoring function, we found that the energy term of the buried surface area of non-hydrogen-bonded hydrophilic atoms is the most important factor contributing to the high energy of the interface regions. In spite of lacking hydrophilic residues, the interface regions were found to be rich in buried non-hydrogen-bonded hydrophilic atoms. Although the calculation results could be affected by the inaccuracy of the scoring function, patch analysis of side-chain energy on the surface of an isolated protein may be helpful in identifying the possible protein-protein interface. A patch was defined as 20 residues surrounding the central residue on the protein surface, and patch energy was calculated as the mean value of the side-chain energy of all residues in the patch. In 12 of the studied monomers, the patch with the highest energy overlaps with the observed interface. The results are more remarkable when only three residues with the highest energy in a patch are averaged to derive the patch energy. All three highest-energy residues of the top energy patch belong to interfacial residues in four of the eight small protomers. We also found that the residue with the highest energy score on the surface of a small protomer is very possibly the key interaction residue.     

Exploring the environmental preference of weak interactions in (alpha/beta)8 barrel proteins

The environmental preference for the occurrence of noncanonical hydrogen bonding and cation-pi interactions, in a data set containing 71 nonredundant (alpha/beta)(8) barrel proteins, with respect to amino acid type, secondary structure, solvent accessibility, and stabilizing residues has been performed. Our analysis reveals some important findings, which include (a) higher contribution of weak interactions mediated by main-chain atoms irrespective of the amino acids involved; (b) domination of the aromatic amino acids among interactions involving side-chain atoms; (c) involvement of strands as the principal secondary structural unit, accommodating cross strand ion pair interaction and clustering of aromatic amino acid residues; (d) significant contribution to weak interactions occur in the solvent exposed areas of the protein; (e) majority of the interactions involve long-range contacts; (f) the preference of Arg is higher than Lys to form cation-pi interaction; and (g) probability of theoretically predicted stabilizing amino acid residues involved in weak interaction is higher for polar amino acids such as Trp, Glu, and Gln. On the whole, the present study reveals that the weak interactions contribute to the global stability of (alpha/beta)(8) TIM-barrel proteins in an environment-specific manner, which can possibly be exploited for protein engineering applications.     

Carbohydrate-carbohydrate binding of ganglioside to integrin alpha(5) modulates alpha(5)beta(1) function

Gangliosides GT1b and GD3, components of keratinocyte membranes, inhibit keratinocyte adhesion to fibronectin. Although ganglioside sialylation is known to be important, the mechanism of inhibition is unknown. Using purified insect recombinant alpha(5) and beta(1) proteins and alpha(5)beta(1) integrin from lysed keratinocyte-derived SCC12 cells, we have shown that GT1b and GD3 inhibit the binding of alpha(5)beta(1) to fibronectin. Co-immunoprecipitation of GT1b and alpha(5)beta(1) from SCC12 cells and direct binding of GT1b and GD3 to affinity-purified alpha(5)beta(1) from SCC12 cells and insect recombinant alpha(5)beta(1), particularly the alpha(5) subunit, further suggest interaction between ganglioside and alpha(5)beta(1). The carbohydrate moieties of integrin appear to be critical since gangliosides are unable to bind deglycosylated forms of alpha(5)beta(1) from SCC12 and insect cells or poorly glycosylated recombinant alpha(5)beta(1) from Escherichia coli cells. The GT1b-alpha(5)beta(1) interaction is inhibited by concanavalin A, suggesting that GT1b binds to mannose structures in alpha(5)beta(1). The preferential binding of GT1b to high mannose rather than reduced mannose ovalbumin further implicates the binding of GT1b to mannose structures. These data provide evidence that highly sialylated gangliosides regulate alpha(5)beta(1)-mediated adhesion of epithelial cells to fibronectin through carbohydrate-carbohydrate interactions between GT1b and the alpha(5) subunit of alpha(5)beta(1) integrin.     

A new arrangement of (beta/alpha)8 barrels in the synthase subunit of PLP synthase

Pyridoxal 5'-phosphate (PLP, vitamin B6), a cofactor in many enzymatic reactions, has two distinct biosynthetic routes, which do not coexist in any organism. Two proteins, known as PdxS and PdxT, together form a PLP synthase in plants, fungi, archaea, and some eubacteria. PLP synthase is a heteromeric glutamine amidotransferase in which PdxT produces ammonia from glutamine and PdxS combines ammonia with five- and three-carbon phosphosugars to form PLP. In the 2.2-A crystal structure, PdxS is a cylindrical dodecamer of subunits having the classic (beta/alpha)8 barrel fold. PdxS subunits form two hexameric rings with the active sites positioned on the inside. The hexamer and dodecamer forms coexist in solution. A novel phosphate-binding site is suggested by bound sulfate. The sulfate and another bound molecule, methyl pentanediol, were used to model the substrate ribulose 5-phosphate, and to propose catalytic roles for residues in the active site. The distribution of conserved surfaces in the PdxS dodecamer was used to predict a docking site for the glutaminase partner, PdxT.