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.     

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.     

Mutations of endo-beta-N-acetylglucosaminidase H active site residueAs sp130 anG glu132: activities and conformations

Endo-beta-N-acetylglucosaminidase H hydrolyzes the beta-(1-4)-glycosidic link of the N,N'-diacetylchitobiose core of high-mannose and hybrid asparagine-linked oligosaccharides. Seven mutants of the active site residues, Asp130 and Glu132, have been prepared, assayed, and crystallized. They include single site mutants of each residue to the corresponding amide, to Ala and to the alternate acidic residue, and to the double amide mutant. The mutants of Asp130 are more active than the corresponding Glu132 mutants, consistent with the assignment of the latter residue as the primary catalytic residue. The amide mutants are more active than the alternate acidic residue mutants, which in turn are more active than the Ala mutants. The structures of the Asn mutant of Asp130 and the double mutant are very similar to that of the wild-type enzyme. Several residues surrounding the mutated residues, including some that form part of the core of the beta-barrel and especially Tyr168 and Tyr244, adopt a very different conformation in the structures of the other two mutants of Asp130 and in the Asp mutant of Glu132. The results show that the residues in the upper layers of the beta-barrel can organize into two very distinct packing arrangements that depend on subtle electrostatic and steric differences and that greatly affect the geometry of the substrate-binding cleft. Consequently, the relative activities of several of the mutants are defined by structural changes, leading to impaired substrate binding, in addition to changes in functionality.     

Three-dimensional structure of (1,4)-beta-D-mannan mannanohydrolase from tomato fruit

The three-dimensional crystal structure of tomato (Lycopersicon esculentum) beta-mannanase 4a (LeMAN4a) has been determined to 1.5 A resolution. The enzyme adopts the (beta/alpha)(8) fold common to the members of glycohydrolase family GH5. The structure is comparable with those of the homologous Trichoderma reesei and Thermomonospora fusca beta-mannanases: There is a conserved three-stranded beta-sheet located near the N terminus that stacks against the central beta-barrel at the end opposite the active site. Three noncanonical beta-helices surround the active site. Similar helices are found in T. reesei but not T. fusca beta-mannanase. By analogy with other beta-mannanases, the catalytic acid/base residue is E204 and the nucleophile residue is E318. The active site cleft of L. esculentum beta-mannanase most closely resembles that of the T. reesei isozyme. A model of substrate binding in LeMAN4a is proposed in which the mannosyl residue occupying the -1 subsite of the enzyme adopts the (1)S(5) skew-boat conformation.     

Comparative model of EutB from coenzyme B12-dependent ethanolamine ammonia-lyase reveals a beta8alpha8, TIM-barrel fold and radical catalytic site structural features

The structure of the EutB protein from Salmonella typhimurium, which contains the active site of the coenzyme B12 (adenosylcobalamin)-dependent enzyme, ethanolamine ammonia-lyase, has been predicted by using structural proteomics techniques of comparative modelling. The 453-residue EutB protein displays no significant sequence identity with proteins of known structure. Therefore, secondary structure prediction and fold recognition algorithms were used to identify templates. Multiple three-dimensional template matching (threading) servers identified predominantly beta8alpha8, TIM-barrel proteins, and in particular, the large subunits of diol dehydratase (PDB: 1eex:A, 1dio:A) and glycerol dehydratase (PDB: 1mmf:A), as templates. Consistent with this identification, the dehydratases are, like ethanolamine ammonia-lyase, Class II coenzyme B12-dependent enzymes. Model building was performed by using MODELLER. Models were evaluated by using different programs, including PROCHECK and VERIFY3D. The results identify a beta8alpha8, TIM-barrel fold for EutB. The beta8alpha8, TIM-barrel fold is consistent with a central role of the alpha/beta-barrel structures in radical catalysis conducted by the coenzyme B12- and S-adenosylmethionine-dependent (radical SAM) enzyme superfamilies. The EutB model and multiple sequence alignment among ethanolamine ammonia-lyase, diol dehydratase, and glycerol dehydratase from different species reveal the following protein structural features: (1) a "cap" loop segment that closes the N-terminal region of the barrel, (2) a common cobalamin cofactor binding topography at the C-terminal region of the barrel, and (3) a beta-barrel-internal guanidinium group from EutB R160 that overlaps the position of the active-site potassium ion found in the dehydratases. R160 is proposed to have a role in substrate binding and radical catalysis.     

