BrEPS: a flexible and automatic protocol to compute enzyme-specific sequence profiles for functional annotation.

BACKGROUND: Models for the simulation of metabolic networks require the accurate prediction of enzyme function. Based on a genomic sequence, enzymatic functions of gene products are today mainly predicted by sequence database searching and operon analysis. Other methods can support these techniques: We have developed an automatic method "BrEPS" that creates highly specific sequence patterns for the functional annotation of enzymes. RESULTS: The enzymes in the UniprotKB are identified and their sequences compared against each other with BLAST. The enzymes are then clustered into a number of trees, where each tree node is associated with a set of EC-numbers. The enzyme sequences in the tree nodes are aligned with ClustalW. The conserved columns of the resulting multiple alignments are used to construct sequence patterns. In the last step, we verify the quality of the patterns by computing their specificity. Patterns with low specificity are omitted and recomputed further down in the tree. The final high-quality patterns can be used for functional annotation. We ran our protocol on a recent Swiss-Prot release and show statistics, as well as a comparison to PRIAM, a probabilistic method that is also specialized on the functional annotation of enzymes. We determine the amount of true positive annotations for five common microorganisms with data from BRENDA and AMENDA serving as standard of truth. BrEPS is almost on par with PRIAM, a fact which we discuss in the context of five manually investigated cases. CONCLUSIONS: Our protocol computes highly specific sequence patterns that can be used to support the functional annotation of enzymes. The main advantages of our method are that it is automatic and unsupervised, and quite fast once the patterns are evaluated. The results show that BrEPS can be a valuable addition to the reconstruction of metabolic networks.

The BRENDA enzyme information system-From a database to an expert system.

Enzymes, representing the largest and by far most complex group of proteins, play an essential role in all processes of life, including metabolism, gene expression, cell division, the immune system, and others. Their function, also connected to most diseases or stress control makes them interesting targets for research and applications in biotechnology, medical treatments, or diagnosis. Their functional parameters and other properties are collected, integrated, and made available to the scientific community in the BRaunschweig ENzyme DAtabase (BRENDA). In the last 30 years BRENDA has developed into one of the most highly used biological databases worldwide. The data contents, the process of data acquisition, data integration and control, the ways to access the data, and visualizations provided by the website are described and discussed.

The open conformation of a Pseudomonas lipase.

BACKGROUND: . The interfacial activation of lipases results primarily from conformational changes in the enzymes which expose the active site and provide a hydrophobic surface for interaction with the lipid substrate. Comparison of the crystallization conditions used and the structures observed for a variety of lipases suggests that the enzyme conformation is dependent on solution conditions. Pseudomonas cepacia lipase (PCL) was crystallized in conditions from which the open, active conformation of the enzyme was expected. Its three-dimensional structure was determined independently in three different laboratories and was compared with the previously reported closed conformations of the closely related lipases from Pseudomonas glumae (PGL) and Chromobacterium viscosum (CVL). These structures provide new insights into the function of this commercially important family of lipases. RESULTS: . The three independent structures of PCL superimpose with only small differences in the mainchain conformations. As expected, the observed conformation reveals a catalytic site exposed to the solvent. Superposition of PCL with the PGL and CVL structures indicates that the rearrangement from the closed to the open conformation involves three loops. The largest movement involves a 40 residue stretch, within which a helical segment moves to afford access to the catalytic site. A hydrophobic cleft that is presumed to be the lipid binding site is formed around the active site. CONCLUSIONS: . The interfacial activation of Pseudomonas lipases involves conformational rearrangements of surface loops and appears to conform to models of activation deduced from the structures of fungal and mammalian lipases. Factors controlling the conformational rearrangement are not understood, but a comparison of crystallization conditions and observed conformation suggests that the conformation of the protein is determined by the solution conditions, perhaps by the dielectric constant.

Effects of carbohydrate depletion on the structure, stability and activity of glucose oxidase from Aspergillus niger.

Glucose oxidase from Aspergillus niger was purified to homogeneity by hydrophobic interaction and ion-exchange chromatography. Approx. 95% of the carbohydrate moiety was cleaved from the protein by incubation of glucose oxidase with endoglycosidase H and alpha-mannosidase. Cleavage of the carbohydrate moiety effected a 24-30% decrease in the molecular weight and a reduction in the number of isoforms of glucose oxidase. No significant changes were observed in the circular dichroism spectra of the deglycosylated enzyme. Other properties, such as thermal stability, pH and temperature optima of glucose oxidase activity and substrate specificity were not affected. However, removal of the carbohydrate moiety marginally affected the kinetics of glucose oxidation and stability at low pH. From these results it appears that the carbohydrate chain of glucose oxidase does not contribute significantly to the structure, stability and activity of glucose oxidase.

