Topological variation in the evolution of new reactions in functionally diverse enzyme superfamilies

Distinct or Overlapping Areas of Mitochondrial Thioredoxin 2 May Be Used for Its Covalent and Strong Non-Covalent Interactions with Protein Ligands

In silico approaches were employed to examine the characteristics of interactions between human mitochondrial thioredoxin 2 (HsTrx2) and its 38 previously identified mitochondrial protein ligands. All interactions appeared driven mainly by electrostatic forces. The statistically significant residues of HsTrx2 for interactions were characterized as "contact hot spots". Since these were identical/adjacent to putative thermodynamic hot spots, an energy network approach identified their neighbors to highlight possible contact interfaces. Three distinct areas for binding emerged: (i) one around the active site for covalent interactions, (ii) another antipodal to the active site for strong non-covalent interactions, and (iii) a third area involved in both kinds of interactions. The contact interfaces of HsTrx2 were projected as respective interfaces for Escherichia coli Trx1 (EcoTrx1), 2, and HsTrx1. Comparison of the interfaces and contact hot spots of HsTrx2 to the contact residues of EcoTx1 and HsTrx1 from existing crystal complexes with protein ligands supported the hypothesis, except for a part of the cleft/groove adjacent to Trp30 preceding the active site. The outcomes of this study raise the possibility for the rational design of selective inhibitors for the interactions of HsTrx2 with specific protein ligands without affecting the entirety of the functions of the Trx system.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Hydrogen-Deuterium Exchange within Adenosine Deaminase, a TIM Barrel Hydrolase, Identifies Networks for Thermal Activation of Catalysis

Proteins are intrinsically flexible macromolecules that undergo internal motions with time scales spanning femtoseconds to milliseconds. These fluctuations are implicated in the optimization of reaction barriers for enzyme catalyzed reactions. Time, temperature, and mutation dependent hydrogen-deuterium exchange coupled to mass spectrometry (HDX-MS) has been previously employed to identify spatially resolved, catalysis-linked dynamical regions of enzymes. We now extend this technique to pursue the correlation of protein flexibility and chemical reactivity within the diverse and widespread TIM barrel proteins, targeting murine adenosine deaminase (mADA) that catalyzes the irreversible deamination of adenosine to inosine and ammonia. Following a structure-function analysis of rate and activation energy for a series of mutations at a second sphere phenylalanine positioned in proximity to the bound substrate, the catalytically impaired Phe61Ala with an elevated activation energy (Ea = 7.5 kcal/mol) and the wild type (WT) mADA (Ea = 5.0 kcal/mol) were selected for HDX-MS experiments. The rate constants and activation energies of HDX for peptide segments are quantified and used to assess mutation-dependent changes in local and distal motions. Analyses reveal that approximately 50% of the protein sequence of Phe61Ala displays significant changes in the temperature dependence of HDX behaviors, with the dominant change being an increase in protein flexibility. Utilizing Phe61Ile, which displays the same activation energy for k_cat as WT, as a control, we were able to further refine the HDX analysis, highlighting the regions of mADA that are altered in a functionally relevant manner. A map is constructed that illustrates the regions of protein that are proposed to be essential for the thermal optimization of active site configurations that dominate reaction barrier crossings in the native enzyme.

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 Structures of L-DOPA Dioxygenase from Streptomyces sclerotialus

