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

Interrogating l-fuconate dehydratase with tartronate and 3-hydroxypyruvate reveals subtle differences within the mandelate racemase-subgroup of the enolase superfamily

Enzymes of the enolase superfamily share a conserved structure and a common partial reaction (i.e., metal-assisted, Brønsted base-catalyzed enol(ate) formation). The architectures of the enolization apparatus at the active sites of the mandelate racemase (MR)-subgroup members MR and l-fuconate dehydratase (FucD) are almost indistinguishable at the structural level. Tartronate and 3-hydroxypyruvate (3-HP) recognize the enolization apparatus and can be used to interrogate the active sites for differences that may not be apparent from structural data. We report a circular dichroism-based assay of FucD activity that monitors the change in ellipticity at 216 nm (Δ[Θ]S-P = 8985 ± 87 deg cm2 mol-1) accompanying the conversion of l-fuconate to 2-keto-3-deoxy-l-fuconate. Tartronate was a linear mixed-type inhibitor of FucD (Ki = 8.4 ± 0.7 mM, αKi = 63 ± 11 mM), binding 18-fold weaker than l-fuconate, compared with 2-fold weaker binding of tartronate by MR relative to mandelate. 3-HP irreversibly inactivated FucD (kinact/KI = 0.018 ± 0.002 M-1s-1) with an efficiency that was ∼4.6 × 103-fold less than that observed with MR. The inactivation arose predominantly from modifications at multiple sites and Tris-HCl, but not l-fuconate, afforded protection against inactivation. Similar to the reaction of 3-HP with MR, 3-HP modified the Brønsted base catalyst (Lys 220) at the active site of FucD, which was facilitated by the Brønsted acid catalyst His 351. Thus, the interactions of tartronate and 3-HP with MR and FucD revealed differences in binding affinity and reactivity that differentiated between the enzymes' enolization apparatuses.

Identification of sequence determinants for the ABHD14 enzymes

Over the course of evolution, enzymes have developed remarkable functional diversity in catalyzing important chemical reactions across various organisms, and understanding how new enzyme functions might have evolved remains an important question in modern enzymology. To systematically annotate functions, based on their protein sequences and available biochemical studies, enzymes with similar catalytic mechanisms have been clustered together into an enzyme superfamily. Typically, enzymes within a superfamily have similar overall three-dimensional structures, conserved catalytic residues, but large variations in substrate recognition sites and residues to accommodate the diverse biochemical reactions that are catalyzed within the superfamily. The serine hydrolases are an excellent example of such an enzyme superfamily. Based on known enzymatic activities and protein sequences, they are split almost equally into the serine proteases and metabolic serine hydrolases. Within the metabolic serine hydrolases, there are two outlying members, ABHD14A and ABHD14B, that have high sequence similarity, but their biological functions remained cryptic till recently. While ABHD14A still lacks any functional annotation to date, we recently showed that ABHD14B functions as a lysine deacetylase in mammals. Given their high sequence similarity, automated databases often wrongly assign ABHD14A and ABHD14B as the same enzyme, and therefore, annotating functions to them in various organisms has been problematic. In this article, we present a bioinformatics study coupled with biochemical experiments, which identifies key sequence determinants for both ABHD14A and ABHD14B, and enable better classification for them. In addition, we map these enzymes on an evolutionary timescale and provide a much-wanted resource for studying these interesting enzymes in different organisms.

Mechanistic and evolutionary insights into alkaline phosphatase superfamily through structure-function studies on Sphingomonas alkaline phosphatase

Alkaline phosphatases (APs), represented by E. coli AP (ECAP), employ an arginine residue to stabilize the phosphoryl group in the active site; whereas, AP from Sphingomonas (SPAP) shows a unique combination of substrate-binding residues; Thr89, Asn110, Lys171, and Arg173. Although such combination has been observed only in SPAP, these residues are present separately in different members of the AP superfamily. Here, we establish the presence of two distinct classes of APs; ECAP-type and SPAP-type. Bioinformatic analyses show that SPAP-type of APs are widely distributed in the bacterial kingdom. The role of active site residues in the catalytic mechanism has been delineated through a set of crystal structures reported here. These structures, representing different stages of the reaction pathway provide wealth of information for the catalytic mechanism. Despite critical differences in the substrate binding residues, SPAP follows a mechanism similar to that of ECAP-type of APs. Structure-based phylogenetic analysis suggests that SPAP and ECAP may have diverged very early during the evolution from a common ancestor. Moreover, it is proposed that the SPAP-type of APs are fundamental members of the AP superfamily and are more closely related to other members of the superfamily as compared to the ECAP-type of APs.

Anti-alpha enolase multi-antibody specificity in human diseases. Clinical significance and molecular mechanisms

Alpha-enolase (Eno) is an ubiquitary glycolytic enzyme playing multiple functions that go well beyond its principal metabolic role of energy supplier during glycolysis. Eno is localized in the cytoplasm, but also expressed on the cell membrane, where it binds plasminogen allowing its activation. Its shorter form, in the nucleus, acts as transcription factor. In inflammatory conditions, Eno undergoes post-translational modifications, such as citrullination, oxidation and phosphorylation. Eno is also an autoantigen in different disorders. In fact, autoantibodies to Eno have been detected in rheumatoid arthritis, lupus nephritis, primary glomerulonephritis, cancer, infections and other disorders, and in many cases they represent specific markers to be utilized in clinical practice. Anti-Eno antibodies in the different clinical conditions are not equal: they differ in isotype and often recognize different epitopes on the enzyme. IgG1 and IgG3 are prevalent in Rheumatoid Arthritis, IgG2 in Lupus nephritis and IgG4 in primary autoimmune glomerulopathy. This review analyzes the characteristics of anti-Eno autoantibodies in autoimmune disorders and cancer, describing their fine specificity and isotype restriction. The post-translational modifications that are target of autoantibodies are also discussed, as they represent the basis for elucidating the molecular mechanisms responsible for epitope generation. Despite an impressive amount of experimental work on anti-Eno antibodies, it is still necessary to validate the use of anti-Eno antibodies as biomarkers of selected diseases and extend the knowledge on the mechanisms of anti-Eno autoantibody production. Strategies that downmodulate the immune response to Eno may represent in the future novel approaches in the treatment of autoimmune disorders.

Enzymology of Alternative Carbohydrate Catabolic Pathways

