Understanding enzyme superfamilies. Chemistry As the fundamental determinant in the evolution of new catalytic activities

Biosynthesis of the antitumor agent chartreusin involves the oxidative rearrangement of an anthracyclic polyketide

Chartreusin is a potent antitumor agent with a mixed polyketide-carbohydrate structure produced by Streptomyces chartreusis. Three type II polyketide synthase (PKS) gene clusters were identified from an S. chartreusis HKI-249 genomic cosmid library, one of which encodes chartreusin (cha) biosynthesis, as confirmed by heterologous expression of the entire cha gene cluster in Streptomyces albus. Molecular analysis of the approximately 37 kb locus and structure elucidation of a linear pathway intermediate from an engineered mutant reveal that the unusual bis-lactone aglycone chartarin is derived from an anthracycline-type polyketide. A revised biosynthetic model involving an oxidative rearrangement is presented.

Role of Protein Dynamics in Reaction Rate Enhancement by Enzymes

An integrated view of protein structure, dynamics, and function is emerging, where proteins are considered as dynamically active assemblies and internal motions are closely linked to function such as enzyme catalysis. Further, the motion of solvent bound to external regions of protein impacts internal motions and, therefore, protein function. Recently, we discovered a network of protein vibrations in enzyme cyclophilin A, coupled to its catalytic activity of peptidyl-prolyl cis-trans isomerization. Detailed studies suggest that this network, extending from surface regions to active site, is a conserved part of enzyme structure and has a role in promoting catalysis. In this report, theoretical investigations of concerted conformational fluctuations occurring on microsecond and longer time scales within the discovered network are presented. Using a new technique, kinetic energy was added to protein vibrational modes corresponding to conformational fluctuations in the network. The results reveal that protein dynamics promotes catalysis by altering transition state barrier crossing behavior of reaction trajectories. An increase in transmission coefficient and number of productive trajectories with increasing amounts of kinetic energy in vibrational modes is observed. Variations in active site enzyme-substrate interactions near transition state are found to be correlated with barrier recrossings. Simulations also showed that energy transferred from first solvation shell to surface residues impacts catalysis through network fluctuations. The detailed characterization of network presented here indicates that protein dynamics plays a role in rate enhancement by enzymes. Therefore, coupled networks in enzymes have wide implications in understanding allostericity and cooperative effects, as well as protein engineering and drug design.

Reconstitution of a defunct glycolytic pathway via recruitment of ambiguous sugar kinases

During a recent investigation of the persistence of substrate ambiguity in contemporary enzymes, we identified three distinct ambiguous sugar kinases embedded within the modern Escherichia coli genome [Miller, B. G., and Raines, R. T. (2004) Biochemistry 43, 6387-6392]. These catalysts are the YajF, YcfX, and NanK polypeptides, all of which possess rudimentary glucokinase activities. Here, we report on the discovery of a fourth bacterial kinase with ambiguous substrate specificity. AlsK phosphorylates the glucose epimer, d-allose, with a k(cat)/K(m) value of 6.5 x 10(4) M(-)(1) s(-)(1). AlsK also phosphorylates d-glucose, with a k(cat)/K(m) value that is 10(5)-fold lower than the k(cat)/K(m) value displayed by native E. coli glucokinase. Overexpression of the alsK gene relieves the auxotrophy of a glucokinase-deficient bacterium, demonstrating that weak enzymatic activities derived from ambiguous catalysts can provide organisms with elaborated metabolic capacities. To explore how ambiguous catalysts are recruited to provide new functions, we placed the glucokinase-deficient bacterium under selection for growth at the expense of glucose. Under these conditions, the bacterium acquires a spontaneous mutation in the putative promoter region of the yajF gene, a locus previously shown to encode a sugar kinase with relaxed substrate specificity. The point mutation regenerates a consensus sigma(70) promoter sequence that leads to a 94-fold increase in the level of yajF expression. This increase provides sufficient glucokinase activity for reconstitution of the defunct glycolytic pathway of the bacterial auxotroph. Our current findings indicate that ambiguous enzymatic activities continue to play an important role in the evolution of new metabolic pathways, and provide insight into the molecular mechanisms that facilitate the recruitment of such catalysts during periods of natural selection.

Molecular machines for protein degradation

One of the most precisely regulated processes in living cells is intracellular protein degradation. The main component of the degradation machinery is the 20S proteasome present in both eukaryotes and prokaryotes. In addition, there exist other proteasome-related protein-degradation machineries, like HslVU in eubacteria. Peptides generated by proteasomes and related systems can be used by the cell, for example, for antigen presentation. However, most of the peptides must be degraded to single amino acids, which are further used in cell metabolism and for the synthesis of new proteins. Tricorn protease and its interacting factors are working downstream of the proteasome and process the peptides into amino acids. Here, we summarise the current state of knowledge about protein-degradation systems, focusing in particular on the proteasome, HslVU, Tricorn protease and its interacting factors and DegP. The structural information about these protein complexes opens new possibilities for identifying, characterising and elucidating the mode of action of natural and synthetic inhibitors, which affects their function. Some of these compounds may find therapeutic applications in contemporary medicine.

In vitro selection and characterization of a stable subdomain of phosphoribosylanthranilate isomerase

The (beta(alpha))8-barrel is the most common enzyme fold and it is capable of catalyzing an enormous diversity of reactions. It follows that this scaffold should be an ideal starting point for engineering novel enzymes by directed evolution. However, experiments to date have utilized in vivo screens or selections and the compatibility of (beta(alpha))8-barrels with in vitro selection methods remains largely untested. We have investigated plasmid display as a suitable in vitro format by engineering a variant of phosphoribosylanthranilate isomerase (PRAI) that carried the FLAG epitope in active-site-forming loop 6. Trial enrichments for binding of mAb M2 (a mAb to FLAG) demonstrated that FLAG-PRAI could be identified from a 10(6)-fold excess of a FLAG-negative competitor in three rounds of in vitro selection. These results suggest PRAI as a useful scaffold for epitope and peptide grafting experiments. Further, we constructed a FLAG-PRAI loop library of approximately 10(7) clones, in which the epitope residues most critical for binding mAb M2 were randomized. Four rounds of selection for antibody binding identified and enriched for a variant in which a single nucleotide insertion produced a truncated (beta(alpha))8-barrel consisting of (beta(alpha))1-5beta6. Biophysical characterization of this clone, trPRAI, demonstrated that it was selected because of a 21-fold increase in mAb M2 affinity compared with full-length FLAG-PRAI. Remarkably, this truncated barrel was found to be soluble, structured, thermostable and monomeric, implying that it represents a genuine subdomain of PRAI and providing further evidence that such subdomains have played an important role in the evolution of the (beta(alpha))8-barrel fold.

Perturbing the Hydrophobic Pocket of Mandelate Racemase To Probe Phenyl Motion during Catalysis

