Evolutionary Aspects of Enzyme Dynamics

The role of evolutionary pressure on the chemical step catalyzed by enzymes is somewhat enigmatic, in part because chemistry is not rate-limiting for many optimized systems. Herein, we present studies that examine various aspects of the evolutionary relationship between protein dynamics and the chemical step in two paradigmatic enzyme families, dihydrofolate reductases and alcohol dehydrogenases. Molecular details of both convergent and divergent evolution are beginning to emerge. The findings suggest that protein dynamics across an entire enzyme can play a role in adaptation to differing physiological conditions. The growing tool kit of kinetics, kinetic isotope effects, molecular biology, biophysics, and bioinformatics provides means to link evolutionary changes in structure-dynamics function to the vibrational and conformational states of each protein.     

Introduction to the Thematic Minireview Series on Enzyme Evolution

In this thematic minireview series, the JBC presents five provocative articles on Enzyme Evolution. The reviews discuss stimulating concepts that include the emergence of primordial catalysts at temperatures that were considerably warmer than present day ones and the impact of the cooling environment on the evolution of catalytic fitness and the preservation of catalysis-promoting conformational dynamics. They also discuss the use of Urzymes or invariant modules in enzyme superfamilies as paradigms for understanding the evolution of catalytic efficiency and specificity, the use of bioinformatics approaches to understand the roles of substrate ambiguity and catalytic promiscuity as drivers of evolution, and the challenges associated with assigning catalytic function as the number of superfamily members grows rapidly.     

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.     

Escherichia coli d-Malate Dehydrogenase,a Generalist Enzyme Active in the Leucine Biosynthesis Pathway

The enzymes of the β-decarboxylating dehydrogenase superfamily catalyze the oxidative decarboxylation of d-malate-based substrates with various specificities. Here, we show that, in addition to its natural function affording bacterial growth on d-malate as a carbon source, the d-malate dehydrogenase of Escherichia coli (EcDmlA) naturally expressed from its chromosomal gene is capable of complementing leucine auxotrophy in a leuB⁻ strain lacking the paralogous isopropylmalate dehydrogenase enzyme. To our knowledge, this is the first example of an enzyme that contributes with a physiologically relevant level of activity to two distinct pathways of the core metabolism while expressed from its chromosomal locus. EcDmlA features relatively high catalytic activity on at least three different substrates (l(+)-tartrate, d-malate, and 3-isopropylmalate). Because of these properties both in vivo and in vitro, EcDmlA may be defined as a generalist enzyme. Phylogenetic analysis highlights an ancient origin of DmlA, indicating that the enzyme has maintained its generalist character throughout evolution. We discuss the implication of these findings for protein evolution.     

Crystal Structure of the Bacillus subtilis Phosphodiesterase PhoD Reveals an Iron and Calcium-containing Active Site

The PhoD family of extra-cytoplasmic phosphodiesterases are among the most commonly occurring bacterial phosphatases. The exemplars for this family are the PhoD protein of Bacillus subtilis and the phospholipase D of Streptomyces chromofuscus. We present the crystal structure of B. subtilis PhoD. PhoD is most closely related to purple acid phosphatases (PAPs) with both types of enzyme containing a tyrosinate-ligated Fe³⁺ ion. However, the PhoD active site diverges from that found in PAPs and uses two Ca²⁺ ions instead of the single extra Fe²⁺, Mn²⁺, or Zn²⁺ ion present in PAPs. The PhoD crystals contain a phosphate molecule that coordinates all three active site metal ions and that is proposed to represent a product complex. A C-terminal helix lies over the active site and controls access to the catalytic center. The structure of PhoD defines a new phosphatase active site architecture based on Fe³⁺ and Ca²⁺ ions.     

An Acyl-covalent Enzyme Intermediate of Lecithin:Retinol Acyltransferase

