Prediction of secondary structure by evolutionary comparison: application to the alpha subunit of tryptophan synthase

The amino acid sequences of the a subunits of tryptophan synthase from ten different microorganisms were aligned by standard procedures. The alpha helices, beta strands and turns of each sequence were predicted separately by two standard prediction algorithms and averaged at homologous sequence positions. Additional evidence for conserved secondary structure was derived from profiles of average hydropathy and chain flexibility values, leading to a joint prediction. There is good agreement between (1) predicted beta strands, maximal hydropathy and minimal flexibility, and (2) predicted loops, great chain flexibility, and protein segments that accept insertions of various lengths in individual sequences. The a subunit is predicted to have eight repeated beta-loop-alpha-loop motifs with an extra N-terminal alpha helix and an intercalated segment of highly conserved residues. This pattern suggests that the territory structure of the a subunit is an eightfold alpha/beta barrel. The distribution of conserved amino acid residues and published data on limited proteolysis, chemical modification, and mutagenesis are consistent with the alpha/beta barrel structure. Both the active site of the a subunit and the combining site for the beta 2 subunit are at the end of the barrel formed by the carboxyl-termini of the beta strands.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

PFIT and PFRIT: bioinformatic algorithms for detecting glycosidase function from structure and sequence

The identification of the enzymes involved in the metabolism of simple and complex carbohydrates presents one bioinformatic challenge in the post-genomic era. Here, we present the PFIT and PFRIT algorithms for identifying those proteins adopting the alpha/beta barrel fold that function as glycosidases. These algorithms are based on the observation that proteins adopting the alpha/beta barrel fold share positions in their tertiary structures having equivalent sets of atomic interactions. These are conserved tertiary interaction positions, which have been implicated in both structure and function. Glycosidases adopting the alpha/beta barrel fold share more conserved tertiary interactions than alpha/beta barrel proteins having other functions. The enrichment pattern of conserved tertiary interactions in the glycosidases is the information that PFIT and PFRIT use to predict whether any given alpha/beta barrel will function as a glycosidase or not. Using as a test set a database of 19 glycosidase and 45 nonglycosidase alpha/beta barrel proteins with low sequence similarity, PFIT and PFRIT can correctly predict glycosidase function for 84% of the proteins known to function as glycosidases. PFIT and PFRIT incorrectly predict glycosidase function for 25% of the nonglycosidases. The program PSI-BLAST can also correctly identify 84% of the 19 glycosidases, however, it incorrectly predicts glycosidase function for 50% of the nonglycosidases (twofold greater than PFIT and PFRIT). Overall, we demonstrate that the structure-based PFIT and PFRIT algorithms are both more selective and sensitive for predicting glycosidase function than the sequence-based PSI-BLAST algorithm.     

Role of hydrophobic clusters and long-range contact networks in the folding of (alpha/beta)8 barrel proteins

Analysis on the three dimensional structures of (alpha/beta)(8) barrel proteins provides ample light to understand the factors that are responsible for directing and maintaining their common fold. In this work, the hydrophobically enriched clusters are identified in 92% of the considered (alpha/beta)(8) barrel proteins. The residue segments with hydrophobic clusters have high thermal stability. Further, these clusters are formed and stabilized through long-range interactions. Specifically, a network of long-range contacts connects adjacent beta-strands of the (alpha/beta)(8) barrel domain and the hydrophobic clusters. The implications of hydrophobic clusters and long-range networks in providing a feasible common mechanism for the folding of (alpha/beta)(8) barrel proteins are proposed.     

Common structural features of the luxF protein and the subunits of bacterial luciferase: evidence for a (beta alpha)8 fold in luciferase.

The amino acid sequence identity and potential structural similarity between the subunits of bacterial luciferase and the recently determined structure of the luxF molecule are examined. The unique beta/alpha barrel fold found in luxF appears to be conserved in part in the luciferase subunits. From secondary structural predictions of both luciferase subunits, and from structural comparisons between the protein product of the luxF gene, NFP, and glycolate oxidase, we propose that it is feasible for both luciferase subunits to adopt a (beta alpha)8 barrel fold with at least 2 excursions from the (beta alpha)8 topology. Amino acids conserved between NFP and the luciferase subunits cluster together in 3 distinct "pockets" of NFP, which are located at hydrophobic interfaces between the beta-strands and alpha-helices. Several tight turns joining the C-termini of beta-strands and the N-termini of alpha-helices are found as key components of these conserved regions. Helix start and end points are easily demarcated in the luciferase subunit protein sequences; the N-cap residues are the most strongly conserved structural features. A partial model of the luciferase beta subunit from Photobacterium leiognathi has been built based on our crystallographically determined structure of luxF at 1.6 A resolution.     

