Theoretical investigation of the reaction mechanism of the dinuclear zinc enzyme dihydroorotase

Cysteine Radical and Glutamate Collaboratively Enable C-H Bond Activation and C-N Bond Cleavage in a Glycyl Radical Enzyme HplG

Hydroxyprolines are abundant in nature and widely utilized by many living organisms. Isomerization of trans-4-hydroxy-d-proline (t4D-HP) to generate 2-amino-4-ketopentanoate has been found to need a glycyl radical enzyme HplG, which catalyzes the cleavage of the C-N bond, while dehydration of trans-4-hydroxy-l-proline involves a homologous enzyme of HplG. Herein, molecular dynamics simulations and quantum mechanics/molecular mechanics (QM/MM) calculations are employed to understand the reaction mechanism of HplG. Two possible reaction pathways of HplG have been explored to decipher the origin of its chemoselectivity. The QM/MM calculations reveal that the isomerization proceeds via an initial hydrogen shift from the Cγ site of t4D-HP to a catalytic cysteine radical, followed by cleavage of the Cδ-N bond in t4D-HP to form a radical intermediate that captures a hydrogen atom from the cysteine. Activation of the Cδ-H bond in t4D-HP to bring about dehydration of t4D-HP possesses an extremely high energy barrier, thus rendering the dehydration pathway implausible in HplG. On the basis of the current calculations, conserved residue Glu429 plays a pivotal role in the isomerization pathway: the hydrogen bonding between it and t4D-HP weakens the hydroxyalkyl Cγ-Hγ bond, and it acts as a proton acceptor to trigger the cleavage of the C-N bond in t4D-HP. Our current QM/MM calculations rationalize the origin of the experimentally observed chemoselectivity of HplG and propose an H-bond-assisted bond activation strategy in radical-containing enzymes. These findings have general implications on radical-mediated enzymatic catalysis and expand our understanding of how nature wisely and selectively activates the C-H bond to modulate catalytic selectivity.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase

Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of β-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-β-amino acids. This enzyme has the potential to be used as a tool in the production of β-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol^-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.

Role of water coordination at zinc binding site and its catalytic pathway of dizinc creatininase: insights from quantum cluster approach

Creatininase is a key enzyme of creatinine-metabolizing pathway in mammals, and has a great potential for diagnostic application. It catalyzes the reversible conversion of creatinine to creatine. Here, we investigated its reaction mechanism with density functional theory in conjunction with the quantum cluster approach. Three reaction pathways in which several possible proton transfers assisted by either His178 or a water ligand to Zn1 (Wat2) or both were considered. DFT calculations reveal, depending on Wat2 coordination mode at Zn1, two competitive ring-opening pathways where His178 playing a central role as a proton shuttle or both His178 and Wat2 serving as a dual catalytic role as a base and an acid, respectively. Three elementary steps were proposed for the reaction: the first involves nucleophilic attack by a bridging hydroxide to the substrate and forms a gem-diolate intermediate, followed by a proton transfer from the gem-diolate to His178 (His178 protonation is a required step for efficient proton transfers). Finally, the second proton transfer from the protonated His178 or Wat2 to the amide of substrate leads to the ring opening. The first proton transfer is the rate-limiting step of the whole reaction, in consistent with previous experimental and computational studies. A detailed understanding of the reaction mechanism of the creatininase enzyme family will also be helpful for developing a biosensor for kidney function.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

How does Mo-dependent perchlorate reductase work in the decomposition of oxyanions?

The mechanisms of Mo-dependent perchlorate reductase (PcrAB)-catalyzed decomposition of perchlorate, bromate, iodate, and nitrate were revealed by density functional calculations.

How does binuclear zinc amidohydrolase FwdA work in the initial step of methanogenesis: From formate to formyl-methanofuran