Structural analysis of the 2.8 A model of Xylose isomerase from Actinoplanes missouriensis

The structure of Xylose isomerase (X.I.) from Actinoplanes missouriensis has been solved to 2.8 Angstroms resolution. Phases were determined from a single Eu3+ derivative and from the noncrystallographic 222 symmetry of the tetrameric molecule. An atomic model was built and subjected to restrained crystallographic refinement. The resulting model is shown to be closely similar to the recently reported X.I.'s structures from three other bacterial sources. Each monomer is found to be composed of an eight-stranded alpha/beta "T.I.M." barrel forming an N-terminal domain of 328 residues followed by a large loop of 66 residues embracing an adjacent subunit. Analysis of intersubunit packing shows that the X.I. tetramer is an assembly of two tight dimers. The beta barrel fits a simple hyperboloid model as other T.I.M. barrels do. The active site, identified as the binding site for the inhibitor xylitol, is located at the carboxyl end of the beta strands in the barrel next to a pair of binding sites for Eu3+ ions, which are assumed to be sites for the divalent ions involved in catalysis. Active sites in the tetramer are oriented towards the interface between dimers. It is suggested that subunit interfaces might stabilize the active site region and this might explain the oligomeric nature of other alpha/beta barrel enzymes.     

Protonation state and substrate binding to B2 metallo-beta-lactamase CphA from Aeromonas hydrofila

The zinc enzymes metallo beta-lactamases counteract the beneficial action of beta-lactam antibiotics against bacterial infections, by hydrolyzing their beta-lactam rings. To understand structure/function relationships on a representative member of this class, the B2 M beta L CphA, we have investigated the H-bond pattern at the Zn enzymatic active site and substrate binding mode by molecular simulation methods. Extensive QM calculations at the DFT-BLYP level on eleven models of the protein active site, along with MD simulations of the protein in aqueous solution, allow us to propose two plausible protonation states for the unbound enzyme, which are probably in equilibrium. Docking procedures along with MD simulations and QM calculations suggest that in the complex between the enzyme and its substrate (biapenem), the latter is stable in only one of the two protonation states, in addition it exhibits two different binding modes, of which only one agrees with previous proposals. In both cases, the substrate is polarized as in aqueous solution. We conclude that addressing mechanistic issues on this class of enzymes requires a careful procedure to assign protonation states and substrate docking modes.     

The crystal structure of human muscle aldolase at 3.0 A resolution

The three-dimensional structure of fructose-1,6-bisphosphate aldolase from human muscle has been determined at 3.0 A resolution by X-ray crystallography. The active protein is a tetramer of 4 identical subunits each of which is composed of an eight-stranded alpha/beta-barrel structure. The lysine residue responsible for Schiff base formation with the substrate is located near the centre of the barrel in the middle of the sixth beta-strand. While the overall topology of the alpha/beta-barrel is very similar to those found in several other enzymes, the distribution of charged residues inside the core of the barrel seems distinct. The quaternary fold of human muscle aldolase uses interfacial regions also involved in the subunit association of other alpha/beta-barrel proteins found in glycolysis, but exploits these regions in a manner not seen previously.     

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.     

Similarity of different β-strands flanked in loops by glycines and prolines from distinct (α/β)8-barrel enzymes: Chance or a homology?