Synthesis of fatty acid retinyl esters determines systemic vitamin A levels and provides substrate for production of visual chromophore (11-cis-retinal) in vertebrates. Lecithin:retinol acyltransferase (LRAT), the main enzyme responsible for retinyl ester formation, catalyzes the transfer of an acyl group from the sn-1 position of phosphatidylcholine to retinol. To delineate the catalytic mechanism of this reaction, we expressed and purified a fully active, soluble form of this enzyme and used it to examine the possible formation of a transient acyl-enzyme intermediate. Detailed mass spectrometry analyses revealed that LRAT undergoes spontaneous, covalent modification upon incubation with a variety of phosphatidylcholine substrates. The addition of an acyl chain occurs at the Cys¹⁶¹ residue, indicating formation of a thioester intermediate. This observation provides the first direct experimental evidence of thioester intermediate formation that constitutes the initial step in the proposed LRAT catalytic reaction. Additionally, we examined the effect of increasing fatty acyl side chain length in phosphatidylcholine on substrate accessibility in this reaction, which provided insights into the function of the single membrane-spanning domain of LRAT. These observations are critical to understanding the catalytic mechanism of LRAT protein family members as well as other lecithin:acyltransferases wherein Cys residues are required for catalysis.     

Structural basis for the acyltransferase activity of lecithin:retinol acyltransferase-like proteins

Lecithin:retinol acyltransferase-like proteins, also referred to as HRAS-like tumor suppressors, comprise a vertebrate subfamily of papain-like or NlpC/P60 thiol proteases that function as phospholipid-metabolizing enzymes. HRAS-like tumor suppressor 3, a representative member of this group, plays a key role in regulating triglyceride accumulation and energy expenditure in adipocytes and therefore constitutes a novel pharmacological target for treatment of metabolic disorders causing obesity. Here, we delineate a catalytic mechanism common to lecithin:retinol acyltransferase-like proteins and provide evidence for their alternative robust lipid-dependent acyltransferase enzymatic activity. We also determined high resolution crystal structures of HRAS-like tumor suppressor 2 and 3 to gain insight into their active site architecture. Based on this structural analysis, two conformational states of the catalytic Cys-113 were identified that differ in reactivity and thus could define the catalytic properties of these two proteins. Finally, these structures provide a model for the topology of these enzymes and allow identification of the protein-lipid bilayer interface. This study contributes to the enzymatic and structural understanding of HRAS-like tumor suppressor enzymes.     

Predicting Enzyme Adsorption to Lignin Films by Calculating Enzyme Surface Hydrophobicity

The inhibitory action of lignin on cellulase cocktails is a major challenge to the biological saccharification of plant cell wall polysaccharides. Although the mechanism remains unclear, hydrophobic interactions between enzymes and lignin are hypothesized to drive adsorption. Here we evaluate the role of hydrophobic interactions in enzyme-lignin binding. The hydrophobicity of the enzyme surface was quantified using an estimation of the clustering of nonpolar atoms, identifying potential interaction sites. The adsorption of enzymes to lignin surfaces, measured using the quartz crystal microbalance, correlates to the hydrophobic cluster scores. Further, these results suggest a minimum hydrophobic cluster size for a protein to preferentially adsorb to lignin. The impact of electrostatic contribution was ruled out by comparing the isoelectric point (pI) values to the adsorption of proteins to lignin surfaces. These results demonstrate the ability to predict enzyme-lignin adsorption and could potentially be used to design improved cellulase cocktails, thus lowering the overall cost of biofuel production.     

Comparative Laboratory Evolution of Ordered and Disordered Enzymes

Intrinsically disordered proteins are ubiquitous in nature. To assess potential evolutionary advantages and disadvantages of structural disorder under controlled laboratory conditions, we directly compared the evolvability of weakly active ordered and disordered variants of dihydrofolate reductase by genetic selection. The circularly permuted Escherichia coli enzyme, which exists as a molten globule in the absence of ligands, and a well folded deletion mutant of the Bacillus stearothermophilus enzyme served as starting points. Both scaffolds evolved at similar rates and to similar extents, reaching near-native activity after three rounds of mutagenesis and selection. Surprisingly, however, the starting structural properties of the two scaffolds changed only marginally during optimization. Although the ordered and disordered proteins accumulated distinct sets of mutations, the changes introduced likely improved catalytic efficiency indirectly in both cases by bolstering the network of dynamic conformational fluctuations that productively couple into the reaction coordinate.     