The initial step of methanogenesis is the fixation of CO₂ to formyl-methanofuran (formyl-MFR) catalyzed by formyl-MFR dehydrogenase, which can be divided into two half reactions. Herein, the second half reaction catalyzed by FwdA (formyl-methanofuran dehydrogenase subunit A), i.e., from formate to formyl-methanofuran, has been investigated using density functional theory and a chemical model based on the X-ray crystal structure. The calculations indicate that, compared with other well-known di-zinc hydrolases, the FwdA reaction employs a reverse mechanism, including the nucleophilic attack of MFR amine on formate carbon leading to a tetrahedral gem-diolate intermediate, two steps of proton transfer from amine to formate moieties assisted by the Asp385, and the CO bond dissociation to form the formyl-MFR product. The second step of proton transfer from the amine moiety to the Asp385 is rate-limiting with an overall barrier of 21.2 kcal/mol. The two zinc ions play an important role in stabilizing the transition states and intermediates, in particular the negative charge at the formate moiety originated from the nucleophilic attack of the MFR amine. The work here appends a crucial piece in the methanogenic mechanistics and advances the understanding of the global carbon cycle.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase

Creatininase catalyzes the conversion of creatinine (a biosensor for kidney function) to creatine via a two-step mechanism: water addition followed by ring opening. Water addition is common to other known cyclic amidohydrolases, but the precise mechanism for ring opening is still under debate. The proton donor in this step is either His178 or a water molecule bound to one of the metal ions, and the roles of His178 and Glu122 are unclear. Here, the two possible reaction pathways have been fully examined by means of combined quantum mechanics/molecular mechanics simulations at the SCC-DFTB/CHARMM22 level of theory. The results indicate that His178 is the main catalytic residue for the whole reaction and explain its role as proton shuttle during the ring-opening step. In the first step, His178 provides electrostatic stabilization to the gem-diolate tetrahedral intermediate. In the second step, His178 abstracts the hydroxyl proton of the intermediate and delivers it to the cyclic amide nitrogen, leading to ring opening. The latter is the rate-limiting step with a free energy barrier of 18.5 kcal/mol, in agreement with the experiment. We find that Glu122 must be protonated during the enzyme reaction, so that it can form a stable hydrogen bond with its neighboring water molecule. Simulations of the E122Q mutant showed that this replacement disrupts the H-bond network formed by three conserved residues (Glu34, Ser78, and Glu122) and water, increasing the energy barrier. Our computational studies provide a comprehensive explanation for previous structural and kinetic observations, including why the H178A mutation causes a complete loss of activity but the E122Q mutation does not.

Evolutionary expansion of the amidohydrolase superfamily in bacteria in response to the synthetic compounds molinate and diuron

The amidohydrolase superfamily has remarkable functional diversity, with considerable structural and functional annotation of known sequences. In microbes, the recent evolution of several members of this family to catalyze the breakdown of environmental xenobiotics is not well understood. An evolutionary transition from binuclear to mononuclear metal ion coordination at the active sites of these enzymes could produce large functional changes such as those observed in nature, but there are few clear examples available to support this hypothesis. To investigate the role of binuclear-mononuclear active-site transitions in the evolution of new function in this superfamily, we have characterized two recently evolved enzymes that catalyze the hydrolysis of the synthetic herbicides molinate (MolA) and phenylurea (PuhB). In this work, the crystal structures, mutagenesis, metal ion analysis, and enzyme kinetics of both MolA and PuhB establish that these enzymes utilize a mononuclear active site. However, bioinformatics and structural comparisons reveal that the closest putative ancestor of these enzymes had a binuclear active site, indicating that a binuclear-mononuclear transition has occurred. These proteins may represent examples of evolution modifying the characteristics of existing catalysts to satisfy new requirements, specifically, metal ion rearrangement leading to large leaps in activity that would not otherwise be possible.

Theoretical study on the chemical mechanism of enoyl-CoA hydratase and the form of inhibitor binding