Mandelate racemase (MR, EC 5.1.2.2) from Pseudomonas putida catalyzes the Mg(2+)-dependent 1,1-proton transfer that interconverts the enantiomers of mandelate. Crystal structures of MR reveal that the phenyl group of all ground-state ligands is located within a hydrophobic cavity, remote from the site of proton abstraction. MR forms numerous electrostatic and H-bonding interactions with the alpha-OH and carboxyl groups of the substrate, suggesting that these polar groups may remain relatively fixed in position during catalysis while the phenyl group is free to move between two binding sites [i.e., the R-pocket and the S-pocket for binding the phenyl group of (R)-mandelate and (S)-mandelate, respectively]. We show that MR binds benzilate (K(i) = 0.67 +/- 0.12 mM) and (S)-cyclohexylphenylglycolate (K(i) = 0.50 +/- 0.03 mM) as competitive inhibitors with affinities similar to that which the enzyme exhibits for the substrate. Therefore, the active site can simultaneously accommodate two phenyl groups, consistent with the existence of an R-pocket and an S-pocket. Wild-type MR exhibits a slightly higher affinity for (S)-mandelate [i.e., K(m)(S)(-)(man) < K(m)(R)(-)(man)] but catalyzes the turnover of (R)-mandelate slightly more rapidly (i.e., k(cat)(R)(-->)(S) > k(cat)(S)(-->)(R)). Upon introduction of steric bulk into the S-pocket using site-directed mutagenesis (i.e., the F52W, Y54W, and F52W/Y54W mutants), this catalytic preference is reversed. Although the catalytic efficiency (k(cat)/K(m)) of all the mutants was reduced (11-280-fold), all mutants exhibited a higher affinity for (R)-mandelate than for (S)-mandelate, and higher turnover numbers with (S)-mandelate as the substrate, relative to those with (R)-mandelate. (R)- and (S)-2-hydroxybutyrate are expected to be less sensitive to the additional steric bulk in the S-pocket. Unlike those for mandelate, the relative binding affinities for these substrate analogues are not reversed. These results are consistent with steric obstruction in the S-pocket and support the hypothesis that the phenyl group of the substrate may move between an R-pocket and an S-pocket during racemization. These conclusions were also supported by modeling of the binary complexes of the wild-type and F52W/Y54W enzymes with the substrate analogues (R)- and (S)-atrolactate, and of wild-type MR with bound benzilate using molecular dynamics simulations.

Functional analysis of the methylmalonyl-CoA epimerase from Caenorhabditis elegans

Methylmalonyl-CoA epimerase (MCE) is an enzyme involved in the propionyl-CoA metabolism that is responsible for the degradation of branched amino acids and odd-chain fatty acids. This pathway typically functions in the reversible conversion of propionyl-CoA to succinyl-CoA. The Caenorhabditis elegans genome contains a single gene encoding MCE (mce-1) corresponding to a 15 kDa protein. This was expressed in Escherichia coli and the enzymatic activity was determined. Analysis of the protein expression pattern at both the tissue and subcellular level by microinjection of green fluorescent protein constructs revealed expression in the pharynx, hypodermis and, most prominently in body wall muscles. The subcellular pattern agrees with predictions of mitochondrial localization. The sequence similarity to an MCE of known structure was high enough to permit a three-dimensional model to be built, suggesting conservation of ligand and metal binding sites. Comparison with corresponding sequences from a variety of organisms shows more than 1/6 of the sequence is completely conserved. Mutants allelic to mce-1 showed no obvious phenotypic alterations, demonstrating that the enzyme is not essential for normal worm development under laboratory conditions. However, survival of the knockout mutants was altered when exposed to stress conditions, with mutants surprisingly showing an increased resistance to oxidative stress.

From cyclohydrolase to oxidoreductase: discovery of nitrile reductase activity in a common fold

The enzyme YkvM from Bacillus subtilis was identified previously along with three other enzymes (YkvJKL) in a bioinformatics search for enzymes involved in the biosynthesis of queuosine, a 7-deazaguanine modified nucleoside found in tRNA(GUN) of Bacteria and Eukarya. Genetic analysis of ykvJKLM mutants in Acinetobacter confirmed that each was essential for queuosine biosynthesis, and the genes were renamed queCDEF. QueF exhibits significant homology to the type I GTP cyclohydrolases characterized by FolE. Given that GTP is the precursor to queuosine and that a cyclohydrolase-like reaction was postulated as the initial step in queuosine biosynthesis, QueF was proposed to be the putative cyclohydrolase-like enzyme responsible for this reaction. We have cloned the queF genes from B. subtilis and Escherichia coli and characterized the recombinant enzymes. Contrary to the predictions based on sequence analysis, we discovered that the enzymes, in fact, catalyze a mechanistically unrelated reaction, the NADPH-dependent reduction of 7-cyano-7-deazaguanineto7-aminomethyl-7-deazaguanine, a late step in the biosynthesis of queuosine. We report here in vitro and in vivo studies that demonstrate this catalytic activity, as well as preliminary biochemical and bioinformatics analysis that provide insight into the structure of this family of enzymes.

Functional analysis of active site residues of the fosfomycin resistance enzyme FosA from Pseudomonas aeruginosa

The metalloglutathione transferase FosA catalyzes the conjugation of glutathione to carbon-1 of the antibiotic fosfomycin, rendering it ineffective as an antibacterial drug. Codon randomization and selection for the ability of resulting clones to confer fosfomycin resistance to Escherichia coli were used to identify residues critical for FosA function. Of the 24 codons chosen for randomization, 16 were found to be essential because only the wild type amino acid was selected. These included ligands to the Mn(2+) and the K(+), residues that furnish hydrogen bonds to fosfomycin, and residues located in a putative glutathione/fosfomycin-binding site. The remaining eight positions randomized were tolerant to substitutions. Site-directed mutagenesis of some of the essential and tolerant amino acids to alanine was performed, and the activity of the purified proteins was determined. Mutation of the residues that are within hydrogen bonding distance to the oxirane or phosphonate oxygens of fosfomycin resulted in variants with very low or no activity. Mutation of Ser(94), which bridges one of the phosphonate oxygens with a potassium ion, resulted in insoluble protein. The Y39A mutation in the putative glutathione-binding site resulted in a 4-fold increase in the apparent K(m) for glutathione. Only two of the amino acids in the substrate-binding site are conserved in the related fosfomycin resistance proteins FosB and FosX, whereas no amino acids in the putative glutathione-binding site are conserved.

Engineering glutathione transferase to a novel glutathione peroxidase mimic with high catalytic efficiency. Incorporation of selenocysteine into a glutathione-binding scaffold using an auxotrophic expression system

Glutathione peroxidase (GPx, EC 1.11.1.9) protects cells against oxidative damage by catalyzing the reduction of hydroperoxides with glutathione (GSH). Several attempts have been made to imitate its function for mechanical study and for its pharmacological development as an antioxidant. By replacing the active site serine 9 with a cysteine and then substituting it with selenocysteine in a cysteine auxotrophic system, catalytically essential residue selenocysteine was bioincorporated into GSH-specific binding scaffold, and thus, glutathione S-transferase (GST, EC 2.5.1.18) from Lucilia cuprina was converted into a selenium-containing enzyme, seleno-LuGST1-1, by genetic engineering. Taking advantage of the important structure similarities between seleno-LuGST1-1 and naturally occurring GPx in the specific GSH binding sites and the geometric conformation for the active selenocysteine in their common GSH binding domain-adopted thioredoxin fold, the as-generated selenoenzyme displayed a significantly high efficiency for catalyzing the reduction of hydrogen peroxide by glutathione, being comparable with those of natural GPxs. The catalytic behaviors of this engineered selenoenzyme were found to be similar to those of naturally occurring GPx. It exhibited pH and temperature-dependent catalytic activity and a typical ping-pong kinetic mechanism. Engineering GST into an efficient GPx-like biocatalyst provided new proof for the previous assumption that both GPx and GST were evolved from a common thioredoxin-like ancestor to accommodate different functions throughout evolution.

Divergence of function in the thioredoxin fold suprafamily: evidence for evolution of peroxiredoxins from a thioredoxin-like ancestor