Many (alpha/beta)8-barrel enzymes contain their conserved sequence regions at or around the beta-strand segments that are often preceded and succeeded by glycines and prolines, respectively. alpha-Amylase is one of these enzymes. Its sequences exhibit a very low degree of similarity, but strong conservation is seen around its beta-strands. These conserved regions were used in the search for similarities with beta-strands of other (alpha/beta)8-barrel enzymes. The analysis revealed an interesting similarity between the segment around the beta 2-strand of alpha-amylase and the one around the beta 4-strand of glycolate oxidase that are flanked in loops by glycines and prolines. The similarity can be further extended on other members of the alpha-amylase and glycolate oxidase subfamilies, i.e., cyclodextrin glycosyltransferase and oligo-1,6-glucosidase, and flavocytochrome b2, respectively. Moreover, the alpha-subunit of tryptophan synthase, the (alpha/beta)8-barrel enzyme belonging to the other subfamily of (alpha/beta)8-barrels, has both investigated strands, beta 2 and beta 4, similar to beta 2 of alpha-amylase and beta 4 of glycolate oxidase. The possibilities of whether this similarity exists only by chance or is a consequence of some processes during the evolution of (alpha/beta)8-barrel proteins are briefly discussed.     

Is the activity-linked electrostatic gradient of bovine Cu, Zn superoxide dismutases conserved in homologous enzymes irrespective of the number and distribution of charges?

Electrostatic potential calculations have been performed on three different Cu, Zn superoxide dismutases (superoxide: superoxide oxidoreductase, EC 1.15 1.1 ), in order to evaluate the degree of conservation of the pattern of the electrostatic interactions between O₂⁻ and the active site recently pointed out in bovine Cu Zn SOD. The three Cu, Zn SODs that have been selected for this study, namely the bovine, ovine, and porcine enzymes, are highly homologous as to reasonably assume identical three-dimensional structure but display large differences in their net charge, as shown by their pI's which span over a wide range pHrange: 8.0 (sheep), 6.5(pig), 5.2(ox). Despite such a large difference in the net protein charge and in the spatial arrangement of electrostatic charges, electrostatic potential calculations show that the electrostatic channel directing the negatively charged substrate toward the positive catalytic site is strictly preserved with the same features for the three proteins. This suggests that the electrostatic funnel for conducting small anions into the active site is a highly conservative property in the evolution of Cu, Zn SOD.     

Crystal structure of New Delhi metallo-β-lactamase reveals molecular basis for antibiotic resistance

β‐Lactams are the most commonly prescribed class of antibiotics and have had an enormous impact on human health. Thus, it is disquieting that an enzyme called New Delhi metallo‐β‐lactamase‐1 (NDM‐1) can confer Enterobacteriaceae with nearly complete resistance to all β‐lactam antibiotics including the carbapenams. We have determined the crystal structure of Klebsiella pneumoniae apo‐NDM‐1 to 2.1‐Å resolution. From the structure, we see that NDM‐1 has an expansive active site with a unique electrostatic profile, which we propose leads to a broader substrate specificity. In addition, NDM‐1 undergoes important conformational changes upon substrate binding. These changes have not been previously observed in metallo‐β‐lactamase enzymes and may have a direct influence on substrate recognition and catalysis.     

Regulation of hematopoietic-specific G-protein Galpha15 and Galpha16 by protein kinase C

Heterotrimeric G proteins mediate cell growth and differentiation by coupling cell surface receptors to intracellular effector enzymes. The G-protein alpha subunit, Galpha(16), and its murine homologue Galpha(15), are expressed specifically in hematopoietic cells and their expression is highly regulated during differentiation of normal and leukemic cells. In this study, we examined the phosphorylation of Galpha(15)/Galpha(16) and its role in receptor and effector coupling. We observed a PMA-stimulated intact cell phosphorylation of Galpha(15) in COS7 cells transfected with Galpha(15) and protein kinase Calpha (PKCalpha), and phosphorylation of endogenous Galpha(16) in HL60 cells. We also showed that peptides derived from the two G-proteins were phosphorylated in vitro using purified brain PKC. Furthermore, we identified the putative phosphorylation site and showed that mutation or deletion of this PKC phosphorylation site inhibited phospholipase C (PLC) activation. The behavior of double mutants with the constitutively active G-protein mutation (QL-mutant) and mutation in the putative phosphorylation site suggests that the phosphorylation site of Galpha(15/16) is essential for receptor-coupled activation of PLC, but not for direct interaction of the G-protein with PLC-beta.     

