Evolution of an enzyme active site: the structure of a new crystal form of muconate lactonizing enzyme compared with mandelate racemase and enolase

Current successes and remaining challenges in protein function prediction

In recent years, improvements in protein function prediction methods have led to increased success in annotating protein sequences. However, the functions of over 30% of protein-coding genes remain unknown for many sequenced genomes. Protein functions vary widely, from catalyzing chemical reactions to binding DNA or RNA or forming structures in the cell, and some types of functions are challenging to predict due to the physical features associated with those functions. Other complications in understanding protein functions arise due to the fact that many proteins have more than one function or very small differences in sequence or structure that correspond to different functions. We will discuss some of the recent developments in predicting protein functions and some of the remaining challenges.

The Birth of Genomic Enzymology: Discovery of the Mechanistically Diverse Enolase Superfamily

Enzymes are categorized into superfamilies by sequence, structural, and mechanistic similarities. The evolutionary implications can be profound. Until the mid-1990s, the approach was fragmented largely due to limited sequence and structural data. However, in 1996, Babbitt et al. published a paper in Biochemistry that demonstrated the potential power of mechanistically diverse superfamilies to identify common ancestry, predict function, and, in some cases, predict specificity. This Perspective describes the findings of the original work and reviews the current understanding of structure and mechanism in the founding family members. The outcomes of the genomic enzymology approach have reached far beyond the functional assignment of members of the enolase superfamily, inspiring the study of superfamilies and the adoption of sequence similarity networks and genome context and yielding fundamental insights into enzyme evolution.

Promoting a more integrated approach to structure and function

The connection between structure and function is one of the fundamental tenets of biology: a biological unit's structure determines its function, and, conversely, its function depends upon its structure. Historically, important advances have been made either when understanding of structure leads to questions about function or when understanding of function raises questions about the structures involved. Consequently, considering the connections between structure and function from a broader perspective might lead to the development of novel hypotheses that move our understanding of the fundamental connections between structure and function forward. Better integration of structure and function is a key component in the broader goal of reintegrating biology within and across scales. Here, we provide examples of how integrating studies of structure and function as well as comparing structure-function relationships across biological scales can lead to scientific advances. We also emphasize the potential of integrating studies of structure and function across scales for bio-inspired design and for improving biology education.

Enzymes, pseudoenzymes, and moonlighting proteins: diversity of function in protein superfamilies

As more genome sequences are elucidated, there is an increasing need for information about the functions of the millions of proteins they encode. The function of a newly sequenced protein is often estimated by sequence alignment with the sequences of proteins with known functions. However, protein superfamilies can contain members that share significant amino acid sequence and structural homology yet catalyze different reactions or act on different substrates. Some homologous proteins differ by having a second or even third function, called moonlighting proteins. More recently, it was found that most protein superfamilies also include pseudoenzymes, a protein, or a domain within a protein, that has a three-dimensional fold that resembles a conventional catalytically active enzyme, but has no catalytic activity. In this review, we discuss several examples of protein families that contain enzymes, pseudoenzymes, and moonlighting proteins. It is becoming clear that pseudoenzymes and moonlighting proteins are widespread in the evolutionary tree, and in many protein families, and they are often very similar in sequence and structure to their monofunctional and catalytically active counterparts. A greater understanding is needed to clarify when similarities and differences in amino acid sequences and structures correspond to similarities and differences in biochemical functions and cellular roles. This information can help improve programs that identify protein functions from sequence or structure and assist in more accurate annotation of sequence and structural databases, as well as in our understanding of the broad diversity of protein functions.

Catalytic properties of the metal ion variants of mandelate racemase reveal alterations in the apparent electrophilicity of the metal cofactor

Mandalate racemase (MR) from Pseudomonas putida requires a divalent metal cation, usually Mg2+, to catalyse the interconversion of the enantiomers of mandelate. Although the active site Mg2+ may be replaced by Mn2+, Co2+, or Ni2+, substitution by these metal ions does not markedly (<10-fold) alter the kinetic parameters Kappm, kappcat, and (kcat/Km)app for the substrates (R)- and (S)-mandelate, and the alternative substrate (S)-trifluorolactate. Viscosity variation experiments with Mn2+-MR showed that the metal ion plays a role in the uniform binding of the transition states for enzyme-substrate association, the chemical step, and enzyme-product dissociation. Surprisingly, the competitive inhibition constants (Ki) for inhibition of each metalloenzyme variant by benzohydroxamate did not vary significantly with the identity of the metal ion unlike the marked variation of the stability constants (K1) observed for M2+·BzH complex formation in solution. A similar trend was observed for the inhibition of the metalloenzyme variants by F-, except for Mg2+-MR, which bound F- tighter than would be predicted based on the stability constants for formation of M2+·F- complexes in solution. Thus, the enzyme modifies the enatic state of the bound metal ion cofactor so that the apparent electrophilicity of Mg2+ is enhanced, while that of Ni2+ is attenuated, resulting in a levelling effect relative to the trends observed for the free metals in solution.

An Alternative Active Site Architecture for O2 Activation in the Ergothioneine Biosynthetic EgtB from Chloracidobacterium thermophilum