Structural principles of alpha/beta barrel proteins: the packing of the interior of the sheet

Alpha/beta barrel structures very similar to that first observed in triose phosphate isomerase are now known to occur in 14 enzymes. To understand the origin of this fold, we analyzed in three of these proteins the geometry of the eight-stranded beta-sheets and the packing of the residues at the center of the barrel. The packing in this region is seen in its simplest form in glycolate oxidase. It consists of 12 residues arranged in three layers. Each layer contains four side chains. The packing of RubisCO and TIM can be understood in terms of distortions of this simple pattern, caused by residues with small side chains at some of the positions inside the barrel. Two classes of packing are found. In one class, to which RubisCO and TIM belong, the central layer is formed by a residue from the first, third, fifth, and seventh strands; the upper and lower layers are formed by residues from the second, fourth, sixth, and eighth strands. In the second class, to which GAO belongs, this is reversed: it is side chains from the even-numbered strands that form the central layer, and side chains from the odd-numbered strands that form the outer layers. Our results suggest that not all proteins with this fold are related by evolution, but that they represent a common favorable solution to the structural problems involved in the creation of a closed beta barrel.     

Crystal structure of imidazole glycerol phosphate synthase: a tunnel through a (beta/alpha)8 barrel joins two active sites

BACKGROUND: Imidazole glycerol phosphate synthase catalyzes a two-step reaction of histidine biosynthesis at the bifurcation point with the purine de novo pathway. The enzyme is a new example of intermediate channeling by glutamine amidotransferases in which ammonia generated by hydrolysis of glutamine is channeled to a second active site where it acts as a nucleophile. In this case, ammonia reacts in a cyclase domain to produce imidazole glycerol phosphate and an intermediate of purine biosynthesis. The enzyme is also a potential target for drug and herbicide development since the histidine pathway does not occur in mammals. RESULTS: The 2.1 A crystal structure of imidazole glycerol phosphate synthase from yeast reveals extensive interaction of the glutaminase and cyclase catalytic domains. At the domain interface, the glutaminase active site points into the bottom of the (beta/alpha)(8) barrel of the cyclase domain. An ammonia tunnel through the (beta/alpha)(8) barrel connects the glutaminase docking site at the bottom to the cyclase active site at the top. A conserved "gate" of four charged residues controls access to the tunnel. CONCLUSIONS: This is the first structure in which all the components of the ubiquitous (beta/alpha)(8) barrel fold, top, bottom, and interior, take part in enzymatic function. Intimate contacts between the barrel domain and the glutaminase active site appear to be poised for crosstalk between catalytic centers in response to substrate binding at the cyclase active site. The structure provides a number of potential sites for inhibitor development in the active sites and in a conserved interdomain cavity.     

Three-dimensional structure of rat acid phosphatase.

The crystal structure of recombinant rat prostatic acid phosphatase was determined to 3 A resolution with protein crystallographic methods. The enzyme subunit is built up of two domains, an alpha/beta domain consisting of a seven-stranded mixed beta-sheet with helices on both sides of the sheet and a smaller alpha domain. Two disulfide bridges between residues 129-340 and 315-319 were found. Electron density at two of the glycosylation sites for parts of the carbohydrate moieties was observed. The dimer of acid phosphatase is formed through two-fold interactions of edge strand 3 from one subunit with strand 3 from the second subunit, thus extending the beta-sheet from seven to 14 strands. Other subunit-subunit interactions involve conserved residues from loops between helices and beta-strands. The fold of the alpha/beta domain is similar to the fold observed in phosphoglycerate mutase. The active site is at the carboxy end of the parallel strands of the alpha/beta domain. There is a strong residual electron density at the phosphate binding site which probably represents a bound chloride ion. Biochemical properties and results from site-directed mutagenesis experiments of acid phosphatase are correlated to the three-dimensional structure.     

Structural predictions of AgfA, the insoluble fimbrial subunit of Salmonella thin aggregative fimbriae.

