Jeremiah D. Farelli

Enzyme Promiscuity: Engine of Evolutionary Innovation*

Catalytic promiscuity and substrate ambiguity are keys to evolvability, which in turn is pivotal to the successful acquisition of novel biological functions. Action on multiple substrates (substrate ambiguity) can be harnessed for performance of functions in the cell that supersede catalysis of a single metabolite. These functions include proofreading, scavenging of nutrients, removal of antimetabolites, balancing of metabolite pools, and establishing system redundancy. In this review, we present examples of enzymes that perform these cellular roles by leveraging substrate ambiguity and then present the structural features that support both specificity and ambiguity. We focus on the phosphatases of the haloalkanoate dehalogenase superfamily and the thioesterases of the hotdog fold superfamily.

Panoramic view of a superfamily of phosphatases through substrate profiling

Significance Here, we examine the activity profile of the haloalkanoic acid dehalogenase (HAD) superfamily by screening a customized library against >200 enzymes from a broad sampling of the superfamily. From this dataset, we can infer the function of nearly 35% of the superfamily. Overall, the superfamily was found to show high substrate ambiguity, with 75% of the superfamily utilizing greater than five substrates. In addition, the HAD members with the least amount of structural accessorization of the Rossmann fold were found to be the most specific, suggesting that elaboration of the core domain may have led to increased substrate range of the superfamily. Large-scale activity profiling of enzyme superfamilies provides information about cellular functions as well as the intrinsic binding capabilities of conserved folds. Herein, the functional space of the ubiquitous haloalkanoate dehalogenase superfamily (HADSF) was revealed by screening a customized substrate library against >200 enzymes from representative prokaryotic species, enabling inferred annotation of ∼35% of the HADSF. An extremely high level of substrate ambiguity was revealed, with the majority of HADSF enzymes using more than five substrates. Substrate profiling allowed assignment of function to previously unannotated enzymes with known structure, uncovered potential new pathways, and identified iso-functional orthologs from evolutionarily distant taxonomic groups. Intriguingly, the HADSF subfamily having the least structural elaboration of the Rossmann fold catalytic domain was the most specific, consistent with the concept that domain insertions drive the evolution of new functions and that the broad specificity observed in HADSF may be a relic of this process.

Divergence of Structure and Function in the Haloacid Dehalogenase Enzyme Superfamily: Bacteroides thetaiotaomicron BT2127 is an Inorganic Pyrophosphatase

The explosion of protein sequence information requires that current strategies for function assignment evolve to complement experimental approaches with computationally based function prediction. This necessitates the development of strategies based on the identification of sequence markers in the form of specificity determinants and a more informed definition of orthologues. Herein, we have undertaken the function assignment of the unknown haloalkanoate dehalogenase superfamily member BT2127 (Uniprot accession code Q8A5 V9) from Bacteroides thetaiotaomicron using an integrated bioinformatics-structure-mechanism approach. The substrate specificity profile and steady-state rate constants of BT2127 (with a k(cat)/K(m) value for pyrophosphate of ~1 × 10(5) M(-1) s(-1)), together with the gene context, support the assigned in vivo function as an inorganic pyrophosphatase. The X-ray structural analysis of wild-type BT2127 and several variants generated by site-directed mutagenesis shows that substrate discrimination is based, in part, on active site space restrictions imposed by the cap domain (specifically by residues Tyr76 and Glu47). Structure-guided site-directed mutagenesis coupled with kinetic analysis of the mutant enzymes identified the residues required for catalysis, substrate binding, and domain-domain association. On the basis of this structure-function analysis, the catalytic residues Asp11, Asp13, Thr113, and Lys147 as well the metal binding residues Asp171, Asn172, and Glu47 were used as markers to confirm BT2127 orthologues identified via sequence searches. This bioinformatic analysis demonstrated that the biological range of BT2127 orthologue is restricted to the phylum Bacteroidetes/Chlorobi. The key structural determinants in the divergence of BT2127 and its closest homologue, β-phosphoglucomutase, control the leaving group size (phosphate vs glucose phosphate) and the position of the Asp acid/base in the open versus closed conformations. HADSF pyrophosphatases represent a third mechanistic and fold type for bacterial pyrophosphatases.

Structure of the Trehalose-6-phosphate Phosphatase from Brugia malayi Reveals Key Design Principles for Anthelmintic Drugs

Parasitic nematodes are responsible for devastating illnesses that plague many of the worlds poorest populations indigenous to the tropical areas of developing nations. Among these diseases is lymphatic filariasis, a major cause of permanent and long-term disability. Proteins essential to nematodes that do not have mammalian counterparts represent targets for therapeutic inhibitor discovery. One promising target is trehalose-6-phosphate phosphatase (T6PP) from Brugia malayi. In the model nematode Caenorhabditis elegans, T6PP is essential for survival due to the toxic effect(s) of the accumulation of trehalose 6-phosphate. T6PP has also been shown to be essential in Mycobacterium tuberculosis. We determined the X-ray crystal structure of T6PP from B. malayi. The protein structure revealed a stabilizing N-terminal MIT-like domain and a catalytic C-terminal C2B-type HAD phosphatase fold. Structure-guided mutagenesis, combined with kinetic analyses using a designed competitive inhibitor, trehalose 6-sulfate, identified five residues important for binding and catalysis. This structure-function analysis along with computational mapping provided the basis for the proposed model of the T6PP-trehalose 6-phosphate complex. The model indicates a substrate-binding mode wherein shape complementarity and van der Waals interactions drive recognition. The mode of binding is in sharp contrast to the homolog sucrose-6-phosphate phosphatase where extensive hydrogen-bond interactions are made to the substrate. Together these results suggest that high-affinity inhibitors will be bi-dentate, taking advantage of substrate-like binding to the phosphoryl-binding pocket while simultaneously utilizing non-native binding to the trehalose pocket. The conservation of the key residues that enforce the shape of the substrate pocket in T6PP enzymes suggest that development of broad-range anthelmintic and antibacterial therapeutics employing this platform may be possible.

