Barrel structures in proteins: automatic identification and classification including a sequence analysis of TIM barrels.

Automated methods for identifying and characterizing regular beta-barrels from coordinate data have been developed to analyze and classify various kinds of barrel structures based on geometric parameters such as the barrel strand number (n) and shear number (S). In total, we find 1,316 barrels in the January 1998 release of Protein Data Bank. Of 1,316 barrels, 1,277 barrels had an even shear number, corresponding to 50 nonhomologous families. The (beta alpha)8 triose phosphate isomerase (TIM) barrel (n = 8, S = 8) fold has the largest number of apparently nonhomologous entries, 16, although the trypsin like antiparallel (n = 6, S = 8) barrels (representing only three families) are the most common with 527 barrels. Of all the protein families that exhibit barrel structures, 68% are found to be various kinds of enzymes, the remainder being binding proteins and transport membrane proteins. In addition, the layers of side chains, which form the cores of barrels with S = n and S = 2n, are also analyzed. More sophisticated methods were developed for detecting TIM barrels specifically, including consideration of the amino acid propensities for the side chains that form the layers. We found that the residues on the outside of the eight stranded parallel beta-barrel, buried by the alpha-helices, are much more hydrophobic than the residues inside the barrel.     

Energetic approach to the folding of alpha/beta barrels

The folding of a polypeptide into a parallel (alpha/beta)8 barrel (which is also called a circularly permuted beta 8 alpha 8 barrel) has been investigated in terms of energy minimization. According to the arrangement of hydrogen bonds between two neighboring beta-strands of the central barrel therein, such an alpha/beta barrel structure can be folded into six different types: (1) left-tilted, left-handed crossover; (2) left-tilted, right-handed crossover; (3) nontilted, left-handed crossover; (4) nontilted, right-handed crossover; (5) right-tilted, left-handed crossover; and (6) right-tilted, right-handed crossover. Here "tilt" refers to the orientational relation of the beta-strands to the axis of the central beta-barrel, and "crossover" to the beta alpha beta folding connection feature of the parallel beta-barrel. It has been found that the right-tilted, right-handed crossover alpha/beta barrel possesses much lower energy than the other five types of alpha/beta barrels, elucidating why the observed alpha/beta barrels in proteins always assume the form of right tilt and right-handed crossover connection. As observed, the beta-strands in the energy-minimized right-tilted, right-handed crossover (alpha/beta)8-barrel are of strong right-handed twist. The value of root-mean-square fits also indicates that the central barrel contained in the lowest energy (alpha/beta)8 structure thus found coincides very well with the observed 8-stranded parallel beta-barrel in triose phosphate isomerase (TIM). Furthermore, an energetic analysis has been made demonstrating why the right-tilt, right-handed crossover barrel is the most stable structure. Our calculations and analysis support the principle that it is possible to account for the main features of frequently occurring folding patterns in proteins by means of conformational energy calculations even for very complicated structures such as (alpha/beta)8 barrels.     

Insights into the fold organization of TIM barrel from interaction energy based structure networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or "sequence conservation" as the basis for their understanding. Recently "interaction energy" based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the "interaction conservation" viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.     

One fold with many functions: the evolutionary relationships between TIM barrel families based on their sequences,structures and functions

The eightfold (betaalpha) barrel structure, first observed in triose-phosphate isomerase, occurs ubiquitously in nature. It is nearly always an enzyme and most often involved in molecular or energy metabolism within the cell. In this review we bring together data on the sequence, structure and function of the proteins known to adopt this fold. We highlight the sequence and functional diversity in the 21 homologous superfamilies, which include 76 different sequence families. In many structures, the barrels are "mixed and matched" with other domains generating additional variety. Global and local structure-based alignments are used to explore the distribution of the associated functional residues on this common structural scaffold. Many of the substrates/co-factors include a phosphate moiety, which is usually but not always bound towards the C-terminal end of the sequence. Some, but not all, of these structures, exhibit a structurally conserved "phosphate binding motif". In contrast metal-ligating residues and catalytic residues are distributed along the sequence. However, we also found striking structural superposition of some of these residues. Lastly we consider the possible evolutionary relationships between these proteins, whose sequences are so diverse that even the most powerful approaches find few relationships, yet whose active sites all cluster at one end of the barrel. This extreme example of the "one fold-many functions" paradigm illustrates the difficulty of assigning function through a structural genomics approach for some folds.     