Structural and Functional Evolution of Positively Selected Sites in Pine Glutathione S-Transferase Enzyme Family

Phylogenetic analyses have identified positive selection as an important driver of protein evolution, both structural and functional. However, the lack of appropriate combined functional and structural assays has generally hindered attempts to elucidate patterns of positively selected sites and their effects on enzyme activity and substrate specificity. In this study we investigated the evolutionary divergence of the glutathione S-transferase (GST) family in Pinus tabuliformis, a pine that is widely distributed from northern to central China, including cold temperate and drought-stressed regions. GSTs play important roles in plant stress tolerance and detoxification. We cloned 44 GST genes from P. tabuliformis and found that 26 of the 44 belong to the largest (Tau) class of GSTs and are differentially expressed across tissues and developmental stages. Substitution models identified five positively selected sites in the Tau GSTs. To examine the functional significance of these positively selected sites, we applied protein structural modeling and site-directed mutagenesis. We found that four of the five positively selected sites significantly affect the enzyme activity and specificity; thus their variation broadens the GST family substrate spectrum. In addition, positive selection has mainly acted on secondary substrate binding sites or sites close to (but not directly at) the primary substrate binding site; thus their variation enables the acquisition of new catalytic functions without compromising the protein primary biochemical properties. Our study sheds light on selective aspects of the functional and structural divergence of the GST family in pine and other organisms.     

Studies of Human 2,4-Dienoyl CoA Reductase Shed New Light on Peroxisomal β-Oxidation of Unsaturated Fatty Acids

Peroxisomes play an essential role in maintaining fatty acid homeostasis. Although mitochondria are also known to participate in the catabolism of fatty acids via β-oxidation, differences exist between the peroxisomal and mitochondrial β-oxidation. Only peroxisomes, but not mitochondrion, can shorten very long chain fatty acids. Here, we describe the crystal structure of a ternary complex of peroxisomal 2,4-dienoyl CoA reductases (pDCR) with hexadienoyl CoA and NADP, as a prototype for comparison with the mitochondrial 2,4-dienoyl CoA reductase (mDCR) to shed light on the differences between the enzymes from the two organelles at the molecular level. Unexpectedly, the structure of pDCR refined to 1.84 Å resolution reveals the absence of the tyrosine-serine pair seen in the active site of mDCR, which together with a lysine and an asparagine have been deemed a hallmark of the SDR family of enzymes. Instead, aspartate hydrogen-bonded to the Cα hydroxyl via a water molecule seems to perturb the water molecule for protonation of the substrate. Our studies provide the first structural evidence for participation of water in the DCR-catalyzed reactions. Biochemical studies and structural analysis suggest that pDCRs can catalyze the shortening of six-carbon-long substrates in vitro. However, the K_m values of pDCR for short chain acyl CoAs are at least 6-fold higher than those for substrates with 10 or more aliphatic carbons. Unlike mDCR, hinge movements permit pDCR to process very long chain polyunsaturated fatty acids.     

Divergent Molecular Evolution of the Mitochondrial Sulfhydryl:Cytochrome c Oxidoreductase Erv in Opisthokonts and Parasitic Protists