Covalent docking predicts substrates for haloalkanoate dehalogenase superfamily phosphatases

Enzyme function prediction remains an important open problem. Though structure-based modeling, such as metabolite docking, can identify substrates of some enzymes, it is ill-suited to reactions that progress through a covalent intermediate. Here we investigated the ability of covalent docking to identify substrates that pass through such a covalent intermediate, focusing particularly on the haloalkanoate dehalogenase superfamily. In retrospective assessments, covalent docking recapitulated substrate binding modes of known cocrystal structures and identified experimental substrates from a set of putative phosphorylated metabolites. In comparison, noncovalent docking of high-energy intermediates yielded nonproductive poses. In prospective predictions against seven enzymes, a substrate was identified for five. For one of those cases, a covalent docking prediction, confirmed by empirical screening, and combined with genomic context analysis, suggested the identity of the enzyme that catalyzes the orphan phosphatase reaction in the riboflavin biosynthetic pathway of Bacteroides.

Substrate deconstruction and the nonadditivity of enzyme recognition.

Predicting substrates for enzymes of unknown function is a major postgenomic challenge. Substrate discovery, like inhibitor discovery, is constrained by our ability to explore chemotypes; it would be expanded by orders of magnitude if reactive sites could be probed with fragments rather than fully elaborated substrates, as is done for inhibitor discovery. To explore the feasibility of this approach, substrates of six enzymes from three different superfamilies were deconstructed into 41 overlapping fragments that were tested for activity or binding. Surprisingly, even those fragments containing the key reactive group had little activity, and most fragments did not bind measurably, until they captured most of the substrate features. Removing a single atom from a recognized substrate could often reduce catalytic recognition by 6 log-orders. To explore recognition at atomic resolution, the structures of three fragment complexes of the β-lactamase substrate cephalothin were determined by X-ray crystallography. Substrate discovery may be difficult to reduce to the fragment level, with implications for function discovery and for the tolerance of enzymes to metabolite promiscuity. Pragmatically, this study supports the development of libraries of fully elaborated metabolites as probes for enzyme function, which currently do not exist.

Structure and Function of Interacting IcmR-IcmQ Domains from a Type IVb Secretion System in Legionella pneumophila

During infection, Legionella pneumophila creates a replication vacuole within eukaryotic cells and this requires a Type IVb secretion system (T4bSS). IcmQ plays a critical role in the translocase and associates with IcmR. In this paper, we show that the N-terminal domain of IcmQ (Qn) mediates self-dimerization, whereas the C-terminal domain with a basic linker promotes membrane association. In addition, the binding of IcmR to IcmQ prevents self-dimerization and also blocks membrane permeabilization. However, IcmR does not completely block membrane binding by IcmQ. We then determined crystal structures of Qn with the interacting region of IcmR. In this complex, each protein forms an alpha-helical hairpin within a parallel four-helix bundle. The amphipathic nature of helices in Qn suggests two possible models for membrane permeabilization by IcmQ. The Rm-Qn structure also suggests how IcmR-like proteins in other L. pneumophila species may interact with their IcmQ partners.

Structure and catalysis in the Escherichia coli hotdog-fold thioesterase paralogs YdiI and YbdB.

Herein, the structural determinants for substrate recognition and catalysis in two hotdog-fold thioesterase paralogs, YbdB and YdiI from Escherichia coli, are identified and analyzed to provide insight into the evolution of biological function in the hotdog-fold enzyme superfamily. The X-ray crystal structures of YbdB and YdiI, in complex with inert substrate analogs, determined in this study revealed the locations of the respective thioester substrate binding sites and the identity of the residues positioned for substrate binding and catalysis. The importance of each of these residues was assessed through amino acid replacements followed by steady-state kinetic analyses of the corresponding site-directed mutants. Transient kinetic and solvent 18O-labeling studies were then carried out to provide insight into the role of Glu63 posited to function as the nucleophile or general base in catalysis. Finally, the structure–function–mechanism profiles of the two paralogs, along with that of a more distant homolog, were compared to identify conserved elements of substrate recognition and catalysis, which define the core traits of the hotdog-fold thioesterase family, as well as structural features that are unique to each thioesterase. Founded on the insight gained from this analysis, we conclude that the promiscuity revealed by in vitro substrate activity determinations, and posited to facilitate the evolution of new biological function, is the product of intrinsic plasticity in substrate binding as well as in the catalytic mechanism.

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously‐folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent‐exposed hydrophobic pocket, a negative replica of the leucine side‐chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ∼3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α‐helix to π‐helix) in order to precisely preserve the wild‐type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side‐chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a proteins fold.

On unsatisfied hydrogen bonds in the N-terminal subdomain of villin headpiece.

Villin headpiece is a small autonomously folding protein that has emerged as a model system for understanding the fundamental tenets governing protein folding. In this communication, we employ NMR and X-ray crystallography to characterize a point mutant, H41F, which retains actin-binding activity, is more thermostable but, interestingly, does not exhibit the partially folded intermediate observed of either wild-type or other similar point mutants.