The pattern of amino acid replacements in alpha/beta-barrels

The determinants of site-to-site variability in the rate of amino acid replacement in alpha/beta-barrel enzyme structures are investigated. Of 125 available alpha/beta-barrel structures, only 25 meet a variety of phylogenetic and statistical criteria necessary to ensure sufficient data for reliable analysis. These 25 enzyme structures (from a wide variety of taxa with diverse lifestyles in diverse habitats) differ greatly in size, number, and topology of domains in addition to the alpha/beta-barrel, quaternary structure, metabolic role, reaction catalyzed, presence of prosthetic groups, regulatory mechanisms, use of cofactors, and catalytic mechanisms. Yet, with the exception of ribulose-1,5-bisphosphate carboxylase, all structures have similar frequency distributions of amino acid replacement rates. Hence, site-specific variability in rates of evolution is largely independent of differences in biology, biochemistry, and molecular structure. A correlation between site-specific rate variation and (1) distance from the active site, (2) solvent accessibility, and (3) treating glycines in unusual main-chain conformations as a separate class, explains approximately half the causal variation. Secondary structure exerts little influence on the pattern and distribution of replacements. Additional domains and subunits, side-chain hydrogen bonds, unusual side-chain rotamers, nonplanar peptide bonds, strained main-chain conformations, and buried hydrophilic-charged residues contribute little to variability among sites because they are rare. Nonlinear models do not improve the fits. In several enzymes, deviations from the typical pattern of replacements suggest the possible action of natural selection. A statistical analysis shows that, in all cases, much of the remaining unexplained variation is not attributable to chance and that other, as yet unidentified, causal relations must exist.     

Structural comparisons of TIM barrel proteins suggest functional and evolutionary relationships between beta-galactosidase and other glycohydrolases

Beta-galactosidase (lacZ) from Escherichia coli is a 464 kDa homotetramer. Each subunit consists of five domains, the third being an alpha/beta barrel that contains most of the active site residues. A comparison is made between each of the domains and a large set of proteins representative of all structures from the protein data bank. Many structures include an alpha/beta barrel. Those that are most similar to the alpha/beta barrel of E. coli beta-galactosidase have similar catalytic residues and belong to the so-called "4/7 superfamily" of glycosyl hydrolases. The structure comparison suggests that beta-amylase should also be included in this family. Of three structure comparison methods tested, the "ProSup" procedure of Zu-Kang and Sippl and the "Superimpose" procedure of Diederichs were slightly superior in discriminating the members of this superfamily, although all procedures were very powerful in identifying related protein structures. Domains 1, 2, and 4 of E. coli beta-galactosidase have topologies related to "jelly-roll barrels" and "immunoglobulin constant" domains. This fold also occurs in the cellulose binding domains (CBDs) of a number of glycosyl hydrolases. The fold of domain 1 of E. coli beta-galactosidase is closely related to some CBDs, and the domain contributes to substrate binding, but in a manner unrelated to cellulose binding by the CBDs. This is typical of domains 1, 2, 4, and 5, which appear to have been recruited to play roles in beta-galactosidase that are unrelated to the functions that such domains provide in other contexts. It is proposed that beta-galactosidase arose from a prototypical single domain alpha/beta barrel with an extended active site cleft. The subsequent incorporation of elements from other domains could then have reduced the size of the active site from a cleft to a pocket to better hydrolyze the disaccharide lactose and, at the same time, to facilitate the production of inducer, allolactose.     

Improving macromolecular electrostatics calculations.

Electrostatic interactions play a key role in many aspects of protein engineering. Consequently, much effort has been put into the design of software for calculating electrostatic fields around macromolecules. We show that optimization of hydrogen bonding networks can improve both the results of pK(a) calculations and the results of electrostatic calculations performed by commonly used programs such as DelPhi. Further optimization can often be achieved by flipping the side chains of asparagine, histidine and glutamine around their chi2, chi2 and chi3 torsion angles, respectively, when this improves the local hydrogen bonding network. These optimizations are applied to some well characterized proteins: BPTI, hen egg white lysozyme and superoxide dismutase. A search for flipped residues in the PDB revealed that significant improvements in electrostatic calculations in or near the active site of enzymes can be expected for about one quarter of all enzymes in the PDB.     

