Structural constraints in protein engineering--the coenzyme specificity of Escherichia coli isocitrate dehydrogenase

Two Different Isocitrate Dehydrogenases from Pseudomonas aeruginosa: Enzymology and Coenzyme-Evolutionary Implications

Pseudomonas aeruginosa PAO1, as an experimental model for Gram-negative bacteria, harbors two NADP+-dependent isocitrate dehydrogenases (NADP-IDHs) that were evolved from its ancient counterpart NAD-IDHs. For a better understanding of PaIDH1 and PaIDH2, we cloned the genes, overexpressed them in Escherichia coli and purified them to homogeneity. PaIDH1 displayed higher affinity to NADP+ and isocitrate, with lower Km values when compared to PaIDH2. Moreover, PaIDH1 possessed higher temperature tolerance (50 °C) and wider pH range tolerance (7.2-8.5) and could be phosphorylated. After treatment with the bifunctional PaIDH kinase/phosphatase (PaIDH K/P), PaIDH1 lost 80% of its enzymatic activity in one hour due to the phosphorylation of Ser115. Small-molecule compounds like glyoxylic acid and oxaloacetate can effectively inhibit the activity of PaIDHs. The mutant PaIDH1-D346I347A353K393 exhibited enhanced affinity for NAD+ while it lost activity towards NADP+, and the Km value (7770.67 μM) of the mutant PaIDH2-L589 I600 for NADP+ was higher than that observed for NAD+ (5824.33 μM), indicating a shift in coenzyme specificity from NADP+ to NAD+ for both PaIDHs. The experiments demonstrated that the mutation did not alter the oligomeric state of either protein. This study provides a foundation for the elucidation of the evolution and function of two NADP-IDHs in the pathogenic bacterium P. aeruginosa.

Biochemical Characterization and Crystal Structure of a Novel NAD+-Dependent Isocitrate Dehydrogenase from Phaeodactylum tricornutum

The marine diatom Phaeodactylum tricornutum originated from a series of secondary symbiotic events and has been used as a model organism for studying diatom biology. A novel type II homodimeric isocitrate dehydrogenase from P. tricornutum (PtIDH1) was expressed, purified, and identified in detail through enzymatic characterization. Kinetic analysis showed that PtIDH1 is NAD⁺-dependent and has no detectable activity with NADP⁺. The catalytic efficiency of PtIDH1 for NAD⁺ is 0.16 μM^-1·s^-1 and 0.09 μM^-1·s^-1 in the presence of Mn²⁺ and Mg²⁺, respectively. Unlike other bacterial homodimeric NAD-IDHs, PtIDH1 activity was allosterically regulated by the isocitrate. Furthermore, the dimeric structure of PtIDH1 was determined at 2.8 Å resolution, and each subunit was resolved into four domains, similar to the eukaryotic homodimeric NADP-IDH in the type II subfamily. Interestingly, a unique and novel C-terminal EF-hand domain was first defined in PtIDH1. Deletion of this domain disrupted the intact dimeric structure and activity. Mutation of the four Ca²⁺-binding sites in the EF-hand significantly reduced the calcium tolerance of PtIDH1. Thus, we suggest that the EF-hand domain could be involved in the dimerization and Ca²⁺-coordination of PtIDH1. The current report, on the first structure of type II eukaryotic NAD-IDH, provides new information for further investigation of the evolution of the IDH family.

Structure of a highly NADP+-specific isocitrate dehydrogenase