Sulfoxide synthases are nonheme iron enzymes that catalyze oxidative carbon-sulfur bond formation between cysteine derivatives and N-α-trimethylhistidine as a key step in the biosynthesis of thiohistidines. The complex catalytic mechanism of this enzyme reaction has emerged as the controversial subject of several biochemical and computational studies. These studies all used the structure of the γ-glutamyl cysteine utilizing sulfoxide synthase, MthEgtB from Mycobacterium thermophilum (EC 1.14.99.50), as a structural basis. To provide an alternative model system, we have solved the crystal structure of CthEgtB from Chloracidobacterium thermophilum (EC 1.14.99.51) that utilizes cysteine as a sulfur donor. This structure reveals a completely different configuration of active site residues that are involved in oxygen binding and activation. Furthermore, comparison of the two EgtB structures enables a classification of all ergothioneine biosynthetic EgtBs into five subtypes, each characterized by unique active-site features. This active site diversity provides an excellent platform to examine the catalytic mechanism of sulfoxide synthases by comparative enzymology, but also raises the question as to why so many different solutions to the same biosynthetic problem have emerged.

A Paradigm for CH Bond Cleavage: Structural and Functional Aspects of Transition State Stabilization by Mandelate Racemase

Mandelate racemase (MR) from Pseudomonas putida catalyzes the Mg²⁺-dependent, 1,1-proton transfer reaction that racemizes (R)- and (S)-mandelate. MR shares a partial reaction (i.e., the metal ion-assisted, Brønsted base-catalyzed proton abstraction of the α-proton of carboxylic acid substrates) and structural features ((β/α)₇β-barrel and N-terminal α + β capping domains) with a vast group of homologous, yet functionally diverse, enzymes in the enolase superfamily. Mechanistic and structural studies have developed this enzyme into a paradigm for understanding how enzymes such as those of the enolase superfamily overcome kinetic and thermodynamic barriers to catalyze the abstraction of an α-proton from a carbon acid substrate with a relatively high pK_a value. Structural studies on MR bound to intermediate/transition state analogues have delineated those structural features that MR uses to stabilize transition states and enhance reaction rates of proton abstraction. Kinetic, site-directed mutagenesis, and structural studies have also revealed that the phenyl ring of the substrate migrates through the hydrophobic cavity within the active site during catalysis and that the Brønsted acid-base catalysts (Lys 166 and His 297) may be utilized as binding determinants for inhibitor recognition. In addition, structural studies on the adduct formed from the irreversible inhibition of MR by 3-hydroxypyruvate revealed that MR can form and deprotonate a Schiff-base with 3-hydroxypyruvate to yield an enol(ate)-aldehyde adduct, suggesting a possible evolutionary link between MR and the Schiff-base forming aldolases. As the archetype of the enolase superfamily, mechanistic and structural studies on MR will continue to enhance our understanding of enzyme catalysis and furnish insights into the evolution of enzyme function.

Sequence-Based Screening for Rare Enzymes: New Insights into the World of AMDases Reveal a Conserved Motif and 58 Novel Enzymes Clustering in Eight Distinct Families

Arylmalonate Decarboxylases (AMDases, EC 4.1.1.76) are very rare and mostly underexplored enzymes. Currently only four known and biochemically characterized representatives exist. However, their ability to decarboxylate α-disubstituted malonic acid derivatives to optically pure products without cofactors makes them attractive and promising candidates for the use as biocatalysts in industrial processes. Until now, AMDases could not be separated from other members of the aspartate/glutamate racemase superfamily based on their gene sequences. Within this work, a search algorithm was developed that enables a reliable prediction of AMDase activity for potential candidates. Based on specific sequence patterns and screening methods 58 novel AMDase candidate genes could be identified in this work. Thereby, AMDases with the conserved sequence pattern of Bordetella bronchiseptica’s prototype appeared to be limited to the classes of Alpha-, Beta-, and Gamma-proteobacteria. Amino acid homologies and comparison of gene surrounding sequences enabled the classification of eight enzyme clusters. Particularly striking is the accumulation of genes coding for different transporters of the tripartite tricarboxylate transporters family, TRAP transporters and ABC transporters as well as genes coding for mandelate racemases/muconate lactonizing enzymes that might be involved in substrate uptake or degradation of AMDase products. Further, three novel AMDases were characterized which showed a high enantiomeric excess (>99%) of the (R)-enantiomer of flurbiprofen. These are the recombinant AmdA and AmdV from Variovorax sp. strains HH01 and HH02, originated from soil, and AmdP from Polymorphum gilvum found by a data base search. Altogether our findings give new insights into the class of AMDases and reveal many previously unknown enzyme candidates with high potential for bioindustrial processes.

Functional evolution of PLP-dependent enzymes based on active-site structural similarities

Families of distantly related proteins typically have very low sequence identity, which hinders evolutionary analysis and functional annotation. Slowly evolving features of proteins, such as an active site, are therefore valuable for annotating putative and distantly related proteins. To date, a complete evolutionary analysis of the functional relationship of an entire enzyme family based on active-site structural similarities has not yet been undertaken. Pyridoxal-5'-phosphate (PLP) dependent enzymes are primordial enzymes that diversified in the last universal ancestor. Using the comparison of protein active site structures (CPASS) software and database, we show that the active site structures of PLP-dependent enzymes can be used to infer evolutionary relationships based on functional similarity. The enzymes successfully clustered together based on substrate specificity, function, and three-dimensional-fold. This study demonstrates the value of using active site structures for functional evolutionary analysis and the effectiveness of CPASS.