Enoyl-CoA hydratase (ECH) catalyzes the second step in the vital β-oxidation pathway of fatty acid metabolism. This enzyme catalyzes the syn-addition of a water molecule across the double bond of 4-(N,N-dimethylamino) cinnamoyl-CoA (DAC-CoA). In this work, the reaction mechanisms of ECH were investigated using the density functional theory (DFT) methods. The different protonation states in which the important residues Glu164 and Glu144 are either neutral or ionized were considered. Four models of the active site were designed based on the X-ray crystal structure of the enzyme. The calculations gave strong support to the proposed mechanism and confirmed that both Glu164 and Glu144 are in a deprotonated state in the reaction mechanism of ECH. In addition, we constructed a model of the active site with the inhibitor acetoacetyl-CoA based on the crystal structure. Caomparison of the calculated energy barriers showed that binding of the keto-enol form of the inhibitor is more reasonable than that of the di-keto form in the inhibition process. Moreover, acetoacetyl-CoA was found to exhibit a keto-enol tautomerism when it acts as an inhibitor in the reaction. The present theoretical results indicated that both residues Glu164 and Glu144 are unprotonated in ECH with the substrate bound, while only Glu164 is unprotonated when the inhibitor binds ECH.

The mononuclear metal center of type-I dihydroorotase from Aquifex aeolicus

BACKGROUND

Dihydroorotase (DHO) is a zinc metalloenzyme, although the number of active site zinc ions has been controversial. E. coli DHO was initially thought to have a mononuclear metal center, but the subsequent X-ray structure clearly showed two zinc ions, α and β, at the catalytic site. Aquifex aeolicus DHO, is a dodecamer comprised of six DHO and six aspartate transcarbamoylase (ATC) subunits. The isolated DHO monomer, which lacks catalytic activity, has an intact α-site and conserved β-site ligands, but the geometry of the second metal binding site is completely disrupted. However, the putative β-site is restored when the complex with ATC is formed and DHO activity is regained. Nevertheless, the X-ray structure of the complex revealed a single zinc ion at the active site. The structure of DHO from the pathogenic organism, S. aureus showed that it also has a single active site metal ion.

RESULTS

Zinc analysis showed that the enzyme has one zinc/DHO subunit and the addition of excess metal ion did not stimulate catalytic activity, nor alter the kinetic parameters. The metal free apoenzyme was inactive, but the full activity was restored upon the addition of one equivalent of Zn2+ or Co2+. Moreover, deletion of the β-site by replacing the His180 and His232 with alanine had no effect on catalysis in the presence or absence of excess zinc. The 2.2 Å structure of the double mutant confirmed that the β-site was eliminated but that the active site remained otherwise intact.

CONCLUSIONS

Thus, kinetically competent A. aeolicus DHO has a mononuclear metal center. In contrast, elimination of the putative second metal binding site in amidohydrolyases with a binuclear metal center, resulted in the abolition of catalytic activity. The number of active site metal ions may be a consideration in the design of inhibitors that selectively target either the mononuclear or binuclear enzymes.

A dominant homolytic O-Cl bond cleavage with low-spin triplet-state Fe(IV)=O formed is revealed in the mechanism of heme-dependent chlorite dismutase