Extradiol dioxygenases are essential biocatalysts for breaking down catechols. The vicinal oxygen chelate (VOC) superfamily contains a large number of extradiol dioxygenases, most of which are found as part of catabolic pathways degrading a variety of natural and human-made aromatic rings. The l-3,4-dihydroxyphenylalanine (L-DOPA) extradiol dioxygenases compose a multitude of pathways that produce various antibacterial or antitumor natural products. The structural features of these dioxygenases are anticipated to be distinct from those of other VOC extradiol dioxygenases. Herein, we identified a new L-DOPA dioxygenase from the thermophilic bacterium Streptomyces sclerotialus (SsDDO) through a sequence and genome context analysis. The activity of SsDDO was kinetically characterized with L-DOPA using an ultraviolet-visible spectrophotometer and an oxygen electrode. The optimal temperature of the assay was 55 °C, at which the K_m and k_cat of SsDDO were 110 ± 10 μM and 2.0 ± 0.1 s^-1, respectively. We determined the de novo crystal structures of SsDDO in the ligand-free form and as a substrate-bound complex, refined to 1.99 and 2.31 Å resolution, respectively. These structures reveal that SsDDO possesses a form IV arrangement of βαβββ modules, the first characterization of this assembly from among the VOC/type I extradiol dioxygenase protein family. Electron paramagnetic resonance spectra of Fe-NO adducts for the resting and substrate-bound enzyme were obtained. This work contributes to our understanding of a growing class of topologically distinct VOC dioxygenases, and the obtained structural features will improve our understanding of the extradiol cleavage reaction within the VOC superfamily.

Genetic, structural, and functional analysis of pathogenic variations causing methylmalonyl-CoA epimerase deficiency

Human methylmalonyl-CoA epimerase (MCEE) catalyzes the interconversion of d-methylmalonyl-CoA and l-methylmalonyl-CoA in propionate catabolism. Autosomal recessive pathogenic variations in MCEE reportedly cause methylmalonic aciduria (MMAuria) in eleven patients. We investigated a cohort of 150 individuals suffering from MMAuria of unknown origin, identifying ten new patients with pathogenic variations in MCEE. Nine patients were homozygous for the known nonsense variation p.Arg47* (c.139C > T), and one for the novel missense variation p.Ile53Arg (c.158T > G). To understand better the molecular basis of MCEE deficiency, we mapped p.Ile53Arg, and two previously described pathogenic variations p.Lys60Gln and p.Arg143Cys, onto our 1.8 Å structure of wild-type (wt) human MCEE. This revealed potential dimeric assembly disruption by p.Ile53Arg, but no clear defects from p.Lys60Gln or p.Arg143Cys. We solved the structure of MCEE-Arg143Cys to 1.9 Å and found significant disruption of two important loop structures, potentially impacting surface features as well as the active-site pocket. Functional analysis of MCEE-Ile53Arg expressed in a bacterial recombinant system as well as patient-derived fibroblasts revealed nearly undetectable soluble protein levels, defective globular protein behavior, and using a newly developed assay, lack of enzymatic activity - consistent with misfolded protein. By contrast, soluble protein levels, unfolding characteristics and activity of MCEE-Lys60Gln were comparable to wt, leaving unclear how this variation may cause disease. MCEE-Arg143Cys was detectable at comparable levels to wt MCEE, but had slightly altered unfolding kinetics and greatly reduced activity. These studies reveal ten new patients with MCEE deficiency and rationalize misfolding and loss of activity as molecular defects in MCEE-type MMAuria.

•

Identification of ten new patients with MCEE-type methylmalonic aciduria.

•

Crystal structures of MCEE wild-type and p.Arg143Cys differ in local conformations.

•

Novel biochemical assay of methylmalonyl-CoA epimerase activity.

Informing Efforts to Develop Nitroreductase for Amine Production

Nitroreductases (NRs) hold promise for converting nitroaromatics to aromatic amines. Nitroaromatic reduction rate increases with Hammett substituent constant for NRs from two different subgroups, confirming substrate identity as a key determinant of reactivity. Amine yields were low, but compounds yielding amines tend to have a large π system and electron withdrawing substituents. Therefore, we also assessed the prospects of varying the enzyme. Several different subgroups of NRs include members able to produce aromatic amines. Comparison of four NR subgroups shows that they provide contrasting substrate binding cavities with distinct constraints on substrate position relative to the flavin. The unique architecture of the NR dimer produces an enormous contact area which we propose provides the stabilization needed to offset the costs of insertion of the active sites between the monomers. Thus, we propose that the functional diversity included in the NR superfamily stems from the chemical versatility of the flavin cofactor in conjunction with a structure that permits tremendous active site variability. These complementary properties make NRs exceptionally promising enzymes for development for biocatalysis in prodrug activation and conversion of nitroaromatics to valuable aromatic amines. We provide a framework for identifying NRs and substrates with the greatest potential to advance.