Isocitrate dehydrogenase catalyzes the first oxidative and decarboxylation steps in the citric acid cycle. It also lies at a crucial bifurcation point between CO2-generating steps in the cycle and carbon-conserving steps in the glyoxylate bypass. Hence, the enzyme is a focus of regulation. The bacterial enzyme is typically dependent on the coenzyme nicotinamide adenine dinucleotide phosphate. The monomeric enzyme from Corynebacterium glutamicum is highly specific towards this coenzyme and the substrate isocitrate while retaining a high overall efficiency. Here, a 1.9 Å resolution crystal structure of the enzyme in complex with its coenzyme and the cofactor Mg2+ is reported. Coenzyme specificity is mediated by interactions with the negatively charged 2'-phosphate group, which is surrounded by the side chains of two arginines, one histidine and, via a water, one lysine residue, forming ion pairs and hydrogen bonds. Comparison with a previous apoenzyme structure indicates that the binding site is essentially preconfigured for coenzyme binding. In a second enzyme molecule in the asymmetric unit negatively charged aspartate and glutamate residues from a symmetry-related enzyme molecule interact with the positively charged arginines, abolishing coenzyme binding. The holoenzyme from C. glutamicum displays a 36° interdomain hinge-opening movement relative to the only previous holoenzyme structure of the monomeric enzyme: that from Azotobacter vinelandii. As a result, the active site is not blocked by the bound coenzyme as in the closed conformation of the latter, but is accessible to the substrate isocitrate. However, the substrate-binding site is disrupted in the open conformation. Hinge points could be pinpointed for the two molecules in the same crystal, which show a 13° hinge-bending movement relative to each other. One of the two pairs of hinge residues is intimately flanked on both sides by the isocitrate-binding site. This suggests that binding of a relatively small substrate (or its competitive inhibitors) in tight proximity to a hinge point could lead to large conformational changes leading to a closed, presumably catalytically active (or inactive), conformation. It is possible that the small-molecule concerted inhibitors glyoxylate and oxaloacetate similarly bind close to the hinge, leading to an inactive conformation of the enzyme.

Crystal structure of 3-isopropylmalate dehydrogenase in complex with NAD(+) and a designed inhibitor

Isopropylmalate dehydrogenase (IPMDH) is the third enzyme specific to leucine biosynthesis in microorganisms and plants, and catalyzes the oxidative decarboxylation of (2R,3S)-3-isopropylmalate to alpha-ketoisocaproate using NAD(+) as an oxidizing agent. In this study, a thia-analogue of the substrate was designed and synthesized as an inhibitor for IPMDH. The analogue showed strong competitive inhibitory activity with K(i)=62nM toward IPMDH derived from Thermus thermophilus. Moreover, the crystal structure of T. thermophilus IPMDH in a ternary complex with NAD(+) and the inhibitor has been determined at 2.8A resolution. The inhibitor exists as a decarboxylated product with an enol/enolate form in the active site. The product interacts with Arg 94, Asn 102, Ser 259, Glu 270, and a water molecule hydrogen-bonding with Arg 132. All interactions between the product and the enzyme were observed in the position associated with keto-enol tautomerization. This result implies that the tautomerization step of the thia-analogue during the IPMDH reaction is involved in the inhibition.

Biochemical characterization of isocitrate dehydrogenase from Methylococcus capsulatus reveals a unique NAD+-dependent homotetrameric enzyme

The gene encoding isocitrate dehydrogenase (IDH) of Methylococcus capsulatus (McIDH) was cloned and overexpressed in Escherichia coli. The purified enzyme was NAD+-dependent with a thermal optimum for activity at 55-60 degrees C and an apparent midpoint melting temperature (Tm) of 70 degrees C. Analytical ultracentrifugation (AUC) revealed a homotetrameric state, and McIDH thus represents the first homotetrameric NAD+-dependent IDH that has been characterized. Based on a structural alignment of McIDH and homotetrameric homoisocitrate dehydrogenase (HDH) from Thermus thermophilus (TtHDH), we identified the clasp-like domain of McIDH as a likely site for tetramerization. McIDH showed moreover, higher sequence identity (48%) to TtHDH than to previously characterized IDHs. Putative NAD+-IDHs with high sequence identity (48-57%) to McIDH were however identified in a variety of bacteria showing that NAD+-dependent IDHs are indeed widespread within the domain, Bacteria. Phylogenetic analysis including these new sequences revealed a close relationship with eukaryal allosterically regulated NAD+-IDH and the subfamily III of IDH was redefined to include bacterial NAD+- and NADP+-dependent IDHs. This apparent relationship suggests that the mitochondrial genes encoding NAD+-IDH are derived from the McIDH-like IDHs.

Mutation-biased adaptation in a protein NK model

Evolutionary trends responsible for systematic differences in genome and proteome composition have been attributed to GC:AT mutation bias in the context of neutral evolution or to selection acting on genome composition. A possibility that has been ignored, presumably because it is part of neither the Modern Synthesis nor the Neutral Theory, is that mutation may impose a directional bias on adaptation. This possibility is explored here with simulations of the effect of a GC:AT bias on amino acid composition during adaptive walks on an abstract protein fitness landscape called an "NK" model. The results indicate that adaptation does not preclude mutation-biased evolution. In the complete absence of neutral evolution, a modest GC:AT bias of realistic magnitude can displace the trajectory of adaptation in a mutationally favored direction, to such a degree that amino acid composition is biased substantially and persistently. Thus, mutational explanations for evolved patterns need not presuppose neutral evolution.