Mia40 and the sulfhydryl:cytochrome c oxidoreductase Erv1/ALR are essential for oxidative protein import into the mitochondrial intermembrane space in yeast and mammals. Although mitochondrial protein import is functionally conserved in the course of evolution, many organisms seem to lack Mia40. Moreover, except for in organello import studies and in silico analyses, nothing is known about the function and properties of protist Erv homologues. Here we compared Erv homologues from yeast, the kinetoplastid parasite Leishmania tarentolae, and the non-related malaria parasite Plasmodium falciparum. Both parasite proteins have altered cysteine motifs, formed intermolecular disulfide bonds in vitro and in vivo, and could not replace Erv1 from yeast despite successful mitochondrial protein import in vivo. To analyze its enzymatic activity, we established the expression and purification of recombinant full-length L. tarentolae Erv and compared the mechanism with related and non-related flavoproteins. Enzyme assays indeed confirmed an electron transferase activity with equine and yeast cytochrome c, suggesting a conservation of the enzymatic activity in different eukaryotic lineages. However, although Erv and non-related flavoproteins are intriguing examples of convergent molecular evolution resulting in similar enzyme properties, the mechanisms of Erv homologues from parasitic protists and opisthokonts differ significantly. In summary, the Erv-mediated reduction of cytochrome c might be highly conserved throughout evolution despite the apparent absence of Mia40 in many eukaryotes. Nevertheless, the knowledge on mitochondrial protein import in yeast and mammals cannot be generally transferred to all other eukaryotes, and the corresponding pathways, components, and mechanisms remain to be analyzed.     

Biochemical and Structural Studies of 6-Carboxy-5,6,7,8-tetrahydropterin Synthase Reveal the Molecular Basis of Catalytic Promiscuity within the Tunnel-fold Superfamily

6-Pyruvoyltetrahydropterin synthase (PTPS) homologs in both mammals and bacteria catalyze distinct reactions using the same 7,8-dihydroneopterin triphosphate substrate. The mammalian enzyme converts 7,8-dihydroneopterin triphosphate to 6-pyruvoyltetrahydropterin, whereas the bacterial enzyme catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin. To understand the basis for the differential activities we determined the crystal structure of a bacterial PTPS homolog in the presence and absence of various ligands. Comparison to mammalian structures revealed that although the active sites are nearly structurally identical, the bacterial enzyme houses a His/Asp dyad that is absent from the mammalian protein. Steady state and time-resolved kinetic analysis of the reaction catalyzed by the bacterial homolog revealed that these residues are responsible for the catalytic divergence. This study demonstrates how small variations in the active site can lead to the emergence of new functions in existing protein folds.     

A Model of the Membrane-bound Cytochrome b5-Cytochrome P450 Complex from NMR and Mutagenesis Data

Microsomal cytochrome b₅ (cytb₅) is a membrane-bound protein that modulates the catalytic activity of its redox partner, cytochrome P4502B4 (cytP450). Here, we report the first structure of full-length rabbit ferric microsomal cytb₅ (16 kDa), incorporated in two different membrane mimetics (detergent micelles and lipid bicelles). Differential line broadening of the cytb₅ NMR resonances and site-directed mutagenesis data were used to characterize the cytb₅ interaction epitope recognized by ferric microsomal cytP450 (56 kDa). Subsequently, a data-driven docking algorithm, HADDOCK (high ambiguity driven biomolecular docking), was used to generate the structure of the complex between cytP4502B4 and cytb₅ using experimentally derived restraints from NMR, mutagenesis, and the double mutant cycle data obtained on the full-length proteins. Our docking and experimental results point to the formation of a dynamic electron transfer complex between the acidic convex surface of cytb₅ and the concave basic proximal surface of cytP4502B4. The majority of the binding energy for the complex is provided by interactions between residues on the C-helix and β-bulge of cytP450 and residues at the end of helix α4 of cytb₅. The structure of the complex allows us to propose an interprotein electron transfer pathway involving the highly conserved Arg-125 on cytP450 serving as a salt bridge between the heme propionates of cytP450 and cytb₅. We have also shown that the addition of a substrate to cytP450 likely strengthens the cytb₅-cytP450 interaction. This study paves the way to obtaining valuable structural, functional, and dynamic information on membrane-bound complexes.     

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²⁺-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.     

Solution NMR Structure and Functional Analysis of the Integral Membrane Protein YgaP from Escherichia coli