Identification of amino acids responsible for the oxygen sensitivity of ferredoxins from Anabaena variabilis using site-directed mutagenesis.

The filamentous cyanobacterium Anabaena variabilis (ATCC 29413) possesses two molybdenum dependent nitrogenase systems, nif1 and nif2. The nif1 system is regulated by a developmental program involving heterocyst differentiation; the nif2 system is expressed in all cells only under anaerobic conditions and the expression is controlled environmentally. The genes fdxH1 and fdxH2, encoding two [2Fe-2S] ferredoxins, are part of the these two distinct and differently regulated nif gene clusters. The sensitivity of both ferredoxins to oxygen was different; the half-life of FdxH2 in air was only approximately 1.5 h, while FdxH1 retained 80% of its nitrogenase activity after 24 h. We used site-directed mutagenesis to identify the role of individual amino acid residues responsible for oxygen sensitivity and found out that the FdxH2 double mutant I76A/V77L was much more resistant to oxygen than the wild-type ferredoxin (FdxH2) and similar to FdxH1. By modelling it was shown that the accessibility of the cavity around the iron-sulfur cluster was responsible for that.

Amino acid similarity coefficients for protein modeling and sequence alignment derived from main-chain folding angles.

A set of "similarity-parameters" was calculated that reflects the influence of the proteinogenic amino acids on the structure of the protein backbone. The parameters were derived from a detailed analysis of the amino acid specific main-chain torsion angle distributions as they are found in proteins (highly resolved protein structures from the Brookhaven Protein Data Bank). The purpose of these parameters is threefold: (1) they should help in estimating the structural effect of an amino acid substitution during the design of new mutants in protein-engineering; (2) in modeling by homology they should mark places in the protein where changes in the folding are expected; and (3) they should form a scoring matrix in protein sequence alignment superior to identity scoring. The usability of the "structure derived correlation matrix (SCM)" for these purposes is assessed and demonstrated for some examples in the paper.

Solution structure of a variant of human pancreatic secretory trypsin inhibitor determined by nuclear magnetic resonance spectroscopy.

The structure of a variant of human pancreatic secretory trypsin inhibitor (PSTI) has been determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and a combination of distance geometry and molecular dynamics simulations. After complete assignment of the 1H signals, the nuclear Overhauser data imply the existence of a rather well-determined tertiary structure stabilized by a central alpha-helix and a short three-stranded beta-sheet. The tertiary structure of the amino terminus and of the loop 11-17 could not be defined by n.m.r. data, suggesting a high flexibility in these areas. As the crystal structures of two complexes of human PSTI variants and that of an uncomplexed variant are also known a comparison of the PSTI tertiary structure in solution and in the crystal is now possible.

Hydrogen bonding and molecular surface shape complementarity as a basis for protein docking.

A geometric docking algorithm based upon correlation analysis for quantification of geometric complementarity between protein molecular surfaces in close interfacial contact has been developed by a detailed optimization of the conformational search of the algorithm. In order to reduce the entire conformation space search required by the method a physico-chemical pre-filter of conformation space has been developed based upon the a priori assumption that two or more intermolecular hydrogen bonds are intrinsic to the mechanism of binding within protein complexes. Donor sites are defined spatially and directionally by the positions of explicitly calculated donor hydrogen atoms, and the vector space within a defined range about the donor atom-hydrogen atom bond vector. Acceptor sites are represented spatially and directionally by the van der Waals molecular surface points having normal vectors within a predefined range of vector space about the acceptor atom covalent bond vector(s). Geometric conditions necessary for the simultaneous hydrogen bonding interaction between both sites of functionally congruent hydrogen bonding site pairs, located on the individual proteins, are then tested on the basis of a transformation invariant parameterization of the site pair spatial and directional properties. Sterically acceptable conformations defined by interaction of functionally, spatially, and directionally compatible site pairs are then refined to a maximum contact of complementary contact surfaces using the simplex method for the angular search and correlation techniques for the translational search. The utility of the spatial and directional properties of hydrogen bonding donor and acceptor sites for the identification of candidate docking conformations is demonstrated by the reliable preliminary reduction of conformation space, the improved geometric ranking of the minimum RMS conformations of some complexes and the overall reduction of CPU time obtained.