Chlorite dismutase (Cld) is a heme-dependent enzyme that catalyzes the decomposition of toxic chlorite (ClO2(-)) into innocuous chloride and O2. In this paper, using the hybrid B3LYP density functional theory (DFT) method including dispersion interactions, the Cld reaction mechanism has been studied with a chemical model constructed on the X-ray crystal structure. The calculations indicate that the reaction proceeds along a stepwise pathway in the doublet state, i.e. a homolytic O-Cl bond cleavage of the substrate leading to an O-Fe(heme) species and a ClO˙ radical, followed by a rebinding O-O bond formation between them. The O-Fe(heme) species is demonstrated to be a low-spin triplet-state Fe(IV)=O diradicaloid. A low-spin singlet-state Fe(IV)=O is much less stable than the former, with an energy difference of 9.2 kcal mol(-1). The O-Cl bond cleavage is rate-limiting with a barrier of 10.6 kcal mol(-1), in good agreement with the experimental reaction rate of 2.0 × 10(5) s(-1). Furthermore, a heterolytic O-Cl bond dissociation in the initial step is shown to be unreachable, which ensures the high efficiency of the Cld enzyme by avoiding the generation of chlorate byproduct observed in the reactions of synthetic Fe porphyrins. Also, the pathways in the quartet and sextet states are unfavorable for the Cld reaction. The present results reveal a detailed mechanism III (defined in the text) including an interesting di-radical intermediate composed of a low-spin triplet-state Fe(IV)=O and a ClO˙ radical. Compared to a competitive heterolytic Cl-O cleavage in synthetic Fe porphyrins, the revelation of the domination of homolysis in Cld indicates not only the high efficiency of enzyme, but also the sensitivity of a heme and the significance of the enzymatic active-site surroundings (the His170 and Arg183 residues in the present case), which gives more insights into heme chemistry.

A colorimetric assay optimization for high-throughput screening of dihydroorotase by detecting ureido groups

Dihydroorotase (DHOase) is the third enzyme in the de novo pyrimidine biosynthesis pathway and is a potential new antibacterial drug target. No target-based high-throughput screening (HTS) assay for this enzyme has been reported to date. Here, we optimized two colorimetric-based enzymatic assays that detect the ureido moiety of the DHOase substrate, carbamyl-aspartate (Ca-asp). Each assay was developed in a 40-μl assay volume using 384-well plates with a different color mix, diacetylmonoxime (DAMO)-thiosemicarbazide (TSC) or DAMO-antipyrine. The sensitivity and color interference of both color mixes were compared in the presence of common HTS buffer additives, including dimethyl sulfoxide, reducing agents, detergents, and bovine serum albumin. DAMO-TSC (Z'-factors 0.7-0.8) was determined to be superior to DAMO-antipyrine (Z'-factors 0.5-0.6) with significantly less variability within replicates. An HTS pilot screening with 29,552 compounds from four structurally diverse libraries confirmed the quality of our newly optimized colorimetric assay with DAMO-TSC. This robust method has no heating requirement, which was the main obstacle to applying previous assays to HTS. More important, this well-optimized HTS assay for DHOase, the first of its kind, should make it possible to screen large-scale compound libraries to develop new inhibitors against any enzymes that produce ureido functional groups.

High-level expression, purification, and characterization of Staphylococcus aureus dihydroorotase (PyrC) as a cleavable His-SUMO fusion

Staphylococcus aureus is a pathogenic bacterium that causes a variety of mild to lethal human diseases. The rapid spread of multidrug-resistant strains makes the discovery of new antimicrobial agents critical. Dihydroorotase (PyrC), the third enzyme in the bacterial pyrimidine biosynthesis pathway, is structurally and mechanistically distinct from its mammalian counterpart. It has been confirmed to be essential in S. aureus making it an attractive antibacterial drug target. No protocol to express and purify S. aureus PyrC (SaPyrC) has been reported. To obtain the SaPyrC enzyme and overcome anticipated solubility problems, the SaPyrC gene was cloned into the pET-SUMO vector. The N-terminal His-SUMO fused SaPyrC was expressed in Escherichia coli BL21 (DE3) with an HRV 3C protease recognition site inserted between the SUMO tag and SaPyrC to allow for improved cleavage by HRV protease. Purification of cleaved protein using HisTrap affinity and gel filtration columns resulted in native SaPyrC with estimated 95% purity and 40% yield. Both His-SUMO tagged and native SaPyrC form dimers, and enzyme characterization studies have shown that the His-SUMO tag affects enzyme activity slightly. Forward and reverse kinetic rate constants for both tagged and native SaPyrC were determined, and pH profiling studies revealed the optimal pH values for forward and reverse reactions.