α/β Barrel evolution and the modular assembly of enzymes: Emerging trends in the flavin oxidase/dehydrogenase family

α/β barrels have an ill-defined origin. Evidence exists which favours their divergent evolution from a common ancestral barrel and convergent evolution to a stable fold. However, recent sequence and structural information for the flavin oxidase/dehydrogenase family of barrel enzymes indicate that sub-families of α/β barrels have evolved divergently. The modular fusion of barrel domains with core structures from other gene families has also contributed to the evolution of related but catalytically distinct enzyme molecules within each sub-family of the flavin oxidases/dehydrogenases. An analysis of the structures and sequences of the flavin oxidases/dehydrogenases has now enabled studies focusing on the evolutionary origins and modular assembly of this important family of proteins to be initiated.     

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.     

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.     

Crystal structure of muconate lactonizing enzyme at 3 A resolution

The crystal structure of muconate lactonizing enzyme has been solved at 3 A resolution, and an unambiguous alpha-carbon backbone chain trace made. The enzyme contains three domains; the central domain is a parallel-stranded alpha-beta barrel, which has previously been reported in six other enzymes, including triose phosphate isomerase and pyruvate kinase. One novel feature of this enzyme is that its alpha-beta barrel has only seven parallel alpha-helices around the central core of eight parallel beta-strands; all other known alpha-beta barrels contain eight such helices. The N-terminal (alpha + beta) and C-terminal domains cover the cleft where the eighth helix would be. The active site of muconate lactonizing enzyme has been found by locating the manganese ion that is essential for catalytic activity, and by binding and locating an inhibitor, alpha-ketoglutarate. The active site lies in a cleft between the N-terminal and barrel domains; when the active sites of muconate lactonizing enzyme and triose phosphate isomerase are superimposed, barrel-strand 1 of triose phosphate isomerase is aligned with barrel-strand 3 of muconate lactonizing enzyme. This implies that structurally homologous active-site residues in the two enzymes are carried on different parts of the primary sequence; the ancestral gene would had to have been transposed during its evolution to the modern proteins, which seems unlikely. Therefore, these two enzymes may be related by convergent, rather than divergent, evolution.     

Universally conserved positions in protein folds: reading evolutionary signals about stability, folding kinetics and function.

Here, we provide an analysis of molecular evolution of five of the most populated protein folds: immunoglobulin fold, oligonucleotide-binding fold, Rossman fold, alpha/beta plait, and TIM barrels. In order to distinguish between "historic", functional and structural reasons for amino acid conservations, we consider proteins that acquire the same fold and have no evident sequence homology. For each fold we identify positions that are conserved within each individual family and coincide when non-homologous proteins are structurally superimposed. As a baseline for statistical assessment we use the conservatism expected based on the solvent accessibility. The analysis is based on a new concept of "conservatism-of-conservatism". This approach allows us to identify the structural features that are stabilized in all proteins having a given fold, despite the fact that actual interactions that provide such stabilization may vary from protein to protein. Comparison with experimental data on thermodynamics, folding kinetics and function of the proteins reveals that such universally conserved clusters correspond to either: (i) super-sites (common location of active site in proteins having common tertiary structures but not function) or (ii) folding nuclei whose stability is an important determinant of folding rate, or both (in the case of Rossman fold). The analysis also helps to clarify the relation between folding and function that is apparent for some folds.     

BDNF and trkB mRNA Expression in Neurons of the Neonatal Mouse Barrel Field Cortex: Normal Development and Plasticity after Cauterizing Facial Vibrissae

