Determinants of cofactor specificity in isocitrate dehydrogenase: structure of an engineered NADP+ --> NAD+ specificity-reversal mutant

Distribution and Functional Analysis of Isocitrate Dehydrogenases across Kinetoplastids

Isocitrate dehydrogenase is an enzyme converting isocitrate to α-ketoglutarate in the canonical tricarboxylic acid (TCA) cycle. There are three different types of isocitrate dehydrogenase documented in eukaryotes. Our study points out the complex evolutionary history of isocitrate dehydrogenases across kinetoplastids, where the common ancestor of Trypanosomatidae and Bodonidae was equipped with two isoforms of the isocitrate dehydrogenase enzyme: the NADP+-dependent isocitrate dehydrogenase 1 with possibly dual localization in the cytosol and mitochondrion and NADP+-dependent mitochondrial isocitrate dehydrogenase 2. In the extant trypanosomatids, isocitrate dehydrogenase 1 is present only in a few species suggesting that it was lost upon separation of Trypanosoma spp. and replaced by the mainly NADP+-dependent cytosolic isocitrate dehydrogenase 3 of bacterial origin in all the derived lineages. In this study, we experimentally demonstrate that the omnipresent isocitrate dehydrogenase 2 has a dual localization in both mitochondrion and cytosol in at least four species that possess only this isoform. The apparent lack of the NAD+-dependent isocitrate dehydrogenase activity in trypanosomatid mitochondrion provides further support to the existence of the noncanonical TCA cycle across trypanosomatids and the bidirectional activity of isocitrate dehydrogenase 3 when operating with NADP+ cofactor instead of NAD+. This observation can be extended to all 17 species analyzed in this study, except for Leishmania mexicana, which showed only low isocitrate dehydrogenase activity in the cytosol. The variability in isocitrate oxidation capacity among species may reflect the distinct metabolic strategies and needs for reduced cofactors in particular environments.

Characterization of the Isocitrate Dehydrogenase Gene Family and Their Response to Drought Stress in Maize

Isocitrate dehydrogenase (IDH) is a key rate-limiting enzyme in the tricarboxylic acid cycle and acts in glutamine synthesis. IDH also participates in plant growth and development and in response to abiotic stresses. We identified 11 maize IDH genes (ZmIDH) and classified these genes into ZmNAD-IDH and ZmNADP-IDH groups based on their different coenzymes (NAD+ or NADP+). The ZmNAD-IDH group was further divided into two subgroups according to their catalytic and non-catalytic subunits, as in Arabidopsis. The ZmIDHs significantly differed in physicochemical properties, gene structure, conserved motifs, and protein tertiary structure. Promoter prediction analysis revealed that the promoters of these ZmIDHs contain cis-acting elements associated with light response, abscisic acid, phytohormones, and abiotic stresses. ZmIDH is predicted to interact with proteins involved in development and stress resistance. Expression analysis of public data revealed that most ZmIDHs are specifically expressed in anthers. Different types of ZmIDHs responded to abiotic stresses with different expression patterns, but all exhibited responses to abiotic stresses to some extent. In addition, analysis of the public sequence from transcription data in an association panel suggested that natural variation in ZmIDH1.4 will be associated with drought tolerance in maize. These results suggested that ZmIDHs respond differently and/or redundantly to abiotic stresses during plant growth and development, and this analysis provides a foundation to understand how ZmIDHs respond to drought stress in maize.

NAD(P)H binding configurations revealed by time-resolved fluorescence and two-photon absorption

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional aspects of NADH and NADPH binding and clarify the biochemical processes that underlie their contrasting intracellular lifetimes.

Biochemical and Phylogenetic Characterization of a Novel NADP+-Specific Isocitrate Dehydrogenase From the Marine Microalga Phaeodactylum tricornutum

Isocitrate dehydrogenase (IDH) family of proteins is classified into three subfamilies, namely, types I, II, and III. Although IDHs are widely distributed in bacteria, archaea, and eukaryotes, all type III IDHs reported to date are found only in prokaryotes. Herein, a novel type III IDH subfamily member from the marine microalga Phaeodactylum tricornutum (PtIDH2) was overexpressed, purified, and characterized in detail for the first time. Relatively few eukaryotic genomes encode this type of IDH and PtIDH2 shares the highest homology with marine bacterial monomeric IDHs, suggesting that PtIDH2 originated through a horizontal gene transfer event between a marine alga and a bacterium. Size-exclusion chromatography revealed that the native PtIDH2 is a homotetramer (∼320 kDa) in solution, comprising four monomeric IDH-like subunits (80 kDa each). Enzymatic characterization showed that PtIDH2 is a bivalent metal ion-dependent enzyme and Mn²⁺ is the optimal activator. The recombinant PtIDH2 protein exhibited maximal activity at 35°C and pH 8.0 in the presence of Mn²⁺. Heat-inactivation analysis revealed that PtIDH2 is a cold-adapted enzyme. Kinetic analysis demonstrated that PtIDH2 is a completely NADP⁺-specific IDH with no detectable NAD⁺-associated catalytic activity. The three putative key NADP⁺-binding residues (His604, Arg615, and Arg664) in PtIDH2 were also evaluated by site-directed mutagenesis. The H⁶⁰⁴L/R⁶¹⁵D/R⁶⁶⁴S triple mutant showed a 3.25-fold preference for NAD⁺ over NADP⁺, implying that the coenzyme specificity of PtIDH2 can be converted from NADP⁺ to NAD⁺ through rational engineering approaches. Additionally, the roles of the conserved residues Ala718 and Leu742 in the thermostability of PtIDH2 were also explored by site-directed mutagenesis. We found that the L⁷⁴²F mutant displayed higher thermostability than wild-type PtIDH2. This study expands the phylogeny of the IDH family and provides new insights into the evolution of IDHs.

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.

Crystal Structures of the Putative Isocitrate Dehydrogenase from Sulfolobus tokodaii Strain 7 in the Apo and NADP+-Bound Forms