The Embden–Meyerhof–Parnas (EMP) and Entner–Doudoroff (ED) pathways are considered the most abundant catabolic pathways found in microorganisms, and ED enzymes have been shown to also be widespread in cyanobacteria, algae and plants. In a large number of organisms, especially common strains used in molecular biology, these pathways account for the catabolism of glucose. The existence of pathways for other carbohydrates that are relevant to biomass utilization has been recognized as new strains have been characterized among thermophilic bacteria and Archaea that are able to transform simple polysaccharides from biomass to more complex and potentially valuable precursors for industrial microbiology. Many of the variants of the ED pathway have the key dehydratase enzyme involved in the oxidation of sugar derived from different families such as the enolase, IlvD/EDD and xylose-isomerase-like superfamilies. There are the variations in structure of proteins that have the same specificity and generally greater-than-expected substrate promiscuity. Typical biomass lignocellulose has an abundance of xylan, and four different pathways have been described, which include the Weimberg and Dahms pathways initially oxidizing xylose to xylono-gamma-lactone/xylonic acid, as well as the major xylose isomerase pathway. The recent realization that xylan constitutes a large proportion of biomass has generated interest in exploiting the compound for value-added precursors, but few chassis microorganisms can grow on xylose. Arabinose is part of lignocellulose biomass and can be metabolized with similar pathways to xylose, as well as an oxidative pathway. Like enzymes in many non-phosphorylative carbohydrate pathways, enzymes involved in L-arabinose pathways from bacteria and Archaea show metabolic and substrate promiscuity. A similar multiplicity of pathways was observed for other biomass-derived sugars such as L-rhamnose and L-fucose, but D-mannose appears to be distinct in that a non-phosphorylative version of the ED pathway has not been reported. Many bacteria and Archaea are able to grow on mannose but, as with other minor sugars, much of the information has been derived from whole cell studies with additional enzyme proteins being incorporated, and so far, only one synthetic pathway has been described. There appears to be a need for further discovery studies to clarify the general ability of many microorganisms to grow on the rarer sugars, as well as evaluation of the many gene copies displayed by marine bacteria.

Biochemical characterization of 2-phosphinomethylmalate synthase from Streptomyces hygroscopicus: A member of the DRE-TIM metallolyase superfamily

2-Phosphinomethylmalate synthase (PMMS) from Streptomyces hygroscopicus catalyzes the first step in the biosynthesis of the herbicide bialophos using 3-phosphinopyruvic acid and acetyl coenzyme A as substrates to form 2-phosphinomethylmalic acid and coenzyme A. PMMS belongs to the Claisen condensation-like (CC-like) subgroup of the DRE-TIM metallolyase superfamily, which uses conserved active site architecture to catalyze a functionally-diverse set of reactions. Analysis of a sequence similarity network for the CC-like subgroup identified PMMS and the related R-citrate synthase in an early-diverging cluster suggesting that this group of sequences are more distinct in relation to other Claisen-condensation subgroup members. To better understand the structure/function landscape of the CC-like subgroup PMMS was recombinantly expressed in Escherichia coli, purified, and characterized with respect to its enzymatic properties. Using oxaloacetate as a substrate analog, the recombinantly-produced enzyme exhibited improved Michaelis constants relative to the previously reported natively-produced enzyme. Results from pH rate profiles and kinetic isotope effects were consistent with results from other members of the CC-like subgroup supporting acid-base chemistry and hydrolysis of the direct Claisen-condensation product as the rate-determining step. Results of site-directed mutagenesis experiments indicate that PMMS uses an active-site architecture similar to homocitrate synthase to select for a dicarboxylic acid substrate.

Crystal structure of MOA in complex with a peptide fragment: A protease caught in flagranti

The Marasmius oreades agglutinin (MOA) is the holotype of an emerging family of fungal chimerolectins and an active Ca2+/Mn2+-dependent protease, which exhibits a unique papain-like fold with special active site features. Here we investigated the functional significance of the structural elements differentiating MOA from other papain-like cysteine proteases. X-ray crystal structures of MOA co-crystallized with two synthetic substrates reveal cleaved peptides bound to the catalytic site, corresponding to the final products of the proteolytic reaction. Anomalous diffraction data on crystals grown in the presence of calcium and manganese, cadmium or zinc resolve the calcium/manganese preference of MOA and elucidate the inhibitory roles of zinc and cadmium towards papain-like cysteine proteases in general. The reported structures, together with activity data of MOA active site variants, point to a conservation of the general proteolysis mechanism established for papain. Ultimately, the findings suggest that papain and the papain-like domain of MOA are the product of convergent evolution.

A strategy for large-scale comparison of evolutionary- and reaction-based classifications of enzyme function

