Analysis and prediction of the location of catalytic residues in enzymes

The catalytic residues of an enzyme are defined as the amino acids directly involved in chemical catalysis. They mainly act as a general acid--base, electrophilic or nucleophilic catalyst or they polarize and stabilize the transition state. An analysis of the structural features of 36 catalytic residues in 17 enzymes of known structure and with defined mechanism is reported. Residues that bind metal ions (Zn2+ and Cu2+) are considered separately. The features examined are: residue type, location in secondary structure, separation between the residues, accessibility to solvent, intra-protein electrostatic interactions, mobility as evaluated from crystallographic temperature factors, polarity of the environment and the sequence conservation between homologous enzymes of residues that were sequentially or spatially close to the catalytic residue. In general the environment of catalytic residues is similar to that of polar side chains that have low accessibility to solvent. Two algorithms have been developed to identify probable catalytic residues. Scanning an alignment of homologous enzyme sequences for peaks of sequence conservation identifies 13 out of the 16 catalytic residues with 50 residues overpredicted. When the conservation of the spatially close residues is used instead, a different set of 13 residues are identified with 47 residues overpredicted. A combination of the two algorithms identifies 11 residues with 36 residues overpredicted.     

Homology of Escherichia coli B glutathione synthetase with dihydrofolate reductase in amino acid sequence and substrate binding site

Glutathione synthetase from Escherichia coli B showed amino acid sequence homology with mammalian and bacterial dihydrofolate reductases over 40 residues, although these two enzymes are different in their reaction mechanisms and ligand requirements. The effects of ligands of dihydrofolate reductase on the reaction of E. coli B glutathione synthetase were examined to find resemblances in catalytic function to dihydrofolate reductase. The E. coli B enzyme was potently inhibited by 7,8-dihydrofolate, methotrexate, and trimethoprim. Methotrexate was studied in detail and proved to bind to an ATP binding site of the E. coli B enzyme with K1 value of 0.1 mM. The homologous portion of the amino acid sequence in dihydrofolate reductases, which corresponds to the portion coded by exon 3 of mammalian dihydrofolate reductase genes, provided a binding site of the adenosine diphosphate moiety of NADPH in the crystal structure of dihydrofolate reductase. These analyses would indicate that the homologous portion of the amino acid sequence of the E. coli B enzyme provides the ATP binding site. This report gives experimental evidence that amino acid sequences related by sequence homology conserve functional similarity even in enzymes which differ in their catalytic mechanisms.     

Using functional domain composition to predict enzyme family classes

According to their main EC (Enzyme Commission) numbers, enzymes are classified into the following 6 main classes: oxidoreductases, transferases, hydrolases, lyases, isomerases, and ligases. A new method has been developed to predict the enzymatic attribute of proteins by introducing the functional domain composition to formulate a given protein sequence. The advantage by doing so is that both the sequence-order-related features and the function-related features are naturally incorporated in the predictor. As a demonstration, the jackknife cross-validation test was performed on a dataset that consists of proteins with only less than 20% sequence identity to each other in order to get rid of any homologous bias. The overall success rate thus obtained was 85% in identifying the enzyme family classes (including the identification of nonenzyme protein sequences as well). The success rate is significantly higher than those obtained by the other methods on such a stringent dataset. This indicates that using the functional domain composition to represent protein samples for statistical prediction is indeed very promising, and will become a powerful tool in bioinformatics and proteomics.     

Crystal structure of a fragment of mouse ubiquitin-activating enzyme

Protein ubiquitination requires the sequential activity of three enzymes: a ubiquitin-activating enzyme (E1), a ubiquitin-conjugating enzyme (E2), and a ubiquitin-ligase (E3). The ubiquitin-transfer machinery is hierarchically organized; for every ubiquitin-activating enzyme, there are several ubiquitin-conjugating enzymes, and most ubiquitin-conjugating enzymes can in turn interact with multiple ubiquitin ligases. Despite the central role of ubiquitin-activating enzyme in this cascade, a crystal structure of a ubiquitin-activating enzyme is not available. The enzyme is thought to consist of an adenylation domain, a catalytic cysteine domain, a four-helix bundle, and possibly, a ubiquitin-like domain. Its adenylation domain can be modeled because it is clearly homologous to the structurally known adenylation domains of the activating enzymes for the small ubiquitin-like modifier (SUMO) and for the protein encoded by the neuronal precursor cell-expressed, developmentally down-regulated gene 8 (NEDD8). Low sequence similarity and vastly different domain lengths make modeling difficult for the catalytic cysteine domain that results from the juxtaposition of two catalytic cysteine half-domains. Here, we present a biochemical and crystallographic characterization of the two half-domains and the crystal structure of the larger, second catalytic cysteine half-domain of mouse ubiquitin-activating enzyme. We show that the domain is organized around a conserved folding motif that is also present in the NEDD8- and SUMO-activating enzymes, and we propose a tentative model for full-length ubiquitin-activating enzyme.     