Using catalytic atom maps to predict the catalytic functions present in enzyme active sites

Catalytic atom maps (CAMs) are minimal models of enzyme active sites. The structures in the Protein Data Bank (PDB) were examined to determine if proteins with CAM-like geometries in their active sites all share the same catalytic function. We combined the CAM-based search protocol with a filter based on the weighted contact number (WCN) of the catalytic residues, a measure of the "crowdedness" of the microenvironment around a protein residue. Using this technique, a CAM based on the Ser-His-Asp catalytic triad of trypsin was able to correctly identify catalytic triads in other enzymes within 0.5 Å rmsd of the CAM with 96% accuracy. A CAM based on the Cys-Arg-(Asp/Glu) active site residues from the tyrosine phosphatase active site achieved 89% accuracy in identifying this type of catalytic functionality. Both of these CAMs were able to identify active sites across different fold types. Finally, the PDB was searched to locate proteins with catalytic functionality similar to that present in the active site of orotidine 5'-monophosphate decarboxylase (ODCase), whose mechanism is not known with certainty. A CAM, based on the conserved Lys-Asp-Lys-Asp tetrad in the ODCase active site, was used to search the PDB for enzymes with similar active sites. The ODCase active site has a geometry similar to that of Schiff base-forming Class I aldolases, with lowest aldolase rmsd to the ODCase CAM at 0.48 Å. The similarity between this CAM and the aldolase active site suggests that ODCase has the correct catalytic functionality present in its active site for the generation of a nucleophilic lysine.

Modeling regionalized volumetric differences in protein-ligand binding cavities

Identifying elements of protein structures that create differences in protein-ligand binding specificity is an essential method for explaining the molecular mechanisms underlying preferential binding. In some cases, influential mechanisms can be visually identified by experts in structural biology, but subtler mechanisms, whose significance may only be apparent from the analysis of many structures, are harder to find. To assist this process, we present a geometric algorithm and two statistical models for identifying significant structural differences in protein-ligand binding cavities. We demonstrate these methods in an analysis of sequentially nonredundant structural representatives of the canonical serine proteases and the enolase superfamily. Here, we observed that statistically significant structural variations identified experimentally established determinants of specificity. We also observed that an analysis of individual regions inside cavities can reveal areas where small differences in shape can correspond to differences in specificity.

Characterization of a novel Agrobacterium tumefaciens galactarolactone cycloisomerase enzyme for direct conversion of D-galactarolactone to 3-deoxy-2-keto-L-threo-hexarate

Microorganisms use different pathways for D-galacturonate catabolism. In the known microbial oxidative pathway, D-galacturonate is oxidized to D-galactarolactone, the lactone hydrolyzed to galactarate, which is further converted to 3-deoxy-2-keto-hexarate and α-ketoglutarate. We have shown recently that Agrobacterium tumefaciens strain C58 contains an uronate dehydrogenase (At Udh) that oxidizes D-galacturonic acid to D-galactarolactone. Here we report identification of a novel enzyme from the same A. tumefaciens strain, which we named Galactarolactone cycloisomerase (At Gci) (E.C. 5.5.1.-), for the direct conversion of the D-galactarolactone to 3-deoxy-2-keto-hexarate. The At Gci enzyme is 378 amino acids long and belongs to the mandelate racemase subgroup in the enolase superfamily. At Gci was heterologously expressed in Escherichia coli, and the purified enzyme was found to exist as an octameric form. It is active both on D-galactarolactone and D-glucarolactone, but does not work on the corresponding linear hexaric acid forms. The details of the reaction mechanism were further studied by NMR and optical rotation demonstrating that the reaction product of At Gci from D-galactaro-1,4-lactone and D-glucaro-1,4-lactone conversion is in both cases the L-threo form of 3-deoxy-2-keto-hexarate.

FINDSITE: a combined evolution/structure-based approach to protein function prediction

A key challenge of the post-genomic era is the identification of the function(s) of all the molecules in a given organism. Here, we review the status of sequence and structure-based approaches to protein function inference and ligand screening that can provide functional insights for a significant fraction of the approximately 50% of ORFs of unassigned function in an average proteome. We then describe FINDSITE, a recently developed algorithm for ligand binding site prediction, ligand screening and molecular function prediction, which is based on binding site conservation across evolutionary distant proteins identified by threading. Importantly, FINDSITE gives comparable results when high-resolution experimental structures as well as predicted protein models are used.

Evolution of enzymatic activities in the enolase superfamily: stereochemically distinct mechanisms in two families of cis,cis-muconate lactonizing enzymes