Cupincin: A Unique Protease Purified from Rice (Oryza sativa L.) Bran Is a New Member of the Cupin Superfamily

Cupin superfamily is one of the most diverse super families. This study reports the purification and characterization of a novel cupin domain containing protease from rice bran for the first time. Hypothetical protein OsI_13867 was identified and named as cupincin. Cupincin was purified to 4.4 folds with a recovery of 4.9%. Cupincin had an optimum pH and temperature of pH 4.0 and 60°C respectively. Cupincin was found to be a homotrimer, consisting of three distinct subunits with apparent molecular masses of 33.45 kDa, 22.35 kDa and 16.67 kDa as determined by MALDI-TOF, whereas it eluted as a single unit with an apparent molecular mass of 135.33 ± 3.52 kDa in analytical gel filtration and migrated as a single band in native page, suggesting its homogeneity. Sequence identity of cupincin was deduced by determining the amino-terminal sequence of the polypeptide chains and by and de novo sequencing. For understanding the hydrolysing mechanism of cupincin, its three-dimensional model was developed. Structural analysis indicated that cupincin contains His313, His326 and Glu318 with zinc ion as the putative active site residues, inhibition of enzyme activity by 1,10-phenanthroline and atomic absorption spectroscopy confirmed the presence of zinc ion. The cleavage specificity of cupincin towards oxidized B-chain of insulin was highly specific; cleaving at the Leu₁₅-Tyr₁₆ position, the specificity was also determined using neurotensin as a substrate, where it cleaved only at the Glu₁-Tyr₂ position. Limited proteolysis of the protease suggests a specific function for cupincin. These results demonstrated cupincin as a completely new protease.

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.

The robustness and innovability of protein folds

Assignment of protein folds to functions indicates that >60% of folds carry out one or two enzymatic functions, while few folds, for example, the TIM-barrel and Rossmann folds, exhibit hundreds. Are there structural features that make a fold amenable to functional innovation (innovability)? Do these features relate to robustness--the ability to readily accumulate sequence changes? We discuss several hypotheses regarding the relationship between the architecture of a protein and its evolutionary potential. We describe how, in a seemingly paradoxical manner, opposite properties, such as high stability and rigidity versus conformational plasticity and structural order versus disorder, promote robustness and/or innovability. We hypothesize that polarity--differentiation and low connectivity between a protein's scaffold and its active-site--is a key prerequisite for innovability.

The evolution of enzyme function in the isomerases

The advent of computational approaches to measure functional similarity between enzymes adds a new dimension to existing evolutionary studies based on sequence and structure. This paper reviews research efforts aiming to understand the evolution of enzyme function in superfamilies, presenting a novel strategy to provide an overview of the evolution of enzymes belonging to an individual EC class, using the isomerases as an exemplar.

Large-scale determination of sequence, structure, and function relationships in cytosolic glutathione transferases across the biosphere

Global networks of the cytosolic glutathione S-transferases illuminate sequence-structure-function relationships across more than 13,000 members of this superfamily, including experimental confirmation of enzymatic activity for 82 members and new crystal structures for 27.