Prediction of protein three-dimensional structures in insertion and deletion regions: a procedure for searching data bases of representative protein fragments using geometric scoring criteria.

The prediction of protein structure in insertion/deletion regions (referred to as indels) is an important part of protein model building by homology. Here we combine cluster analysis with data base search procedures. Initially, data bases of representative protein fragments are constructed using two different clustering algorithms. In the HCAPD (hierarchical clustering after preliminary division) approach, all protein fragments are divided into classes with similar anchor region structures (a protein fragment consists of two anchoring regions and a central region). Within these classes the fragments are further clustered using a hierarchical cluster algorithm. The DCANN (deterministic clustering by assignment of all nearest neighbours) approach is a variant of the k-nearest neighbours cluster algorithm. Only geometric scoring criteria are used for data base searching. The main advantage of a non-redundant data base is the ability to provide structurally different fragments during the search process, which leads to an improvement in structure prediction. Both methods have been tested on 71 insertions and 74 deletions with lengths between one and eight residues.

Crystal structure of a ternary complex of D-2-hydroxyisocaproate dehydrogenase from Lactobacillus casei, NAD+ and 2-oxoisocaproate at 1.9 A resolution.

D-2-hydroxyisocaproate dehydrogenase (D-HicDH) from Lactobacillus casei is a homodimer with 333 amino acids and a molecular mass of 37 kDa per subunit. The enzyme belongs to the protein family of NAD+-dependent D-2-hydroxycarboxylate dehydrogenases and within this family to the subgroup of D-lactate dehydrogenases (D-LDHs). Compared with other D-LDHs D-HicDH is characterized by a very low specificity regarding size and chemical constitution of the accepted D-2-hydroxycarboxylates. Hexagonal crystals of recombinant D-HicDH in the presence of NAD+ and 2-oxoisocaproate (4-methyl-2-oxopentanoate) were grown with ammonium sulphate as precipitating agent. The structure of these crystals was solved by molecular replacement and refined to a final R-factor of 19.6% for all measured X-ray reflections in the resolution range (infinity to 1.86 A). Both NAD+ and 2-oxoisocaproate were identified in the electron density map; binding of the latter in the active site, however, competes with a sulphate ion, which is also defined by electron density. Additionally the final model contains 182 water molecules and a second sulphate ion. The binding of both an in vitro substrate and the natural cosubstrate in the active site provides substantial insight into the catalytic mechanism and allows us to assess previously published active site models for this enzyme family, in particular the two most controversial points, the role of the conserved Arg234 and substrate binding. Furthermore the overall topology and details of the D-HicDH structure are described, discussed against the background of homologous structures and compared with one closely and one distantly related protein.

Crystal structure of L-2-hydroxyisocaproate dehydrogenase from Lactobacillus confusus at 2.2 A resolution. An example of strong asymmetry between subunits.

L-2-Hydroxyisocaproate dehydrogenase (L-HicDH) from Lactobacillus confusus, a homotetramer with a molecular mass of 33 kDa per subunit, belongs to the protein family of the NAD(+)-dependent L-2-hydroxycarboxylate dehydrogenases. L-HicDH was crystallized with ammonium sulphate as precipitant in the presence of NAD+. The crystals belong to the trigonal space group P3(2)21, with a = 135.9 A and c = 205.9 A, and diffract X-rays to 2.2 A resolution. The crystal structure was solved by Patterson search and molecular replacement techniques and refined to an R-value of 21.4% (2.2 to 8 A). The final structure model contains one NAD+ molecule and one sulphate ion per subunit, with 309 water molecules. An unusual feature of this crystal structure is the deviation of the protein subunits from non-crystallographic symmetry, which is so strong that it can be detected globally by self-rotation calculations in reciprocal space. This asymmetry is especially pronounced in the environment of the active site; it is reflected also in the nicotinamide conformation of NAD+ and allows some conclusions to be drawn about the catalytic mechanism. In this context, an "inner active site loop" is identified as a structural element of fundamental functional importance. Furthermore, with knowledge of the crystal structure of L-HicDH the differences in substrate specificity between L-HicDH and the L-lactate dehydrogenases can be partly explained.