The mechanistically diverse enolase superfamily is a paradigm for elucidating Nature's strategies for divergent evolution of enzyme function. Each of the different reactions catalyzed by members of the superfamily is initiated by abstraction of the alpha-proton of a carboxylate substrate that is coordinated to an essential Mg(2+). The muconate lactonizing enzyme (MLE) from Pseudomonas putida, a member of a family that catalyzes the syn-cycloisomerization of cis,cis-muconate to (4S)-muconolactone in the beta-ketoadipate pathway, has provided critical insights into the structural bases for evolution of function within the superfamily. A second, divergent family of homologous MLEs that catalyzes anti-cycloisomerization has been identified. Structures of members of both families liganded with the common (4S)-muconolactone product (syn, Pseudomonas fluorescens, gi 70731221 ; anti, Mycobacterium smegmatis, gi 118470554 ) document that the conserved Lys at the end of the second beta-strand in the (beta/alpha)(7)beta-barrel domain serves as the acid catalyst in both reactions. The different stereochemical courses (syn and anti) result from different structural strategies for determining substrate specificity: although the distal carboxylate group of the cis,cis-muconate substrate attacks the same face of the proximal double bond, opposite faces of the resulting enolate anion intermediate are presented to the conserved Lys acid catalyst. The discovery of two families of homologous, but stereochemically distinct, MLEs likely provides an example of "pseudoconvergent" evolution of the same function from different homologous progenitors within the enolase superfamily, in which different spatial arrangements of active site functional groups and substrate specificity determinants support catalysis of the same reaction.

Prediction and assignment of function for a divergent N-succinyl amino acid racemase

The protein databases contain many proteins with unknown function. A computational approach for predicting ligand specificity that requires only the sequence of the unknown protein would be valuable for directing experiment-based assignment of function. We focused on a family of unknown proteins in the mechanistically diverse enolase superfamily and used two approaches to assign function: (i) enzymatic assays using libraries of potential substrates, and (ii) in silico docking of the same libraries using a homology model based on the most similar (35% sequence identity) characterized protein. The results matched closely; an experimentally determined structure confirmed the predicted structure of the substrate-liganded complex. We assigned the N-succinyl arginine/lysine racemase function to the family, correcting the annotation (L-Ala-D/L-Glu epimerase) based on the function of the most similar characterized homolog. These studies establish that ligand docking to a homology model can facilitate functional assignment of unknown proteins by restricting the identities of the possible substrates that must be experimentally tested.

Novel enzymes through design and evolution

The generation of enzymes with new catalytic activities remains a major challenge. So far, several different strategies have been developed to tackle this problem, including site-directed mutagenesis, random mutagenesis (directed evolution), antibody catalysis, computational redesign, and de novo methods. Using these techniques, a broad array of novel enzymes has been created (aldolases, decarboxylases, dehydratases, isomerases, oxidases, reductases, and others), although their low efficiencies (10 to 100 M(-1) s(-l)) compared to those of the best natural enzymes (10(6) to 10(8) M(-1) s(-1)) remains a significant concern. Whereas rational design might be the most promising and versatile approach to generating new activities, directed evolution seems to be the best way to optimize the catalytic properties of novel enzymes. Indeed, impressive successes in enzyme engineering have resulted from a combination of rational and random design.

Evolutionary migration of a post-translationally modified active-site residue in the proton-pumping heme-copper oxygen reductases

In the respiratory chains of aerobic organisms, oxygen reductase members of the heme-copper superfamily couple the reduction of O2 to proton pumping, generating an electrochemical gradient. There are three distinct families of heme-copper oxygen reductases: A, B, and C types. The A- and B-type oxygen reductases have an active-site tyrosine that forms a unique cross-linked histidine-tyrosine cofactor. In the C-type oxygen reductases (also called cbb3 oxidases), an analogous active-site tyrosine has recently been predicted by molecular modeling to be located within a different transmembrane helix in comparison to the A- and B-type oxygen reductases. In this work, Fourier-transform mass spectrometry is used to show that the predicted tyrosine forms a histidine-tyrosine cross-linked cofactor in the active site of the C-type oxygen reductases. This is the first known example of the evolutionary migration of a post-translationally modified active-site residue. It also verifies the presence of a unique cofactor in all three families of proton-pumping respiratory oxidases, demonstrating that these enzymes likely share a common reaction mechanism and that the histidine-tyrosine cofactor may be a required component for proton pumping.

FAST-NMR: functional annotation screening technology using NMR spectroscopy

An abundance of protein structures emerging from structural genomics and the Protein Structure Initiative (PSI) are not amenable to ready functional assignment because of a lack of sequence and structural homology to proteins of known function. We describe a high-throughput NMR methodology (FAST-NMR) to annotate the biological function of novel proteins through the structural and sequence analysis of protein-ligand interactions. This is based on basic tenets of biochemistry where proteins with similar functions will have similar active sites and exhibit similar ligand binding interactions, despite global differences in sequence and structure. Protein-ligand interactions are determined through a tiered NMR screen using a library composed of compounds with known biological activity. A rapid co-structure is determined by combining the experimental identification of the ligand binding site from NMR chemical shift perturbations with the protein-ligand docking program AutoDock. Our CPASS (Comparison of Protein Active Site Structures) software and database are then used to compare this active site with proteins of known function. The methodology is demonstrated using unannotated protein SAV1430 from Staphylococcus aureus.