Development of the central somatosensory system is profoundly modulated by the sensory periphery. Cauterization of facial whiskers alters the segregation pattern of barrels in rodents only during a few days just after birth (critical period). Although a molecular basis of the segregation of barrel neurons and the critical period for the anatomical plasticity observed in layer IV barrel neuron is not clear yet, the accumulating evidence suggests that neurotrophins modulate synaptic connections including central nervous system. In this study, we showed by in situ hybridization that mouse barrel side neurons express brain-derived neurotrophic factor (BDNF) mRNA and both catalytic and non-catalytic forms of trkB mRNA. Cautery of row C vibrissae on the right side of the face within 24 h after birth (post natal day 0, PND0) reduced the expression of BDNF and trkB mRNA from the division region between the contralateral row C barrels at PND7. The vibrissae in row A, C, and E were cauterized at PND0 followed by quantitative RT-PCR for BDNF and trkB mRNA with total RNA isolated from the barrel region at PND7. The result showed that BDNF, but not trkB, mRNA was increased several-fold in the contralateral barrel region. These data suggest that the expression of BDNF mRNA is differentially regulated between injured barrels and actively innervated barrels. The differential expression of the mRNA encoding neurotrophins and their receptors may be important in regulating the injury-dependent re-segregation of barrels.     

Alternative Splice Variants in TIM Barrel Proteins from Human Genome Correlate with the Structural and Evolutionary Modularity of this Versatile Protein Fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.     

In vivo fragment complementation of a (beta/alpha)(8) barrel protein: generation of variability by recombination

The high representation of the TIM barrel as a scaffold for enzymatic proteins makes it an interesting model for protein engineering. Based on previous reports of folding mechanisms of TIM barrels that suggest an independent folding unit formed by six (beta/alpha) subunits, we interrupted the gene of phosphoribosylanthranilate isomerase (PRAI) from Escherichia coli at three different positions to yield fragments with different combinations of (beta/alpha) subunits. When these constructions were expressed as polycistrons in a TrpF-E. coli strain, complementation of the function only occurred with fragments beta1-alpha4 and beta5-alpha8, demonstrating that (beta/alpha)(4) subunits are stable enough to survive in vivo conditions and to assemble to yield a functional enzyme. The expression of these fragments in a separated plasmid/phagemid system to complement the function gave a slower complementation in the TrpF-E. coli strain; this was overcome by introducing extra secondary elements to the structure that reinforce their interaction.     

Common evolutionary origin of swapped-hairpin and double-psi beta barrels

The core of swapped-hairpin and double-psi beta barrels is formed by duplication of a conserved betaalphabeta element, suggesting a common evolutionary origin. The path connecting the two folds is unclear as the two barrels are not interconvertible by a simple topological modification, such as circular permutation. We have identified a protein family whose sequence properties are intermediate to the two folds. The structure of one of these proteins, Pyrococcus horikoshii PhS018, is also built by duplication of the conserved betaalphabeta element but shows yet a third topology, which we name the RIFT barrel. This topology is widespread in the structure database and spans three folds of the SCOP classification, including the middle domain of EF-Tu and the N domain of F1-ATPase. We propose that swapped-hairpin beta barrels arose from an ancestral RIFT barrel by strand invasion and double-psi beta barrels by a strand swap. We group the three barrel types into a metafold, the cradle-loop barrels.     

Dipole moment in TIM alpha/beta fold proteins

TIM proteins of alpha/beta barrel fold from alpha/beta class as given in SCOP database were taken for dipole moment analysis. In all, 32 structures were analyzed for their dipole moment contributions. Representative structures from 20 super families in the alpha/beta fold, with different enzyme functions and 12 protein domains of TIM family in TIM super family were considered. The active sites of these proteins are located on the C-terminal side of the beta-strands. The molecules of same alpha/beta fold, but differing in their functionality also showed a common electrostatic field pattern along the barrel axis and had the dipole moment along the barrel axis and towards C-terminal end of the beta-strands. However, it is observed from our calculations that the dipole moment direction is possibly a consequence of the structural fold, with distribution of charges playing a modulatory role, and does not contribute to the location of active site. We show here that apart from the commonly held view as proposed by Hol et al [Hol W G L, van Duijnen PT and Berendsen H J C (1978) Nature (London), 273, 443-446] of the role of the alpha helical dipole moment, the beta-sheets in the barrel can also have a considerable dipole moment contribution. Taken together with our dipole moment analysis on integral membrane proteins [Vasanthi G and Krishnaswamy S (2002) Indian J Biochem Biophys 39, 93-100], this suggests the need to examine the role of dipole moment in the case of especially beta sheets forming barrels.     

