Redesigning the substrate specificity of an enzyme: isocitrate dehydrogenase

Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently

The amino acid sequence of a protein affects both its structure and its function. Thus, the ability to modify the sequence, and hence the structure and activity, of individual proteins in a systematic way, opens up many opportunities, both scientifically and (as we focus on here) for exploitation in biocatalysis. Modern methods of synthetic biology, whereby increasingly large sequences of DNA can be synthesised de novo, allow an unprecedented ability to engineer proteins with novel functions. However, the number of possible proteins is far too large to test individually, so we need means for navigating the 'search space' of possible protein sequences efficiently and reliably in order to find desirable activities and other properties. Enzymologists distinguish binding (Kd) and catalytic (kcat) steps. In a similar way, judicious strategies have blended design (for binding, specificity and active site modelling) with the more empirical methods of classical directed evolution (DE) for improving kcat (where natural evolution rarely seeks the highest values), especially with regard to residues distant from the active site and where the functional linkages underpinning enzyme dynamics are both unknown and hard to predict. Epistasis (where the 'best' amino acid at one site depends on that or those at others) is a notable feature of directed evolution. The aim of this review is to highlight some of the approaches that are being developed to allow us to use directed evolution to improve enzyme properties, often dramatically. We note that directed evolution differs in a number of ways from natural evolution, including in particular the available mechanisms and the likely selection pressures. Thus, we stress the opportunities afforded by techniques that enable one to map sequence to (structure and) activity in silico, as an effective means of modelling and exploring protein landscapes. Because known landscapes may be assessed and reasoned about as a whole, simultaneously, this offers opportunities for protein improvement not readily available to natural evolution on rapid timescales. Intelligent landscape navigation, informed by sequence-activity relationships and coupled to the emerging methods of synthetic biology, offers scope for the development of novel biocatalysts that are both highly active and robust.

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.

Enhancement of the latent 3-isopropylmalate dehydrogenase activity of promiscuous homoisocitrate dehydrogenase by directed evolution

HICDH (homoisocitrate dehydrogenase), which is involved in lysine biosynthesis through α-aminoadipate, is a paralogue of IPMDH [3-IPM (3-isopropylmalate) dehydrogenase], which is involved in leucine biosynthesis. TtHICDH (Thermus thermophilus HICDH) can recognize isocitrate, as well as homoisocitrate, as the substrate, and also shows IPMDH activity, although at a considerably decreased rate. In the present study, the promiscuous TtHICDH was evolved into an enzyme showing distinct IPMDH activity by directed evolution using a DNA-shuffling technique. Through five repeats of DNA shuffling/screening, variants that allowed Escherichia coli C600 (leuB−) to grow on a minimal medium in 2 days were obtained. One of the variants LR5–1, with eight amino acid replacements, was found to possess a 65-fold increased kcat/Km value for 3-IPM, compared with TtHICDH. Introduction of a single back-replacement H15Y change caused a further increase in the kcat/Km value and a partial recovery of the decreased thermotolerance of LR5–1. Site-directed mutagenesis revealed that most of the amino acid replacements found in LR5–1 effectively increased IPMDH activity; replacements around the substrate-binding site contributed to the improved recognition for 3-IPM, and other replacements at sites away from the substrate-binding site enhanced the turnover number for the IPMDH reaction. The crystal structure of LR5–1 was determined at 2.4 Å resolution and revealed that helix α4 was displaced in a manner suitable for recognition of the hydrophobic γ-moiety of 3-IPM. On the basis of the crystal structure, possible reasons for enhancement of the turnover number are discussed.

Isocitrate dehydrogenase-1 mutations: a fundamentally new understanding of diffuse glioma?