α-Methylacyl-CoA Racemase from Mycobacterium tuberculosis

Alpha-methylacyl-CoA racemase (Amacr) catalyzes the racemization of alpha-methyl-branched CoA esters. Sequence comparisons have shown that this enzyme is a member of the family III CoA transferases. The mammalian Amacr is involved in bile acid synthesis and branched-chain fatty acid degradation. In human, mutated variants of Amacr have been shown to be associated with disease states. Amino acid sequence alignment of Amacrs and its homologues from various species revealed 26 conserved protic residues, assumed to be potential candidates as catalytic residues. Amacr from Mycobacterium tuberculosis (MCR) was taken as a representative of the racemases. To determine their importance for efficient catalysis, each of these 26 protic residues of MCR was mutated into an alanine, respectively, and the mutated variants were overexpressed in Escherichia coli. It was found that four variants (R91A, H126A, D156A, and E241A) were properly folded but had much decreased catalytic efficiency. Apparently, Arg91, His126, Asp156, and Glu241 are important catalytic residues of MCR. The importance of these residues for catalysis can be rationalized by the 1.8 A resolution crystal structure of MCR, which shows that the catalytic site is at the interface between the large and small domain of two different subunits of the dimeric enzyme. This crystal structure is the first structure of a complete enzyme of the bile acid synthesis pathway. It shows that MCR has unique structural features, not seen in the structures of the sequence related formyl-CoA transferases, suggesting that the family III CoA transferases can be subdivided in at least two classes, being racemases and CoA transferases.

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.

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.

Automated selection of positions determining functional specificity of proteins by comparative analysis of orthologous groups in protein families

The increasing volume of genomic data opens new possibilities for analysis of protein function. We introduce a method for automated selection of residues that determine the functional specificity of proteins with a common general function (the specificity-determining positions [SDP] prediction method). Such residues are assumed to be conserved within groups of orthologs (that may be assumed to have the same specificity) and to vary between paralogs. Thus, considering a multiple sequence alignment of a protein family divided into orthologous groups, one can select positions where the distribution of amino acids correlates with this division. Unlike previously published techniques, the introduced method directly takes into account nonuniformity of amino acid substitution frequencies. In addition, it does not require setting arbitrary thresholds. Instead, a formal procedure for threshold selection using the Bernoulli estimator is implemented. We tested the SDP prediction method on the LacI family of bacterial transcription factors and a sample of bacterial water and glycerol transporters belonging to the major intrinsic protein (MIP) family. In both cases, the comparison with available experimental and structural data strongly supported our predictions.

Structural Basis for Catalytic Racemization and Substrate Specificity of an N-Acylamino Acid Racemase Homologue from Deinococcus radiodurans

N-acylamino acid racemase (NAAAR) catalyzes the racemization of N-acylamino acids and can be used in concert with an aminoacylase to produce enantiopure alpha-amino acids, a process that has potential industrial applications. Here we have cloned and characterized an NAAAR homologue from a radiation-resistant ancient bacterium, Deinococcus radiodurans. The expressed NAAAR racemized various substrates at an optimal temperature of 60 degrees C and had Km values of 24.8 mM and 12.3 mM for N-acetyl-D-methionine and N-acetyl-L-methionine, respectively. The crystal structure of NAAAR was solved to 1.3 A resolution using multiwavelength anomalous dispersion (MAD) methods. The structure consists of a homooctamer in which each subunit has an architecture characteristic of enolases with a capping domain and a (beta/alpha)7 beta barrel domain. The NAAAR.Mg2+ and NAAAR.N-acetyl-L-glutamine.Mg2+ structures were also determined, allowing us to define the Lys170-Asp195-Glu220-Asp245-Lys269 framework for catalyzing 1,1-proton exchange of N-acylamino acids. Four subsites enclosing the substrate are identified: catalytic site, metal-binding site, side-chain-binding region, and a flexible lid region. The high conservation of catalytic and metal-binding sites in different enolases reflects the essentiality of a common catalytic platform, allowing these enzymes to robustly abstract alpha-protons of various carboxylate substrates efficiently. The other subsites involved in substrate recognition are less conserved, suggesting that divergent evolution has led to functionally distinct enzymes.

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.

A domain for editing by an archaebacterial tRNA synthetase

The rules of the genetic code are established by aminoacylations of transfer RNAs by aminoacyl tRNA synthetases. New codon assignments, and the introduction of new kinds of amino acids, are blocked by vigorous tRNA-dependent editing reactions occurring at hydrolytic sites embedded within specialized domains in the synthetases. For some synthetases, these domains were present at the time of the last common ancestor and were fixed in evolution through all three of the kingdoms of life. Significantly, a well characterized domain for editing found in bacterial and eukaryotic threonyl- and all alanyl-tRNA synthetases is missing from archaebacterial threonine enzymes. Here we show that the archaebacterial Methanosarcina mazei ThrRS efficiently misactivates serine, but does not fuse serine to tRNA. Consistent with this observation, the enzyme cleared serine that was linked to threonine-specific tRNAs. M. mazei and most other archaebacterial ThrRSs have a domain, N2(A), fused to the N terminus and not found in bacterial or eukaryotic orthologs. Mutations at conserved residues in this domain led to an inability to clear threonine-specific tRNA mischarged with serine. Thus, these results demonstrate a domain for editing that is distinct from all others, is restricted to just one branch of the tree of life, and was most likely added to archaebacterial ThrRSs after the eukaryote/archaebacteria split.

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.