The cytosolic glutathione transferase (cytGST) superfamily comprises more than 13,000 nonredundant sequences found throughout the biosphere. Their key roles in metabolism and defense against oxidative damage have led to thousands of studies over several decades. Despite this attention, little is known about the physiological reactions they catalyze and most of the substrates used to assay cytGSTs are synthetic compounds. A deeper understanding of relationships across the superfamily could provide new clues about their functions. To establish a foundation for expanded classification of cytGSTs, we generated similarity-based subgroupings for the entire superfamily. Using the resulting sequence similarity networks, we chose targets that broadly covered unknown functions and report here experimental results confirming GST-like activity for 82 of them, along with 37 new 3D structures determined for 27 targets. These new data, along with experimentally known GST reactions and structures reported in the literature, were painted onto the networks to generate a global view of their sequence-structure-function relationships. The results show how proteins of both known and unknown function relate to each other across the entire superfamily and reveal that the great majority of cytGSTs have not been experimentally characterized or annotated by canonical class. A mapping of taxonomic classes across the superfamily indicates that many taxa are represented in each subgroup and highlights challenges for classification of superfamily sequences into functionally relevant classes. Experimental determination of disulfide bond reductase activity in many diverse subgroups illustrate a theme common for many reaction types. Finally, sequence comparison between an enzyme that catalyzes a reductive dechlorination reaction relevant to bioremediation efforts with some of its closest homologs reveals differences among them likely to be associated with evolution of this unusual reaction. Interactive versions of the networks, associated with functional and other types of information, can be downloaded from the Structure-Function Linkage Database (SFLD; http://sfld.rbvi.ucsf.edu).

Cytosolic glutathione transferases (cytGSTs) are a large and diverse superfamily of enzymes that have important roles in metabolism and defense against oxidative damage. They have been studied for several decades but because of the synthetic nature of the chemicals used to test these proteins to determine if they have cytGST activity, little is known about the physiological reactions and roles of cytGSTs. In this large, collaborative study, we constructed networks where more than 13,000 cytGST sequences were grouped by sequence similarity and then used these networks to prioritize new targets for experimental characterization in relatively unexplored regions of the superfamily. We report here experimental results confirming GST-like activity for 82 of them, along with 37 new three-dimensional molecular structures determined for 27 targets. These new data, along with experimental data previously reported in the literature, were painted onto the networks to generate a global view of their sequence-structure-function relationships. The results show how proteins of both known and unknown function relate to each other across the entire superfamily and illuminate the complex ways in which their variations in sequence and structure affect our ability to predict unknown functional properties.

Enzyme informatics

Over the last 50 years, sequencing, structural biology and bioinformatics have completely revolutionised biomolecular science, with millions of sequences and tens of thousands of three dimensional structures becoming available. The bioinformatics of enzymes is well served by, mostly free, online databases. BRENDA describes the chemistry, substrate specificity, kinetics, preparation and biological sources of enzymes, while KEGG is valuable for understanding enzymes and metabolic pathways. EzCatDB, SFLD and MACiE are key repositories for data on the chemical mechanisms by which enzymes operate. At the current rate of genome sequencing and manual annotation, human curation will never finish the functional annotation of the ever-expanding list of known enzymes. Hence there is an increasing need for automated annotation, though it is not yet widespread for enzyme data. In contrast, functional ontologies such as the Gene Ontology already profit from automation. Despite our growing understanding of enzyme structure and dynamics, we are only beginning to be able to design novel enzymes. One can now begin to trace the functional evolution of enzymes using phylogenetics. The ability of enzymes to perform secondary functions, albeit relatively inefficiently, gives clues as to how enzyme function evolves. Substrate promiscuity in enzymes is one example of imperfect specificity in protein-ligand interactions. Similarly, most drugs bind to more than one protein target. This may sometimes result in helpful polypharmacology as a drug modulates plural targets, but also often leads to adverse side-effects. Many chemoinformatics approaches can be used to model the interactions between druglike molecules and proteins in silico. We can even use quantum chemical techniques like DFT and QM/MM to compute the structural and energetic course of enzyme catalysed chemical reaction mechanisms, including a full description of bond making and breaking.

Evolutionary and structural analyses of mammalian haloacid dehalogenase-type phosphatases AUM and chronophin provide insight into the basis of their different substrate specificities

Mammalian haloacid dehalogenase (HAD)-type phosphatases are an emerging family of phosphatases with important functions in physiology and disease, yet little is known about the basis of their substrate specificity. Here, we characterize a previously unexplored HAD family member (gene annotation, phosphoglycolate phosphatase), which we termed AUM, for aspartate-based, ubiquitous, Mg(2+)-dependent phosphatase. AUM is a tyrosine-specific paralog of the serine/threonine-specific protein and pyridoxal 5'-phosphate-directed HAD phosphatase chronophin. Comparative evolutionary and biochemical analyses reveal that a single, differently conserved residue in the cap domain of either AUM or chronophin is crucial for phosphatase specificity. We have solved the x-ray crystal structure of the AUM cap fused to the catalytic core of chronophin to 2.65 Å resolution and present a detailed view of the catalytic clefts of AUM and chronophin that explains their substrate preferences. Our findings identify a small number of cap domain residues that encode the different substrate specificities of AUM and chronophin.

Brain miffed by macrophage migration inhibitory factor

Macrophage migration inhibitory factor (MIF) is a cytokine which also exhibits enzymatic properties like oxidoreductase and tautomerase. MIF plays a pivotal role in innate and acquired immunity as well as in the neuroendocrine axis. Since it is involved in the pathogenesis of acute and chronic inflammation, neoangiogenesis, and cancer, MIF and its signaling components are considered suitable targets for therapeutic intervention in several fields of medicine. In neurodegenerative and neurooncological diseases, MIF is a highly relevant, but still a hardly investigated mediator. MIF operates via intracellular protein-protein interaction as well as in CD74/CXCR2/CXCR4 receptor-mediated pathways to regulate essential cellular systems such as redox balance, HIF-1, and p53-mediated senescence and apoptosis as well as multiple signaling pathways. Acting as an endogenous glucocorticoid antagonist, MIF thus represents a relevant resistance gene in brain tumor therapies. Alongside this dual action, a functional homolog-annotated D-dopachrome tautomerase/MIF-2 has been uncovered utilizing the same cell surface receptor signaling cascade as MIF. Here we review MIF actions with respect to redox regulation in apoptosis and in tumor growth as well as its extracellular function with a focus on its potential role in brain diseases. We consider the possibility of MIF targeting in neurodegenerative processes and brain tumors by novel MIF-neutralizing approaches.

Divergence and convergence in enzyme evolution

Comparative analysis of the sequences of enzymes encoded in a variety of prokaryotic and eukaryotic genomes reveals convergence and divergence at several levels. Functional convergence can be inferred when structurally distinct and hence non-homologous enzymes show the ability to catalyze the same biochemical reaction. In contrast, as a result of functional diversification, many structurally similar enzyme molecules act on substantially distinct substrates and catalyze diverse biochemical reactions. Here, we present updates on the ATP-grasp, alkaline phosphatase, cupin, HD hydrolase, and N-terminal nucleophile (Ntn) hydrolase enzyme superfamilies and discuss the patterns of sequence and structural conservation and diversity within these superfamilies. Typically, enzymes within a superfamily possess common sequence motifs and key active site residues, as well as (predicted) reaction mechanisms. These observations suggest that the strained conformation (the entatic state) of the active site, which is responsible for the substrate binding and formation of the transition complex, tends to be conserved within enzyme superfamilies. The subsequent fate of the transition complex is not necessarily conserved and depends on the details of the structures of the enzyme and the substrate. This variability of reaction outcomes limits the ability of sequence analysis to predict the exact enzymatic activities of newly sequenced gene products. Nevertheless, sequence-based (super)family assignments and generic functional predictions, even if imprecise, provide valuable leads for experimental studies and remain the best approach to the functional annotation of uncharacterized proteins from new genomes.