The thioredoxin fold is found in proteins that serve a wide variety of functions. Among these are peroxiredoxins, which catalyze the reduction of hydrogen peroxide and alkyl peroxides. Although the common structural fold shared by thioredoxins and peroxiredoxins suggests the possibility that they have evolved from a common progenitor, it has been difficult to examine this hypothesis in depth because pairwise sequence identities between proteins in these two superfamilies are statistically insignificant. Using the Shotgun program, we have found that sequences of reductases involved in maturation of cytochromes in certain bacteria bridge the sequences of thioredoxins and peroxiredoxins. Analysis of motifs found in a divergent set of thioredoxins, cytochrome maturation proteins, and peroxiredoxins provides further support for an evolutionary relationship between these proteins. Within the conserved motifs are specific residues that are characteristic of individual protein classes, and therefore are likely to be involved in the specific functions of those classes. We have used this information, in combination with existing structural and functional information, to gain new insight into the structure-function relationships in these proteins and to construct a model for the emergence of peroxiredoxins from a thioredoxin-like ancestor.

Divergent evolution in the enolase superfamily: the interplay of mechanism and specificity

The members of the mechanistically diverse enolase superfamily catalyze different overall reactions. Each shares a partial reaction in which an active site base abstracts the alpha-proton of the carboxylate substrate to generate an enolate anion intermediate that is stabilized by coordination to the essential Mg(2+) ion; the intermediates are then directed to different products in the different active sites. In this minireview, our current understanding of structure/function relationships in the divergent members of the superfamily is reviewed, and the use of this knowledge for our future studies is proposed.

Mucosal human papillomaviruses encode four different E5 proteins whose chemistry and phylogeny correlate with malignant or benign growth

We performed a phylogenetic study of the E2-L2 region of human mucosal papillomaviruses (PVs) and of the proteins therein encoded. Hitherto, proteins codified in this region were known as E5 proteins. We show that many of these proteins could be spurious translations, according to phylogenetic and chemical coherence criteria between similar protein sequences. We show that there are four separate families of E5 proteins, with different characteristics of phylogeny, chemistry, and rate of evolution. For the sake of clarity, we propose a change in the present nomenclature. E5alpha is present in groups A5, A6, A7, A9, and A11, PVs highly associated with malignant carcinomas of the cervix and penis. E5beta is present in groups A2, A3, A4, and A12, i.e., viruses associated with certain warts. E5gamma is present in group A10, and E5delta is encoded in groups A1, A8, and A10, which are associated with benign transformations. The phylogenetic relationships between mucosal human PVs are the same when considering the oncoproteins E6 and E7 and the E5 proteins and differ from the phylogeny estimated for the structural proteins L1 and L2. Besides, the protein divergence rate is higher in early proteins than in late proteins, increasing in the order L1 < L2 < E6 approximately E7 < E5. Moreover, the same proteins have diverged more rapidly in viruses associated with malignant transformations than in viruses associated with benign transformations. The E5 proteins display, therefore, evolutionary characteristics similar to those of the E6 and E7 oncoproteins. This could reflect a differential involvement of the E5 types in the transformation processes.

The 1.3 A crystal structure of human mitochondrial Delta3-Delta2-enoyl-CoA isomerase shows a novel mode of binding for the fatty acyl group

The crystal structure of Delta3-Delta2-enoyl-CoA isomerase from human mitochondria (hmEci), complexed with the substrate analogue octanoyl-CoA, has been refined at 1.3 A resolution. This enzyme takes part in the beta-oxidation of unsaturated fatty acids by converting both cis-3 and trans-3-enoyl-CoA esters (with variable length of the acyl group) to trans-2-enoyl-CoA. hmEci belongs to the hydratase/isomerase (crotonase) superfamily. Most of the enzymes belonging to this superfamily are hexamers, but hmEci is shown to be a trimer. The mode of binding of the ligand, octanoyl-CoA, shows that the omega-end of the acyl group binds in a hydrophobic tunnel formed by residues of the loop preceding helix H4 as well as by side-chains of the kinked helix H9. From the structure of the complex it can be seen that Glu136 is the only catalytic residue. The importance of Glu136 for catalysis is confirmed by mutagenesis studies. A cavity analysis shows the presence of two large, adjacent empty hydrophobic cavities near the active site, which are shaped by side-chains of helices H1, H2, H3 and H4. The structure comparison of hmEci with structures of other superfamily members, in particular of rat mitochondrial hydratase (crotonase) and yeast peroxisomal enoyl-CoA isomerase, highlights the variable mode of binding of the fatty acid moiety in this superfamily.

Evolution of function in the crotonase superfamily: (3S)-methylglutaconyl-CoA hydratase from Pseudomonas putida

Members of the enoyl-CoA hydratase (crotonase) superfamily catalyze different overall reactions that utilize a common catalytic strategy delivered by a shared structural scaffold; the substrates are usually acyl esters of coenzyme A, and the intermediates are usually thioester enolate anions stabilized by a conserved oxyanion hole. In many bacterial genomes, orthologous members that contain homologues of acid/base catalyst Glu164 but not of Glu144 in rat mitochondrial crotonase are encoded by operons of which the functions have not been assigned. Focusing on the orthologues from Pseudomonas aeruginosa and P. putida, we have determined that these operons encode enzymes in leucine catabolism with the unknown enzyme assigned as (3S)-methylglutaconyl-CoA hydratase (MGCH), which catalyzes the syn-hydration of (E)-3-methylglutaconyl-CoA to (3S)-hydroxymethylglutaryl-CoA. The discovery that bacterial MGCHs catalyze hydration of enoyl-CoAs utilizing a single active-site residue contrasts with the paradigm crotonases as well as with the recently identified mammalian MGCHs that use homologues of both Glu144 and Glu164 in crotonase. Substrate analogues lacking a gamma-carboxylate have been shown to be competitive inhibitors of the enzyme, and installation of a glutamate for the "missing" homologue of Glu144 fails to introduce hydratase activity with the substrate analogues. Thus, bacterial MGCHs may provide an example of opportunistic evolution in which a carboxylate group of the substrate functionally replaces one of the active site glutamate residues in the reactions catalyzed by crotonases and the eukaryotic MGCHs.

Deciphering enzymes. Genetic selection as a probe of structure and mechanism

The efficient engineering of enzymes with novel activities remains an ongoing challenge. Towards this end, genetic selection techniques provide a method for finding rare solutions to catalytic problems that requires only a limited foreknowledge of structure-function relationships. We have used genetic selections to extensively probe the structure and mechanism of chorismate mutases. The insights gained from these investigations will aid future enzyme design efforts.

Finding evolutionary relations beyond superfamilies: fold-based superfamilies

Superfamily classifications are based variably on similarity of sequences, global folds, local structures, or functions. We have examined the possibility of defining superfamilies purely from the viewpoint of the global fold/function relationship. For this purpose, we first classified protein domains according to the beta-sheet topology. We then introduced the concept of kinship relations among the classified beta-sheet topology by assuming that the major elementary event leading to creation of a new beta-sheet topology is either an addition or deletion of one beta-strand at the edge of an existing beta-sheet during the molecular evolution. Based on this kinship relation, a network of protein domains was constructed so that the distance between a pair of domains represents the number of evolutionary events that lead one from the other domain. We then mapped on it all known domains with a specific core chemical function (here taken, as an example, that involving ATP or its analogs). Careful analyses revealed that the domains are found distributed on the network as >20 mutually disjointed clusters. The proteins in each cluster are defined to form a fold-based superfamily. The results indicate that >20 ATP-binding protein superfamilies have been invented independently in the process of molecular evolution, and the conservative evolutionary diffusion of global folds and functions is the origin of the relationship between them.

Directed evolution of protein enzymes using nonhomologous random recombination