The structural determination of phosphosulfolactate synthase from Methanococcus jannaschii at 1.7-A resolution: an enolase that is not an enolase

Members of the enolase mechanistically diverse superfamily catalyze a wide variety of chemical reactions that are related by a common mechanistic feature, the abstraction of a proton adjacent to a carboxylate group. Recent investigations into the function and mechanism of the phosphosulfolactate synthase encoded by the ComA gene in Methanococcus jannaschii have suggested that ComA, which catalyzes the stereospecific Michael addition of sulfite to phosphoenolpyruvate to form phosphosulfolactate, may be a member of the enolase superfamily. The ComA-catalyzed reaction, the first step in the coenzyme M biosynthetic pathway, likely proceeds via a Mg2+ ion-stabilized enolate intermediate in a manner similar to that observed for members of the enolase superfamily. ComA, however, has no significant sequence similarity to any known enolase. Here we report the x-ray crystal structure of ComA to 1.7-A resolution. The overall fold for ComA is an (alpha/beta)8 barrel that assembles with two other ComA molecules to form a trimer in which three active sites are created at the subunit interfaces. From the positions of two ordered sulfate ions in the active site, a model for the binding of phosphoenolpyruvate and sulfite is proposed. Despite its mechanistic similarity to the enolase superfamily, the overall structure and active site architecture of ComA are unlike any member of the enolase superfamily, which suggests that ComA is not a member of the enolase superfamily but instead acquired an enolase-type mechanism through convergent evolution.

The structure of Pseudomonas P51 Cl-muconate lactonizing enzyme: co-evolution of structure and dynamics with the dehalogenation function

Bacterial muconate lactonizing enzymes (MLEs) catalyze the conversion of cis,cis-muconate as a part of the beta-ketoadipate pathway, and some MLEs are also able to dehalogenate chlorinated muconates (Cl-MLEs). The basis for the Cl-MLEs dehalogenating activity is still unclear. To further elucidate the differences between MLEs and Cl-MLEs, we have solved the structure of Pseudomonas P51 Cl-MLE at 1.95 A resolution. Comparison of Pseudomonas MLE and Cl-MLE structures reveals the presence of a large cavity in the Cl-MLEs. The cavity may be related to conformational changes on substrate binding in Cl-MLEs, at Gly52. Site-directed mutagenesis on Pseudomonas MLE core positions to the equivalent Cl-MLE residues showed that the variant Thr52Gly was rather inactive, whereas the Thr52Gly-Phe103Ser variant had regained part of the activity. These residues form a hydrogen bond in the Cl-MLEs. The Cl-MLE structure, as a result of the Thr-to-Gly change, is more flexible than MLE: As a mobile loop closes over the active site, a conformational change at Gly52 is observed in Cl-MLEs. The loose packing and structural motions in Cl-MLE may be required for the rotation of the lactone ring in the active site necessary for the dehalogenating activity of Cl-MLEs. Furthermore, we also suggest that differences in the active site mobile loop sequence between MLEs and Cl-MLEs result in lower active site polarity in Cl-MLEs, possibly affecting catalysis. These changes could result in slower product release from Cl-MLEs and make it a better enzyme for dehalogenation of substrate.

Gene cluster of Arthrobacter ilicis Ru61a involved in the degradation of quinaldine to anthranilate: characterization and functional expression of the quinaldine 4-oxidase qoxLMS genes

A genetic analysis of the anthranilate pathway of quinaldine degradation was performed. A 23-kb region of DNA from Arthrobacter ilicis Rü61a was cloned into the cosmid pVK100. Although Escherichia coli clones containing the recombinant cosmid did not transform quinaldine, cosmids harboring the 23-kb region, or a 10.8-kb stretch of this region, conferred to Pseudomonas putida KT2440 the ability to cometabolically convert quinaldine to anthranilate. The 10.8-kb fragment thus contains the genes coding for quinaldine 4-oxidase (Qox), 1H-4-oxoquinaldine 3-monooxygenase, 1H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase, and N-acetylanthranilate amidase. The qoxLMS genes coding for the molybdopterin cytosine dinucleotide-(MCD-), FeSI-, FeSII-, and FAD-containing Qox were inserted into the expression vector pJB653, generating pKP1. Qox is the first MCD-containing enzyme to be synthesized in a catalytically fully competent form by a heterologous host, P. putida KT2440 pKP1; the catalytic properties and the UV-visible and EPR spectra of Qox purified from P. putida KT2440 pKP1 were essentially like those of wild-type Qox. This provides a starting point for the construction of protein variants of Qox by site-directed mutagenesis. Downstream of the qoxLMS genes, a putative gene whose deduced amino acid sequence showed 37% similarity to the cofactor-inserting chaperone XdhC was located. Additional open reading frames identified on the 23-kb segment may encode further enzymes (a glutamyl tRNA synthetase, an esterase, two short-chain dehydrogenases/reductases, an ATPase belonging to the AAA family, a 2-hydroxyhepta-2,4-diene-1,7-dioate isomerase/5-oxopent-3-ene-1,2,5-tricarboxylate decarboxylase-like protein, and an enzyme of the mandelate racemase group) and hypothetical proteins involved in transcriptional regulation, and metabolite transport.