Isocitrate dehydrogenase is a catabolic enzyme that acts during the third step of the tricarboxylic acid cycle. The hypothetical protein ST2166 from the archaeon Sulfolobus tokodaii was isolated and crystallized. It shares high primary structure homology with prokaryotic NADP⁺-dependent IDHs, suggesting that these enzymes share a common enzymatic mechanism. The crystal structure of ST2166 was determined at 2.0 Å resolution in the apo form, and then the structure of the crystal soaked with NADP⁺ was also determined at 2.4 Å resolution, which contained NADP⁺ bound at the putative active site. Comparisons between the structures of apo and NADP⁺-bound forms and NADP-IDHs from other prokaryotes suggest that prokaryotic NADP-IDHs recognize their cofactors using conserved Lys335, Tyr336, and Arg386 in ST2166 at the opening cleft before the domain closure.

Protein Engineering for Nicotinamide Coenzyme Specificity in Oxidoreductases: Attempts and Challenges

Oxidoreductases are ubiquitous enzymes that catalyze an extensive range of chemical reactions with great specificity, efficiency, and selectivity. Most oxidoreductases are nicotinamide cofactor-dependent enzymes with a strong preference for NADP or NAD. Because these coenzymes differ in stability, bioavailability and costs, the enzyme preference for a specific coenzyme is an important issue for practical applications. Different approaches for the manipulation of coenzyme specificity have been reported, with different degrees of success. Here we present various attempts for the switching of nicotinamide coenzyme preference in oxidoreductases by protein engineering. This review covers 103 enzyme engineering studies from 82 articles and evaluates the accomplishments in terms of coenzyme specificity and catalytic efficiency compared to wild type enzymes of different classes. We analyzed different protein engineering strategies and related them with the degree of success in inverting the cofactor specificity. In general, catalytic activity is compromised when coenzyme specificity is reversed, however when switching from NAD to NADP, better results are obtained. In most of the cases, rational strategies were used, predominantly with loop exchange generating the best results. In general, the tendency of removing acidic residues and incorporating basic residues is the strategy of choice when trying to change specificity from NAD to NADP, and vice versa. Computational strategies and algorithms are also covered as helpful tools to guide protein engineering strategies. This mini review aims to give a general introduction to the topic, giving an overview of tools and information to work in protein engineering for the reversal of coenzyme specificity.

A receptor dependent-4D QSAR approach to predict the activity of mutated enzymes

Screening and selection tools to obtain focused libraries play a key role in successfully engineering enzymes of desired qualities. The quality of screening depends on efficient assays; however, a focused library generated with a priori information plays a major role in effectively identifying the right enzyme. As a proof of concept, for the first time, receptor dependent - 4D Quantitative Structure Activity Relationship (RD-4D-QSAR) has been implemented to predict kinetic properties of an enzyme. The novelty of this study is that the mutated enzymes also form a part of the training data set. The mutations were modeled in a serine protease and molecular dynamics simulations were conducted to derive enzyme-substrate (E-S) conformations. The E-S conformations were enclosed in a high resolution grid consisting of 156,250 grid points that stores interaction energies to generate QSAR models to predict the enzyme activity. The QSAR predictions showed similar results as reported in the kinetic studies with >80% specificity and >50% sensitivity revealing that the top ranked models unambiguously differentiated enzymes with high and low activity. The interaction energy descriptors of the best QSAR model were used to identify residues responsible for enzymatic activity and substrate specificity.

Switch in Cofactor Specificity of a Baeyer-Villiger Monooxygenase

Baeyer-Villiger monooxygenases (BVMOs) catalyze the oxidation of ketones to esters or lactones by using molecular oxygen and a cofactor. Type I BVMOs display a strong preference for NADPH. However, for industrial purposes NADH is the preferred cofactor, as it is ten times cheaper and more stable. Thus, we created a variant of the cyclohexanone monooxygenase from Acinetobacter sp. NCIMB 9871 (CHMO_Acineto ); this used NADH 4200-fold better than NADPH. By combining structure analysis, sequence alignment, and literature data, 21 residues in proximity of the cofactor were identified and targeted for mutagenesis. Two combinatorial variants bearing three or four mutations showed higher conversions of cyclohexanone with NADH (79 %) compared to NADPH (58 %) as well as specificity. The structural reasons for this switch in cofactor specificity of a type I BVMO are especially a hydrogen-bond network coordinating the two hydroxy groups of NADH through direct interactions and bridging water molecules.

Fluorescence-Activated Cell Sorting of Human l-asparaginase Mutant Libraries for Detecting Enzyme Variants with Enhanced Activity

Immunogenicity is one of the most common complications occurring during therapy making use of protein drugs of nonhuman origin. A notable example of such a case is bacterial l-asparaginases (L-ASNases) used for the treatment of acute lymphoblastic leukemia (ALL). The replacement of the bacterial enzymes by human ones is thought to set the basis for a major improvement of antileukemic therapy. Recently, we solved the crystal structure of a human enzyme possessing L-ASNase activity, designated hASNase-3. This enzyme is expressed as an inactive precursor protein and post-translationally undergoes intramolecular processing leading to the generation of two subunits which remain noncovalently, yet tightly associated and constitute the catalytically active form of the enzyme. We discovered that this intramolecular processing can be drastically and selectively accelerated by the free amino acid glycine. In the present study, we report on the molecular engineering of hASNase-3 aiming at the improvement of its catalytic properties. We created a fluorescence-activated cell sorting (FACS)-based high-throughput screening system for the characterization of rationally designed mutant libraries, capitalizing on the finding that free glycine promotes autoproteolytic cleavage, which activates the mutant proteins expressed in an E. coli strain devoid of aspartate biosynthesis. Successive screening rounds led to the isolation of catalytically improved variants showing up to 6-fold better catalytic efficiency as compared to the wild-type enzyme. Our work establishes a powerful strategy for further exploitation of the human asparaginase sequence space to facilitate the identification of in vitro-evolved enzyme species that will lay the basis for improved ALL therapy.

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second β-strand (β2-Asp/Glu). We identified a canonical motif, defined by the β2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose's ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint--geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the β2-Asp/Glu motif.