We recently reported the development of nonhomologous random recombination (NRR) as a method for nucleic acid diversification and applied NRR to the evolution of DNA aptamers. Here, we describe a modified method, protein NRR, that enables proteins to access diversity previously difficult or impossible to generate. We investigated the structural plasticity of protein folds and the ability of helical motifs to function in different contexts by applying protein NRR and in vivo selection to the evolution of chorismate mutase (CM) enzymes. Functional CM mutants evolved using protein NRR contained many insertions, deletions, and rearrangements. The distribution of these changes was not random but clustered in certain regions of the protein. Topologically rearranged but functional enzymes also emerged from these studies, indicating that multiple connectivities can accommodate a functional CM active site and demonstrating the ability to generate new domain connectivities through protein NRR. Protein NRR was also used to randomly recombine CM and fumarase, an unrelated but also alpha-helical protein. Whereas the resulting library contained fumarase fragments in many contexts before functional selection, library members surviving selection for CM activity invariably contained a CM core with fumarase sequences found only at the termini or in one loop. These results imply that internal helical fragments cannot be swapped between these proteins without the loss of nearly all CM activity. Our findings suggest that protein NRR will be useful in probing the functional requirements of enzymes and in the creation of new protein topologies.

Superfamily active site templates

We show that three-dimensional signatures consisting of only a few functionally important residues can be diagnostic of membership in superfamilies of enzymes. Using the enolase superfamily as a model system, we demonstrate that such a signature, or template, can identify superfamily members in structural databases with high sensitivity and specificity. This is remarkable because superfamilies can be highly diverse, with members catalyzing many different overall reactions; the unifying principle can be a conserved partial reaction or chemical capability. Our definition of a superfamily thus hinges on the disposition of residues involved in a conserved function, rather than on fold similarity alone. A clear advantage of basing structure searches on such active site templates rather than on fold similarity is the specificity with which superfamilies with distinct functional characteristics can be identified within a large set of proteins with the same fold, such as the (beta/alpha)8 barrels. Preliminary results are presented for an additional group of enzymes with a different fold, the haloacid dehalogenase superfamily, suggesting that this approach may be generally useful for assigning reading frames of unknown function to specific superfamilies and thereby allowing inference of some of their functional properties.

Understanding the importance of protein structure to nature's routes for divergent evolution in TIM barrel enzymes

It is widely agreed that new enzymes evolve from existing ones through the duplication of genes encoding existing enzymes followed by sequence divergence. While evolution is an inherently random process, studies of divergently related enzymes have shown that the evolution of new enzymes follows one of three general routes in which the substrate specificity, reaction mechanism, or active site architecture of the progenitor enzyme is reused in the new enzyme. Recent developments in structural biology relating to divergently related (beta/alpha)8 enzymes have brought new insight into these processes and have revealed that conserved structural elements play an important role in divergent evolution. These studies have shown that, although evolution occurs as a series of random mutations, stable folds such as the (beta/alpha)8 barrel and structural features of the active sites of enzymes are frequently reused in evolution and adapted for new catalytic purposes.

Hydrophobic Nature of the Active Site of Mandelate Racemase

Mandelate racemase (EC 5.1.2.2) from Pseudomonas putida catalyzes the interconversion of the two enantiomers of mandelic acid with remarkable proficiency, stabilizing the altered substrate in the transition state by approximately 26 kcal/mol. We have used a series of substrate analogues (glycolates) and intermediate analogues (hydroxamates) to evaluate the contribution of the hydrophobic cavity within the enzyme's active site to ligand binding. Free energy changes accompanying binding of glycolate derivatives correlated well with the hydrophobic substituent constant pi and the van der Waals surface areas of the ligands. The observed dependence of the apparent binding free energy on surface area of the ligand was -30 +/- 5 cal mol(-1) A(-2) at 25 degrees C. Free energy changes accompanying binding of hydroxamate derivatives also correlated well with pi values and the van der Waals surface areas of the ligands, giving a slightly greater free energy dependence equal to -41 +/- 3 cal mol(-1) A(-2) at 25 degrees C. Surprisingly, mandelate racemase exhibited a binding affinity for the intermediate analogue benzohydroxamate that was 2 orders of magnitude greater than that predicted solely on the basis of hydrophobic interactions. This suggests that there are additional specific interactions that stabilize the altered substrate in the transition state. Mandelate racemase was competitively inhibited by (R,S)-1-naphthylglycolate (apparent K(i) = 1.9 +/- 0.1 mM) and (R,S)-2-naphthylglycolate (apparent K(i) = 0.52 +/- 0.03 mM), demonstrating the plasticity of the hydrophobic pocket. Both (R)- (K(m) = 0.46 +/- 0.06 mM, k(cat) = 33 +/- 1 s(-1)) and (S)-2-naphthylglycolate (K(m) = 0.41 +/- 0.03 mM, k(cat) = 25 +/- 1 s(-1)) were shown to be alternative substrates for mandelate racemase. These kinetic results demonstrate that no major steric restrictions are imposed on the binding of this bulkier substrate in the ground state but that steric factors appear to impair transition state/intermediate stabilization. 2-Naphthohydroxamate was identified as a competitive inhibitor of mandelate racemase, binding with an affinity (K(i) = 57 +/- 18 microM) that was reduced relative to that observed for benzohydroxamate and that was in accord with the approximately 10-fold reduction in the value of k(cat)/K(m) for the racemization of 2-naphthylglycolate. These findings indicate that, for mandelate racemase, steric constraints within the hydrophobic cavity of the enzyme-intermediate complex are more stringent than those in the enzyme-substrate complex.

MosA, a protein implicated in rhizopine biosynthesis in Sinorhizobium meliloti L5-30, is a dihydrodipicolinate synthase

MosA is a gene product encoded on a pSym megaplasmid of Sinorhizobium meliloti L5-30. The gene is part of an operon reported to be essential for the synthesis of the rhizopine 3-O-methyl-scyllo-inosamine. MosA has been assigned the function of an O-methyltransferase. However, the reported sequence of this protein is very much like that of dihydrodipicolinate synthase (DHDPS), except for a 40 amino acid residue C-terminal domain. This similarity contradicts accepted ideas regarding structure-function relationships of enzymes. We have cloned and overexpressed the recombinant gene in Escherichia coli, and discovered that the reported sequence contains an error resulting in a frame-shift. The correct sequence contains a new stop codon, truncating the C-terminal 41 amino acid residues of the reported sequence. The expressed protein, bearing an N-terminal polyhistidine tag, catalyzes the condensation of pyruvate and aspartate beta-semialdehyde efficiently, suggesting that this activity is not a side-reaction, but an activity for which this enzyme has evolved. Electro-spray mass spectrometry experiments and inhibition by L-lysine are consistent with the enzyme being a DHDPS. E.coli AT997, a mutant host normally requiring exogenous diaminopimelate for growth, could be complemented by transformation with a plasmid bearing the gene encoding MosA. A role for this enzyme in rhizopine synthesis cannot be ruled out, but is called into question.

Reactions of 4-Oxalocrotonate Tautomerase and YwhB with 3-Halopropiolates: Analysis and Implications