The N-terminal β-sheet of the hyperthermophilic endoglucanase from Pyrococcus horikoshii is critical for thermostability

The β-1,4-endoglucanase (EC 3.2.1.4) from the hyperthermophilic archaeon Pyrococcus horikoshii (EGPh) has strong hydrolyzing activity toward crystalline cellulose. When EGPh is used in combination with β-glucosidase (EC 3.2.1.21), cellulose is completely hydrolyzed to glucose at high temperature, suggesting great potential for EGPh in bioethanol industrial applications. The crystal structure of EGPh shows a triosephosphate isomerase (TIM) (β/α)(8)-barrel fold with an N-terminal antiparallel β-sheet at the opposite side of the active site and a very short C-terminal sequence outside of the barrel structure. We describe here the function of the peripheral sequences outside of the TIM barrel core structure. Sequential deletions were performed from both N and C termini. The activity, thermostability, and pH stability of the expressed mutants were assessed and compared to the wild-type EGPh enzyme. Our results demonstrate that the TIM barrel core is essential for enzyme activity and that the N-terminal β-sheet is critical for enzyme thermostability. Bioinformatics analyses identified potential key residues which may contribute to enzyme hyperthermostability.     

Toward the detection and validation of repeats in protein structure

We present a method called DAVROS to detect, localize, and validate repeating motifs in protein structure allowing for insertions and deletions. DAVROS uses the score matrix from a structural alignment program (SAP) to search for repeating motifs using an algorithm based on concepts from signal processing and the statistical properties of the alignments. The method was tested against a nonredundant Protein Data Bank, and each chain was assigned a score. For the top 50 chains ranked by score, 70% contain repeating motifs detected without error. These represent 14 types of fold covering alpha, beta, and alphabeta protein classes. A second data set comprising protein chains in different sequence families for triosephosphate isomerase (TIM) barrel, leucine-rich repeat (LRR), trefoil, and alpha-alpha barrel folds was used to assess the ability of DAVROS to detect all motifs within a specific fold. For the second test set, the percentage of motifs detected was highest for the LRR chains (88.7%) and least for the TIM barrels (60%). This variability results from the regularity of the LRR motif compared to the alphabeta units of the TIM barrel, which generally have many more indels. These reduce the strength of the repeat signal in the SAP matrix, making repeat detection more difficult.     

Structural basis for the function of pyridoxine 5'-phosphate synthase

BACKGROUND: Pyridoxal 5'-phosphate is the active form of vitamin B(6) that acts as an essential, ubiquitous coenzyme in amino acid metabolism. In Escherichia coli, the pathway of the de novo biosynthesis of vitamin B(6) results in the formation of pyridoxine 5'-phosphate (PNP), which can be regarded as the first synthesized B(6) vitamer. PNP synthase (commonly referred to as PdxJ) is a homooctameric enzyme that catalyzes the final step in this pathway, a complex intramolecular condensation reaction between 1-deoxy-D-xylulose-5'-phosphate and 1-amino-acetone-3-phosphate. RESULTS: The crystal structure of E. coli PNP synthase was solved by single isomorphous replacement with anomalous scattering and refined at a resolution of 2.0 A. The monomer of PNP synthase consists of one compact domain that adopts the abundant TIM barrel fold. Intersubunit contacts are mediated by three additional helices, respective to the classical TIM barrel helices, generating a tetramer of symmetric dimers with 422 symmetry. In the shared active sites of the active dimers, Arg20 is directly involved in substrate binding of the partner monomer. Furthermore, the structure of PNP synthase with its physiological products, PNP and P(i), was determined at 2.3 A resolution, which provides insight into the dynamic action of the enzyme and allows us to identify amino acids critical for enzymatic function. CONCLUSION: The high-resolution structures of the free enzyme and the enzyme-product complex of E. coli PNP synthase suggest essentials of the enzymatic mechanism. The main catalytic features are active site closure upon substrate binding by rearrangement of one C-terminal loop of the TIM barrel, charge-charge stabilization of the protonated Schiff-base intermediate, the presence of two phosphate binding sites, and a water channel that penetrates the beta barrel and allows the release of water molecules in the closed state. All related PNP synthases are predicted to fold into a similar TIM barrel pattern and have comparable active site architecture. Thus, a common mechanism can be anticipated.     