Biochemical Characterization and Complete Conversion of Coenzyme Specificity of Isocitrate Dehydrogenase from Bifidobacterium longum

Bifidobacterium longum is a very important gram-positive non-pathogenic bacterium in the human gastrointestinal tract for keeping the digestive and immune system healthy. Isocitrate dehydrogenase (IDH) from B. longum (BlIDH), a novel member in Type II subfamily, was overexpressed, purified and biochemically characterized in detail. The active form of BlIDH was an 83-kDa homodimer. Kinetic analysis showed BlIDH was a NADP⁺-dependent IDH (NADP-IDH), with a 567- and 193-fold preference for NADP⁺ over NAD⁺ in the presence of Mg(2+) and Mn(2+), respectively. The maximal activity for BlIDH occurred at 60 °C (with Mn(2+)) and 65 °C (with Mg(2+)), and pH 7.5 (with Mn(2+)) and pH 8.0 (with Mg(2+)). Heat-inactivation profiles revealed that BlIDH retained 50% of maximal activity after incubation at 45 °C for 20 min with either Mn(2+) or Mg(2+). Furthermore, the coenzyme specificity of BlIDH can be completely reversed from NADP⁺ to NAD⁺ by a factor of 2387 by replacing six residues. This current work, the first report on the coenzyme specificity conversion of Type II NADP-IDHs, would provide better insight into the evolution of NADP⁺ use by the IDH family.

Mutations in adenine-binding pockets enhance catalytic properties of NAD(P)H-dependent enzymes

NAD(P)H-dependent enzymes are ubiquitous in metabolism and cellular processes and are also of great interest for pharmaceutical and industrial applications. Here, we present a structure-guided enzyme engineering strategy for improving catalytic properties of NAD(P)H-dependent enzymes toward native or native-like reactions using mutations to the enzyme's adenine-binding pocket, distal to the site of catalysis. Screening single-site saturation mutagenesis libraries identified mutations that increased catalytic efficiency up to 10-fold in 7 out of 10 enzymes. The enzymes improved in this study represent three different cofactor-binding folds (Rossmann, DHQS-like, and FAD/NAD binding) and utilize both NADH and NADPH. Structural and biochemical analyses show that the improved activities are accompanied by minimal changes in other properties (cooperativity, thermostability, pH optimum, uncoupling), and initial tests on two enzymes (ScADH6 and EcFucO) show improved functionality in Escherichia coli.

In Silico Constraint-Based Strain Optimization Methods: the Quest for Optimal Cell Factories

Shifting from chemical to biotechnological processes is one of the cornerstones of 21st century industry. The production of a great range of chemicals via biotechnological means is a key challenge on the way toward a bio-based economy. However, this shift is occurring at a pace slower than initially expected. The development of efficient cell factories that allow for competitive production yields is of paramount importance for this leap to happen. Constraint-based models of metabolism, together with in silico strain design algorithms, promise to reveal insights into the best genetic design strategies, a step further toward achieving that goal. In this work, a thorough analysis of the main in silico constraint-based strain design strategies and algorithms is presented, their application in real-world case studies is analyzed, and a path for the future is discussed.

Isocitrate dehydrogenase from Streptococcus mutans: biochemical properties and evaluation of a putative phosphorylation site at Ser102

Isocitrate deyhdrogenase (IDH) is a reversible enzyme in the tricarboxylic acid cycle that catalyzes the NAD(P)⁺-dependent oxidative decarboxylation of isocitrate to α-ketoglutarate (αKG) and the NAD(P)H/CO₂-dependent reductive carboxylation of αKG to isocitrate. The IDH gene from Streptococcus mutans was fused with the icd gene promoter from Escherichia coli to initiate its expression in the glutamate auxotrophic strain E. coli Δicd::kan^r of which the icd gene has been replaced by kanamycin resistance gene. The expression of S. mutans IDH (SmIDH) may restore the wild-type phenotype of the icd-defective strain on minimal medium without glutamate. The molecular weight of SmIDH was estimated to be 70 kDa by gel filtration chromatography, suggesting a homodimeric structure. SmIDH was divalent cation-dependent and Mn²⁺ was found to be the most effective cation. The optimal pH of SmIDH was 7.8 and the maximum activity was around 45°C. SmIDH was completely NAD⁺ dependent and its apparent K_m for NAD⁺ was 137 μM. In order to evaluate the role of the putative phosphorylation site at Ser102 in catalysis, two “stably phosphorylated” mutants were constructed by converting Ser102 into Glu102 or Asp102 in SmIDH to mimick a constitutively phosphorylated state. Meanwhile, the functional roles of another four amino acids (threonine, glycine, alanine and tyrosine) containing variant size of side chains were investigated. The replacement of Asp102 or Glu102 totally inactivated the enzyme, while the S102T, S102G, S102A and S102Y mutants decreased the affinity to isocitrate and only retained 16.0%, 2.8%, 3.3% and 1.1% of the original activity, respectively. These results reveal that Ser102 plays important role in substrate binding and is required for the enzyme function. Also, Ser102 in SmIDH is a potential phosphorylation site, indicating that the ancient NAD-dependent IDHs might be the underlying origin of “phosphorylation mechanism” used by their bacterial NADP-dependent homologs.

Functional relevance of dynamic properties of Dimeric NADP-dependent Isocitrate Dehydrogenases

Background

Isocitrate Dehydrogenases (IDHs) are important enzymes present in all living cells. Three subfamilies of functionally dimeric IDHs (subfamilies I, II, III) are known. Subfamily I are well-studied bacterial IDHs, like that of Escherischia coli. Subfamily II has predominantly eukaryotic members, but it also has several bacterial members, many being pathogens or endosymbionts. subfamily III IDHs are NAD-dependent.