The unusually stable and multifunctional, thin aggregative fimbriae common to all Salmonella spp. are principally polymers of the fimbrin subunit, AgfA. AgfA of Salmonella enteritidis consists of two domains: a protease-sensitive, 22 amino acid residue N-terminal region and a protease-resistant, 109 residue C-terminal core. The unusual amino acid sequence of the AgfA core region comprises two-, five- and tenfold internal sequence homology patterns reflected in five conserved, 18-residue tandem repeats. These repeats have the consensus sequence, Sx5QxGx2NxAx3Q and are linked together by four or five residues, (x)xAx2. The predicted secondary structure for this unusual arrangement of tandem repeats in AgfA indicates mainly extended conformation with the beta strands linked by four to six residues. Candidate proteins of known structure with motifs of alternating beta strands and short loops were selected from folds described in SCOP as a source of coordinates for AgfA model construction. Three all-beta class motifs selected from the Serratia marcescens metalloprotease, myelin P2 protein or vitelline membrane outer protein I were used for initial AgfA homology build-up procedures ultimately resulting in three structural models; beta barrel, beta prism and parallel beta helix. The beta barrel model is a compact, albeit irregular structure, with the beta strands arranged in two antiparallel beta sheet faces. The beta prism model does not reflect the 5 or 10-fold symmetry of the AgfA primary sequence. However, the favored, parallel beta helix model is a compact coil of ten helically arranged beta strands forming two parallel beta sheet faces. This arrangement predicts a regular, potentially stable, C-terminal core region consistent with the observed tandem repeat sequences, protease-resistance and strong tendency of this fimbrin to oligomerize and aggregate. Positional conservation of amino acid residues in AgfA and the Escherichia coli AgfA homologue, CsgA, provides strong support for this model. The parallel beta helix model of AgfA offers an interesting solution to a multifunctional fimbrin molecular surface having solvent exposed areas, regions for major and minor subunit interactions as well as fiber-fiber interactions common to many bacterial fimbriae.     

A six-stranded double-psi beta barrel is shared by several protein superfamilies.

BACKGROUND: Six-stranded beta barrels with a pseudo-twofold axis are found in several proteins. One group comprises a Greek-key structure with all strands antiparallel; an example is the N-terminal domain of ferredoxin reductase. Others involve parallel strands forming two psi structures (the double-psi beta barrel). A recently discovered example of the latter class is aspartate-alpha-decarboxylase (ADC) from Escherichia coli, a pyruvoyl-dependent tetrameric enzyme involved in the synthesis of pantothenate. RESULTS: Visual inspection and automated database searches identified the six-stranded double-psi beta barrel in ADC, Rhodobacter sphaeroides dimethylsulfoxide (DMSO) reductase, E. coli formate dehydrogenase H (FDHH), the plant defense protein barwin, Humicola insolens endoglucanase V (EGV) and, with a circular permutation, in the aspartic proteinases. Structure-based sequence alignments revealed several interactions including hydrophobic contacts or sidechain-mainchain hydrogen bonds that position the middle beta strand under a psi loop, which may significantly contribute to stabilizing the fold. The identification of key interactions allowed the filtering of weak sequence similarities to some of these proteins, which had been detected by sequence database searches. This led to the prediction of the double-psi beta-barrel domain in several families of proteins in eukaryotes and archaea. CONCLUSIONS: The structure comparison and clustering study of double-psi beta barrels suggests that there could be a common homodimeric ancestor to ADC, FDHH and DMSO reductase, and also to barwin and EGV. There are other protein families with unknown structure that are likely to adopt the same fold. In the known structures, the protein active sites cluster around the psi loop, indicating that its rigidity, protrusion and free mainchain functional groups may be well suited to providing a framework for catalysis.     

Protein-segment universe exhibiting transitions at intermediate segment length in conformational subspaces

Background

Lipocalins are widely distributed in nature and are found in bacteria, plants, arthropoda and vertebra. In hematophagous arthropods, they are implicated in the successful accomplishment of the blood meal, interfering with platelet aggregation, blood coagulation and inflammation and in the transmission of disease parasites such as Trypanosoma cruzi and Borrelia burgdorferi. The pairwise sequence identity is low among this family, often below 30%, despite a well conserved tertiary structure. Under the 30% identity threshold, alignment methods do not correctly assign and align proteins. The only safe way to assign a sequence to that family is by experimental determination. However, these procedures are long and costly and cannot always be applied. A way to circumvent the experimental approach is sequence and structure analyze. To further help in that task, the residues implicated in the stabilisation of the lipocalin fold were determined. This was done by analyzing the conserved interactions for ten lipocalins having a maximum pairwise identity of 28% and various functions.

Results