Electrostatic fields near the active site of human aldose reductase: 1. New inhibitors and vibrational stark effect measurements

Vibrational Stark effect spectroscopy was used to measure electrostatic fields in the hydrophobic region of the active site of human aldose reductase (hALR2). A new hALR2 inhibitor was designed and synthesized that contains a nitrile probe with a Stark tuning rate of 0.77 cm-1/(MV/cm). Mutations to amino acid residues in the vicinity of the nitrile functional group were selected based on electrostatics calculations, possible complications from hydrogen bonds near the nitrile, and comparison with the active site of human aldehyde reductase, whose structure is very similar. Changes in the absorption energy of the nitrile probe when bound to those mutated proteins were then used to quantify perturbations to the protein's electrostatic field. Electrostatic field changes as large as -10 MV/cm were observed. Measured electrostatic fields were compared to predictions based on continuum electrostatics calculations, revealing that substantial modifications to the calculation strategy are necessary. The effects of hydrogen bonding of amino acid side chains to the nitrile probe are considered, and applications of vibrational Stark effect spectroscopy to investigations of ligand binding and biological function are discussed.     

Electrostatics in the active site of an alpha-amylase.

The importance of electrostatics in catalysis has been emphasized in the literature for a large number of enzymes. We examined this hypothesis for the Bacillus licheniformis alpha-amylase by constructing site-directed mutants that were predicted to change the pKa values of the catalytic residues and thus change the pH-activity profile of the enzyme. To change the pKa of the catalytic residues in the active site, we constructed mutations that altered the hydrogen bonding network, mutations that changed the solvent accessibility, and mutations that altered the net charge of the molecule. The results show that changing the hydrogen bonding network near an active site residue or changing the solvent accessibility of an active site residue will very likely result in an enzyme with drastically reduced activity. The differences in the pH-activity profiles for these mutants were modest. pH-activity profiles of mutants which change the net charge on the molecule were significantly different from the wild-type pH-activity profile. The differences were, however, difficult to correlate with the electrostatic field changes calculated. In several cases we observed that pH-activity profiles shifted in the opposite direction compared to the shift predicted from electrostatic calculations. This strongly suggests that electrostatic effects cannot be solely responsible for the pH-activity profile of the B. licheniformis alpha-amylase.     

Three dimensional structure of porcine pancreatic alpha-amylase at 2.9 A resolution. Role of calcium in structure and activity.

The crystal structure of porcine pancreatic alpha-amylase (PPA) has been solved at 2.9 A resolution by X-ray crystallographic methods. The enzyme contains three domains. The larger, in the N-terminal part, consists of 330 amino acid residues. This central domain has the typical parallel-stranded alpha-beta barrel structure (alpha beta)8, already found in a number of other enzymes like triose phosphate isomerase and pyruvate kinase. The C-terminal domain forms a distinct globular unit where the chain folds into an eight-stranded antiparallel beta-barrel. The third domain lies between a beta-strand and a alpha-helix of the central domain, in a position similar to those found for domain B in triose phosphate isomerase and pyruvate kinase. It is essentially composed of antiparallel beta-sheets. The active site is located in a cleft within the N-terminal central domain, at the carboxy-end of the beta-strands of the (alpha beta)8 barrel. Binding of various substrate analogues to the enzyme suggests that the amino acid residues involved in the catalytic reaction are a pair of aspartic acids. A number of other residues surround the substrate and seem to participate in its binding via hydrogen bonds and hydrophobic interactions. The 'essential' calcium ion has been located near the active site region and between two domains, each of them providing two calcium ligands. On the basis of sequence comparisons this calcium binding site is suggested to be a common structural feature of all alpha-amylases. It represents a new type of calcium-protein interaction pattern.(ABSTRACT TRUNCATED AT 250 WORDS)