The eukaryotic-like subfamily II IDH from pathogenic bacteria such as Mycobacterium tuberculosis IDH1 are expected to have regulation similar to that of bacteria which use the glyoxylate bypass to survive starvation. Yet they are structurally different from IDHs of subfamily I, such as the E. coli IDH.

Results

We have used phylogeny, structural comparisons and molecular dynamics simulations to highlight the similarity and differences between NADP-dependent dimeric IDHs with an emphasis on regulation. Our phylogenetic study indicates that an additional subfamily (IV) may also be present. Variation in sequence and structure in an aligned region may indicate functional importance concerning regulation in bacterial subfamily I IDHs.

Correlation in movement of prominent loops seen from molecular dynamics may explain the adaptability and diversity of the predominantly eukaryotic subfamily II IDHs.

Conclusion

This study discusses possible regulatory mechanisms operating in various IDHs and implications for regulation of eukaryotic-like bacterial IDHs such as that of M. tuberculosis, which may provide avenues for intervention in disease.

Induced fit and the catalytic mechanism of isocitrate dehydrogenase

NADP(+) dependent isocitrate dehydrogenase (IDH; EC 1.1.1.42) belongs to a large family of α-hydroxyacid oxidative β-decarboxylases that catalyze similar three-step reactions, with dehydrogenation to an oxaloacid intermediate preceding β-decarboxylation to an enol intermediate followed by tautomerization to the final α-ketone product. A comprehensive view of the induced fit needed for catalysis is revealed on comparing the first "fully closed" crystal structures of a pseudo-Michaelis complex of wild-type Escherichia coli IDH (EcoIDH) and the "fully closed" reaction product complex of the K100M mutant with previously obtained "quasi-closed" and "open" conformations. Conserved catalytic residues, binding the nicotinamide ring of NADP(+) and the metal-bound substrate, move as rigid bodies during domain closure by a hinge motion that spans the central β-sheet in each monomer. Interactions established between Thr105 and Ser113, which flank the "phosphorylation loop", and the nicotinamide mononucleotide moiety of NADP(+) establish productive coenzyme binding. Electrostatic interactions of a Lys100-Leu103-Asn115-Glu336 tetrad play a pivotal role in assembling a catalytically competent active site. As predicted, Lys230* is positioned to deprotonate/reprotonate the α-hydroxyl in both reaction steps and Tyr160 moves into position to protonate C3 following β-decarboxylation. A proton relay from the catalytic triad Tyr160-Asp307-Lys230* connects the α-hydroxyl of isocitrate to the bulk solvent to complete the picture of the catalytic mechanism.

The complex structures of isocitrate dehydrogenase from Clostridium thermocellum and Desulfotalea psychrophila suggest a new active site locking mechanism

Isocitrate dehydrogenase (IDH) catalyzes the oxidative NAD(P)⁺-dependent decarboxylation of isocitrate into α-ketoglutarate and CO₂ and is present in organisms spanning the biological range of temperature. We have solved two crystal structures of the thermophilic Clostridium thermocellum IDH (CtIDH), a native open apo CtIDH to 2.35 Å and a quaternary complex of CtIDH with NADP⁺, isocitrate and Mg²⁺ to 2.5 Å. To compare to these a quaternary complex structure of the psychrophilic Desulfotalea psychrophila IDH (DpIDH) was also resolved to 1.93 Å. CtIDH and DpIDH showed similar global thermal stabilities with melting temperatures of 67.9 and 66.9 °C, respectively. CtIDH represents a typical thermophilic enzyme, with a large number of ionic interactions and hydrogen bonds per residue combined with stabilization of the N and C termini. CtIDH had a higher activity temperature optimum, and showed greater affinity for the substrates with an active site that was less thermolabile compared to DpIDH. The uncompensated negative surface charge and the enlarged methionine cluster in the hinge region both of which are important for cold activity in DpIDH, were absent in CtIDH. These structural comparisons revealed that prokaryotic IDHs in subfamily II have a unique locking mechanism involving Arg310, Asp251′ and Arg255 (CtIDH). These interactions lock the large domain to the small domain and direct NADP⁺ into the correct orientation, which together are important for NADP⁺ selectivity.

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.

Enzymatic and whole cell catalysis: finding new strategies for old processes

The use of enzymes and whole bacterial cells has allowed the production of a plethora of compounds that have been used for centuries in foods and beverages. However, only recently we have been able to master techniques that allow the design and development of new biocatalysts with high stability and productivity. Rational redesign and directed evolution have lead to engineered enzymes with new characteristics whilst the understanding of adaptation mechanisms in bacterial cells has allowed their use under new operational conditions. Bacteria able to thrive under the most extreme conditions have also provided new and extraordinary catalytic processes. In this review, the new tools available for the improvement of biocatalysts are presented and discussed.

Different biochemical mechanisms ensure network-wide balancing of reducing equivalents in microbial metabolism

To sustain growth, the catabolic formation of the redox equivalent NADPH must be balanced with the anabolic demand. The mechanisms that ensure such network-wide balancing, however, are presently not understood. Based on 13C-detected intracellular fluxes, metabolite concentrations, and cofactor specificities for all relevant central metabolic enzymes, we have quantified catabolic NADPH production in Agrobacterium tumefaciens, Bacillus subtilis, Escherichia coli, Paracoccus versutus, Pseudomonas fluorescens, Rhodobacter sphaeroides, Sinorhizobium meliloti, and Zymomonas mobilis. For six species, the estimated NADPH production from glucose catabolism exceeded the requirements for biomass synthesis. Exceptions were P. fluorescens, with balanced rates, and E. coli, with insufficient catabolic production, in which about one-third of the NADPH is supplied via the membrane-bound transhydrogenase PntAB. P. versutus and B. subtilis were the only species that appear to rely on transhydrogenases for balancing NADPH overproduction during growth on glucose. In the other four species, the main but not exclusive redox-balancing mechanism appears to be the dual cofactor specificities of several catabolic enzymes and/or the existence of isoenzymes with distinct cofactor specificities, in particular glucose 6-phosphate dehydrogenase. An unexpected key finding for all species, except E. coli and B. subtilis, was the lack of cofactor specificity in the oxidative pentose phosphate pathway, which contrasts with the textbook view of the pentose phosphate pathway dehydrogenases as being NADP+ dependent.