The formate complex of the cytochrome bo quinol oxidase of Escherichia coli exhibits a 'g = 12' EPR feature analogous to that of 'slow' cytochrome oxidase.

The cytochrome bo quinol oxidase of Escherichia coli is homologous in sequence and in structure to cytochrome aa3 type cytochrome oxidase in subunit I, which contains the catalytic core. The cytochrome bo enzyme forms a formate complex which exhibits 'g = 12' and 'g = 2.9' EPR signals at X band; similar signals have previously been observed only in association with the 'slow' and formate-ligand states of cytochrome oxidase. These signals arise from transitions within integral spin multiples identified with the homologous heme-copper binuclear catalytic centers in both enzymes.     

Formation of a Hybrid Enzyme between Echinoderm Arginine Kinase and Mammalian Creatine Kinase

MOST intracellular enzymes are composed of more than one subunit. Where this is so, specificity in the interaction between subunits must be as important as specificity for substrates and the integrity of the catalytic structures of the molecule. All these features of enzyme proteins arise from their specific primary structures and may be modified during evolution as a result of changes in amino-acid sequence. The problem is posed as to whether the chemical structures within a subunit enabling it to combine with others are conserved during evolution or, on the other hand, whether they are modified during the divergence of homologous enzymes and isoenzymes.     

Molecular characterization of xynX, a gene encoding a multidomain xylanase with a thermostabilizing domain from Clostridium thermocellum

A Clostridium thermocellum gene, xynX, coding for a xylanase was cloned and the complete nucleotide sequence was determined. The xylanase gene of Clostridium thermocellum consists of an ORF of 3261 nucleotide encoding a xylanase (XynX) of 1087 amino acid residues (116 kDa). Sequence analysis of XynX showed a multidomain structure that consisted of four different domains: an N-terminal thermostabilizing domain homologous to sequences found in several thermophilic enzymes, a catalytic domain homologous to family 10 glycosyl hydrolases, a duplicated cellulose-binding domain (CBD) homologous to family IX CBDs, and a triplicated S-layer homologous domain. A deletion mutant of xynX having only the catalytic region produced a mutant enzyme XynX-C which retained catalytic activity but lost thermostability. In terms of half-life at 70 °C, the thermostability of XynX-C was about six times lower than that of the other mutant enzyme, XynX-TC, produced by a mutant containing both the thermostabilizing domain and the catalytic domain. The optimum temperature of XynX-C was about 5–10 °C lower than that of XynX-TC. Received: 12 January 2000 / Received revision: 24 April 2000 / Accepted: 1 May 2000     

Histidine-272 of isopenicillin N synthase of Cephalosporium acremonium, which is possibly involved in iron binding, is essential for its catalytic activity

Abstract Amino acid sequence alignment of the Cephalosporium acremonium isopenicillin N synthase (cIPNS) to similar non-heme Fe²⁺-containing enzymes from 28 different sources (bacterial, fungal, plant and animals) revealed a homologous region of high sequence conservation containing an invariant histidine residue at position 272 in cIPNS. The importance of this histidine residue in cIPNS was investigated through site-directed mutagenesis by replacing the histidine residue with leucine. The mutated gene was verified by DNA sequence analysis and expressed in Escherichia coli . When analyzed by denaturing gel electrophoresis and immunoblotting, the mutant cIPNS had identical mobility as that of the wild-type enzyme. Enzyme studies on the mutant enzyme showed loss of enzymatic activity indicating that His272 is essential for the catalytic function of cIPNS, possibly as a ligand for iron binding.     

Structure and evolution of a group of related aminoacyl-tRNA synthetases