The crystal structure of a hyperthermostable subfamily II isocitrate dehydrogenase from Thermotoga maritima

Isocitrate dehydrogenase (IDH) from the hyperthermophile Thermotoga maritima (TmIDH) catalyses NADP+- and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. It belongs to the beta-decarboxylating dehydrogenase family and is the only hyperthermostable IDH identified within subfamily II. Furthermore, it is the only IDH that has been characterized as both dimeric and tetrameric in solution. We solved the crystal structure of the dimeric apo form of TmIDH at 2.2 A. The R-factor of the refined model was 18.5% (R(free) 22.4%). The conformation of the TmIDH structure was open and showed a domain rotation of 25-30 degrees compared with closed IDHs. The separate domains were found to be homologous to those of the mesophilic mammalian IDHs of subfamily II and were subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that stabilization of the N- and C-termini via long-range electrostatic interactions were important for the higher thermostability of TmIDH. Moreover, the number of intra- and intersubunit ion pairs was higher and the ionic networks were larger compared with the mesophilic IDHs. Other factors likely to confer higher stability in TmIDH were a less hydrophobic and more charged accessible surface, a more hydrophobic subunit interface, more hydrogen bonds per residue and a few loop deletions. The residues responsible for the binding of isocitrate and NADP+ were found to be highly conserved between TmIDH and the mammalian IDHs and it is likely that the reaction mechanism is the same.

Thr373, Asp375, and Lys260 are in the coenzyme site of porcine NADP-dependent isocitrate dehydrogenase

Thr(373), Lys(374), Asp(375), and Lys(260) were chosen as site-directed mutagenesis targets within porcine NADP-dependent isocitrate dehydrogenase based on structurally corrected sequence alignment among prokaryotic and eukaryotic NADP-isocitrate dehydrogenases. Wild-type and all mutant enzymes were expressed in Escherichia coli and purified to homogeneity. These mutations do not alter the secondary structure or dimerization state of the mutants. The D375N and K260Q mutants exhibit, respectively, a 15- and 28-fold increase in K(m) for NADP, along with marked decreases in V(max) as compared to wild-type enzyme. In contrast, replacing Lys(374), which was previously proposed to contribute to apparent coenzyme affinity, does not change the enzyme's kinetic parameters. T373S exhibits similar kinetic parameters to those of wild-type while T373A and T373V mutations reduce the V(max) values of the resulting enzymes to 1 and 20%, respectively of that of wild-type. We conclude that a hydroxyl group at position 373 is required for effective enzyme function and that Asp(375) and Lys(260) are critical amino acids contributing to coenzyme affinity as well as catalysis by porcine NADP-isocitrate dehydrogenase.

Complete reversal of coenzyme specificity of isocitrate dehydrogenase from Haloferax volcanii

Haloferax volcanii Ds-threo-isocitrate dehydrogenase (ICDH) was highly expressed in bacteria as inclusion bodies. The recombinant enzyme was refolded, purified and characterized, and was found to be NADP-dependent like the wild-type protein. Sequence alignment of several isocitrate dehydrogenases from evolutionarily divergent organisms including H. volcanii revealed that the amino acid residues involved in coenzyme specificity are highly conserved. Our objective was to switch the coenzyme specificity of halophilic ICDH by altering these conserved amino acids. We were able to switch coenzyme specificity from NADP+ to NAD+ by changing five amino acids by site-directed mutagenesis (Arg291, Lys343, Tyr344, Val350 and Tyr390). The five mutants of ICDH were overexpressed in Escherichia coli as inclusion bodies and each recombinant ICDH protein was refolded and purified, and its kinetic parameters were determined. Coenzyme specificity did not switch until all five amino acids were substituted.

Identification of a novel trifunctional homoisocitrate dehydrogenase and modulation of the broad substrate specificity through site-directed mutagenesis