Role of residues in the adenosine binding site of NAD of the Ascaris suum malic enzyme

Ascaris suum mitochondrial malic enzyme catalyzes the divalent metal ion dependent conversion of l-malate to pyruvate and CO(2), with concomitant reduction of NAD(P) to NAD(P)H. In this study, some of the residues that form the adenosine binding site of NAD were mutated to determine their role in binding of the cofactor and/or catalysis. D361, which is completely conserved among species, is located in the dinucleotide-binding Rossmann fold and makes a salt bridge with R370, which is also highly conserved. D361 was mutated to E, A and N. R370 was mutated to K and A. D361E and A mutant enzymes were inactive, likely a result of the increase in the volume in the case of the D361E mutant enzyme that caused clashes with the surrounding residues, and loss of the ionic interaction between D361 and R370, for D361A. Although the K(m) for the substrates and isotope effect values did not show significant changes for the D361N mutant enzyme, V/E(t) decreased by 1400-fold. Data suggested the nonproductive binding of the cofactor, giving a low fraction of active enzyme. The R370K mutant enzyme did not show any significant changes in the kinetic parameters, while the R370A mutant enzyme gave a slight change in V/E(t), contrary to expectations. Overall, results suggest that the salt bridge between D361 and R370 is important for maintaining the productive conformation of the NAD binding site. Mutation of residues involved leads to nonproductive binding of NAD. The interaction stabilizes one of the Rossmann fold loops that NAD binds. Mutation of H377 to lysine, which is conserved in NADP-specific malic enzymes and proposed to be a cofactor specificity determinant, did not cause a shift in cofactor specificity of the Ascaris malic enzyme from NAD to NADP. However, it is confirmed that this residue is an important second layer residue that affects the packing of the first layer residues that directly interact with the cofactor.

Structure and quantum chemical analysis of NAD+-dependent isocitrate dehydrogenase: hydride transfer and co-factor specificity

The crystal structure of Acidithiobacillus thiooxidans isocitrate dehydrogenase complexed with NAD+ and citrate has been solved to a resolution of 1.9 A. The protein fold of this NAD+-dependent enzyme shares a high similarity with those of NADP+-dependent bacterial ICDHs. The NAD+ and the citrate are clearly identified in the active site cleft with a well-defined electron density. Asp-357 is the direct cofactor-specificity determinant that interacts with 2'-OH and 3'-OH of the adenosine ribose. The adenosine ribose takes a C2'-endo puckering conformation as previously reported for an NAD+-specific isopropylmalate dehydrogenase. The nicotinamide moiety of NAD+ has the amide NH2 group oriented in cis conformation with respect to the C4 carbon of the nicotinamide ring, slanted toward the bound citrate molecule with a dihedral angle of -21 degrees . The semi-empirical molecular orbital calculation suggests that the pro-R hydrogen atom at C4 of NADH would bear the largest negative charge when the amide NH2 group is in such conformation, suggesting that the amide group has a catalytically significant role in stabilizing the transition state as NADH is being formed during the hydride transfer catalysis.

Empirical fitness landscapes reveal accessible evolutionary paths

When attempting to understand evolution, we traditionally rely on analysing evolutionary outcomes, despite the fact that unseen intermediates determine its course. A handful of recent studies has begun to explore these intermediate evolutionary forms, which can be reconstructed in the laboratory. With this first view on empirical evolutionary landscapes, we can now finally start asking why particular evolutionary paths are taken.

Amino acid residues that determine functional specificity of NADP- and NAD-dependent isocitrate and isopropylmalate dehydrogenases

Isocitrate and isopropylmalalte dehydrogenases are homologous enzymes important for the cell metabolism. They oxidize their substrates using NAD or NADP as cofactors. Thus, they have two specificities, towards the substrate and the cofactor, appearing in three combinations. Although many three-dimensional (3D) structures are resolved, identification of amino acids determining these specificities remains a challenge. We present computational identification and analysis of specificity-determining positions (SDPs). Besides many experimentally proven SDPs, we predict new SDPs, for example, four substrate-specific positions (103Leu, 105Thr, 337Ala, and 341Thr in IDH from E. coli) that contact the cofactor and may play a role in the recognition process.

Preparative use of isolated CYP102 monooxygenases—A critical appraisal

Isolated P450 monooxygenases have for long been neglected catalysts in enzyme technology. This is surprising as they display a remarkable substrate specificity catalyzing reactions, which represent a challenge for classic organic chemistry. On the other hand, many P450 monooxygenases are membrane bound, depend on rather complicated electron transfer systems and require expensive cofactors such as NAD(P)H. Their activities are low, and stability leaves much to be desired. The use of bacterial P450 monooxygenases from CYP102 family allows overcoming some of these handicaps. They are soluble and their turnovers are high, presumably because their N-terminal heme monooxygenase and their C-terminal diflavin reductase domain are covalently linked. In recent years, protein engineering approaches have been successfully used to turn CYP102 monooxgenases into powerful biocatalysts.

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.

Structural plasticity of an aminoacyl-tRNA synthetase active site

Recently, tRNA aminoacyl-tRNA synthetase pairs have been evolved that allow one to genetically encode a large array of unnatural amino acids in both prokaryotic and eukaryotic organisms. We have determined the crystal structures of two substrate-bound Methanococcus jannaschii tyrosyl aminoacyl-tRNA synthetases that charge the unnatural amino acids p-bromophenylalanine and 3-(2-naphthyl)alanine (NpAla). A comparison of these structures with the substrate-bound WT synthetase, as well as a mutant synthetase that charges p-acetylphenylalanine, shows that altered specificity is due to both side-chain and backbone rearrangements within the active site that modify hydrogen bonds and packing interactions with substrate, as well as disrupt the alpha8-helix, which spans the WT active site. The high degree of structural plasticity that is observed in these aminoacyl-tRNA synthetases is rarely found in other mutant enzymes with altered specificities and provides an explanation for the surprising adaptability of the genetic code to novel amino acids.

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.