A yeast nuclear gene, designated MSK1, has been selected from a yeast genomic library by transformation of a respiratory deficient mutant impaired in acylation of mitochondrial lysine tRNA. This gene confers a respiratory competent phenotype and restores the mutant's ability to acylate the mitochondrial lysine tRNA. The amino acid sequence of the protein encoded by MSK1 is homologous to yeast cytoplasmic lysyl-tRNA synthetase and to the product of the herC gene, which has recently been suggested to code for the Escherichia coli enzyme. These observations indicate that MSK1 codes for the lysyl-tRNA synthetase of yeast mitochondria. Several regions of high primary sequence conservation have been identified in the bacterial and yeast lysyl-tRNA synthetases. These domains are also present in the aspartyl- and asparaginyl-tRNA synthetases, further confirming the notion that all three present-day enzymes originated from a common ancestral gene. The most conserved domain, located near the carboxyl terminal ends of this group of synthetases is characterized by a cluster of glycines and is also highly homologous to the carboxyl-terminal region of the E. coli ammonia-dependent asparagine synthetase. A catalytic function of the carboxyl terminal domain is indicated by in vitro mutagenesis of the yeast mitochondrial lysyl-tRNA synthetase. Replacement of any one of three glycine residues by alanine and in one case by aspartic acid completely suppresses the activity of the enzymes, as evidenced by the inability of the mutant genes to complement an msk1 mutant, even when present in high copy. Other mutations result in partial loss of activity. Only one glycine replacement affects the stability of the protein in vivo. The observed presence of a homologous domain in asparagine synthetase, which, like the aminoacyl-tRNA synthetases, catalyzes the formation of an aminoacyladenylate, suggests that the glycine-rich sequence is part of a catalytic site involved in binding of ATP and of the aminoacyladenylate intermediate.     

Ligand binding and structural changes associated with allostery in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an octameric enzyme composed of four each of regulatory IDH1 and catalytic IDH2 subunits that share 42% sequence identity. IDH2 contains catalytic isocitrate/Mg2+ and NAD+ binding sites whereas IDH1 contains homologous binding sites, respectively, for cooperative binding of isocitrate and for allosteric binding of AMP. Ligand binding is highly ordered in vitro, and IDH exhibits the unusual property of half-site binding for all ligands. The structures of IDH solved in the absence or presence of ligands have shown: (a) a heterodimer to be the basic structural/functional unit of the enzyme, (b) the organization of heterodimers to form tetramer and octamer structures, (c) structural differences that may underlie cooperative and allosteric regulatory mechanisms, and (d) the possibility for formation of a disulfide bond that could reduce catalytic activity. In vivo analyses of mutant enzymes have elucidated the physiological importance of catalytic activity and allosteric regulation of this tricarboxylic acid cycle enzyme. Other studies have established the importance of a disulfide bond in regulation of IDH activity in vivo, as well as contributions of this bond to the property of half-site ligand binding exhibited by the wild-type enzyme.     

Nucleotide and deduced amino acid sequence of the Rhodobacter sphaeroides gene encoding form II ribulose-1,5-bisphosphate carboxylase/oxygenase and comparison with other deduced form I and II sequences

The nucleotide sequence for the Rhodobacter sphaeroides form II ribulose 1,5-bisphosphate carboxylase/oxygenase was determined. The deduced product is highly homologous with the form II-like enzyme of Rhodospirillum rubrum , but appears to be more distantly related to the large subunit of the L₈S₈ enzyme found in autotrophic bacteria, cyanobacteria and higher plants. Several regions highly conserved among L₈S₈ and L_X enzymes correspond with regions previously implicated in catalytic activity and subunit interactions. An imperfect palindrome and a stem loop structure were identified in the 5' and 3' flanking sequences, respectively, of R. sphaeroides rbpL .     

Studies into factors contributing to substrate specificity of membrane-bound 3-ketoacyl-CoA synthases.