It was determined that two hydrophobic clusters of residues are conserved by analysing the ten lipocalin structures and sequences. One cluster is internal to the barrel, involving all strands and the 3₁₀ helix. The other is external, involving four strands and the helix lying parallel to the barrel surface. These clusters are also present in RaHBP2, a unusual "outlier" lipocalin from tick Rhipicephalus appendiculatus. This information was used to assess assignment of LIR2 a protein from Ixodes ricinus and to build a 3D model that helps to predict function. FTIR data support the lipocalin fold for this protein.

Conclusion

By sequence and structural analyzes, two conserved clusters of hydrophobic residues in interactions have been identified in lipocalins. Since the residues implicated are not conserved for function, they should provide the minimal subset necessary to confer the lipocalin fold. This information has been used to assign LIR2 to lipocalins and to investigate its structure/function relationship. This study could be applied to other protein families with low pairwise similarity, such as the structurally related fatty acid binding proteins or avidins.     

The central cavity from the (alpha/alpha)6 barrel structure of Anabaena sp. CH1 N-acetyl-D-glucosamine 2-epimerase contains two key histidine residues for reversible conversion

N-acetyl-D-glucosamine 2-epimerase (GlcNAc 2-epimerase) catalyzes the reversible epimerization between N-acetyl-D-glucosamine (GlcNAc) and N-acetyl-D-mannosamine (ManNAc). We report here the 2.0 A resolution crystal structure of the GlcNAc 2-epimerase from Anabaena sp. CH1. The structure demonstrates an (alpha/alpha)(6) barrel fold, which shows structural homology with porcine GlcNAc 2-epimerase as well as a number of glycoside hydrolase enzymes and other sugar-metabolizing enzymes. One side of the barrel structure consists of short loops involved in dimer interactions. The other side of the barrel structure is comprised of long loops containing six short beta-sheets, which enclose a putative central active-site pocket. Site-directed mutagenesis of conserved residues near the N-terminal region of the inner alpha helices shows that R57, H239, E308, and H372 are strictly required for activity. E242 and R375 are also essential in catalysis. Based on the structure and kinetic analysis, H239 and H372 may serve as the key active site acid/base catalysts. These results suggest that the (alpha/alpha)(6) barrel represents a steady fold for presenting active-site residues in a cleft at the N-terminal ends of the inner alpha helices, thus forming a fine-tuned catalytic site in GlcNAc 2-epimerase.     

Stabilizing interactions in the dimer interface of alpha-subunit in Escherichia coli RNA polymerase: a graph spectral and point mutation study

The formation of alpha(2) dimer in Escherichia coli core RNA polymerase (RNAP) is thought to be the first step toward the assembly of the functional enzyme. A large number of evidences indicate that the alpha-subunit dimerizes through its N-terminal domain (NTD). The crystal structures of the alpha-subunit NTD and that of a homologous Thermus aquaticus core RNAP are known. To identify the stabilizing interactions in the dimer interface of the alpha-NTD of E. coli RNAP, we identified side-chain clusters by using the crystal structure coordinates of E. coli alpha-NTD. A graph spectral algorithm was used to identify side-chain clusters. This algorithm considers the global nonbonded side-chain interactions of the residues for the clustering procedure and is unique in identifying residues that make the largest number of interactions among the residues that form clusters in a very quantitative way. By using this algorithm, a nine-residue cluster consisting of polar and hydrophobic residues was identified in the subunit interface adjacent to the hydrophobic core. The residues forming the cluster are relatively rigid regions of the interface, as measured by the thermal factors of the residues. Most of the cluster residues in the E. coli enzyme were topologically and sequentially conserved in the T. aquaticus RNAP crystal structure. Residues 35F and 46I were predicted to be important in the stability of the alpha-dimer interface, with 35F forming the center of the cluster. The predictions were tested by isolating single-point mutants alpha-F35A and alpha-I46S on the dimer interface, which were found to disrupt dimerization. Thus, the identified cluster at the edge of the dimer interface seems to be a vital component in stabilizing the alpha-NTD.     

Crystal structure of amylomaltase from thermus aquaticus, a glycosyltransferase catalysing the production of large cyclic glucans