Identification and analysis of catalytic TIM barrel domains in seven further glycoside hydrolase families

Fold recognition results allocate catalytic triose phosphate isomerase (TIM) barrels to seven previously unassigned glycoside hydrolase (GH) families, numbers 29, 44, 50, 71, 84, 85 and 89, enabling prediction of catalytic residues. Modelling of GH family 50 suggests that it may be the common evolutionary ancestor of families 42 and 14. TIM barrels now comprise the catalytic domains of more than half of the assigned GH families, and catalyse a much larger variety of GH reactions than any other catalytic domain architecture. Only 327 GH sequences still have no structurally identified catalytic domain.     

Kinetic and structural properties of triosephosphate isomerase from Helicobacter pylori

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate and D-glyceraldehyde-3-phosphate in the glycolysis-gluconeogenesis metabolism pathway. The Helicobacter pylori TIM gene (HpTIM) was cloned, and HpTIM was expressed and purified. The enzymatic activity of HpTIM for the substrate GAP was determined (K(m) = 3.46 +/- 0.23 mM and k(cat) = 8.8 x 10(4) min(-1)). The crystal structure of HpTIM was determined by molecular replacement at 2.3 A resolution. The overall structure of HpTIM was (beta/alpha)beta(beta/alpha)(6), which resembles the common TIM barrel fold, (beta/alpha)(8); however, a helix is missing after the second beta-strand. The conformation of loop 6 and binding of phosphate ion suggest that the determined structure of HpTIM was in the "closed" state. A highly conserved Arg-Asp salt bridge in the "DX(D/N)G" motif of most TIMs is absent in HpTIM because the sequence of this motif is "(211)SVDG(214)." To determine the significance of this salt bridge to HpTIM, four mutants, including K183S, K183A, D213Q, and D213A, were constructed and characterized. The results suggest that this conserved salt bridge is not essential for the enzymatic activity of HpTIM; however, it might contribute to the conformational stability of HpTIM.     

What distinguishes GroEL substrates from other Escherichia coli proteins?

Experimental studies and theoretical considerations have shown that only a small subset of Escherichia coli proteins fold in vivo with the help of the GroE chaperone system. These proteins, termed GroE substrates, have been divided into three classes: (a) proteins that can fold independently, but are found to associate with GroEL; (b) proteins that require GroE when the cell is under stress; and (c) 'obligatory' proteins that require GroE assistance even under normal conditions. It remains unclear, however, why some proteins need GroE and others do not. Here, we review experimental and computational studies that addressed this question by comparing the sequences and structural, biophysical and evolutionary properties of GroE substrates with those of nonsubstrates. In general, obligatory substrates are found to have lower folding propensities and be more aggregation prone. GroE substrates are also more conserved than other proteins and tend to utilize more optimal codons, but this latter feature is less apparent for obligatory substrates. There is no evidence, however, for any specific sequence signatures although there is a tendency for sequence periodicity. Our review shows that reliable sequence- or structure-based predictions of GroE dependency remain a challenge. We suggest that the different classes of GroE substrates be studied separately and that proper control test sets (e.g. TIM barrel proteins that need GroE for folding versus TIM barrels that fold independently) be used more extensively in such studies.