Expanding protein universe and its origin from the biological Big Bang

The bottom-up approach to understanding the evolution of organisms is by studying molecular evolution. With the large number of protein structures identified in the past decades, we have discovered peculiar patterns that nature imprints on protein structural space in the course of evolution. In particular, we have discovered that the universe of protein structures is organized hierarchically into a scale-free network. By understanding the cause of these patterns, we attempt to glance at the very origin of life.

Using orthologous and paralogous proteins to identify specificity-determining residues in bacterial transcription factors

Concepts of orthology and paralogy are become increasingly important as whole-genome comparison allows their identification in complete genomes. Functional specificity of proteins is assumed to be conserved among orthologs and is different among paralogs. We used this assumption to identify residues which determine specificity of protein-DNA and protein-ligand recognition. Finding such residues is crucial for understanding mechanisms of molecular recognition and for rational protein and drug design. Assuming conservation of specificity among orthologs and different specificity of paralogs, we identify residues that correlate with this grouping by specificity. The method is taking advantage of complete genomes to find multiple orthologs and paralogs. The central part of this method is a procedure to compute statistical significance of the predictions. The procedure is based on a simple statistical model of protein evolution. When applied to a large family of bacterial transcription factors, our method identified 12 residues that are presumed to determine the protein-DNA and protein-ligand recognition specificity. Structural analysis of the proteins and available experimental results strongly support our predictions. Our results suggest new experiments aimed at rational re-design of specificity in bacterial transcription factors by a minimal number of mutations.

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.

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.

The structure of 3-methylaspartase from Clostridium tetanomorphum functions via the common enolase chemical step

Methylaspartate ammonia-lyase (3-methylaspartase, MAL; EC ) catalyzes the reversible anti elimination of ammonia from L-threo-(2S,3S)-3-methylaspartic acid to give mesaconic acid. This reaction lies on the main catabolic pathway for glutamate in Clostridium tetanomorphum. MAL requires monovalent and divalent cation cofactors for full catalytic activity. The enzyme has attracted interest because of its potential use as a biocatalyst. The structure of C. tetanomorphum MAL has been solved to 1.9-A resolution by the single-wavelength anomalous diffraction method. A divalent metal ion complex of the protein has also been determined. MAL is a homodimer with each monomer consisting of two domains. One is an alpha/beta-barrel, and the other smaller domain is mainly beta-strands. The smaller domain partially occludes the C terminus of the barrel and forms a large cleft. The structure identifies MAL as belonging to the enolase superfamily of enzymes. The metal ion site is located in a large cleft between the domains. Potential active site residues have been identified based on a combination of their proximity to a metal ion site, molecular modeling, and sequence homology. In common with all members of the enolase superfamily, the carboxylic acid of the substrate is co-ordinated by the metal ions, and a proton adjacent to a carboxylic acid group of the substrate is abstracted by a base. In MAL, it appears that Lys(331) removes the alpha-proton of methylaspartic acid. This motif is the defining mechanistic characteristic of the enolase superfamily of which all have a common fold. The degree of structural conservation is remarkable given only four residues are absolutely conserved.

Proton transfer at carbon

The viability of living systems requires that C--H bonds of biological molecules be stable in water, but that there also be a mechanism for shortening the timescale for their heterolytic cleavage through enzymatic catalysis of a variety of catabolic and metabolic reactions. An understanding of the mechanism of enzymatic catalysis of proton transfer at carbon requires the integration of results of studies to determine the structure of the enzyme-substrate complex with model studies on the mechanism for the non-enzymatic reaction in water, and the effect of the local protein environment on the stability of the transition state for this reaction. A common theme is the importance of electrostatic interactions in providing stabilization of bound carbanion intermediates of enzyme-catalyzed proton-transfer reactions.

Evolution of function in protein superfamilies, from a structural perspective

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.

The TIM-barrel fold: a versatile framework for efficient enzymes

Recent studies on triosephosphate isomerase (TIM)-barrel enzymes highlight the remarkable versatility of the TIM-barrel scaffold. At least 15 distinct enzyme families use this framework to generate the appropriate active site geometry, always at the C-terminal end of the eight parallel beta-strands of the barrel. Sequence and structure comparisons now suggest that many of the TIM-barrel enzymes are evolutionarily related. Common structural properties of TIM-barrel enzymes are discussed.

Methodologies for target selection in structural genomics

As the number of complete genomes that have been sequenced keeps growing, unknown areas of the protein space are revealed and new horizons open up. Most of this information will be fully appreciated only when the structural information about the encoded proteins becomes available. The goal of structural genomics is to direct large-scale efforts of protein structure determination, so as to increase the impact of these efforts. This review focuses on current approaches in structural genomics aimed at selecting representative proteins as targets for structure determination. We will discuss the concept of representative structures/folds, the current methodologies for identifying those proteins, and computational techniques for identifying proteins which are expected to adopt new structural folds.