A gene encoding homoisocitrate dehydrogenase (HICDH) of Deinococcus radiodurans was cloned, sequenced, and overexpressed in Escherichia coli. The amino acid sequence was 63% identical to HICDH from Thermus thermophilus (Tth-HICDH). Similar to Tth-HICDH, purified, recombinant Dra-HICDH was a tetramer, required K+ and Mn2+ for activity, and used NAD+ as a coenzyme. However, unlike Tth-HICDH, which has a 20-fold preference for isocitrate over homoisocitrate, Dra-HICDH preferred homoisocitrate to isocitrate by 1.5-fold. Moreover, it catalyzed the oxidation of 3-isopropylmalate, albeit at approximately 0.1% the rate seen with homoisocitrate and isocitrate. Saturation mutagenesis of Dra-HICDH Arg87 was next performed because an orthologous Arg85 to valine mutation in Tth-HICDH results in loss of activity toward isocitrate, but in retention of activity toward homoisocitrate. Unexpectedly, the Arg85Val variant became able to catalyze the oxidation of 3-isopropylmalate. Screening of the saturation mutagenesis library identified two variants, Arg87Val and Arg87Thr, that were able to catalyze the oxidation of homoisocitrate, but not isocitrate or 3-isopropylmalate. Deletion of Dra-HICDH Ala80, a residue missing from Tth-HICDH and predicted to reside at the entrance of alpha-helix Arg87, resulted in alterations in substrate specificity that rendered Dra-HICDH similar to Tth-HICDH; i.e., a 4-fold preference for isocitrate over homoisocitrate and inability to catalyze the oxidation of 3-isopropylmalate. Seemingly minor changes in primary sequence result in changes in substrate specificity of beta-decarboxylating dehydrogenases.

Isocitrate dehydrogenase from the hyperthermophile Aeropyrum pernix: X-ray structure analysis of a ternary enzyme-substrate complex and thermal stability

Isocitrate dehydrogenase from Aeropyrum pernix (ApIDH) is a homodimeric enzyme that belongs to the beta-decarboxylating dehydrogenase family and is the most thermostable IDH identified. It catalyzes the NADP+ and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. We have solved the crystal structures of a native ApIDH at 2.2 A, a pseudo-native ApIDH at 2.1 A, and of ApIDH in complex with NADP+, Ca2+ and d-isocitrate at 2.3 A. The pseudo-native ApIDH is in complex with etheno-NADP+ which was located at the surface instead of in the active site revealing a novel adenine-nucleotide binding site in ApIDH. The native and the pseudo-native ApIDHs were found in an open conformation, whereas one of the subunits of the ternary complex was closed upon substrate binding. The closed subunit showed a domain rotation of 19 degrees compared to the open subunit. The binding of isocitrate in the closed subunit was identical with that of the binary complex of porcine mitochondrial IDH, whereas the binding of NADP+ was similar to that of the ternary complex of IDH from Escherichiacoli. The reaction mechanism is likely to be conserved in the different IDHs. A proton relay chain involving at least five solvent molecules, the 5'-phosphate group of the nicotinamide-ribose and a coupled lysine-tyrosine pair in the active site, is postulated as essential in both the initial and the final steps of the catalytic reaction of IDH. ApIDH was found to be highly homologous to the mesophilic IDHs and was subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that a disulfide bond at the N terminus and a seven-membered inter-domain ionic network at the surface are major determinants for the higher thermostability of ApIDH compared to EcIDH. Furthermore, the total number of ion pairs was dramatically higher in ApIDH compared to the mesophilic IDHs if a cutoff of 4.2 A was used. A calculated net charge of only +1 compared to -19 and -25 in EcIDH and BsIDH, respectively, suggested a high degree of electrostatic optimization, which is known to be an important determinant for increased thermostability.

Directed evolution relieves product inhibition and confers in vivo function to a rationally designed tyrosine aminotransferase