We are interested in constructing a model for the substrate-binding site of fatty acid elongase-1 3-ketoacyl CoA synthase (FAE1 KCS), the enzyme responsible for production of very long chain fatty acids of plant seed oils. Arabidopsis thaliana and Brassica napus FAE1 KCS enzymes are highly homologous but the seed oil content of these plants suggests that their substrate specificities differ with respect to acyl chain length. We used in vivo and in vitro assays of Saccharomyces cerevisiae-expressed FAE1 KCSs to demonstrate that the B. napus FAE1 KCS enzyme favors longer chain acyl substrates than the A. thaliana enzyme. Domains/residues responsible for substrate specificity were investigated by determining catalytic activity and substrate specificity of chimeric enzymes of A. thaliana and B. napus FAE1 KCS. The N-terminal region, excluding the transmembrane domain, was shown to be involved in substrate specificity. One chimeric enzyme that included A. thaliana sequence from the N terminus to residue 114 and B. napus sequence from residue 115 to the C terminus had substrate specificity similar to that of A. thaliana FAE1 KCS. However, a K92R substitution in this chimeric enzyme changed the specificity to that of the B. napus enzyme without loss of catalytic activity. Thus, this study was successful in identifying a domain involved in determining substrate specificity in FAE1 KCS and in engineering an enzyme with novel activity.     

Effects on specific antibodies on the catalytic activity of L-asparaginase from Serratia marcescens and Escherichia coli.

Rabbit antisera against highly purified L-asparaginase from Serratia marcescens and from Escherichia coli showed up to 60% inhibition of the catalytic amidohydrolysis of L-asparagine when combined with the homologous enzyme. This inhibition was diminished somewhat against the heterologous enzyme. Kinetic studies in the presence of these antisera showed an increased Kmapp for both homologous and heterologous enzymes using L-asparagine as substrate. In contrast, kinetic studies employing the poor substrate, L-glutamine, showed activation attributable to specific antibodies. This was seen in lower Kmapp values and up to twofold increases in the Vmax over the normal rabbit serum controls. The high degree of cross-inhibition (approximately 80%) and the low degree of cross-reactivity in the quantitative precipitin test (approximately 34%) suggest that these two enzymes possess structural similarities located mainly in the regions of the catalytic sites.     

Yeast mitochondrial and cytoplasmic valyl-tRNA synthetases

Yeast mitochondrial tRNA synthetase has been partially purified and chromatographic, catalytic and antigenic properties have been compared to the cytoplasmic homologous enzyme from yeast. No significant differences could be observed between the two enzymes with respect to their behaviour during ammonium sulfate precipitation or in chromatographic separation on DEAE cellulose, hydroxylapatite and Sephadex G 200. The Km of the two enzymes toward tRNAs from yeast mitochondria, yeast cytoplasm or E. coli are pratically identical. The antigenic properties of the two enzymes are very similar; antisera against either the mitochondria or the cytoplasmic enzyme lead to the inhibition of their catalytic properties. The mitochondrial ValRS is formed by a single polypeptide chain whose molecular weight is 125,000 daltons, a value very close to that of the yeast cytoplasmic enzyme.     

A highly specific monomeric isocitrate dehydrogenase from Corynebacterium glutamicum

The monomeric isocitrate dehydrogenase (IDH) of Corynebacterium glutamicum is compared to the topologically distinct dimeric IDH of Escherichia coli. Both IDHs have evolved to efficiently catalyze identical reactions with similar pH optimum as well as striking specificity toward NADP and isocitrate. However, the monomeric IDH is 10-fold more active (calculated as kcat/Km.isocitrate/Km.NADP) and 7-fold more NADP-specific than the dimeric enzyme, favoring NADP over NAD by a factor of 50,000. Such an extraordinary coenzyme specificity is not rivaled by any other characterized dehydrogenases. In addition, the monomeric enzyme is 10-fold more specific for isocitrate. The spectacular substrate specificity may be predominantly attributed to the isocitrate-assisted stabilization of catalytic complex during hydride transfer. No significant overall sequence identity is found between the monomeric and dimeric enzymes. However, structure-based alignment leads to the identification of three regions in the monomeric enzyme that match closely the three motifs located in the central region of dimeric IDHs and the homologous isopropylmalate dehydrogenases. The role of Lys253 as catalytic residue has been demonstrated by site-directed mutagenesis. Our results suggest that monomeric and dimeric forms of IDHs are functionally and structurally homologous.     

Prediction of catalytic residues in enzymes based on known tertiary structure,stability profile,and sequence conservation

The catalytic or functionally important residues of a protein are known to exist in evolutionarily constrained regions. However, the patterns of residue conservation alone are sometimes not very informative, depending on the homologous sequences available for a given query protein. Here, we present an integrated method to locate the catalytic residues in an enzyme from its sequence and structure. Mutations of functional residues usually decrease the activity, but concurrently often increase stability. Also, catalytic residues tend to occupy partially buried sites in holes or clefts on the molecular surface. After confirming these general tendencies by carrying out statistical analyses on 49 representative enzymes, these data together with amino acid conservation were evaluated. This novel method exhibited better sensitivity in the prediction accuracy than traditional methods that consider only the residue conservation. We applied it to some so-called "hypothetical" proteins, with known structures but undefined functions. The relationships among the catalytic, conserved, and destabilizing residues in enzymatic proteins are discussed.     