4-Oxalocrotonate tautomerase (4-OT) and YwhB, a 4-OT homologue found in Bacillus subtilis, exhibit a low level hydratase activity that converts trans-3-haloacrylates to acetaldehyde, presumably through a malonate semialdehyde intermediate. The mechanism for the initial transformation of the 3-haloacrylate to malonate semialdehyde involves Pro-1 as well as an arginine, two residues that play critical roles in the 4-OT-catalyzed isomerization reaction and the YwhB-catalyzed tautomerization reaction. These residues are also critical for the trans-3-chloroacrylic acid dehalogenase (CaaD)-catalyzed conversion of trans-3-haloacrylates to malonate semialdehyde. Recently, 3-bromo- and 3-chloropropiolate, the acetylene analogues of 3-haloacrylates, were characterized as potent irreversible inhibitors of CaaD due to the covalent modification of the catalytic proline. In view of these observations, an investigation of the behavior of 4-OT and YwhB with the 3-halopropiolates was undertaken. The results show that these compounds are potent irreversible inhibitors of 4-OT and YwhB with Pro-1 being the sole site of covalent modification by 3-bromopropiolate. The inactivation process could involve the enzyme-catalyzed addition of water to the 3-halopropiolate yielding an acyl halide, which would inactivate the enzyme or be initiated by the nucleophilic attack of Pro-1 at the C-3 position of the 3-halopropiolate in a Michael type reaction. The presence of the halogen along with Arg-11 could facilitate both reactions with the latter causing the polarization of the alpha,beta-unsaturated acids. The 3-halopropiolates are the first identified inhibitors of YwhB and confirm the importance of Pro-1 in its mechanism. In addition, the results set the stage for the use of these compounds as mechanistic probes of the primary as well as low level activities of 4-OT and YwhB.

Steady-State Kinetics and Molecular Evolution of Escherichia coli MenD [(1R,6R)-2-Succinyl-6-hydroxy-2,4-cyclohexadiene-1-carboxylate Synthase], an Anomalous Thiamin Diphosphate-Dependent Decarboxylase−Carboligase

(1R,6R)-2-Succinyl-6-hydroxy-2,4-cyclohexadiene-1-carboxylate (SHCHC) synthase, or MenD, catalyzes the thiamin diphosphate- (ThDP-) dependent decarboxylation of 2-oxoglutarate, the subsequent addition of the resulting succinyl-ThDP moiety to isochorismate, and the delta-elimination of pyruvate to yield SHCHC, pyruvate, and carbon dioxide. The enzyme is part of a superfamily of ThDP-dependent 2-oxo acid decarboxylases that includes pyruvate decarboxylase, benzoylformate decarboxylase, and acetohydroxy acid synthase, among others. However, this is the only enzyme known to catalyze a Stetter-like 1,4-addition of a ThDP adduct to the beta-carbon of an unsaturated carboxylate. Herein we report properties of the MenD protein from Escherichia coli, including the results of the first steady-state kinetic studies of the SHCHC synthase reaction. The protein is a dimer and shows cooperativity with respect to both substrates. The enzyme prefers divalent manganese as its metal ion cofactor and shows no dependence on FAD. MenD, required for biosynthesis of menaquinone and phylloquinone, is found in the genomes of a wide range of bacteria, as well as that of the archaeon Halobacterium sp. NRC-1 and the eukaryote Arabidopsis thaliana. Sequence alignments with other members of the superfamily are used to predict amino acid residues likely to be important in the binding and activation of ThDP. A site-directed mutant that replaces the conserved glutamic acid residue (E55), predicted to interact with N1' of the aminopyrimidine ring, with glutamine was generated, with catastrophic results for catalysis. There is no evidence for the release of succinate semialdehyde as a product; therefore, EC 4.1.1.71 should not be used for this enzyme.

A perspective on enzyme catalysis

The seminal hypotheses proposed over the years for enzymatic catalysis are scrutinized. The historical record is explored from both biochemical and theoretical perspectives. Particular attention is given to the impact of molecular motions within the protein on the enzyme's catalytic properties. A case study for the enzyme dihydrofolate reductase provides evidence for coupled networks of predominantly conserved residues that influence the protein structure and motion. Such coupled networks have important implications for the origin and evolution of enzymes, as well as for protein engineering.

Catalysing new reactions during evolution: economy of residues and mechanism

The diversity of function in some enzyme superfamilies shows that during evolution, enzymes have evolved to catalyse different reactions on the same structure scaffold. In this analysis, we examine in detail how enzymes can modify their chemistry, through a comparison of the catalytic residues and mechanisms in 27 pairs of homologous enzymes of totally different functions. We find that evolution is very economical. Enzymes retain structurally conserved residues to aid catalysis, including residues that bind catalytic metal ions and modulate cofactor chemistry. We examine the conservation of residue type and residue function in these structurally conserved residue pairs. Additionally, enzymes often retain common mechanistic steps catalyzed by structurally conserved residues. We have examined these steps in the context of their overall reactions.

Transcarboxylase 12S crystal structure: hexamer assembly and substrate binding to a multienzyme core

Transcarboxylase from Propionibacterium shermanii is a 1.2 MDa multienzyme complex that couples two carboxylation reactions, transferring CO(2)(-) from methylmalonyl-CoA to pyruvate, yielding propionyl-CoA and oxaloacetate. The 1.9 A resolution crystal structure of the central 12S hexameric core, which catalyzes the first carboxylation reaction, has been solved bound to its substrate methylmalonyl-CoA. Overall, the structure reveals two stacked trimers related by 2-fold symmetry, and a domain duplication in the monomer. In the active site, the labile carboxylate group of methylmalonyl-CoA is stabilized by interaction with the N-termini of two alpha-helices. The 12S domains are structurally similar to the crotonase/isomerase superfamily, although only domain 1 of each 12S monomer binds ligand. The 12S reaction is similar to that of human propionyl-CoA carboxylase, whose beta-subunit has 50% sequence identity with 12S. A homology model of the propionyl-CoA carboxylase beta-subunit, based on this 12S crystal structure, provides new insight into the propionyl-CoA carboxylase mechanism, its oligomeric structure and the molecular basis of mutations responsible for enzyme deficiency in propionic acidemia.

Bioorganic chemistry à la baguette: studies on molecular recognition in biological systems using rigid-rod molecules

Initial studies using rigid-rod molecules or "baguettes" to address bioorganic topics of current scientific concern are reported. It is illustrated how transmembrane oligo(p-phenylene)s as representative model rods can be tuned to recognize lipid bilayer membranes either by their thickness or polarization. The construction of otherwise problematic hydrogen-bonded chains along transmembrane rods yields "proton wires," which act by a mechanism that is central in bioenergetics but poorly explored by means of synthetic models. Another example focuses on multivalent ligands assembling rigid-rod cell-surface receptors into transmembrane dynamic arene arrays. The potassium transport mediated by these ligand-receptor complexes provides experimental support for the potential biological importances of the controversial cation-pi mechanism. More complex supramolecular architecture is portrayed in the first artificial beta-barrels. It is shown how programmed assembly of toroidal rigid-rod supramolecules in detergent-free water permits control of diameter of the chemical nature of their interior. Reversed rigid-rod beta-barrels are assembled to function as self-assembled ionophores, ion channel models, and transmembrane nanopores. The potential of future intratoroidal chemistry is exemplified by encapsulation and planarization of beta-carotene in water and the construction of transmembrane B-DNA at the center of a second-sphere host-guest complex à al baguette.

The biochemistry of peroxisomal beta-oxidation in the yeast Saccharomyces cerevisiae

Peroxisomal fatty acid degradation in the yeast Saccharomyces cerevisiae requires an array of beta-oxidation enzyme activities as well as a set of auxiliary activities to provide the beta-oxidation machinery with the proper substrates. The corresponding classical and auxiliary enzymes of beta-oxidation have been completely characterized, many at the structural level with the identification of catalytic residues. Import of fatty acids from the growth medium involves passive diffusion in combination with an active, protein-mediated component that includes acyl-CoA ligases, illustrating the intimate linkage between fatty acid import and activation. The main factors involved in protein import into peroxisomes are also known, but only one peroxisomal metabolite transporter has been characterized in detail, Ant1p, which exchanges intraperoxisomal AMP with cytosolic ATP. The other known transporter is Pxa1p-Pxa2p, which bears similarity to the human adrenoleukodystrophy protein ALDP. The major players in the regulation of fatty acid-induced gene expression are Pip2p and Oaf1p, which unite to form a transcription factor that binds to oleate response elements in the promoter regions of genes encoding peroxisomal proteins. Adr1p, a transcription factor, binding upstream activating sequence 1, also regulates key genes involved in beta-oxidation. The development of new, postgenomic-era tools allows for the characterization of the entire transcriptome involved in beta-oxidation and will facilitate the identification of novel proteins as well as the characterization of protein families involved in this process.