The solution NMR structure of the α-helical integral membrane protein YgaP from Escherichia coli in mixed 1,2-diheptanoyl-sn-glycerol-3-phosphocholine/1-myristoyl-2-hydroxy-sn-glycero-3-phospho-(1′-rac-glycerol) micelles is presented. In these micelles, YgaP forms a homodimer with the two transmembrane helices being the dimer interface, whereas the N-terminal cytoplasmic domain includes a rhodanese-fold in accordance to its sequence homology to the rhodanese family of sulfurtransferases. The enzymatic sulfur transfer activity of full-length YgaP as well as of the N-terminal rhodanese domain only was investigated performing a series of titrations with sodium thiosulfate and potassium cyanide monitored by NMR and EPR. The data indicate the thiosulfate concentration-dependent addition of several sulfur atoms to the catalytic Cys-63, which process can be reversed by the addition of potassium cyanide. The catalytic reaction induces thereby conformational changes within the rhodanese domain, as well as on the transmembrane α-helices of YgaP. These results provide insights into a potential mechanism of YgaP during the catalytic thiosulfate activity in vivo.     

A Triple Mutant in the Ω-loop of TEM-1 β-Lactamase Changes the Substrate Profile via a Large Conformational Change and an Altered General Base for Catalysis

β-Lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. TEM-1 is a prevalent plasmid-encoded β-lactamase in Gram-negative bacteria that efficiently catalyzes the hydrolysis of penicillins and early cephalosporins but not oxyimino-cephalosporins. A previous random mutagenesis study identified a W165Y/E166Y/P167G triple mutant that displays greatly altered substrate specificity with increased activity for the oxyimino-cephalosporin, ceftazidime, and decreased activity toward all other β-lactams tested. Surprisingly, this mutant lacks the conserved Glu-166 residue critical for enzyme function. Ceftazidime contains a large, bulky side chain that does not fit optimally in the wild-type TEM-1 active site. Therefore, it was hypothesized that the substitutions in the mutant expand the binding site in the enzyme. To investigate structural changes and address whether there is an enlargement in the active site, the crystal structure of the triple mutant was solved to 1.44 Å. The structure reveals a large conformational change of the active site Ω-loop structure to create additional space for the ceftazidime side chain. The position of the hydroxyl group of Tyr-166 and an observed shift in the pH profile of the triple mutant suggests that Tyr-166 participates in the hydrolytic mechanism of the enzyme. These findings indicate that the highly conserved Glu-166 residue can be substituted in the mechanism of serine β-lactamases. The results reveal that the robustness of the overall β-lactamase fold coupled with the plasticity of an active site loop facilitates the evolution of enzyme specificity and mechanism.     

Inference of functional properties from large-scale analysis of enzyme superfamilies

As increasingly large amounts of data from genome and other sequencing projects become available, new approaches are needed to determine the functions of the proteins these genes encode. We show how large-scale computational analysis can help to address this challenge by linking functional information to sequence and structural similarities using protein similarity networks. Network analyses using three functionally diverse enzyme superfamilies illustrate the use of these approaches for facile updating and comparison of available structures for a large superfamily, for creation of functional hypotheses for metagenomic sequences, and to summarize the limits of our functional knowledge about even well studied superfamilies.     

Blm10 protein promotes proteasomal substrate turnover by an active gating mechanism

For optimal proteolytic function, the central core of the proteasome (core particle (CP) or 20S) has to associate with activators. We investigated the impact of the yeast activator Blm10 on proteasomal peptide and protein degradation. We found enhanced degradation of peptide substrates in the presence of Blm10 and demonstrated that Blm10 has the capacity to accelerate proteasomal turnover of the unstructured protein tau-441 in vitro. Mechanistically, proteasome activation requires the opening of a closed gate, which allows passage of unfolded proteins into the catalytic chamber. Our data indicate that gate opening by Blm10 is achieved via engagement of its C-terminal segment with the CP. Crucial for this activity is a conserved C-terminal YYX motif, with the penultimate tyrosine playing a preeminent role. Thus, Blm10 utilizes a gate opening strategy analogous to the proteasomal ATPases HbYX-dependent mechanism. Because gating incompetent Blm10 C-terminal point mutants confers a loss of function phenotype, we propose that the cellular function of Blm10 is based on CP association and activation to promote the degradation of proteasome substrates.