Crystallization and preliminary characterization of crystals of D-2-hydroxyisocaproate dehydrogenase from Lactobacillus casei.

D-2-hydroxyisocaproate dehydrogenase (D-HicDH) from Lactobacillus casei is a homodimeric enzyme with a molecular mass of 74.6 kDa. It catalyzes the reduction of a wide range of 2-ketocarboxylic acids to D-2-hydroxycarboxylic acids using NADH as co-substrate. The enzyme has been crystallized by vapor diffusion using ammonium sulfate as precipitant. The crystals belong to hexagonal space group type P6(3)22 with a = b = 134.1 A, c = 124.1 A and diffract X-rays to 3.0 A resolution. Packing considerations show that there are either one or two D-HicDH monomers in the asymmetric unit.

X-ray structure determination and comparison of two crystal forms of a variant (Asn115Arg) of the alkaline protease from Bacillus alcalophilus refined at 1.85 A resolution.

The X-ray structure determination, refinement and comparison of two crystal forms of a variant (Asn115Arg) of the alkaline protease from Bacillus alcalophilus is described. Under identical conditions crystals were obtained in the orthorhombic space group P2(1)2(1)2(1) (form I) and the rhombohedral space group R32 (form II). For both space groups the structures of the protease were solved by molecular replacement and refined at 1.85 A resolution. The final R-factors are 17.9% and 17.1% for form I and form II, respectively. The root-mean-square deviation between the two forms is 0.48 A and 0.86 A for main-chain and side-chain atoms, respectively. Due to differences in crystal lattice contacts and packing, the structures of the two crystal forms differ in intermolecular interaction affecting the local conformation of three flexible polypeptide sequences (Ser50-Glu55, Ser99-Gly102, Gly258-Ser259) at the surface of the protein. While the two overall structures are very similar, the differences are significantly larger than the errors inherent in the structure determination. As expected, the differences in the temperature factors in form I and II are correlated with the solvent accessibility of the corresponding amino acid residues. In form II, two symmetry-related substrate binding sites face each other, forming a tight intermolecular interaction. Some residues contributing to this intermolecular interaction are also found to be involved in the formation of the complex between subtilisin Carlsberg and the proteinaceous inhibitor eglin C. This demonstrates that the two symmetry-related molecules interact with each other at the same molecular surface area that is used for binding of substrates and inhibitors.

Three-dimensional structure of the complexes between bovine chymotrypsinogen A and two recombinant variants of human pancreatic secretory trypsin inhibitor (Kazal-type).

Variants of the human pancreatic secretory trypsin inhibitor (PSTI) have been created during a protein design project to generate a high-affinity inhibitor with respect to some serine proteases other than trypsin. Two modified versions of human PSTI with high affinity for chymotrypsin were crystallized as a complex with chymotrypsinogen. Both crystallize isomorphously in space group P4(1)2(1)2 with lattice constants a = 84.4 A, c = 86.7 A and diffract to 2.3 A resolution. The structure was solved by molecular replacement. The final R-value after refinement with 8.0 to 2.3 A resolution data was 19.5% for both complexes after inclusion of about 50 bound water molecules. The overall three-dimensional structure of PSTI is similar to the structure of porcine PSTI in the trypsinogen complex (1TGS). Small differences in the relative orientation of the binding loop and the core of the inhibitors indicate flexible adaptation to the proteases. The chymotrypsinogen part of the complex is similar to chymotrypsin. After refolding induced by binding of the inhibitor the root-mean-square difference of the active site residues A186 to A195 and A217 to A222 compared to chymotrypsin was 0.26 A.

Crystallization and preliminary X-ray diffraction studies of a deglycosylated glucose oxidase from Aspergillus niger.

Deglycosylation was shown to be an important prerequisite step for the crystallization of glucose oxidase from Aspergillus niger. Whereas the glycosylated enzyme could not be crystallized, crystals of the deglycosylated enzyme suitable for X-ray diffraction analysis were reproducibly grown in the presence of 1.6 M-ammonium sulphate and octanetriol at pH 5.3 to 5.6. The crystals belong to the space group P3(1)21 or P3(2)21 with refined lattice constants of a = 66.5 A and c = 214.4 A, indicating a cell content of one monomer per asymmetric unit of the crystal. Crystals diffract to at least 2.5 A resolution. Cleavage of 95% of its carbohydrate moiety affected the kinetics of glucose oxidation, stability at low pH and some electrophoretic properties of glucose oxidase, such as molecular mass and the number of isoelectric forms. However, other properties, such as thermal stability, pH and temperature optima of activity were not affected.