Conservation of electrostatic properties within enzyme families and superfamilies

Electrostatic interactions play a key role in enzyme catalytic function. At long range, electrostatics steer the incoming ligand/substrate to the active site, and at short distances, electrostatics provide the specific local interactions for catalysis. In cases in which electrostatics determine enzyme function, orthologs should share the electrostatic properties to maintain function. Often, electrostatic potential maps are employed to depict how conserved surface electrostatics preserve function. We expand on previous efforts to explain conservation of function, using novel electrostatic sequence and structure analyses of four enzyme families and one enzyme superfamily. We show that the spatial charge distribution is conserved within each family and superfamily. Conversely, phylogenetic analysis of key electrostatic residues provide the evolutionary origins of functionality.

Functional analysis of BamHI DNA cytosine-N4 methyltransferase

We show that the kinetic mechanism of the DNA (cytosine-N(4)-)-methyltransferase M.BamHI, which modifies the underlined cytosine (GGATCC), differs from cytosine C(5) methyltransferases, and is similar to that observed with adenine N(6) methyltransferases. This suggests that the obligate order of ternary complex assembly and disassembly depends on the type of methylation reaction. In contrast, the single-turnover rate of catalysis for M.BamHI (0.10s(-1)) is closer to the DNA (cytosine-C(5)-)-methyltransferases (0.14s(-1)) than the DNA (adenine-N(6)-)-methyltransferases (>200s(-1)). The nucleotide flipping transition dominates the single-turnover constant for adenine N(6) methyltransferases, and, since the disruption of the guanine-cytosine base-pair is essential for both types of cytosine DNA methyltransferases, this transition may be a common, rate-limiting step for methylation for these two enzyme subclasses. The similar overall rate of catalysis by M.BamHI and other DNA methyltransferases is consistent with a common rate-limiting catalytic step of product dissociation. Our analyses of M.BamHI provide functional insights into the relationship between the three different classes of DNA methyltransferases that complement both prior structural and evolutionary insights.

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.

Organization of the multifunctional enzyme type 1: interaction between N- and C-terminal domains is required for the hydratase-1/isomerase activity

Rat peroxisomal multifunctional enzyme type 1 (perMFE-1) is a monomeric protein of beta-oxidation. We have defined five functional domains (A, B, C, D and E) in the perMFE-1 based on comparison of the amino acid sequence with homologous proteins from databases and structural data of the hydratase-1/isomerases (H1/I) and (3 S )-hydroxyacyl-CoA dehydrogenases (HAD). Domain A (residues 1-190) comprises the H1/I fold and catalyses both 2-enoyl-CoA hydratase-1 and Delta(3)-Delta(2)-enoyl-CoA isomerase reactions. Domain B (residues 191-280) links domain A to the (3 S )-dehydrogenase region, which includes both domain C (residues 281-474) and domain D (residues 480-583). Domains C and D carry features of the dinucleotide-binding and the dimerization domains of monofunctional HADs respectively. Domain E (residues 584-722) has sequence similarity to domain D of the perMFE-1, which suggests that it has evolved via partial gene duplication. Experiments with engineered perMFE-1 variants demonstrate that the H1/I competence of domain A requires stabilizing interactions with domains D and E. The variant His-perMFE (residues 288-479)Delta, in which the domain C is deleted, is stable and has hydratase-1 activity. It is proposed that the extreme C-terminal domain E in perMFE-1 serves the following three functions: (i) participation in the folding of the N-terminus into a functionally competent H1/I fold, (ii) stabilization of the dehydrogenation domains by interaction with the domain D and (iii) the targeting of the perMFE-1 to peroxisomes via its C-terminal tripeptide.

Sequence and structural differences between enzyme and nonenzyme homologs

To improve our understanding of the evolution of novel functions, we performed a sequence, structural, and functional analysis of homologous enzymes and nonenzymes of known three-dimensional structure. In most examples identified, the nonenzyme is derived from an ancestral catalytic precursor (as opposed to the reverse evolutionary scenario, nonenzyme to enzyme), and the active site pocket has been disrupted in some way, owing to the substitution of critical catalytic residues and/or steric interactions that impede substrate binding and catalysis. Pairwise sequence identity is typically insignificant, and almost one-half of the enzyme and nonenzyme pairs do not share any similarity in function. Heterooligomeric enzymes comprising homologous subunits in which one chain is catalytically inactive and enzyme polypeptides that contain internal catalytic and noncatalytic duplications of an ancient enzyme domain are also discussed.

Transmutation of human glutathione transferase A2-2 with peroxidase activity into an efficient steroid isomerase

A major goal in protein engineering is the tailor-making of enzymes for specified chemical reactions. Successful attempts have frequently been based on directed molecular evolution involving libraries of random mutants in which variants with desired properties were identified. For the engineering of enzymes with novel functions, it would be of great value if the necessary changes of the active site could be predicted and implemented. Such attempts based on the comparison of similar structures with different substrate selectivities have previously met with limited success. However, the present work shows that the knowledge-based redesign restricted to substrate-binding residues in human glutathione transferase A2-2 can introduce high steroid double-bond isomerase activity into the enzyme originally characterized by glutathione peroxidase activity. Both the catalytic center activity (k(cat)) and catalytic efficiency (k(cat)/K(m)) match the values of the naturally evolved glutathione transferase A3-3, the most active steroid isomerase known in human tissues. The substrate selectivity of the mutated glutathione transferase was changed 7000-fold by five point mutations. This example demonstrates the functional plasticity of the glutathione transferase scaffold as well as the potential of rational active-site directed mutagenesis as a complement to DNA shuffling and other stochastic methods for the redesign of proteins with novel functions.

Plasticity of enzyme active sites

The expectation is that any similarity in reaction chemistry shared by enzyme homologues is mediated by common functional groups conserved through evolution. However, detailed enzyme studies have revealed the flexibility of many active sites, in that different functional groups, unconserved with respect to position in the primary sequence, mediate the same mechanistic role. Nevertheless, the catalytic atoms might be spatially equivalent. More rarely, the active sites have completely different locations in the protein scaffold. This variability could result from: (1) the hopping of functional groups from one position to another to optimize catalysis; (2) the independent specialization of a low-activity primordial enzyme in different phylogenetic lineages; (3) functional convergence after evolutionary divergence; or (4) circular permutation events.

A multidomain TIGR/olfactomedin protein family with conserved structural similarity in the N-terminal region and conserved motifs in the C-terminal region

Based on the similarity between the TIGR (trabecular-meshwork inducible glucocorticoid response) (also known as myocilin) and olfactomedin protein families identified throughout the length of the TIGR protein, we have identified more distantly related proteins to determine the elements essential to the function/structure of the TIGR and olfactomedin proteins. Using a sequence walk method and the Shotgun program, we have identified a family including 31 olfactomedin domain-containing sequences. Multiple sequence alignments and secondary structure analyses were used to identify conserved sequence elements. Pairwise identity in the olfactomedin domain ranges from 8 to 64%, with an average pairwise identity of 24%. The N-terminal regions of the proteins fall into two subgroups, one including the TIGR and olfactomedin families and another group of apparently unrelated domains. The TIGR and olfactomedin sequences display conserved motifs including a residual leucine zipper region and maintain a similar secondary structure throughout the N-terminal region. The correlation between conserved elements and disease-associated mutations and apparent polymorphisms in human TIGR was also examined to evaluate the apparent importance of conserved residues to the function/structure of TIGR. Several residues have been identified as essential to the function and/or structure of the human TIGR protein based on their degree of conservation across the family and their implication in the pathogenesis of primary open-angle glaucoma. Additionally, we have identified a group of chitinase sequences containing several of the highly conserved motifs present in the C-terminal region of the olfactomedin domain-containing sequences.