The discovery of somatic mutations in the gene encoding isocitrate dehydrogenase-1 (IDH1) in glioblastomas was remarkable because the enzyme was not previously identified with any known oncogenic pathway. IDH1 is mutated in up to 75% of grade II and grade III diffuse gliomas. Apart from acute myeloid leukaemia, other tumour types do not carry IDH1 mutations. Mutations in a homologous gene, IDH2, have also been identified, although they are much rarer. Although TP53 mutations and 1p/19q codeletions are mutually exclusive in gliomas, in both of these genotypes IDH1 mutations are common. IDH1 and IDH2 mutations are early events in the development of gliomas. Moreover, IDH1 and IDH2 mutations are a major prognostic marker for overall and progression-free survival in grade II-IV gliomas. Mutated IDH1 has an altered catalytic activity that results in the accumulation of 2-hydroxyglutarate. Molecularly, IDH1 and IDH2 mutations are heterozygous, affect only a single codon, and rarely occur together. Because IDH1 does not belong to a traditional oncogenic pathway and is specifically and commonly mutated in gliomas, the altered enzymatic activity of IDH1 may provide a fundamentally new understanding of diffuse glioma.

Chromosome stability and gene loss in cockroach endosymbionts

Insect endosymbiont genomes reflect massive gene loss. Two Blattabacterium genomes display colinearity and similar gene contents, despite high orthologous gene divergence, reflecting over 140 million years of independent evolution in separate cockroach lineages. We speculate that distant homologs may replace the functions of some eliminated genes through broadened substrate specificity.

Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate

Engineering enzymes to degrade anthropogenic compounds efficiently is challenging. We obtained Rhodococcus rhodochrous haloalkane dehalogenase mutants with up to 32-fold higher activity than wild type toward the toxic, recalcitrant anthropogenic compound 1,2,3-trichloropropane (TCP) using a new strategy. We identified key residues in access tunnels connecting the buried active site with bulk solvent by rational design and randomized them by directed evolution. The most active mutant has large aromatic residues at two out of three randomized positions and two positions modified by site-directed mutagenesis. These changes apparently enhance activity with TCP by decreasing accessibility of the active site for water molecules, thereby promoting activated complex formation. Kinetic analyses confirmed that the mutations improved carbon-halogen bond cleavage and shifted the rate-limiting step to the release of products. Engineering access tunnels by combining computer-assisted protein design with directed evolution may be a valuable strategy for refining catalytic properties of enzymes with buried active sites.

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.

Crystal structure of tetrameric homoisocitrate dehydrogenase from an extreme thermophile, Thermus thermophilus: involvement of hydrophobic dimer-dimer interaction in extremely high thermotolerance

The crystal structure of homoisocitrate dehydrogenase involved in lysine biosynthesis from Thermus thermophilus (TtHICDH) was determined at 1.85-A resolution. Arg85, which was shown to be a determinant for substrate specificity in our previous study, is positioned close to the putative substrate binding site and interacts with Glu122. Glu122 is highly conserved in the equivalent position in the primary sequence of ICDH and archaeal 3-isopropylmalate dehydrogenase (IPMDH) but interacts with main- and side-chain atoms in the same domain in those paralogs. In addition, a conserved Tyr residue (Tyr125 in TtHICDH) which extends its side chain toward a substrate and thus has a catalytic function in the related beta-decarboxylating dehydrogenases, is flipped out of the substrate-binding site. These results suggest the possibility that the conformation of the region containing Glu122-Tyr125 is changed upon substrate binding in TtHICDH. The crystal structure of TtHICDH also reveals that the arm region is involved in tetramer formation via hydrophobic interactions and might be responsible for the high thermotolerance. Mutation of Val135, located in the dimer-dimer interface and involved in the hydrophobic interaction, to Met alters the enzyme to a dimer (probably due to steric perturbation) and markedly decreases the thermal inactivation temperature. Both the crystal structure and the mutation analysis indicate that tetramer formation is involved in the extremely high thermotolerance of TtHICDH.

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.

Characterization of homoisocitrate dehydrogenase involved in lysine biosynthesis of an extremely thermophilic bacterium, Thermus thermophilus HB27, and evolutionary implication of beta-decarboxylating dehydrogenase