Three-dimensional structure of a recombinant variant of human pancreatic secretory trypsin inhibitor (Kazal type).

A modified version of the human pancreatic trypsin inhibitor (PSTI), generated in a protein-design project, has been crystallized in spacegroup P4(3) with lattice constants a = 40.15 A, c = 33.91 A. The structure has been solved by molecular replacement. Refinement of the structure by simulated annealing and conventional restrained least-squares yielded for 8.0 to 2.3 A data a final R-value of 19.1%. Differences to the known structures of porcine PSTI complexed with trypsinogen and modified human PSTI complexed with chymotrypsinogen occur at the flexible N-terminal part of the molecule. These differences are influenced by crystal packing, as are low temperature factors for the binding loop. The geometry of the binding loop is similar to the complexed structures.

Crystallization and preliminary X-ray data of bromoperoxidase from Streptomyces aureofaciens ATCC 10762.

Bromoperoxidase from Streptomyces aureofaciens ATCC 10762, a non-haem haloperoxidase, has been crystallized using the hanging drop method. Preliminary X-ray diffraction studies show that the crystals belong to the cubic space group P2(1)3 with a = 123.4 A. The asymmetric unit contains a dimer of Mr = 60,200. The crystals diffract to at least 2.3 A resolution and are suitable for crystallographic structure analysis.

Crystallization and preliminary X-ray diffraction studies of a deglycosylated glucose oxidase from Penicillium amagasakiense.

The dimeric glucose oxidase from Penicillium amagasakiense was deglycosylated, purified and crystallized as a complex with its coenzyme FAD. Deglycosylation and purification to isoelectric homogeneity were shown to be an important prerequisite step to obtain crystals suitable for X-ray investigations. Crystals of the deglycosylated enzyme were reproducibly grown using ammonium sulfate as precipitant at pH 7.4 to 7.5. Crystals diffract to at least 2.0 A resolution and belong to the orthorhombic space group P2(1)2(1)2(1), with refined lattice constants of a = 59.3 A, b = 136.3 A and c = 156.7 A. Assuming two monomers (approximately 135 kDa) per asymmetric unit the Vm value is 2.3 A3/Da.

Crystal structure of a bacterial lipase from Chromobacterium viscosum ATCC 6918 refined at 1.6 angstroms resolution.

The crystal structure of a lipase from the bacterium Chromobacterium viscosum ATCC 6918 (CVL) has been determined by isomorphous replacement and refined at 1.6 angstroms resolution to an R-factor of 17.8%. The lipase has the overall topology of an alpha/beta type protein, which was also found for previously determined lipase structures. The catalytic triad of the active center consists of the residues Ser87, Asp263 and His285. These residues are not exposed to the solvent, but a narrow channel connects them with the molecular surface. This conformation is very similar to the previously reported closed conformation of Pseudomonas glumae lipase (PGL), but superposition of the two lipase structures reveals several conformational differences. r.m.s. deviations greater than 2 angstroms are found for the C alpha-atoms of the polypeptide chains from His15 to Asp28, from Leu49 to Ser54 and from Lys128 to Gln158. Compared to the PGL structure in the CVL structure, three alpha-helical fragments are shorter, one beta-strand is longer and an additional antiparallel beta-sheet is found. In contrast to PGL, CVL displays an oxyanion hole, which is stabilized by the amide nitrogen atoms of Leu17 and Gln88, and a cis-peptide bond between Gln291 and Leu292. CVL contains a Ca2+, like the PGL, which is coordinated by four oxygen atoms from the protein and two water molecules.

Crystallization and preliminary X-ray diffraction studies of a NADH oxidase from Thermus thermophilus HB8.

The thermophile NADH oxidase from Thermus thermophilus, cloned and expressed in Escherichia coli, has been purified to homogeneity and crystallized. Three different crystal forms were found to be suitable for X-ray diffraction analysis. Crystals of the tetragonal form, grown in the presence of 25% polyethylene glycol 4000 and 0.25 M-NaCl at pH 6.6, were chosen for further analysis. These crystals belong to the space group P4(1)(3)2(1)2 with refined lattice constants of a = 94.8 A and c = 49.0 A, indicating a cell content of one monomer per asymmetric unit of the crystal. The crystals diffract to a resolution of 2.2 A.