The mechanism of cyclization in chromophore maturation of green fluorescent protein: a theoretical study

An intriguing aspect of the green fluorescent protein (GFP) is the autocatalytic post-translational modification that results in the formation of its chromophore. Numerous experimental and theoretical studies indicate that cyclization is the first and the most important step in the maturation process. In this work, two proposed mechanisms for the cyclization were investigated by using the hybrid density functional theory method B3LYP. Cluster models corresponding to the two mechanisms proposed by Wachter et al. [J. Biol. Chem. 2005, 280, 26248-26255] are constructed on the basis of the X-ray crystal structure (PDB entry 2AWJ) and corresponding reaction path potential energy profiles for the two cyclization mechanisms are presented. Our results suggest that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 alpha-carbon. Moreover, Arg96 fulfills the role of stabilizing the enolate moiety, and Glu222 plays the role of a general base. The formation of the cyclized product is found to be 16.0 and 18.6 kcal/mol endothermic with respect to the two models, which is in agreement with experimental observation.

Phosphate mono- and diesterase activities of the trinuclear zinc enzyme nuclease P1--insights from quantum chemical calculations

Nuclease P1 is a trinuclear zinc enzyme that catalyzes the hydrolysis of single-stranded DNA and RNA. Density functional calculations are used to elucidate the reaction mechanism of this enzyme with a model of the active site designed on the basis of the X-ray crystal structure. 2-Tetrahydrofuranyl phosphate and methyl 2-tetrahydrofuranyl phosphate substrates are used to explore the phosphomonoesterase and phosphodiesterase activities of this enzyme, respectively. The calculations reveal that for both activities, a bridging hydroxide performs an in-line attack on the phosphorus center, resulting in inversion of the configuration. Simultaneously, the P-O bond is cleaved, and Zn2 stabilizes the negative charge of the leaving alkoxide anion and assists its departure. All three zinc ions, together with Arg48, provide electrostatic stabilization to the penta-coordinated transition state, thereby lowering the reaction barrier.

Quantum chemical modeling of enzymatic reactions: the case of histone lysine methyltransferase

Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N-terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p).

Reaction mechanism of cis-3-chloroacrylic acid dehalogenase: a theoretical study

The reaction mechanism of cis-3-chloroacrylic acid dehalogenase (cis-CaaD) is studied using the B3LYP density functional theory method. This enzyme catalyzes the hydrolytic dehalogenation of cis-3-chloroacrylic acid to yield malonate semialdehyde and HCl. The uncatalyzed reaction is first considered, and excellent agreement is found between the calculated barrier and the measured rate constant. The enzymatic reaction is then studied with an active site model consisting of 159 atoms. The results suggest an alternative mechanism for cis-CaaD catalysis and different roles for some active site residues in this mechanism.

Target selection and annotation for the structural genomics of the amidohydrolase and enolase superfamilies

To study the substrate specificity of enzymes, we use the amidohydrolase and enolase superfamilies as model systems; members of these superfamilies share a common TIM barrel fold and catalyze a wide range of chemical reactions. Here, we describe a collaboration between the Enzyme Specificity Consortium (ENSPEC) and the New York SGX Research Center for Structural Genomics (NYSGXRC) that aims to maximize the structural coverage of the amidohydrolase and enolase superfamilies. Using sequence- and structure-based protein comparisons, we first selected 535 target proteins from a variety of genomes for high-throughput structure determination by X-ray crystallography; 63 of these targets were not previously annotated as superfamily members. To date, 20 unique amidohydrolase and 41 unique enolase structures have been determined, increasing the fraction of sequences in the two superfamilies that can be modeled based on at least 30% sequence identity from 45% to 73%. We present case studies of proteins related to uronate isomerase (an amidohydrolase superfamily member) and mandelate racemase (an enolase superfamily member), to illustrate how this structure-focused approach can be used to generate hypotheses about sequence-structure-function relationships.