Understanding nature's strategies for enzyme-catalyzed racemization and epimerization

Epimerases and racemases are enzymes that catalyze the inversion of stereochemistry in biological molecules. In this article, three distinct examples are used to illustrate the wide range of chemical strategies employed during catalysis, and the diverse set of ancestors from which these enzymes have evolved. Glutamate racemase is an example of an enzyme that operates at an "activated" stereocenter (bearing a relatively acidic proton) and employs a nonstereospecific deprotonation/reprotonation mechanism. UDP-N-Acetylglucosamine 2-epimerase acts at an "unactivated" stereocenter and uses a mechanism involving a nonstereospecific elimination/addition of UDP. L-Ribulose phosphate 4-epimerase also acts at an unactivated stereocenter and uses a nonstereospecific retroaldol/aldol mechanism.

Formation and stability of enolates of acetamide and acetate anion: an Eigen plot for proton transfer at alpha-carbonyl carbon

Second-order rate constants were determined in D(2)O for deprotonation of acetamide, N,N-dimethylacetamide, and acetate anion by deuterioxide ion and for deprotonation of acetamide by quinuclidine. The values of k(B) = 4.8 x 10(-8) M(-1) s(-1) for deprotonation of acetamide by quinuclidine (pK(BH) = 11.5) and k(BH) = 2-5 x 10(9) M(-1) s(-1) for the encounter-limited reverse protonation of the enolate by protonated quinuclidine give pK(a)(C) = 28.4 for ionization of acetamide as a carbon acid. The limiting value of k(HOH) = 1 x 10(11) s(-1) for protonation of the enolate of acetate anion by solvent water and k(HO) = 3.5 x 10(-9) M(-1) s(-1) for deprotonation of acetate anion by HO(-) give pK(a)(C) approximately 33.5 for acetate anion. The change in the rate-limiting step from chemical proton transfer to solvent reorganization results in a downward break in the slope of the plot of log k(HO) against carbon acid pK(a) for deprotonation of a wide range of neutral alpha-carbonyl carbon acids by hydroxide ion, from -0.40 to -1.0. Good estimates are reported for the stabilization of the carbonyl group relative to the enol tautomer by electron donation from alpha-SEt, alpha-OMe, alpha-NH(2), and alpha-O(-) substituents. The alpha-NH(2) and alpha-OMe groups show similar stabilizing interactions with the carbonyl group, while the interaction of alpha-O(-) is only 3.4 kcal/mol more stabilizing than for alpha-OH. We propose that destabilization of the enolate intermediates of enzymatic reactions results in an increasing recruitment of metal ions by the enzyme to provide electrophilic catalysis of enolate formation.

Kinetics and thermodynamics of mandelate racemase catalysis