Mutations in catB, the gene encoding muconate cycloisomerase, activate transcription of the distal ben genes and contribute to a complex regulatory circuit in Acinetobacter sp. strain ADP1

Mutants of the bacterium Acinetobacter sp. strain ADP1 were selected to grow on benzoate without the BenM transcriptional activator. In the wild type, BenM responds to benzoate and cis,cis-muconate to activate expression of the benABCDE operon, which is involved in benzoate catabolism. This operon encodes enzymes that convert benzoate to catechol, a compound subsequently degraded by cat gene-encoded enzymes. In this report, four spontaneous mutants were found to carry catB mutations that enabled BenM-independent growth on benzoate. catB encodes muconate cycloisomerase, an enzyme required for benzoate catabolism. Its substrate, cis,cis-muconate, is enzymatically produced from catechol by the catA-encoded catechol 1,2-dioxygenase. Muconate cycloisomerase was purified to homogeneity from the wild type and the catB mutants. Each purified enzyme was active, although there were differences in the catalytic properties of the wild type and variant muconate cycloisomerases. Strains with a chromosomal benA::lacZ transcriptional fusion were constructed and used to investigate how catB mutations affect growth on benzoate. All of the catB mutations increased cis,cis-muconate-activated ben gene expression in strains lacking BenM. A model is presented in which the catB mutations reduce muconate cycloisomerase activity during growth on benzoate, thereby increasing intracellular cis, cis-muconate concentrations. This, in turn, may allow CatM, an activator similar to BenM in sequence and function, to activate ben gene transcription. CatM normally responds to cis,cis-muconate to activate cat gene expression. Consistent with this model, muconate cylcoisomerase specific activities in cell extracts of benzoate-grown catB mutants were low relative to that of the wild type. Moreover, the catechol 1,2-dioxygenase activities of the mutants were elevated, which may result from CatM responding to the altered intracellular levels of cis,cis-muconate and increasing catA expression. Collectively, these results support the important role of metabolite concentrations in controlling benzoate degradation via a complex transcriptional regulatory circuit.

Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion

The atomic structures of two proteins in the histidine biosynthesis pathway consist of beta/alpha barrels with a twofold repeat pattern. It is likely that these proteins evolved by twofold gene duplication and gene fusion from a common half-barrel ancestor. These ancestral domains are not visible as independent domains in the extant proteins but can be inferred from a combination of sequence and structural analysis. The detection of subdomain structures may be useful in efforts to search genome sequences for functionally and structurally related proteins.

Directed evolution of a (beta alpha)8-barrel enzyme to catalyze related reactions in two different metabolic pathways

Enzymes participating in different metabolic pathways often have similar catalytic mechanisms and structures, suggesting their evolution from a common ancestral precursor enzyme. We sought to create a precursor-like enzyme for N'-[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) isomerase (HisA; EC ) and phosphoribosylanthranilate (PRA) isomerase (TrpF; EC ), which catalyze similar reactions in the biosynthesis of the amino acids histidine and tryptophan and have a similar (betaalpha)(8)-barrel structure. Using random mutagenesis and selection, we generated several HisA variants that catalyze the TrpF reaction both in vivo and in vitro, and one of these variants retained significant HisA activity. A more detailed analysis revealed that a single amino acid exchange could establish TrpF activity on the HisA scaffold. These findings suggest that HisA and TrpF may have evolved from an ancestral enzyme of broader substrate specificity and underscore that (betaalpha)(8)-barrel enzymes are very suitable for the design of new catalytic activities.

Directed evolution of new catalytic activity using the alpha/beta-barrel scaffold

In biological systems, enzymes catalyse the efficient synthesis of complex molecules under benign conditions, but widespread industrial use of these biocatalysts depends crucially on the development of new enzymes with useful catalytic functions. The evolution of enzymes in biological systems often involves the acquisition of new catalytic or binding properties by an existing protein scaffold. Here we mimic this strategy using the most common fold in enzymes, the alpha/beta-barrel, as the scaffold. By combining an existing binding site for structural elements of phosphoribosylanthranilate with a catalytic template required for isomerase activity, we are able to evolve phosphoribosylanthranilate isomerase activity from the scaffold of indole-3-glycerol-phosphate synthase. We find that targeting the catalytic template for in vitro mutagenesis and recombination, followed by in vivo selection, results in a new phosphoribosylanthranilate isomerase that has catalytic properties similar to those of the natural enzyme, with an even higher specificity constant. Our demonstration of divergent evolution and the widespread occurrence of the alpha/beta-barrel suggest that this scaffold may be a fold of choice for the directed evolution of new biocatalysts.

Evolution of protein function, from a structural perspective

The recent growth in structural data, and ensuing analyses, have revealed the structural and functional versatility of protein families. With respect to enzymes, local active-site mutations, variations in surface loops and recruitment of additional domains accommodate the diverse substrate specificities and catalytic activities observed within several superfamilies. Conversely, some functions have more than one structural solution, having evolved independently several times during evolution. Combined with the existence of multi-functional genes, which have arisen by gene recruitment, these phenomena must be considered in the process of genome annotation.