Structural characterization of a p-acetylphenylalanyl aminoacyl-tRNA synthetase

It has been recently shown that orthogonal tRNA/aminoacyl-tRNA synthetase pairs can be evolved to allow genetic incorporation of unnatural amino acids into proteins in both prokaryotes and eukaryotes. Here we describe the crystal structure of an evolved aminoacyl-tRNA synthetase that charges the unnatural amino acid p-acetylphenylalanine. Molecular recognition is due to altered hydrogen bonding and packing interactions with bound substrate that result from changes in both side-chain and backbone conformation.

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.

The coenzyme specificity of Candida tenuis xylose reductase (AKR2B5) explored by site-directed mutagenesis and X-ray crystallography

CtXR (xylose reductase from the yeast Candida tenuis; AKR2B5) can utilize NADPH or NADH as co-substrate for the reduction of D-xylose into xylitol, NADPH being preferred approx. 33-fold. X-ray structures of CtXR bound to NADP+ and NAD+ have revealed two different protein conformations capable of accommodating the presence or absence of the coenzyme 2'-phosphate group. Here we have used site-directed mutagenesis to replace interactions specific to the enzyme-NADP+ complex with the aim of engineering the co-substrate-dependent conformational switch towards improved NADH selectivity. Purified single-site mutants K274R (Lys274-->Arg), K274M, K274G, S275A, N276D, R280H and the double mutant K274R-N276D were characterized by steady-state kinetic analysis of enzymic D-xylose reductions with NADH and NADPH at 25 degrees C (pH 7.0). The results reveal between 2- and 193-fold increases in NADH versus NADPH selectivity in the mutants, compared with the wild-type, with only modest alterations of the original NADH-linked xylose specificity and catalytic-centre activity. Catalytic reaction profile analysis demonstrated that all mutations produced parallel effects of similar magnitude on ground-state binding of coenzyme and transition state stabilization. The crystal structure of the double mutant showing the best improvement of coenzyme selectivity versus wild-type and exhibiting a 5-fold preference for NADH over NADPH was determined in a binary complex with NAD+ at 2.2 A resolution.

Strain engineering for stereoselective bioreduction of dicarbonyl compounds by yeast reductases

Pure chiral molecules are needed in the pharmaceutical and chemical industry as intermediates for the production of drugs or fine chemicals. Microorganisms represent an attractive alternative to chemical synthesis since they have the potential to generate single stereoisomers in high enantiomeric excess (ee). The baker's yeast Saccharomyces cerevisiae can notably reduce dicarbonyl compounds (in particular alpha- and beta-diketones and keto esters) to chiral alcohols with high ee. However, products are formed at a low rate. Moreover, large amounts of co-substrate are required for the regeneration of NADPH that is the preferred co-factor in almost all the known dicarbonyl reductions. Traditionally, better ee, reduction rate and product titre have been achieved via process engineering. The advent of recombinant DNA technology provides an alternative strategy to improve productivity and yield by strain engineering. This review discusses two aspects of strain engineering: (i) the generation of strains with higher reductase activity towards dicarbonyl compounds and (ii) the optimisation of co-substrate utilisation for NADPH cofactor regeneration.

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.

The selective cause of an ancient adaptation

Phylogenetic analysis reveals that the use of nicotinamide adenine dinucleotide phosphate (NADP) by prokaryotic isocitrate dehydrogenase (IDH) arose around the time eukaryotic mitochondria first appeared, about 3.5 billion years ago. We replaced the wild-type gene that encodes the NADP-dependent IDH of Escherichia coli with an engineered gene that possesses the ancestral NAD-dependent phenotype. The engineered enzyme is disfavored during competition for acetate. The selection intensifies in genetic backgrounds where other sources of reduced NADP have been removed. A survey of sequenced prokaryotic genomes reveals that those genomes that encode isocitrate lyase, which is essential for growth on acetate, always have an NADP-dependent IDH. Those with only an NAD-dependent IDH never have isocitrate lyase. Hence, the NADP dependence of prokaryotic IDH is an ancient adaptation to anabolic demand for reduced NADP during growth on acetate.

Complete reversal of coenzyme specificity of xylitol dehydrogenase and increase of thermostability by the introduction of structural zinc

Pichia stipitis NAD(+)-dependent xylitol dehydrogenase (XDH), a medium-chain dehydrogenase/reductase, is one of the key enzymes in ethanol fermentation from xylose. For the construction of an efficient biomass-ethanol conversion system, we focused on the two areas of XDH, 1) change of coenzyme specificity from NAD(+) to NADP(+) and 2) thermostabilization by introducing an additional zinc atom. Site-directed mutagenesis was used to examine the roles of Asp(207), Ile(208), Phe(209), and Asn(211) in the discrimination between NAD(+) and NADP(+). Single mutants (D207A, I208R, F209S, and N211R) improved 5 approximately 48-fold in catalytic efficiency (k(cat)/K(m)) with NADP(+) compared with the wild type but retained substantial activity with NAD(+). The double mutants (D207A/I208R and D207A/F209S) improved by 3 orders of magnitude in k(cat)/K(m) with NADP(+), but they still preferred NAD(+) to NADP(+). The triple mutant (D207A/I208R/F209S) and quadruple mutant (D207A/I208R/F209S/N211R) showed more than 4500-fold higher values in k(cat)/K(m) with NADP(+) than the wild-type enzyme, reaching values comparable with k(cat)/K(m) with NAD(+) of the wild-type enzyme. Because most NADP(+)-dependent XDH mutants constructed in this study decreased the thermostability compared with the wild-type enzyme, we attempted to improve the thermostability of XDH mutants by the introduction of an additional zinc atom. The introduction of three cysteine residues in wild-type XDH gave an additional zinc-binding site and improved the thermostability. The introduction of this mutation in D207A/I208R/F209S and D207A/I208R/F209S/N211R mutants increased the thermostability and further increased the catalytic activity with NADP(+).