Amylomaltase is involved in the metabolism of starch, one of the most important polysaccharides in nature. A unique feature of amylomaltase is its ability to catalyze the formation of cyclic amylose. In contrast to the well studied cyclodextrin glucanotransferases (CGTases), which synthesize cycloamylose with a ring size (degree of polymerization or DP) of 6-8, the amylomaltase from Thermus aquaticus produces cycloamyloses with a DP of 22 and higher. The crystal structure of amylomaltase from Thermus aquaticus was determined to 2.0 A resolution. It is a member of the alpha-amylase superfamily of enzymes, whose core structure consists of a (beta, alpha)(8) barrel. In amylomaltase, the 8-fold symmetry of this barrel is disrupted by several insertions between the barrel strands. The largest insertions are between the third and fifth barrel strands, where two insertions form subdomain B1, as well as between the second and third barrel strands, forming the alpha-helical subdomain B2. Whereas part of subdomain B1 is also present in other enzyme structures of the alpha-amylase superfamily, subdomain B2 is unique to amylomaltase. Remarkably, the C-terminal domain C, which is present in all related enzymes of the alpha-amylase family, is missing in amylomaltase. Amylomaltase shows a similar arrangement of the catalytic side-chains (two Asp residues and one Glu residue) as in previously characterized members of the alpha-amylase superfamily, indicating similar mechanisms of the glycosyl transfer reaction. In amylomaltase, a conserved loop of around eight amino acid residues is partially shielding the active center. This loop, which is well conserved among other amylomaltases, may sterically hinder the formation of small cyclic products.     

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.     

Solution structure of the SEA domain from the murine homologue of ovarian cancer antigen CA125 (MUC16)

Human CA125, encoded by the MUC16 gene, is an ovarian cancer antigen widely used for a serum assay. Its extracellular region consists of tandem repeats of SEA domains. In this study we determined the three-dimensional structure of the SEA domain from the murine MUC16 homologue using multidimensional NMR spectroscopy. The domain forms a unique alpha/beta sandwich fold composed of two alpha helices and four antiparallel beta strands and has a characteristic turn named the TY-turn between alpha1 and alpha2. The internal mobility of the main chain is low throughout the domain. The residues that form the hydrophobic core and the TY-turn are fully conserved in all SEA domain sequences, indicating that the fold is common in the family. Interestingly, no other residues are conserved throughout the family. Thus, the sequence alignment of the SEA domain family was refined on the basis of the three-dimensional structure, which allowed us to classify the SEA domains into several subfamilies. The residues on the surface differ between these subfamilies, suggesting that each subfamily has a different function. In the MUC16 SEA domains, the conserved surface residues, Asn-10, Thr-12, Arg-63, Asp-75, Asp-112, Ser-115, and Phe-117, are clustered on the beta sheet surface, which may be functionally important. The putative epitope (residues 58-77) for anti-MUC16 antibodies is located around the beta2 and beta3 strands. On the other hand the tissue tumor marker MUC1 has a SEA domain belonging to another subfamily, and its GSVVV motif for proteolytic cleavage is located in the short loop connecting beta2 and beta3.     

Structure of the nuclear factor ALY: insights into post-transcriptional regulatory and mRNA nuclear export processes

ALY is a ubiquitously expressed nuclear protein which interacts with proteins such as TAP that are involved in export of mRNA from the nucleus to the cytoplasm, as well as with proteins that bind the T cell receptor alpha gene enhancer. ALY has also been shown to bind mRNA and to co-localize in the nucleus with components of a multiprotein postsplicing complex that is deposited 20-24 nucleotides upstream of exon-exon junctions. ALY has a conserved RNA binding domain (RBD) flanked by Gly-Arg rich N-terminal and C-terminal sequences. We determined the solution structure of the RBD homology region in ALY by nuclear magnetic resonance methods. The RBD motif in ALY has a characteristic beta(1)alpha(1)beta(2)-beta(3)alpha(2)beta(4) fold, consisting of a beta sheet composed of four antiparallel beta strands and two alpha helices that pack on one face of the beta sheet. As in other RBD structures, the beta sheet has an exposed face with hydrophobic and charged residues that could modulate interactions with other molecules. The loop that connects beta strands 2 and 3 is in intermediate motion in the NMR time scale, which is also characteristic of other RBDs. This loop presents side chains close to the exposed surface of the beta sheet and is a primary candidate site for intermolecular interactions. The structure of the conserved RBD of ALY provides insight into the nature of interactions involving this multifunctional protein.     

Structure of L-xylulose-5-Phosphate 3-epimerase (UlaE) from the anaerobic L-ascorbate utilization pathway of Escherichia coli: identification of a novel phosphate binding motif within a TIM barrel fold