Determining the molecular function of enzymes discovered by genome sequencing represents a primary foundation for understanding many aspects of biology. Historically, classification of enzyme reactions has used the enzyme nomenclature system developed to describe the overall reactions performed by biochemically characterized enzymes, irrespective of their associated sequences. In contrast, functional classification and assignment for the millions of protein sequences of unknown function now available is largely done in two computational steps, first by similarity-based assignment of newly obtained sequences to homologous groups, followed by transferring to them the known functions of similar biochemically characterized homologs. Due to the fundamental differences in their etiologies and practice, `how' these chemistry- and evolution-centric functional classification systems relate to each other has been difficult to explore on a large scale. To investigate this issue in a new way, we integrated two published ontologies that had previously described each of these classification systems independently. The resulting infrastructure was then used to compare the functional assignments obtained from each classification system for the well-studied and functionally diverse enolase superfamily. Mapping these function assignments to protein structure and reaction similarity networks shows a profound and complex disconnect between the homology- and chemistry-based classification systems. This conclusion mirrors previous observations suggesting that except for closely related sequences, facile annotation transfer from small numbers of characterized enzymes to the huge number uncharacterized homologs to which they are related is problematic. Our extension of these comparisons to large enzyme superfamilies in a computationally intelligent manner provides a foundation for new directions in protein function prediction for the huge proportion of sequences of unknown function represented in major databases. Interactive sequence, reaction, substrate and product similarity networks computed for this work for the enolase and two other superfamilies are freely available for download from the Structure Function Linkage Database Archive (http://sfld.rbvi.ucsf.edu).

Structural comparisons of phosphoenolpyruvate carboxykinases reveal the evolutionary trajectories of these phosphodiester energy conversion enzymes

Inorganic pyrophosphate (PP_i) consists of two phosphate molecules and can act as an energy and phosphate donor in cellular reactions, similar to ATP. Several kinases use PP_i as a substrate, and these kinases have recently been suggested to have evolved from ATP-dependent functional homologs, which have significant amino acid sequence similarity to PP_i-utilizing enzymes. In contrast, phosphoenolpyruvate carboxykinase (PEPCK) can be divided into three types according to the phosphate donor (ATP, GTP, or PPi), and the amino acid sequence similarity of these PEPCKs is too low to confirm that they share a common ancestor. Here we solved the crystal structure of a PP_i-PEPCK homolog from the bacterium Actinomyces israelii at 2.6 Å resolution and compared it with previously reported structures from ATP- and GTP-specific PEPCKs to assess the degrees of similarities and divergences among these PEPCKs. These comparisons revealed that they share a tertiary structure with significant value and that amino acid residues directly contributing to substrate recognition, except for those that recognize purine moieties, are conserved. Furthermore, the order of secondary structural elements between PP_i-, ATP-, and GTP-specific PEPCKs was strictly conserved. The structure-based comparisons of the three PEPCK types provide key insights into the structural basis of PP_i specificity and suggest that all of these PEPCKs are derived from a common ancestor.

Exploring the sequence, function, and evolutionary space of protein superfamilies using sequence similarity networks and phylogenetic reconstructions

Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions by enabling the observation and investigation of complex sequence-structure-function and evolutionary relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) and phylogenetic reconstructions to map the functional divergence and evolutionary history of protein superfamilies. We exemplify this approach using the nitroreductase (NTR) flavoenzyme superfamily, demonstrating that SSN investigations can provide a rapid and effective means to classify groups of proteins, expose sequence similarity relationships across the global scale of a protein superfamily, and efficiently support detailed phylogenetic analyses. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes, their evolution, and their associated physiological roles.

New imine-reducing enzymes from β-hydroxyacid dehydrogenases by single amino acid substitutions

We report the exploration of the evolutionary relationship between imine reductases (IREDs) and other dehydrogenases. This approach is informed by the sequence similarity between these enzyme families and the recently described promiscuous activity of IREDs for the highly reactive carbonyl compound 2,2,2-trifluoroacetophenone. Using the structure of the R-selective IRED from Streptosporangium roseum (R-IRED-Sr) as a model, β-hydroxyacid dehydrogenases (βHADs) were identified as the dehydrogenases most similar to IREDs. To understand how active site differences in IREDs and βHADs enable the reduction of predominantly C = N or C = O bonds respectively, we substituted amino acid residues in βHADs with the corresponding residues from the R-IRED-Sr and were able to increase the promiscuous activity of βHADs for C = N functions by a single amino acid substitution. Variants βHADAt_K170D and βHADAt_K170F lost mainly their keto acid reduction activity and gained the ability to catalyze the reduction of imines. Moreover, the product enantiomeric purity for a bulky imine substrate could be increased from 23% ee (R-IRED-Sr) to 97% ee (βHADAt_K170D/F_F231A) outcompeting already described IRED selectivity.

A platform for chemical modification of mandelate racemase: characterization of the C92S/C264S and γ-thialysine 166 variants

Mandelate racemase (MR) serves as a paradigm for our understanding of enzyme-catalyzed deprotonation of a carbon acid substrate. To facilitate structure-function studies on MR using non-natural amino acid substitutions, we engineered the Cys92Ser/Cys264Ser variant (dmMR) as a platform for introducing Cys residues at specific locations for subsequent covalent modification. While the highly reactive thiol of Cys furnishes a site for chemical modification, site-specificity requires that other Cys residues be non-reactive or replaced by a non-reactive amino acid, especially if chemical modification is conducted under denaturing conditions. The catalytic efficiency of dmMR is reduced only ~2-fold relative to wild-type MR, making dmMR a viable platform for the site-specific introduction of Cys. As an example, the inactive Lys166Cys variant of dmMR was treated with ethylenimine under denaturing conditions to replace the Brønsted acid-base catalyst Lys 166 with the non-natural amino acid γ-thialysine. Comparison of the pH-activity profiles of dmMR and the active γ-thialysine variant revealed a reduction in the pKa for the side chain amino group of ~0.4 units for the latter variant. Unlike wild-type MR for which diffusion is partially rate-limiting, dmMR and the γ-thialysine variant showed no dependence on the solvent viscosity suggesting that the chemical step is fully rate-limiting.

Revealing Unexplored Sequence-Function Space Using Sequence Similarity Networks

The rapidly expanding number of protein sequences found in public databases can improve our understanding of how protein functions evolve. However, our current knowledge of protein function likely represents a small fraction of the diverse repertoire that exists in nature. Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions through the observation and investigation of the complex sequence-structure-function relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) to identify previously unexplored sequence and function space. We exemplify this approach using the nitroreductase (NTR) superfamily. We demonstrate that SSN investigations can provide a rapid and effective means to classify groups of proteins, therefore exposing experimentally unexplored sequences that may exhibit novel functionality. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes and their associated physiological roles.

Functional assignment of multiple catabolic pathways for D-apiose

Colocation of the genes encoding ABC, TRAP, and TCT transport systems and catabolic pathways for the transported ligand provides a strategy for discovering novel microbial enzymes and pathways. We screened solute-binding proteins (SBPs) for ABC transport systems and identified three that bind D-apiose, a branched pentose in the cell walls of higher plants. Guided by sequence similarity networks (SSNs) and genome neighborhood networks (GNNs), the identities of the SBPs enabled the discovery of four catabolic pathways for D-apiose with eleven previously unknown reactions. The new enzymes include D-apionate oxidoisomerase, which catalyzes hydroxymethyl group migration, as well as 3-oxo-isoapionate-4-phosphate decarboxylase and 3-oxo-isoapionate-4-phosphate transcarboxylase/hydrolase, which are RuBisCO-like proteins (RLPs). The web tools for generating SSNs and GNNs are publicly accessible ( http://efi.igb.illinois.edu/efi-est/ ), so similar 'genomic enzymology' strategies for discovering novel pathways can be used by the community.

A continuous assay for l-talarate/galactarate dehydratase using circular dichroism

l-Talarate/galactarate dehydratase (TGD) is a member of the enolase superfamily of enzymes and catalyzes the dehydration of either meso-galactarate or l-talarate to form 5-keto-4-deoxy-d-glucarate (5-KDG). To facilitate study of this enzyme and other galactarate dehydratases, a continuous circular dichroism-based assay has been developed. Using recombinant enzyme from Salmonella typhimurium (StTGD), the rates of StTGD-catalyzed conversion of m-galactarate to 5-KDG were determined by following the change in ellipticity at 323 nm. The apparent molar ellipticity ([θ]₃₂₃) for the 5-KDG formed was determined to be 202 ± 2 deg cm² dmol^-1, which was used to convert observed rates (Δθ/Δt) into concentration-dependent rates (Δc/Δt). The kinetic parameters K_m, k_cat, and k_cat/K_m were 0.38 ± 0.05 mM, 4.8 ± 0.1 s^-1, and 1.3 (±0.2) × 10⁴ M^-1s^-1, respectively. These values are in excellent agreement with those published previously [Yew, W.S. et al. (2007) Biochemistry46, 9564-9577] using a coupled assay system. To demonstrate the utility of the assay, the inhibition constant (K_i = 10.7 ± 0.4 mM) was determined for the competitive inhibitor tartronate. The continuous CD-based assay offers a practical and efficient alternative method to the coupled assay that requires access to 5-KDG aldolase, and to the labor-intensive, fixed-time assays.

Evolutionary and molecular foundations of multiple contemporary functions of the nitroreductase superfamily

Insight regarding how diverse enzymatic functions and reactions have evolved from ancestral scaffolds is fundamental to understanding chemical and evolutionary biology, and for the exploitation of enzymes for biotechnology. We undertook an extensive computational analysis using a unique and comprehensive combination of tools that include large-scale phylogenetic reconstruction to determine the sequence, structural, and functional relationships of the functionally diverse flavin mononucleotide-dependent nitroreductase (NTR) superfamily (>24,000 sequences from all domains of life, 54 structures, and >10 enzymatic functions). Our results suggest an evolutionary model in which contemporary subgroups of the superfamily have diverged in a radial manner from a minimal flavin-binding scaffold. We identified the structural design principle for this divergence: Insertions at key positions in the minimal scaffold that, combined with the fixation of key residues, have led to functional specialization. These results will aid future efforts to delineate the emergence of functional diversity in enzyme superfamilies, provide clues for functional inference for superfamily members of unknown function, and facilitate rational redesign of the NTR scaffold.

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.

Structural and functional innovations in the real-time evolution of new (βα)8 barrel enzymes

New genes can arise by duplication and divergence, but there is a fundamental gap in our understanding of the relationship between these genes, the evolving proteins they encode, and the fitness of the organism. Here we used crystallography, NMR dynamics, kinetics, and mass spectrometry to explain the molecular innovations that arose during a previous real-time evolution experiment. In that experiment, the (βα)₈ barrel enzyme HisA was under selection for two functions (HisA and TrpF), resulting in duplication and divergence of the hisA gene to encode TrpF specialists, HisA specialists, and bifunctional generalists. We found that selection affects enzyme structure and dynamics, and thus substrate preference, simultaneously and sequentially. Bifunctionality is associated with two distinct sets of loop conformations, each essential for one function. We observed two mechanisms for functional specialization: structural stabilization of each loop conformation and substrate-specific adaptation of the active site. Intracellular enzyme performance, calculated as the product of catalytic efficiency and relative expression level, was not linearly related to fitness. Instead, we observed thresholds for each activity above which further improvements in catalytic efficiency had little if any effect on growth rate. Overall, we have shown how beneficial substitutions selected during real-time evolution can lead to manifold changes in enzyme function and bacterial fitness. This work emphasizes the speed at which adaptive evolution can yield enzymes with sufficiently high activities such that they no longer limit the growth of their host organism, and confirms the (βα)₈ barrel as an inherently evolvable protein scaffold.

The interdigitating loop of the enolase superfamily as a specificity binding determinant or 'flying buttress'

BACKGROUND

Enzymes of the enolase superfamily (ENS) are mechanistically diverse, yet share a common partial reaction (abstraction of the α-proton from a carboxylate substrate). While the catalytic machinery responsible for the deprotonation reaction has been conserved, divergent evolution has led to numerous ENS members that catalyze different overall reactions. This rich functional diversity has made the ENS an excellent model system for developing the approaches necessary to validate enzyme function. However, enzymes of the ENS also share a common bidomain structure ((β/α)₇β-barrel domain and α+β capping domain) which makes validation of function from structural information challenging.

SCOPE OF THE REVIEW

This review presents a comparative survey of the structural data obtained over the past decade for enzymes from all seven subgroups that comprise the ENS.

MAJOR CONCLUSIONS

Of the seven ENS subgroups (enolase, mandelate racemase (MR), muconate lactonizing enzyme, β-methylaspartate ammonia lyase, d-glucarate dehydratase, d-mannonate dehydratase (ManD), and galactarate dehydratase 2), only enzymes of the MR and ManD subgroups exhibit an additional feature of structural complexity-an interdigitating loop. This loop emanates from one protomer of a homodimeric pair and penetrates into the adjacent, symmetry-related protomer to either contribute a binding determinant to the active site of the adjacent protomer, or act as a 'flying buttress' to support residues of the active site.

GENERAL SIGNIFICANCE

The analysis presented in this review suggests that the interdigitating loop is the only gross structural element that permits functional distinction between ENS subgroups at the tertiary level of protein structure.

Identification and characterization of two new 5-keto-4-deoxy-D-Glucarate Dehydratases/Decarboxylases

Background

Hexuronic acids such as D-galacturonic acid and D-glucuronic acid can be utilized via different pathways within the metabolism of microorganisms. One representative, the oxidative pathway, generates α-keto-glutarate as the direct link entering towards the citric acid cycle. The penultimate enzyme, keto-deoxy glucarate dehydratase/decarboxylase, catalyses the dehydration and decarboxylation of keto-deoxy glucarate to α-keto-glutarate semialdehyde. This enzymatic reaction can be tracked continuously by applying a pH-shift assay.

Results

Two new keto-deoxy glucarate dehydratases/decarboxylases (EC 4.2.1.41) from Comamonas testosteroni KF-1 and Polaromonas naphthalenivorans CJ2 were identified and expressed in an active form using Escherichia coli ArcticExpress(DE3). Subsequent characterization concerning K_m, k_cat and thermal stability was conducted in comparison with the known keto-deoxy glucarate dehydratase/decarboxylase from Acinetobacter baylyi ADP1. The kinetic constants determined for A. baylyi were K_m 1.0 mM, k_cat 4.5 s⁻¹, for C. testosteroni K_m 1.1 mM, k_cat 3.1 s⁻¹, and for P. naphthalenivorans K_m 1.1 mM, k_cat 1.7 s⁻¹. The two new enzymes had a slightly lower catalytic activity (increased K_m and a decreased k_cat) but showed a higher thermal stability than that of A. baylyi. The developed pH-shift assay, using potassium phosphate and bromothymol blue as the pH indicator, enables a direct measurement. The use of crude extracts did not interfere with the assay and was tested for wild-type landscapes for all three enzymes.

Conclusions

By establishing a pH-shift assay, an easy measurement method for keto-deoxy glucarate dehydratase/decarboxylase could be developed. It can be used for measurements of the purified enzymes or using crude extracts. Therefore, it is especially suitable as the method of choice within an engineering approach for further optimization of these enzymes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-016-0308-3) contains supplementary material, which is available to authorized users.

Large-Scale Analysis Exploring Evolution of Catalytic Machineries and Mechanisms in Enzyme Superfamilies

Enzymes, as biological catalysts, form the basis of all forms of life. How these proteins have evolved their functions remains a fundamental question in biology. Over 100 years of detailed biochemistry studies, combined with the large volumes of sequence and protein structural data now available, means that we are able to perform large-scale analyses to address this question. Using a range of computational tools and resources, we have compiled information on all experimentally annotated changes in enzyme function within 379 structurally defined protein domain superfamilies, linking the changes observed in functions during evolution to changes in reaction chemistry. Many superfamilies show changes in function at some level, although one function often dominates one superfamily. We use quantitative measures of changes in reaction chemistry to reveal the various types of chemical changes occurring during evolution and to exemplify these by detailed examples. Additionally, we use structural information of the enzymes active site to examine how different superfamilies have changed their catalytic machinery during evolution. Some superfamilies have changed the reactions they perform without changing catalytic machinery. In others, large changes of enzyme function, in terms of both overall chemistry and substrate specificity, have been brought about by significant changes in catalytic machinery. Interestingly, in some superfamilies, relatives perform similar functions but with different catalytic machineries. This analysis highlights characteristics of functional evolution across a wide range of superfamilies, providing insights that will be useful in predicting the function of uncharacterised sequences and the design of new synthetic enzymes.

•

Examining how enzyme function evolves using sequence, structure, and reaction mechanism data.

•

Quantifying changes in reaction mechanisms reveals how function has diverged in many superfamilies.

•

Homologous domains frequently use different catalytic residues, which sometimes perform the same enzyme chemistry.

•

This large-scale analysis has significance in protein function prediction and enzyme design.

An additional role for the Brønsted acid-base catalysts of mandelate racemase in transition state stabilization

Mandelate racemase (MR) catalyzes the interconversion of the enantiomers of mandelate and serves as a paradigm for understanding the enzyme-catalyzed abstraction of an α-proton from a carbon acid substrate with a high pKa. The enzyme utilizes a two-base mechanism with Lys 166 and His 297 acting as Brønsted acid and base catalysts, respectively, in the R → S reaction direction. In the S → R reaction direction, their roles are reversed. Using isothermal titration calorimetry (ITC), MR is shown to bind the intermediate/transition state (TS) analogue inhibitor benzohydroxamate (BzH) in an entropy-driven process with a value of ΔCp equal to -358 ± 3 cal mol(-1) K(-1), consistent with an increased number of hydrophobic interactions. However, MR binds BzH with an affinity that is ∼2 orders of magnitude greater than that predicted solely on the basis of hydrophobic interactions [St. Maurice, M., and Bearne, S. L. (2004) Biochemistry 43, 2524], suggesting that additional specific interactions contribute to binding. To test the hypothesis that cation-π/NH-π interactions between the side chains of Lys 166 and His 297 and the aromatic ring and/or the hydroxamate/hydroximate moiety of BzH contribute to the binding of BzH, site-directed mutagenesis was used to generate the MR variants K166M, K166C, H297N, and K166M/H297N and their binding affinity for various ligands determined using ITC. Comparison of the binding affinities of these MR variants with the intermediate/TS analogues BzH and cyclohexanecarbohydroxamate revealed that cation-π/NH-π interactions between His 297 and the hydroxamate/hydroximate moiety and the phenyl ring of BzH contribute approximately 0.26 and 0.91 kcal/mol to binding, respectively, while interactions with Lys 166 contribute approximately 1.74 and 1.74 kcal/mol, respectively. Similarly, comparison of the binding affinities of these mutants with substrate analogues revealed that Lys 166 contributes >2.93 kcal/mol to the binding of (R)-atrolactate, and His 297 contributes 2.46 kcal/mol to the binding of (S)-atrolactate. These results are consistent with Lys 166 and His 297 playing dual roles in catalysis: they act as Brønsted acid-base catalysts, and they stabilize both the enolate moiety and phenyl ring of the altered substrate in the TS.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background

In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to “Domains of Unknown Function” (DUF) or “Uncharacterized Protein Family” (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function.

Results

We applied a ‘computational structural genomics’ approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low- confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659.

Conclusions

This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still ‘non-trivial’ with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners.

Reviewers

This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0069-2) contains supplementary material, which is available to authorized users.

Inactivation of Mandelate Racemase by 3-Hydroxypyruvate Reveals a Potential Mechanistic Link between Enzyme Superfamilies

Mandelate racemase (MR), a member of the enolase superfamily, catalyzes the Mg(2+)-dependent interconversion of the enantiomers of mandelate. Several α-keto acids are modest competitive inhibitors of MR [e.g., mesoxalate (Ki = 1.8 ± 0.3 mM) and 3-fluoropyruvate (Ki = 1.3 ± 0.1 mM)], but, surprisingly, 3-hydroxypyruvate (3-HP) is an irreversible, time-dependent inhibitor (kinact/KI = 83 ± 8 M(-1) s(-1)). Protection from inactivation by the competitive inhibitor benzohydroxamate, trypsinolysis and electrospray ionization tandem mass spectrometry analyses, and X-ray crystallographic studies reveal that 3-HP undergoes Schiff-base formation with Lys 166 at the active site, followed by formation of an aldehyde/enol(ate) adduct. Such a reaction is unprecedented in the enolase superfamily and may be a relic of an activity possessed by a promiscuous progenitor enzyme. The ability of MR to form and deprotonate a Schiff-base intermediate furnishes a previously unrecognized mechanistic link to other α/β-barrel enzymes utilizing Schiff-base chemistry and is in accord with the sequence- and structure-based hypothesis that members of the metal-dependent enolase superfamily and the Schiff-base-forming N-acetylneuraminate lyase superfamily and aldolases share a common ancestor.

Prediction of substrate specificity and preliminary kinetic characterization of the hypothetical protein PVX_123945 from Plasmodium vivax

Members of the haloacid dehalogenase (HAD) superfamily are emerging as an important group of enzymes by virtue of their role in diverse chemical reactions. In different Plasmodium species their number varies from 16 to 21. One of the HAD superfamily members, PVX_123945, a hypothetical protein from Plasmodium vivax, was selected for examining its substrate specificity. Based on distant homology searches and structure comparisons, it was predicted to be a phosphatase. Thirty-eight metabolites were screened to identify potential substrates. Further, to validate the prediction, biochemical and kinetic studies were carried out that showed that the protein was a monomer with high catalytic efficiency for β-glycerophosphate followed by pyridoxal 5'-phosphate. The enzyme also exhibited moderate catalytic efficiencies for α-glycerophosphate, xanthosine 5'-monophosphate and adenosine 5'-monophosphate. It also hydrolyzed the artificial substrate p-nitrophenyl phosphate (pNPP). Mg(2+) was the most preferred divalent cation and phosphate inhibited the enzyme activity. The study is the first attempt at understanding the substrate specificity of a hypothetical protein belonging to HAD superfamily from the malarial parasite P. vivax.

Using the structure-function linkage database to characterize functional domains in enzymes

The Structure-Function Linkage Database (SFLD; http://sfld.rbvi.ucsf.edu/) is a Web-accessible database designed to link enzyme sequence, structure, and functional information. This unit describes the protocols by which a user may query the database to predict the function of uncharacterized enzymes and to correct misannotated functional assignments. The information in this unit is especially useful in helping a user discriminate functional capabilities of a sequence that is only distantly related to characterized sequences in publicly available databases.

Identification of Functionally Related Enzymes by Learning-to-Rank Methods

Enzyme sequences and structures are routinely used in the biological sciences as queries to search for functionally related enzymes in online databases. To this end, one usually departs from some notion of similarity, comparing two enzymes by looking for correspondences in their sequences, structures or surfaces. For a given query, the search operation results in a ranking of the enzymes in the database, from very similar to dissimilar enzymes, while information about the biological function of annotated database enzymes is ignored. In this work, we show that rankings of that kind can be substantially improved by applying kernel-based learning algorithms. This approach enables the detection of statistical dependencies between similarities of the active cleft and the biological function of annotated enzymes. This is in contrast to search-based approaches, which do not take annotated training data into account. Similarity measures based on the active cleft are known to outperform sequence-based or structure-based measures under certain conditions. We consider the Enzyme Commission (EC) classification hierarchy for obtaining annotated enzymes during the training phase. The results of a set of sizeable experiments indicate a consistent and significant improvement for a set of similarity measures that exploit information about small cavities in the surface of enzymes.

One origin for metallo-β-lactamase activity, or two? An investigation assessing a diverse set of reconstructed ancestral sequences based on a sample of phylogenetic trees

Bacteria use metallo-β-lactamase enzymes to hydrolyse lactam rings found in many antibiotics, rendering them ineffective. Metallo-β-lactamase activity is thought to be polyphyletic, having arisen on more than one occasion within a single functionally diverse homologous superfamily. Since discovery of multiple origins of enzymatic activity conferring antibiotic resistance has broad implications for the continued clinical use of antibiotics, we test the hypothesis of polyphyly further; if lactamase function has arisen twice independently, the most recent common ancestor (MRCA) is not expected to possess lactam-hydrolysing activity. Two major problems present themselves. Firstly, even with a perfectly known phylogeny, ancestral sequence reconstruction is error prone. Secondly, the phylogeny is not known, and in fact reconstructing a single, unambiguous phylogeny for the superfamily has proven impossible. To obtain a more statistical view of the strength of evidence for or against MRCA lactamase function, we reconstructed a sample of 98 MRCAs of the metallo-β-lactamases, each based on a different tree in a bootstrap sample of reconstructed phylogenies. InterPro sequence signatures and homology modelling were then used to assess our sample of MRCAs for lactamase functionality. Only 5 % of these models conform to our criteria for metallo-β-lactamase functionality, suggesting that the ancestor was unlikely to have been a metallo-β-lactamase. On the other hand, given that ancestral proteins may have had metallo-β-lactamase functionality with variation in sequence and structural properties compared with extant enzymes, our criteria are conservative, estimating a lower bound of evidence for metallo-β-lactamase functionality but not an upper bound.

New insights about enzyme evolution from large scale studies of sequence and structure relationships

Understanding how enzymes have evolved offers clues about their structure-function relationships and mechanisms. Here, we describe evolution of functionally diverse enzyme superfamilies, each representing a large set of sequences that evolved from a common ancestor and that retain conserved features of their structures and active sites. Using several examples, we describe the different structural strategies nature has used to evolve new reaction and substrate specificities in each unique superfamily. The results provide insight about enzyme evolution that is not easily obtained from studies of one or only a few enzymes.

Identification of the in vivo function of the high-efficiency D-mannonate dehydratase in Caulobacter crescentus NA1000 from the enolase superfamily

The d-mannonate dehydratase (ManD) subgroup of the enolase superfamily contains members with varying catalytic activities (high-efficiency, low-efficiency, or no activity) that dehydrate d-mannonate and/or d-gluconate to 2-keto-3-deoxy-d-gluconate [Wichelecki, D. J., et al. (2014) Biochemistry53, 2722–2731]. Despite extensive in vitro characterization, the in vivo physiological role of a ManD has yet to be established. In this study, we report the in vivo functional characterization of a high-efficiency ManD from Caulobacter crescentus NA1000 (UniProt entry B8GZZ7) by in vivo discovery of its essential role in d-glucuronate metabolism. This in vivo functional annotation may be extended to ∼50 additional proteins.

The natural history of biocatalytic mechanisms

Phylogenomic analysis of the occurrence and abundance of protein domains in proteomes has recently showed that the α/β architecture is probably the oldest fold design. This holds important implications for the origins of biochemistry. Here we explore structure-function relationships addressing the use of chemical mechanisms by ancestral enzymes. We test the hypothesis that the oldest folds used the most mechanisms. We start by tracing biocatalytic mechanisms operating in metabolic enzymes along a phylogenetic timeline of the first appearance of homologous superfamilies of protein domain structures from CATH. A total of 335 enzyme reactions were retrieved from MACiE and were mapped over fold age. We define a mechanistic step type as one of the 51 mechanistic annotations given in MACiE, and each step of each of the 335 mechanisms was described using one or more of these annotations. We find that the first two folds, the P-loop containing nucleotide triphosphate hydrolase and the NAD(P)-binding Rossmann-like homologous superfamilies, were α/β architectures responsible for introducing 35% (18/51) of the known mechanistic step types. We find that these two oldest structures in the phylogenomic analysis of protein domains introduced many mechanistic step types that were later combinatorially spread in catalytic history. The most common mechanistic step types included fundamental building blocks of enzyme chemistry: "Proton transfer," "Bimolecular nucleophilic addition," "Bimolecular nucleophilic substitution," and "Unimolecular elimination by the conjugate base." They were associated with the most ancestral fold structure typical of P-loop containing nucleotide triphosphate hydrolases. Over half of the mechanistic step types were introduced in the evolutionary timeline before the appearance of structures specific to diversified organisms, during a period of architectural diversification. The other half unfolded gradually after organismal diversification and during a period that spanned ∼2 billion years of evolutionary history.

Discovery of function in the enolase superfamily: D-mannonate and d-gluconate dehydratases in the D-mannonate dehydratase subgroup

The continued increase in the size of the protein sequence databases as a result of advances in genome sequencing technology is overwhelming the ability to perform experimental characterization of function. Consequently, functions are assigned to the vast majority of proteins via automated, homology-based methods, with the result that as many as 50% are incorrectly annotated or unannotated (Schnoes et al. PLoS Comput. Biol.2009, 5 (12), e100060520011109). This manuscript describes a study of the d-mannonate dehydratase (ManD) subgroup of the enolase superfamily (ENS) to investigate how function diverges as sequence diverges. Previously, one member of the subgroup had been experimentally characterized as ManD [dehydration of d-mannonate to 2-keto-3-deoxy-d-mannonate (equivalently, 2-keto-3-deoxy-d-gluconate)]. In this study, 42 additional members were characterized to sample sequence–function space in the ManD subgroup. These were found to differ in both catalytic efficiency and substrate specificity: (1) high efficiency (k_cat/K_M = 10³ to 10⁴ M^–1 s^–1) for dehydration of d-mannonate, (2) low efficiency (k_cat/K_M = 10¹ to 10² M^–1 s^–1) for dehydration of d-mannonate and/or D-gluconate, and 3) no-activity with either d-mannonate or d-gluconate (or any other acid sugar tested). Thus, the ManD subgroup is not isofunctional and includes d-gluconate dehydratases (GlcDs) that are divergent from the GlcDs that have been characterized in the mandelate racemase subgroup of the ENS (Lamble et al. FEBS Lett.2004, 576, 133–13615474024) (Ahmed et al. Biochem. J.2005, 390, 529–54015869466). These observations signal caution for functional assignment based on sequence homology and lay the foundation for the studies of the physiological functions of the GlcDs and the promiscuous ManDs/GlcDs.

A two-electron-shell game: intermediates of the extradiol-cleaving catechol dioxygenases

Extradiol-cleaving catechol dioxygenases function by binding both the organic substrate and O2 at a divalent metal center in the active site. They have proven to be a particularly versatile group of enzymes with which to study the O2 activation process. Here, recent studies of homoprotocatechuate 2,3-dioxygenase are summarized, showing how nature can utilize the enzyme structure and the properties of the metal and the substrate to select among many possible chemical paths to achieve both specificity and efficiency. Possible intermediates in the mechanism have been trapped by swapping active-site metals, introducing active-site amino acid substituted variants, and using substrates with different electron-donating capacities. Although each of these intermediates could form part of a viable reaction pathway, kinetic measurements significantly limit the likely candidates. Structural, kinetic, spectroscopic, and computational analyses of the various intermediates shed light on how catalytic efficiency can be achieved.

Discovery of new enzymes and metabolic pathways by using structure and genome context

Assigning valid functions to proteins identified in genome projects is challenging: overprediction and database annotation errors are the principal concerns. We and others are developing computation-guided strategies for functional discovery with 'metabolite docking' to experimentally derived or homology-based three-dimensional structures. Bacterial metabolic pathways often are encoded by 'genome neighbourhoods' (gene clusters and/or operons), which can provide important clues for functional assignment. We recently demonstrated the synergy of docking and pathway context by 'predicting' the intermediates in the glycolytic pathway in Escherichia coli. Metabolite docking to multiple binding proteins and enzymes in the same pathway increases the reliability of in silico predictions of substrate specificities because the pathway intermediates are structurally similar. Here we report that structure-guided approaches for predicting the substrate specificities of several enzymes encoded by a bacterial gene cluster allowed the correct prediction of the in vitro activity of a structurally characterized enzyme of unknown function (PDB 2PMQ), 2-epimerization of trans-4-hydroxy-L-proline betaine (tHyp-B) and cis-4-hydroxy-D-proline betaine (cHyp-B), and also the correct identification of the catabolic pathway in which Hyp-B 2-epimerase participates. The substrate-liganded pose predicted by virtual library screening (docking) was confirmed experimentally. The enzymatic activities in the predicted pathway were confirmed by in vitro assays and genetic analyses; the intermediates were identified by metabolomics; and repression of the genes encoding the pathway by high salt concentrations was established by transcriptomics, confirming the osmolyte role of tHyp-B. This study establishes the utility of structure-guided functional predictions to enable the discovery of new metabolic pathways.

Evolution of a designed retro-aldolase leads to complete active site remodeling

Evolutionary advances are often fueled by unanticipated innovation. Directed evolution of a computationally designed enzyme suggests that dramatic molecular changes can also drive the optimization of primitive protein active sites. The specific activity of an artificial retro-aldolase was boosted >4,400 fold by random mutagenesis and screening, affording catalytic efficiencies approaching those of natural enzymes. However, structural and mechanistic studies reveal that the engineered catalytic apparatus, consisting of a reactive lysine and an ordered water molecule, was unexpectedly abandoned in favor of a new lysine residue in a substrate binding pocket created during the optimization process. Structures of the initial in silico design, a mechanistically promiscuous intermediate, and one of the most evolved variants highlight the importance of loop mobility and supporting functional groups in the emergence of the new catalytic center. Such internal competition between alternative reactive sites may have characterized the early evolution of many natural enzymes.

Evolutionary history and biotechnological future of carboxylases

Carbon dioxide (CO2) is a potent greenhouse gas whose presence in the atmosphere is a critical factor for global warming. At the same time atmospheric CO2 is also a cheap and readily available carbon source that can in principle be used to synthesize value-added products. However, as uncatalyzed chemical CO2-fixation reactions usually require quite harsh conditions to functionalize the CO2 molecule, not many processes have been developed that make use of CO2. In contrast to synthetical chemistry, Nature provides a multitude of different carboxylating enzymes whose carboxylating principle(s) might be exploited in biotechnology. This review focuses on the biochemical features of carboxylases, highlights possible evolutionary scenarios for the emergence of their reactivity, and discusses current, as well as potential future applications of carboxylases in organic synthesis, biotechnology and synthetic biology.

Conformational plasticity of an enzyme during catalysis: intricate coupling between cyclophilin A dynamics and substrate turnover

Enzyme catalysis is central to almost all biochemical processes, speeding up rates of reactions to biological relevant timescales. Enzymes make use of a large ensemble of conformations in recognizing their substrates and stabilizing the transition states, due to the inherent dynamical nature of biomolecules. The exact role of these diverse enzyme conformations and the interplay between enzyme conformational dynamics and catalysis is, according to the literature, not well understood. Here, we use molecular dynamics simulations to study human cyclophilin A (CypA), in order to understand the role of enzyme motions in the catalytic mechanism and recognition. Cyclophilin A is a tractable model system to study using classical simulation methods, because catalysis does not involve bond formation or breakage. We show that the conformational dynamics of active site residues of substrate-bound CypA is inherent in the substrate-free enzyme. CypA interacts with its substrate via conformational selection as the configurations of the substrate changes during catalysis. We also show that, in addition to tight intermolecular hydrophobic interactions between CypA and the substrate, an intricate enzyme-substrate intermolecular hydrogen-bonding network is extremely sensitive to the configuration of the substrate. These enzyme-substrate intermolecular interactions are loosely formed when the substrate is in the reactant and product states and become well formed and reluctant to break when the substrate is in the transition state. Our results clearly suggest coupling among enzyme-substrate intermolecular interactions, the dynamics of the enzyme, and the chemical step. This study provides further insights into the mechanism of peptidyl-prolyl cis/trans isomerases and the general interplay between enzyme conformational dynamics and catalysis.

ECOH: an enzyme commission number predictor using mutual information and a support vector machine

MOTIVATION

The enzyme nomenclature system, commonly known as the enzyme commission (EC) number, plays a key role in classifying and predicting enzymatic reactions. However, numerous reactions have been described in various pathways that do not have an official EC number, and the reactions are not expected to have an EC number assigned because of a lack of articles published on enzyme assays. To predict the EC number of a non-classified enzymatic reaction, we focus on the structural similarity of its substrate and product to the substrate and product of reactions that have been classified.

RESULTS

We propose a new method to assign EC numbers using a maximum common substructure algorithm, mutual information and a support vector machine, termed the Enzyme COmmission numbers Handler (ECOH). A jack-knife test shows that the sensitivity, precision and accuracy of the method in predicting the first three digits of the official EC number (i.e. the EC sub-subclass) are 86.1%, 87.4% and 99.8%, respectively. We furthermore demonstrate that, by examining the ranking in the candidate lists of EC sub-subclasses generated by the algorithm, the method can successfully predict the classification of 85 enzymatic reactions that fall into multiple EC sub-subclasses. The better performance of the ECOH as compared with existing methods and its flexibility in predicting EC numbers make it useful for predicting enzyme function.

AVAILABILITY

ECOH is freely available via the Internet at http://www.bioinfo.sk.ritsumei.ac.jp/apps/ecoh/. This program only works on 32-bit Windows.

Is EC class predictable from reaction mechanism?

BACKGROUND

We investigate the relationships between the EC (Enzyme Commission) class, the associated chemical reaction, and the reaction mechanism by building predictive models using Support Vector Machine (SVM), Random Forest (RF) and k-Nearest Neighbours (kNN). We consider two ways of encoding the reaction mechanism in descriptors, and also three approaches that encode only the overall chemical reaction. Both cross-validation and also an external test set are used.

RESULTS

The three descriptor sets encoding overall chemical transformation perform better than the two descriptions of mechanism. SVM and RF models perform comparably well; kNN is less successful. Oxidoreductases and hydrolases are relatively well predicted by all types of descriptor; isomerases are well predicted by overall reaction descriptors but not by mechanistic ones.

CONCLUSIONS

Our results suggest that pairs of similar enzyme reactions tend to proceed by different mechanisms. Oxidoreductases, hydrolases, and to some extent isomerases and ligases, have clear chemical signatures, making them easier to predict than transferases and lyases. We find evidence that isomerases as a class are notably mechanistically diverse and that their one shared property, of substrate and product being isomers, can arise in various unrelated ways.The performance of the different machine learning algorithms is in line with many cheminformatics applications, with SVM and RF being roughly equally effective. kNN is less successful, given the role that non-local information plays in successful classification. We note also that, despite a lack of clarity in the literature, EC number prediction is not a single problem; the challenge of predicting protein function from available sequence data is quite different from assigning an EC classification from a cheminformatics representation of a reaction.

Structure of mandelate racemase with bound intermediate analogues benzohydroxamate and cupferron

Mandelate racemase (MR, EC 5.1.2.2) from Pseudomonas putida catalyzes the Mg(2+)-dependent interconversion of the enantiomers of mandelate, stabilizing the altered substrate in the transition state by 26 kcal/mol relative to the substrate in the ground state. To understand the origins of this binding discrimination, we determined the X-ray crystal structures of wild-type MR complexed with two analogues of the putative aci-carboxylate intermediate, benzohydroxamate and Cupferron, to 2.2-Å resolution. Benzohydroxamate is shown to be a reasonable mimic of the transition state and/or intermediate because its binding affinity for 21 MR variants correlates well with changes in the free energy of transition state stabilization afforded by these variants. Both benzohydroxamate and Cupferron chelate the active site divalent metal ion and are bound in a conformation with the phenyl ring coplanar with the hydroxamate and diazeniumdiolate moieties, respectively. Structural overlays of MR complexed with benzohydroxamate, Cupferron, and the ground state analogue (S)-atrolactate reveal that the para carbon of the substrate phenyl ring moves by 0.8-1.2 Å between the ground state and intermediate state, consistent with the proposal that the phenyl ring moves during MR catalysis while the polar groups remain relatively fixed. Although the overall protein structure of MR with bound intermediate analogues is very similar to that of MR with bound (S)-atrolactate, the intermediate-Mg(2+) distance becomes shorter, suggesting a tighter complex with the catalytic Mg(2+). In addition, Tyr 54 moves closer to the phenyl ring of the bound intermediate analogues, contributing to an overall constriction of the active site cavity. However, site-directed mutagenesis experiments revealed that the role of Tyr 54 in MR catalysis is relatively minor, suggesting that alterations in enzyme structure that contribute to discrimination between the altered substrate in the transition state and the ground state by this proficient enzyme are extremely subtle.

The Enzyme Function Initiative

The Enzyme Function Initiative (EFI) was recently established to address the challenge of assigning reliable functions to enzymes discovered in bacterial genome projects; in this Current Topic, we review the structure and operations of the EFI. The EFI includes the Superfamily/Genome, Protein, Structure, Computation, and Data/Dissemination Cores that provide the infrastructure for reliably predicting the in vitro functions of unknown enzymes. The initial targets for functional assignment are selected from five functionally diverse superfamilies (amidohydrolase, enolase, glutathione transferase, haloalkanoic acid dehalogenase, and isoprenoid synthase), with five superfamily specific Bridging Projects experimentally testing the predicted in vitro enzymatic activities. The EFI also includes the Microbiology Core that evaluates the in vivo context of in vitro enzymatic functions and confirms the functional predictions of the EFI. The deliverables of the EFI to the scientific community include (1) development of a large-scale, multidisciplinary sequence/structure-based strategy for functional assignment of unknown enzymes discovered in genome projects (target selection, protein production, structure determination, computation, experimental enzymology, microbiology, and structure-based annotation), (2) dissemination of the strategy to the community via publications, collaborations, workshops, and symposia, (3) computational and bioinformatic tools for using the strategy, (4) provision of experimental protocols and/or reagents for enzyme production and characterization, and (5) dissemination of data via the EFI's Website, http://enzymefunction.org. The realization of multidisciplinary strategies for functional assignment will begin to define the full metabolic diversity that exists in nature and will impact basic biochemical and evolutionary understanding, as well as a wide range of applications of central importance to industrial, medicinal, and pharmaceutical efforts.