Crystal structure of the monomeric isocitrate dehydrogenase in the presence of NADP+: insight into the cofactor recognition, catalysis, and evolution

NADP+-dependent monomeric isocitrate dehydrogenase (IDH) from the nitrogen-fixing bacterium Azotobacter vinelandii (AvIDH) is one of members of the beta-decarboxylating dehydrogenase family and catalyzes the dehydration and decarboxylation of isocitrate to yield 2-oxoglutrate and CO2 in the Krebs cycle. We solved the crystal structure of the AvIDH in complex with cofactor NADP+ (AvIDH-NADP+ complex). The final refined model shows the closed form that has never been detected in any previously solved structures of beta-decarboxylating dehydrogenases. The structure also reveals all of the residues that interact with NADP+. The structure-based sequence alignment reveals that these residues were not conserved in any other dimeric NADP+-dependent IDHs. Therefore the NADP+ specificity of the monomeric and dimeric IDHs was independently acquired through the evolutional process. The AvIDH was known to show an exceptionally high turnover rate. The structure of the AvIDH-NADP+ complex indicates that one loop, which is not present in the Escherichia coli IDHs, reliably stabilizes the conformation of the nicotinamide mononucleotide of the bound NADP+ by forming a few hydrogen bonds, and such interactions are considered to be important for the monomeric enzyme to initiate the hydride transfer reaction immediately. Finally, the structure of the AvIDH is compared with that of other dimeric NADP-IDHs. Several structural features demonstrate that the monomeric IDHs are structurally more related to the eukaryotic dimeric IDHs than to the bacterial dimeric IDHs.

Homologous binding sites in yeast isocitrate dehydrogenase for cofactor (NAD+) and allosteric activator (AMP)

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an allosterically regulated octameric enzyme composed of two types of homologous subunits designated IDH1 and IDH2. Based on sequence comparisons and structural models, both subunits are predicted to have adenine nucleotide binding sites. This was tested by alanine replacement of residues in putative sites in each subunit. Targets included adjacent aspartate/isoleucine residues implicated as important for determining cofactor specificity in related dehydrogenases and a residue in each IDH subunit in a position occupied by histidine in other cofactor binding sites. The primary kinetic effects of D286A/I287A and of H281A replacements in IDH2 were found to be a dramatic reduction in apparent affinity of the holoenzyme for NAD(+) and a concomitant reduction in V(max). Ligand binding assays also showed that the H281A mutant enzyme fails to bind NAD(+) under conditions that are saturating for the wild-type enzyme. In contrast, the primary effect of corresponding D279A/D280A and of R274A replacements in IDH1 is a reduction in holoenzyme binding of AMP, with concomitant alterations in kinetic and isocitrate binding properties normally associated with activation by this allosteric effector. These results suggest that the nucleotide cofactor binding site is primarily contributed by the IDH2 subunit, whereas the homologous nucleotide binding site in IDH1 has evolved for regulatory binding of AMP. These results are consistent with previous studies demonstrating that the catalytic isocitrate binding sites are comprised of residues primarily contributed by IDH2, whereas sites for regulatory binding of isocitrate are contributed by analogous residues of IDH1. In this study, we also demonstrate that a prerequisite for holoenzyme binding of NAD(+) is binding of isocitrate/Mg(2+) at the IDH2 catalytic site. This is comparable to the dependence of AMP binding upon binding of isocitrate at the IDH1 regulatory site.

Structure of the monomeric isocitrate dehydrogenase: evidence of a protein monomerization by a domain duplication

NADP(+)-dependent isocitrate dehydrogenase is a member of the beta-decarboxylating dehydrogenase family and catalyzes the oxidative decarboxylation reaction from 2R,3S-isocitrate to yield 2-oxoglutarate and CO(2) in the Krebs cycle. Although most prokaryotic NADP(+)-dependent isocitrate dehydrogenases (IDHs) are homodimeric enzymes, the monomeric IDH with a molecular weight of 80-100 kDa has been found in a few species of bacteria. The 1.95 A crystal structure of the monomeric IDH revealed that it consists of two distinct domains, and its folding topology is related to the dimeric IDH. The structure of the large domain repeats a motif observed in the dimeric IDH. Such a fusional structure by domain duplication enables a single polypeptide chain to form a structure at the catalytic site that is homologous to the dimeric IDH, the catalytic site of which is located at the interface of two identical subunits.

Conserved positions for ribose recognition: importance of water bridging interactions among ATP, ADP and FAD-protein complexes

Analysis of the spatial arrangement of protein and water atoms that form polar interactions with ribose has been performed for a structurally non-redundant dataset of ATP, ADP and FAD-protein complexes. The 26 ligand-protein structures were separated into two groups corresponding to the most populated furanose ring conformations (N and S-domains). Four conserved positions were found for S-domain protein-ligand complexes and five for N-domain complexes. Multiple protein folds and secondary structural elements were represented at a single conserved position. The following novel points were revealed: (i) Two complementary positions sometimes combine to describe a putative atomic spatial location for a specific conserved binding spot. (ii) More than one third of the interactions scored were water-mediated. Thus, conserved spatial positions rich in water atoms are a significant feature of ribose-protein complexes.

Biochemical and molecular characterization of the NAD(+)-dependent isocitrate dehydrogenase from the chemolithotroph Acidithiobacillus thiooxidans