Mandelate racemase (EC 5.1.2.2) from Pseudomonas putida catalyzes the interconversion of the two enantiomers of mandelic acid with remarkable proficiency, producing a rate enhancement exceeding 15 orders of magnitude. The rates of the forward and reverse reactions catalyzed by the wild-type enzyme and by a sluggish mutant (N197A) have been studied in the absence and presence of several viscosogenic agents. A partial dependence on relative solvent viscosity was observed for values of kcat and kcat/Km for the wild-type enzyme in sucrose-containing solutions. The value of kcat for the sluggish mutant was unaffected by varying solvent viscosity. However, sucrose did have a slight activating effect on mutant enzyme efficiency. In the presence of the polymeric viscosogens poly(ethylene glycol) and Ficoll, no effect on kcat or kcat/Km for the wild-type enzyme was observed. These results are consistent with both substrate binding and product dissociation being partially rate-determining in both directions. The viscosity variation method was used to estimate the rate constants comprising the steady-state expressions for kcat and kcat/Km. The rate constant for the conversion of bound (R)-mandelate to bound (S)-mandelate (k2) was found to be 889 +/- 40 s(-1) compared with a value of 654 +/- 58 s(-1) for kcat in the same direction. From the temperature dependence of Km (shown to equal K(S)), k2, and the rate constant for the uncatalyzed reaction [Bearne, S. L., and Wolfenden, R. (1997) Biochemistry 36, 1646-1656], we estimated the enthalpic and entropic changes associated with substrate binding (DeltaH = -8.9 +/- 0.8 kcal/mol, TDeltaS = -4.8 +/- 0.8 kcal/mol), the activation barrier for conversion of bound substrate to bound product (DeltaH# = +15.4 +/- 0.4 kcal/mol, TDeltaS# = +2.0 +/- 0.1 kcal/mol), and transition state stabilization (DeltaH(tx) = -22.9 +/- 0.8 kcal/mol, TDeltaS(tx) = +1.8 +/- 0.8 kcal/mol) during mandelate racemase-catalyzed racemization of (R)-mandelate at 25 degrees C. Although the high proficiency of mandelate racemase is achieved principally by enthalpic reduction, there is also a favorable and significant entropic contribution.

Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies

The protein sequence and structure databases are now sufficiently representative that strategies nature uses to evolve new catalytic functions can be identified. Groups of divergently related enzymes whose members catalyze different reactions but share a common partial reaction, intermediate, or transition state (mechanistically diverse superfamilies) have been discovered, including the enolase, amidohydrolase, thiyl radical, crotonase, vicinal-oxygen-chelate, and Fe-dependent oxidase superfamilies. Other groups of divergently related enzymes whose members catalyze different overall reactions that do not share a common mechanistic strategy (functionally distinct suprafamilies) have also been identified: (a) functionally distinct suprafamilies whose members catalyze successive transformations in the tryptophan and histidine biosynthetic pathways and (b) functionally distinct suprafamilies whose members catalyze different reactions in different metabolic pathways. An understanding of the structural bases for the catalytic diversity observed in super- and suprafamilies may provide the basis for discovering the functions of proteins and enzymes in new genomes as well as provide guidance for in vitro evolution/engineering of new enzymes.

Evolution of enzymatic activities in the enolase superfamily: crystal structures of the L-Ala-D/L-Glu epimerases from Escherichia coli and Bacillus subtilis

The members of the enolase superfamily catalyze different overall reactions, yet share a partial reaction that involves Mg(2+)-assisted enolization of the substrate carboxylate anion. The fate of the resulting enolate intermediate is determined by the active site of each enzyme. Several members of this superfamily have been structurally characterized to permit an understanding of the evolutionary strategy for using a common structural motif to catalyze different overall reactions. In the preceding paper, two new members of the superfamily were identified that catalyze the epimerization of the glutamate residue in L-Ala-D/L-Glu. These enzymes belong to the muconate lactonizing enzyme subgroup of the enolase superfamily, and their sequences are only 31% identical. The structure of YcjG, the epimerase from Escherichia coli, was determined by MAD phasing using both the SeMet-labeled protein and a heavy atom derivative. The structure of YkfB, the epimerase from Bacillus subtilis, was determined by molecular replacement using the muconate lactonizing enzyme as a search model. In this paper, we report the three-dimensional structures of these enzymes and compare them to the structure of o-succinylbenzoate synthase, another member of the muconate lactonizing enzyme subgroup.

Role of active site binding interactions in 4-chlorobenzoyl-coenzyme A dehalogenase catalysis

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA (4-HBA-CoA) via a multistep mechanism involving initial attack of Asp145 on C(4) of the substrate benzoyl ring to form a Meisenheimer intermediate (EMc), followed by expulsion of the chloride ion to form an arylated enzyme intermediate (EAr) and then ester hydrolysis in the EAr to form product. This study examines the role of binding interactions in dehalogenase catalysis. The enzyme and substrate groups positioned for favorable binding interaction were identified from the X-ray crystal structure of the enzyme-4-HBA-3'-dephospho-CoA complex. These groups were individually modified (via site-directed mutagenesis or chemical synthesis) for the purpose of disrupting the binding interaction. The changes in the Gibbs free energy of the enzyme-substrate complex (DeltaDeltaG(ES)) and enzyme-transition state complex (DeltaDeltaG) brought about by the modification were measured. Cases where DeltaDeltaG exceeds DeltaDeltaG(ES) are indicative of binding interactions used for catalysis. On the basis of this analysis, we show that the H-bond interactions between the Gly114 and Phe64 backbone amide NHs and the substrate benzoyl C=O group contribute an additional 3.1 kcal/mol of stabilization at the rate-limiting transition state. The binding interactions between the enzyme and the substrate CoA nucleotide moiety also intensify in the rate-limiting transition state, reducing the energy barrier to catalysis by an additional 3.3 kcal/mol. Together, these binding interactions contribute approximately 10(6) to the k(cat)/K(m).

Evolution of enzymatic activities in the enolase superfamily: functional assignment of unknown proteins in Bacillus subtilis and Escherichia coli as L-Ala-D/L-Glu epimerases

The members of the mechanistically diverse enolase superfamily catalyze different overall reactions by using a common catalytic strategy and structural scaffold. In the muconate lactonizing enzyme (MLE) subgroup of the superfamily, abstraction of a proton adjacent to a carboxylate group initiates reactions, including cycloisomerization (MLE), dehydration [o-succinylbenzoate synthase (OSBS)], and 1,1-proton transfer (catalyzed by an OSBS that also catalyzes a promiscuous N-acylamino acid racemase reaction). The realization that a member of the MLE subgroup could catalyze a 1,1-proton transfer reaction, albeit poorly, led to a search for other enzymes which might catalyze a 1,1-proton transfer as their physiological reaction. YcjG from Escherichia coli and YkfB from Bacillus subtilis, proteins of previously unknown function, were discovered to be L-Ala-D/L-Glu epimerases, although they also catalyze the epimerization of other dipeptides. The values of k(cat)/K(M) for L-Ala-D/L-Glu for both proteins are approximately 10(4) M(-1) s(-1). The genomic context and the substrate specificity of both YcjG and YkfB suggest roles in the metabolism of the murein peptide, of which L-Ala-D-Glu is a component. Homologues possessing L-Ala-D/L-Glu epimerase activity have been identified in at least two other organisms.

The structure of L-ribulose-5-phosphate 4-epimerase: an aldolase-like platform for epimerization

The structure of L-ribulose-5-phosphate 4-epimerase from E. coli has been solved to 2.4 A resolution using X-ray diffraction data. The structure is homo-tetrameric and displays C(4) symmetry. Each subunit has a single domain comprised of a central beta-sheet flanked on either side by layers of alpha-helices. The active site is identified by the position of the catalytic zinc residue and is located at the interface between two adjacent subunits. A remarkable feature of the structure is that it shows a very close resemblance to that of L-fuculose-1-phosphate aldolase. This is consistent with the notion that both enzymes belong to a superfamily of epimerases/aldolases that catalyze carbon-carbon bond cleavage reactions via a metal-stabilized enolate intermediate. Detailed inspection of the epimerase structure, however, indicates that despite the close overall structural similarity to class II aldolases, the enzyme has evolved distinct active site features that promote its particular chemistry.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.

Elementary steps in the acquisition of Mn2+ by the fosfomycin resistance protein (FosA)

The fosfomycin resistance protein, FosA, catalyzes the Mn(2+)-dependent addition of glutathione to the antibiotic fosfomycin, (1R,2S)-epoxypropylphosphonic acid, rendering the antibiotic inactive. The enzyme is a homodimer of 16 kDa subunits, each of which contains a single mononuclear metal site. Stopped-flow absorbance/fluorescence spectrometry provides evidence suggesting a complex kinetic mechanism for the acquisition of Mn(2+) by apoFosA. The binding of Mn(H(2)O)(6)(2+) to apoFosA alters the UV absorption and intrinsic fluorescence characteristics of the protein sufficiently to provide sensitive spectroscopic probes of metal binding. The acquisition of metal is shown to be a multistep process involving rapid preequilibrium formation of an initial complex with release of approximately two protons (k(obsd) > or = 800 s(-1)). The initial complex either rapidly dissociates or forms an intermediate coordination complex (k > 300 s(-1)) with rapid isomerization (k > or = 20 s(-1)) to a set of tight protein-metal complexes. The observed bimolecular rate constant for formation of the intermediate coordination complex is 3 x 10(5) M(-1) s(-1). The release of Mn(2+) from the protein is slow (k approximately 10(-2) s(-1)). The kinetic results suggest a more complex chelate effect than is typically observed for metal binding to simple multidentate ligands. Although the addition of the substrate, fosfomycin, has no appreciable effect on the association kinetics of enzyme and metal, it significantly decreases the dissociation rate, suggesting that the substrate interacts directly with the metal center.

Involvement of coenzyme A esters and two new enzymes, an enoyl-CoA hydratase and a CoA-transferase, in the hydration of crotonobetaine to L-carnitine by Escherichia coli

Two proteins (CaiB and CaiD) were found to catalyze the reversible biotransformation of crotonobetaine to L-carnitine in Escherichia coli in the presence of a cosubstrate (e.g., gamma-butyrobetainyl-CoA or crotonobetainyl-CoA). CaiB (45 kDa) and CaiD (27 kDa) were purified in two steps to electrophoretic homogeneity from overexpression strains. CaiB was identified as crotonobetainyl-CoA:carnitine CoA-transferase by MALDI-TOF mass spectrometry and enzymatic assays. The enzyme exhibits high cosubstrate specificity to CoA derivatives of trimethylammonium compounds. In particular, the N-terminus of CaiB shows significant identity with other CoA-transferases (e.g., FldA from Clostridium sporogenes, Frc from Oxalobacter formigenes, and BbsE from Thauera aromatica) and CoA-hydrolases (e.g., BaiF from Eubacterium sp.). CaiD was shown to be a crotonobetainyl-CoA hydratase using MALDI-TOF mass spectrometry and enzymatic assays. Besides crotonobetainyl-CoA CaiD is also able to hydrate crotonyl-CoA with a significantly lower Vmax (factor of 10(3)) but not crotonobetaine. The substrate specificity of CaiD and its homology to the crotonase confirm this enzyme as a new member of the crotonase superfamily. Concluding these results, it was verified that hydration of crotonobetaine to L-carnitine proceeds at the CoA level in two steps: the CaiD catalyzed hydration of crotonobetainyl-CoA to L-carnitinyl-CoA, followed by a CoA transfer from L-carnitinyl-CoA to crotonobetaine, catalyzed by CaiB. When gamma-butyrobetainyl-CoA was used as a cosubstrate (CoA donor), the first reaction is the CoA transfer. The optimal ratios of CaiB and CaiD during this hydration reaction, determined to be 4:1 when crotonobetainyl-CoA was used as cosubstrate and 5:1 when gamma-butyrobetainyl-CoA was used as cosubstrate, are different from that found for in vivo conditions (1:3).