Three catabolic enzymes, UlaD, UlaE, and UlaF, are involved in a pathway leading to fermentation of l-ascorbate under anaerobic conditions. UlaD catalyzes a beta-keto acid decarboxylation reaction to produce L-xylulose-5-phosphate, which undergoes successive epimerization reactions with UlaE (L-xylulose-5-phosphate 3-epimerase) and UlaF (L-ribulose-5-phosphate 4-epimerase), yielding D-xylulose-5-phosphate, an intermediate in the pentose phosphate pathway. We describe here crystallographic studies of UlaE from Escherichia coli O157:H7 that complete the structural characterization of this pathway. UlaE has a triosephosphate isomerase (TIM) barrel fold and forms dimers. The active site is located at the C-terminal ends of the parallel beta-strands. The enzyme binds Zn(2+), which is coordinated by Glu155, Asp185, His211, and Glu251. We identified a phosphate-binding site formed by residues from the beta1/alpha1 loop and alpha3' helix in the N-terminal region. This site differs from the well-characterized phosphate-binding motif found in several TIM barrel superfamilies that is located at strands beta7 and beta8. The intrinsic flexibility of the active site region is reflected by two different conformations of loops forming part of the substrate-binding site. Based on computational docking of the L-xylulose 5-phosphate substrate to UlaE and structural similarities of the active site of this enzyme to the active sites of other epimerases, a metal-dependent epimerization mechanism for UlaE is proposed, and Glu155 and Glu251 are implicated as catalytic residues. Mutation and activity measurements for structurally equivalent residues in related epimerases supported this mechanistic proposal.     

Solution structure of a conserved domain of antizyme: a protein regulator of polyamines

Antizyme and its isoforms are members of an unusual yet broadly conserved family of proteins, with roles in regulating polyamine levels within cells. Antizyme has the ability to bind and inhibit the enzyme ornithine decarboxylase (ODC), targeting it for degradation at the proteasome; antizyme is also known to affect the transport of polyamines and interact with the antizyme inhibitor protein (AZI), as well as the cell-cycle protein cyclin D1. In the present work, NMR methods were used to determine the solution structure of a stable, folded domain of mammalian antizyme isoform-1 (AZ-1), consisting of amino acid residues 87-227. The protein was found to contain eight beta strands and two alpha helices, with the strands forming a mixed parallel and antiparallel beta sheet. At the level of primary sequence, antizyme is not similar to any protein of known structure, and results show that antizyme exhibits a novel arrangement of its strands and helices. Interestingly, however, the fold of antizyme is similar to that found in a family of acetyl transferases, as well as translation initiation factor IF3, despite a lack of functional relatedness between these proteins. Structural results, combined with amino acid sequence comparisons, were used to identify conserved features among the various homologues of antizyme and their isoforms. Conserved surface residues, including a cluster of acidic amino acids, were found to be located on a single face of antizyme, suggesting this surface is a possible site of interaction with target proteins such as ODC. This structural model provides an essential framework for an improved future understanding of how the different parts of antizyme play their roles in polyamine regulation.     

TOUCHSTONE II: a new approach to ab initio protein structure prediction

We have developed a new combined approach for ab initio protein structure prediction. The protein conformation is described as a lattice chain connecting C(alpha) atoms, with attached C(beta) atoms and side-chain centers of mass. The model force field includes various short-range and long-range knowledge-based potentials derived from a statistical analysis of the regularities of protein structures. The combination of these energy terms is optimized through the maximization of correlation for 30 x 60,000 decoys between the root mean square deviation (RMSD) to native and energies, as well as the energy gap between native and the decoy ensemble. To accelerate the conformational search, a newly developed parallel hyperbolic sampling algorithm with a composite movement set is used in the Monte Carlo simulation processes. We exploit this strategy to successfully fold 41/100 small proteins (36 approximately 120 residues) with predicted structures having a RMSD from native below 6.5 A in the top five cluster centroids. To fold larger-size proteins as well as to improve the folding yield of small proteins, we incorporate into the basic force field side-chain contact predictions from our threading program PROSPECTOR where homologous proteins were excluded from the data base. With these threading-based restraints, the program can fold 83/125 test proteins (36 approximately 174 residues) with structures having a RMSD to native below 6.5 A in the top five cluster centroids. This shows the significant improvement of folding by using predicted tertiary restraints, especially when the accuracy of side-chain contact prediction is >20%. For native fold selection, we introduce quantities dependent on the cluster density and the combination of energy and free energy, which show a higher discriminative power to select the native structure than the previously used cluster energy or cluster size, and which can be used in native structure identification in blind simulations. These procedures are readily automated and are being implemented on a genomic scale.