An isocitrate dehydrogenase (ICDH) with an unique coenzyme specificity from Acidithiobacillus thiooxidans was purified and characterized, and its gene was cloned. The native enzyme was homodimeric with a subunit of M(r) 45000 and showed a 78-fold preference for NAD(+) over NADP(+). The cloned ICDH gene (icd) was expressed in an icd-deficient strain of Escherichia coli EB106; the activity was found in the cell extract. The gene encodes a 429-amino acid polypeptide and is located between open reading frames encoding a putative aconitase gene (upstream of icd) and a putative succinyl-CoA synthase beta-subunit gene (downstream of icd). A. thiooxidans ICDH showed high sequence similarity to bacterial NADP(+)-dependent ICDH rather than eukaryotic NAD(+)-dependent ICDH, but the NAD(+)-preference of the enzyme was suggested due to residues conserved in the coenzyme binding site of the NAD(+)-dependent decarboxylating dehydrogenase.

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.

Analysis of homodimeric protein interfaces by graph-spectral methods

The quaternary structures impart structural and functional credibility to proteins. In a multi-subunit protein, it is important to understand the factors that drive the association or dissociation of the subunits. It is a well known fact that both hydrophobic and charged interactions contribute to the stability of the protein interface. The interface residues are also known to be highly conserved. Though they are buried in the oligomer, these residues are either exposed or partially exposed in the monomer. It is felt that a systematic and objective method of identifying interface clusters and their analysis can significantly contribute to the identification of a residue or a collection of residues important for oligomerization. Recently, we have applied the techniques of graph-spectral methods to a variety of problems related to protein structure and folding. A major advantage of this methodology is that the problem is viewed from a global protein topology point of view rather than localized regions of the protein structure. In the present investigation, we have applied the methods of graph-spectral analysis to identify side chain clusters at the interface and the centers of these clusters in a set of homodimeric proteins. These clusters are analyzed in terms of properties such as amino acid composition, accessibility to solvent and conservation of residues. Interesting results such as participation of charged and aromatic residues like arginine, glutamic acid, histidine, phenylalanine and tyrosine, consistent with earlier investigations, have emerged from these analyses. Important additional information is that the residues involved are a part of a cluster(s) and that they are sequentially distant residues which have come closer to each other in the three-dimensional structure of the protein. These residues can easily be detected using our graph-spectral algorithm. This method has also been used to identify important residues ('hot spots') in dimerization and also to detect dimerization sites on the monomer. The residues predicted using the present algorithm have correlated well with the experiments indicating the efficacy of this method in predicting residues involved in dimer stability.

Comparison of isocitrate dehydrogenase from three hyperthermophiles reveals differences in thermostability, cofactor specificity, oligomeric state, and phylogenetic affiliation

With the aim of gaining insight into the molecular and phylogenetic relationships of isocitrate dehydrogenase (IDH) from hyperthermophiles, we carried out a comparative study of putative IDHs identified in the genomes of the eubacterium Thermotoga maritima and the archaea Aeropyrum pernix and Pyrococcus furiosus. An optimum for activity at 90 degrees C or above was found for each IDH. PfIDH and ApIDH were the most thermostable with a melting temperature of 103.7 and 109.9 degrees C, respectively, compared with 98.3 and 98.5 degrees C for TmIDH and AfIDH, respectively. Analytical ultracentrifugation revealed a tetrameric oligomeric state for TmIDH and a homodimeric state for ApIDH and PfIDH. TmIDH and ApIDH were NADP-dependent (K(m)((NADP)) of 55.2 and 44.4 microm, respectively) whereas PfIDH was NAD-dependent (K(m)((NAD)) of 68.3 microm). These data document that TmIDH represents a novel tetrameric NADP-dependent form of IDH and that PfIDH is a homodimeric NAD-dependent IDH not previously found among the archaea. The homodimeric NADP-IDH present in A. pernix is the most common form of IDH known so far. The evolutionary relationships of ApIDH, PfIDH, and TmIDH with all of the available amino acid sequences of di- and multimeric IDHs are described and discussed.

Engineering novel biocatalytic routes for production of semisynthetic opiate drugs

The morphine alkaloids and their semisynthetic derivatives provide a diverse range of important pharmaceutical drugs. Current production of semisynthetic opiate drugs is by chemical means from naturally occurring morphine, codeine and thebaine. Although various microbial transformations of morphine alkaloids have been identified since the 1960s, more recently there has been considerable effort devoted to engineering biocatalytic routes for producing these important compounds. Such biocatalytic routes are attractive, as they would provide an alternative to the chemical production processes which suffer from limited supply of precursors, often low yields and toxic wastes. The biotransformation of morphine and codeine to the potent analgesic hydromorphone and the mild analgesic/antitussive hydrocodone, respectively, by recombinant Escherichia coli has been demonstrated and the problems encountered when engineering such a system will be discussed.

Probing the determinants of coenzyme specificity in ferredoxin-NADP+ reductase by site-directed mutagenesis

On the basis of sequence and three-dimensional structure comparison between Anabaena PCC7119 ferredoxin-NADP(+) reductase (FNR) and other reductases from its structurally related family that bind either NADP(+)/H or NAD(+)/H, a set of amino acid residues that might determine the FNR coenzyme specificity can be assigned. These residues include Thr-155, Ser-223, Arg-224, Arg-233 and Tyr-235. Systematic replacement of these amino acids was done to identify which of them are the main determinants of coenzyme specificity. Our data indicate that all of the residues interacting with the 2'-phosphate of NADP(+)/H in Anabaena FNR are not involved to the same extent in determining coenzyme specificity and affinity. Thus, it is found that Ser-223 and Tyr-235 are important for determining NADP(+)/H specificity and orientation with respect to the protein, whereas Arg-224 and Arg-233 provide only secondary interactions in Anabaena FNR. The analysis of the T155G FNR form also indicates that the determinants of coenzyme specificity are not only situated in the 2'-phosphate NADP(+)/H interacting region but that other regions of the protein must be involved. These regions, although not interacting directly with the coenzyme, must produce specific structural arrangements of the backbone chain that determine coenzyme specificity. The loop formed by residues 261-268 in Anabaena FNR must be one of these regions.