The Escherichia coli aspartate (AATase) and tyrosine (TATase) aminotransferases share 43% sequence identity and 72% similarity, but AATase has only 0.08% and 0.01% of the TATase activities (k(cat)/K(m)) for tyrosine and phenylalanine, respectively. Approximately 5% of TATase activity was introduced into the AATase framework earlier both by rational design (six mutations, termed HEX) and by directed evolution (9-17 mutations). The enzymes realized from the latter procedure complement tyrosine auxotrophy in TATase deficient E. coli. HEX complements even more poorly than does wild-type AATase, even though the (k(cat)/K(m)) value for tyrosine exhibited by HEX is similar to those of the enzymes found from directed evolution. HEX, however, is characterized by very low values of K(m) and K(D) for dicarboxylic ligands, and by a particularly slow release for oxaloacetate, the product of the reaction with aspartate and a TCA cycle intermediate. These observations suggest that HEX exists largely as an enzyme-product complex in vivo. HEX was therefore subjected to a single round of directed evolution with selection for complementation of tyrosine auxotrophy. A variant with a single amino acid substitution, A293D, exhibited substantially improved TATase function in vivo. The A293D mutation alleviates the tight binding to dicarboxylic ligands as K(m)s for aspartate and alpha-ketoglutarate are >20-fold higher in the HEX + A293D construct compared to HEX. This mutation also increased k(cat)/K(m)(Tyr) threefold. A second mutation, I73V, elicited smaller but similar effects. Both residues are in close proximity to Arg292 and the mutations may function to modulate the arginine switch mechanism responsible for dual substrate recognition in TATases and HEX.

Implication by site-directed mutagenesis of Arg314 and Tyr316 in the coenzyme site of pig mitochondrial NADP-dependent isocitrate dehydrogenase

Sequence alignment of pig mitochondrial NADP-dependent isocitrate dehydrogenase with eukaryotic (human, rat, and yeast) and Escherichia coli isocitrate dehydrogenases reveals that Tyr316 is completely conserved and is equivalent to the E. coli Tyr345, which interacts with the 2'-phosphate of NADP in the crystal structure [Hurley et al., Biochemistry 30 (1991) 8671-8678]. Lys321 is also completely conserved in the five isocitrate dehydrogenases. Either an arginine or lysine residue is found among the enzymes from other species at the position corresponding to the pig enzyme Arg314. While Arg323 is not conserved among all species, its proximity to the coenzyme site makes it a good candidate for investigation. The importance of these four amino acids to the function of pig mitochondrial NADP-isocitrate dehydrogenase was studied by site-directed mutagenesis. Mutants (R314Q, Y316F, Y316L, K321Q, and R323Q) were generated by a megaprimer polymerase chain reaction method. Wild-type and mutant enzymes were expressed in E. coli and purified to homogeneity. All mutant and wild-type enzymes exhibited comparable molecular weights indicative of the dimeric enzyme. Mutations do not cause an appreciable change in enzyme secondary structure as revealed by circular dichroism measurements. The kinetic parameters (V(max) and K(M) values) of K321Q and R323Q are similar to those of wild-type, indicating that Lys321 and Arg323 are not involved in enzyme function. R314Q exhibits a 10-fold increase in K(M) for NADP as compared to that of wild-type, while they have comparable V(max) values. These results suggest that Arg314 contributes to the affinity between the enzyme and NADP. The hydroxyl group of Tyr316 is not required for enzyme function since Y316F exhibits similar kinetic parameters to those of wild-type. Y316L shows a 4-fold increase in K(M) for NADP and a decrease in V(max) as compared to wild-type, suggesting that the aromatic ring of the Tyr of isocitrate dehydrogenase contributes to the affinity for coenzyme, as well as to catalysis. The K(i) for NAD of R314Q, Y316F, and Y316L is comparable to that of wild-type, indicating that the Arg314 and Tyr316 may be located near the 2'-phosphate of enzyme-bound NADP.

Enzyme engineering: rational redesign versus directed evolution

Enzyme engineering is undergoing the most profound and exciting transformation in its history (1). It promises unprecedented expansion in the scope and applications of modified or improved enzymes with desired physical and catalytic properties. Two complementary strategies are currently available: rational redesign (2-5) and directed evolution (6,7). Although both approaches have met with great success, each has limitations. In this article, the perspectives for these enzyme-engineering strategies are discussed briefly.

Functional prediction: identification of protein orthologs and paralogs

Orthologs typically retain the same function in the course of evolution. Using beta-decarboxylating dehydrogenase family as a model, we demonstrate that orthologs can be confidently identified. The strategy is based on our recent findings that substitutions of only a few amino acid residues in these enzymes are sufficient to exchange substrate and coenzyme specificities. Hence, the few major specificity determinants can serve as reliable markers for determining orthologous or paralogous relationships. The power of this approach has been demonstrated by correcting similarity-based functional misassignment and discovering new genes and related pathways, and should be broadly applicable to other enzyme families.