Although the presence of an enzyme that catalyzes beta-decarboxylating dehydrogenation of homoisocitrate to synthesize 2-oxoadipate has been postulated in the lysine biosynthesis pathway through alpha-aminoadipate (AAA), the enzyme has not yet been analyzed at all, because no gene encoding the enzyme has been identified until recently. A gene encoding a protein with a significant amino acid sequence identity to both isocitrate dehydrogenase and 3-isopropylmalate dehydrogenase was cloned from Thermus thermophilus HB27. The gene product produced in recombinant Escherichia coli cells demonstrated homoisocitrate dehydrogenase (HICDH) activity. A knockout mutant of the gene showed an AAA-auxotrophic phenotype, indicating that the gene product is involved in lysine biosynthesis through AAA. We therefore named this gene hicdh. HICDH, the gene product, did not catalyze the conversion of 3-isopropylmalate to 2-oxoisocaproate, a leucine biosynthetic reaction, but it did recognize isocitrate, a related compound in the tricarboxylic acid cycle, as well as homoisocitrate as a substrate. It is of interest that HICDH catalyzes the reaction with isocitrate about 20 times more efficiently than the reaction with the putative native substrate, homoisocitrate. The broad specificity and possible dual function suggest that this enzyme represents a key link in the evolution of the pathways utilizing citrate derivatives. Site-directed mutagenesis study reveals that replacement of Arg(85) with Val in HICDH causes complete loss of activity with isocitrate but significant activity with 3-isopropylmalate and retains activity with homoisocitrate. These results indicate that Arg(85) is a key residue for both substrate specificity and evolution of beta-decarboxylating dehydrogenases.

The application and usefulness of the ratio k(cat)/K(M)

The ratio of two important constituents in enzyme action, k(cat) and K(M), has become of value and provides insight into enzymatic mechanisms and the functional effects of enzyme mutations. It is timely to examine how this ratio is used and where it can be effectively applied. It has been called on some occasions the "specificity constant" and on other occasions, the "performance constant," and it is of interest to examine which then is the most accurate and useful way to utilize this ratio.

Directed evolution of enzymes and pathways for industrial biocatalysis

Directed evolution has become a powerful tool for developing enzyme and whole cell based biocatalysts. Significant recent advances include the creation of novel enzyme functions and the development of several new efficient directed evolution methods. The combination of directed evolution and rational design promises to accelerate the development of biocatalysts for applications in the pharmaceutical, chemical and food industries.

Use of indirect site-directed mutagenesis to alter the substrate specificity of methylamine dehydrogenase

Methylamine dehydrogenase (MADH) is a tryptophan tryptophylquinone-dependent enzyme that catalyzes the oxidative deamination of primary amines. Native MADH exhibits a strong preference for methylamine over longer carbon chain amines. Residue alphaPhe(55) controls this substrate specificity. When alphaPhe(55) is replaced with Ala, this preference is reversed with alphaF55A MADH preferring long-chain amines with at least seven carbons (Zhu, Z., Sun, D., and Davidson, V. L. (2000) Biochemistry 39, 11184-11186). To further modulate the substrate specificity of MADH, the side-chain of alphaPhe(55) was repositioned by site-directed mutagenesis of residue betaIle(107). This residue makes close contact with alphaPhe(55) and restricts its movement. betaI107V MADH exhibits a strong preference for propylamine, and betaI107N MADH exhibits a preference for 1-aminopentane. Thus, it has been possible to create forms of MADH that exhibit a preference for amines with carbon chain lengths of one, three, five, or seven carbons. The ability to discriminate between amines of different chain length was essentially abolished by an alphaF55I mutation. Molecular modeling studies with the known crystal structure are described that provide an explanation for these results. These results provide an example of a design-based approach to protein engineering in which site-directed mutagenesis on one residue can be used to reposition another residue to specifically alter enzyme specificity.

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.