Genome-wide enzyme annotation with precision control: catalytic families (CatFam) databases

In this article, we present a new method termed CatFam (Catalytic Families) to automatically infer the functions of catalytic proteins, which account for 20-40% of all proteins in living organisms and play a critical role in a variety of biological processes. CatFam is a sequence-based method that generates sequence profiles to represent and infer protein catalytic functions. CatFam generates profiles through a stepwise procedure that carefully controls profile quality and employs nonenzymes as negative samples to establish profile-specific thresholds associated with a predefined nominal false-positive rate (FPR) of predictions. The adjustable FPR allows for fine precision control of each profile and enables the generation of profile databases that meet different needs: function annotation with high precision and hypothesis generation with moderate precision but better recall. Multiple tests of CatFam databases (generated with distinct nominal FPRs) against enzyme and nonenzyme datasets show that the method's predictions have consistently high precision and recall. For example, a 1% FPR database predicts protein catalytic functions for a dataset of enzymes and nonenzymes with 98.6% precision and 95.0% recall. Comparisons of CatFam databases against other established profile-based methods for the functional annotation of 13 bacterial genomes indicate that CatFam consistently achieves higher precision and (in most cases) higher recall, and that (on average) CatFam provides 21.9% additional catalytic functions not inferred by the other similarly reliable methods. These results strongly suggest that the proposed method provides a valuable contribution to the automated prediction of protein catalytic functions. The CatFam databases and the database search program are freely available at http://www.bhsai.org/downloads/catfam.tar.gz.     

L-Glutamate dehydrogenase from the antarctic fish Chaenocephalus aceratus. Primary structure, function and thermodynamic characterisation: relationship with cold adaptation

In order to study the molecular mechanisms of enzyme cold adaptation, direct amino acid sequence, catalytic features, thermal stability and thermodynamics of the reaction and of heat inactivation of L-glutamate dehydrogenase (GDH) from the liver of the Antarctic fish Chaenocephalus aceratus (suborder Notothenioidei, family Channichthyidae) were investigated. The enzyme shows dual coenzyme specificity, is inhibited by GTP and the forward reaction is activated by ADP and ATP. The complete primary structure of C. aceratus GDH has been established; it is the first amino acid sequence of a fish GDH to be described. In comparison with homologous mesophilic enzymes, the amino acid substitutions suggest a less compact molecular structure with a reduced number of salt bridges. Functional characterisation indicates efficient compensation of Q(10), achieved by increased k(cat) and modulation of S(0.5), which produce a catalytic efficiency at low temperature very similar to that of bovine GDH at its physiological temperature. The structural and functional characteristics are indicative of a high extent of protein flexibility. This property seems to find correspondence in the heat inactivation of Antarctic and bovine enzymes, which are inactivated at very similar temperature, but with different thermodynamics.     

Crystal structure of a thwarted mismatch glycosylase DNA repair complex

The bacterial mismatch-specific uracil-DNA glycosylase (MUG) and eukaryotic thymine-DNA glycosylase (TDG) enzymes form a homologous family of DNA glycosylases that initiate base-excision repair of G:U/T mismatches. Despite low sequence homology, the MUG/TDG enzymes are structurally related to the uracil-DNA glycosylase enzymes, but have a very different mechanism for substrate recognition. We have now determined the crystal structure of the Escherichia coli MUG enzyme complexed with an oligonucleotide containing a non-hydrolysable deoxyuridine analogue mismatched with guanine, providing the first structure of an intact substrate-nucleotide productively bound to a hydrolytic DNA glycosylase. The structure of this complex explains the preference for G:U over G:T mispairs, and reveals an essentially non-specific pyrimidine-binding pocket that allows MUG/TDG enzymes to excise the alkylated base, 3, N(4)-ethenocytosine. Together with structures for the free enzyme and for an abasic-DNA product complex, the MUG-substrate analogue complex reveals the conformational changes accompanying the catalytic cycle of substrate binding, base excision and product release.