Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase

Characterization and designing of an SAM riboswitch to establish a high-throughput screening platform for SAM overproduction in Saccharomyces cerevisiae

S-adenosyl- l-methionine (SAM) is a high-value compound widely used in the treatment of various diseases. SAM can be produced through fermentation, but further enhancing the microbial production of SAM requires novel high-throughput screening methods for rapid detection and screening of mutant libraries. In this work, an SAM-OFF riboswitch capable of responding to the SAM concentration was obtained and a high-throughput platform for screening SAM overproducers was established. SAM synthase was engineered by semirational design and directed evolution, which resulted in the SAM2S203F,W164R,T251S,Y285F,S365R mutant with almost twice higher catalytic activity than the parental enzyme. The best mutant was then introduced into Saccharomyces cerevisiae BY4741, and the resulting strain BSM8 produced a sevenfold higher SAM titer in shake-flask fermentation, reaching 1.25 g L-1 . This work provides a reference for designing biosensors to dynamically detect metabolite concentrations for high-throughput screening and the construction of effective microbial cell factories.

Protein cost minimization promotes the emergence of coenzyme redundancy

Metabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes [e.g., NAD(H) and NADP(H)] associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling and structural and sequence analysis indicates that coenzyme redundancy may not be essential for metabolism but could rather constitute an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate and acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules [e.g., NAD(H) and NADP(H)] that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here, we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy and that coenzymes typically thought to exist in a single pool (e.g., coenzyme A [CoA]) may exist in more than one form (e.g., dephospho-CoA).

Biochemical Characterization of Recombinant Isocitrate Dehydrogenase and Its Putative Role in the Physiology of an Acidophilic Micrarchaeon

Despite several discoveries in recent years, the physiology of acidophilic Micrarchaeota, such as "Candidatus Micrarchaeum harzensis A_DKE", remains largely enigmatic, as they highly express numerous genes encoding hypothetical proteins. Due to a lacking genetic system, it is difficult to elucidate the biological function of the corresponding proteins and heterologous expression is required. In order to prove the viability of this approach, A_DKE's isocitrate dehydrogenase (MhIDH) was recombinantly produced in Escherichia coli and purified to electrophoretic homogeneity for biochemical characterization. MhIDH showed optimal activity around pH 8 and appeared to be specific for NADP+ yet promiscuous regarding divalent cations as cofactors. Kinetic studies showed KM-values of 53.03 ± 5.63 µM and 1.94 ± 0.12 mM and kcat-values of 38.48 ± 1.62 and 43.99 ± 1.46 s-1 resulting in kcat/KM-values of 725 ± 107.62 and 22.69 ± 2.15 mM-1 s-1 for DL-isocitrate and NADP+, respectively. MhIDH's exceptionally low affinity for NADP+, potentially limiting its reaction rate, can likely be attributed to the presence of a proline residue in the NADP+ binding pocket, which might cause a decrease in hydrogen bonding of the cofactor and a distortion of local secondary structure.

Revolutionizing enzyme engineering through artificial intelligence and machine learning

The combinatorial space of an enzyme sequence has astronomical possibilities and exploring it with contemporary experimental techniques is arduous and often ineffective. Multi-target objectives such as concomitantly achieving improved selectivity, solubility and activity of an enzyme have narrow plausibility under approaches of restricted mutagenesis and combinatorial search. Traditional enzyme engineering approaches have a limited scope for complex optimization due to the requirement of a priori knowledge or experimental burden of screening huge protein libraries. The recent surge in high-throughput experimental methods including Next Generation Sequencing and automated screening has flooded the field of molecular biology with big-data, which requires us to re-think our concurrent approaches towards enzyme engineering. Artificial Intelligence (AI) and Machine Learning (ML) have great potential to revolutionize smart enzyme engineering without the explicit need for a complete understanding of the underlying molecular system. Here, we portray the role and position of AI techniques in the field of enzyme engineering along with their scope and limitations. In addition, we explain how the traditional approaches of directed evolution and rational design can be extended through AI tools. Recent successful examples of AI-assisted enzyme engineering projects and their deviation from traditional approaches are highlighted. A comprehensive picture of current challenges and future avenues for AI in enzyme engineering are also discussed.

Innovations in CAZyme gene diversity and its modification for biorefinery applications

For sustainable growth, concept of biorefineries as recourse to the "fossil derived" energy source is important. Here, the Carbohydrate Active enZymes (CAZymes) play decisive role in generation of biofuels and related sugar-based products utilizing lignocellulose as a carbon source. Given their industrial significance, extensive studies on the evolution of CAZymes have been carried out. Various bacterial and fungal organisms have been scrutinized for the development of CAZymes, where advance techniques for strain enhancement such as CRISPR and analysis of specific expression systems have been deployed. Specific Omic-based techniques along with protein engineering have been adopted to unearth novel CAZymes and improve applicability of existing enzymes. In-Silico computational research and functional annotation of new CAZymes to synergy experiments are being carried out to devise cocktails of enzymes for use in biorefineries. Thus, with the establishment of these technologies, increased diversity of CAZymes with broad span of functions and applications is seen.

The fitness challenge of studying molecular adaptation

Advances in bioinformatics and high-throughput genetic analysis increasingly allow us to predict the genetic basis of adaptive traits. These predictions can be tested and confirmed, but the molecular-level changes - i.e. the molecular adaptation - that link genetic differences to organism fitness remain generally unknown. In recent years, a series of studies have started to unpick the mechanisms of adaptation at the molecular level. In particular, this work has examined how changes in protein function, activity, and regulation cause improved organismal fitness. Key to addressing molecular adaptations is identifying systems and designing experiments that integrate changes in the genome, protein chemistry (molecular phenotype), and fitness. Knowledge of the molecular changes underpinning adaptations allow new insight into the constraints on, and repeatability of adaptations, and of the basis of non-additive interactions between adaptive mutations. Here we critically discuss a series of studies that examine the molecular-level adaptations that connect genetic changes and fitness.

Improving Catalytic Efficiency and Changing Substrate Spectrum for Asymmetric Biocatalytic Reductive Amination

With the advantages of biocatalytic method, enzymes have been excavated for the synthesis of chiral amino acids by the reductive amination of ketones, offering a promising way of producing pharmaceutical intermediates. In this work, a robust phenylalanine dehydrogenase (PheDH) with wide substrate spectrum and high catalytic efficiency was constructed through rational design and active-site-targeted, site-specific mutagenesis by using the parent enzyme from Bacillus halodurans. Active sites with bonding substrate and amino acid residues surrounding the substrate binding pocket, 49L-50G-51G, 74M,77K, 122G-123T-124D-125M, 275N, 305L and 308V of the PheDH, were identified. Noticeably, the new mutant PheDH (E113D-N276L) showed approximately 6.06-fold increment of k_cat/Km in the oxidative deamination and more than 1.58-fold in the reductive amination compared to that of the wide type. Meanwhile, the PheDHs exhibit high capacity of accepting benzylic and aliphatic ketone substrates. The broad specificity, high catalytic efficiency and selectivity, along with excellent thermal stability, render these broad-spectrum enzymes ideal targets for further development with potential diagnostic reagent and pharmaceutical compounds applications.

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 genetically encoded tool for manipulation of NADP+/NADPH in living cells

NADH and NADPH are redox coenzymes broadly required for energy metabolism, biosynthesis and detoxification. Despite detailed knowledge of specific enzymes and pathways that utilize these coenzymes, a holistic understanding of the regulation and compartmentalization of NADH and NADPH-dependent pathways is lacking, in part because of a lack of tools with which to investigate them in living cells. We previously reported the use of the naturally occurring Lactobacillus brevis H₂O-forming NADH oxidase (LbNOX) as a genetic tool for manipulation of the NAD⁺/NADH ratio in human cells. Here we present TPNOX (triphosphopyridine nucleotide oxidase), a rationally designed and engineered mutant of LbNOX that is strictly specific towards NADPH. We characterize the effects of TPNOX expression on cellular metabolism and use it in combination with LbNOX to show how the redox states of mitochondrial NADPH and NADH pools are connected.

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.

Engineering Nucleotide Specificity of Succinyl-CoA Synthetase in Blastocystis: The Emerging Role of Gatekeeper Residues

Charged, solvent-exposed residues at the entrance to the substrate binding site (gatekeeper residues) produce electrostatic dipole interactions with approaching substrates, and control their access by a novel mechanism called “electrostatic gatekeeper effect”. This proof-of-concept study demonstrates that the nucleotide specificity can be engineered by altering the electrostatic properties of the gatekeeper residues outside the binding site. Using Blastocystis succinyl-CoA synthetase (SCS, EC 6.2.1.5), we demonstrated that the gatekeeper mutant (ED) resulted in ATP-specific SCS to show high GTP specificity. Moreover, nucleotide binding site mutant (LF) had no effect on GTP specificity and remained ATP-specific. However, via combination of the gatekeeper mutant with the nucleotide binding site mutant (ED+LF), a complete reversal of nucleotide specificity was obtained with GTP, but no detectable activity was obtained with ATP. This striking result of the combined mutant (ED+LF) was due to two changes; negatively charged gatekeeper residues (ED) favored GTP access, and nucleotide binding site residues (LF) altered ATP binding, which was consistent with the hypothesis of the “electrostatic gatekeeper effect”. These results were further supported by molecular modeling and simulation studies. Hence, it is imperative to extend the strategy of the gatekeeper effect in a different range of crucial enzymes (synthetases, kinases, and transferases) to engineer substrate specificity for various industrial applications and substrate-based drug design.

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.

A computational strategy for altering an enzyme in its cofactor preference to NAD(H) and/or NADP(H)

Coenzyme engineering, especially for altered coenzyme specificity, has been a research hotspot for more than a decade. In the present study, a novel computational strategy that enhances the hydrogen-bond interaction between an enzyme and a coenzyme was developed and utilized to alter the coenzyme preference. This novel computational strategy only required the structure of the target enzyme. No other homologous enzymes were needed to achieve alteration in the coenzyme preference of a certain enzyme. Using our novel strategy, Gox2181 was reconstructed from exhibiting complete NADPH preference to exhibiting dual cofactor specificity for NADH and NADPH. Structure-guided Gox2181 mutants were designed in silico and molecular dynamics simulations were performed to evaluate the strength of hydrogen-bond interactions between the enzyme and the coenzyme NADPH. Three Gox2181 mutants displaying high structure stability and structural compatibility to NADH/NADPH were chosen for experimental confirmation. Among the three Gox2181 mutants, Gox2181-Q20R&D43S showed the highest enzymatic activity by utilizing NADPH as its coenzyme, which was even better than the wild-type enzyme. In addition, isothermal titration calorimetry analysis further verified that Gox2181-Q20R&D43S was able to interact with NADPH but the wild-type enzyme could not. This novel computational strategy represents an insightful approach for altering the cofactor preference of target enzymes.

Engineering of global regulator cAMP receptor protein (CRP) in Escherichia coli for improved lycopene production

Transcriptional engineering has received significant attention for improving strains by modulating the behavior of transcription factors, which could be used to reprogram a series of gene transcriptions and enable multiple simultaneous modifications at the genomic level. In this study, engineering of the cAMP receptor protein (CRP) was explored with the aim of subtly balancing entire pathway networks and potentially improving lycopene production without significant genetic intervention in other pathways. Amino acid mutations were introduced to CRP by error-prone PCR, and three variants (mcrp26, mcrp159 and mcrp424) with increased lycopene productivity were screened. Combinations of three point mutations were then created via site-directed mutagenesis. The best mutant gene (mcrp26) was integrated into the genome of E. coli BW25113-BIE to replace the wild-type crp gene (MT-1), which resulted in a higher lycopene production (18.49mg/g DCW) compared to the original strain (WT). The mutant strain MT-1 was further investigated in a 10-L bench-top fermentor with a lycopene yield of 128mg/l at 20h, approximately 25% higher than WT. DNA microarray analyses showed that 396 genes (229 up-regulated and 167 down-regulated) were differentially expressed in the mutant MT-1 compared to WT. Finally, the introduction of the mutant crp gene (mcrp26) increased β-carotene production in E. coli. This is the first report of improving the phenotype for metabolite overproduction in E. coli using a CRP engineering strategy.

A novel serine hydroxymethyltransferase from Arthrobacter nicotianae: characterization and improving catalytic efficiency by rational design

BACKGROUND

Serine hydroxymethyltransferase (SHMT) is the key enzyme in L-serine enzymatic production, suggesting the importance of obtaining a SHMT with high activity.

RESULTS

Here, a novel SHMT gene, glyA, was obtained through degenerate oligonucleotide-primed PCR and encoded a novel SHMT with 54.3% similarity to the known SHMT from Escherichia coli. The obtained protein AnSHMT showed the optimal activity at 40 °C and pH 7.5, and was more stable in weakly alkali conditions (pH 6.5-8.5) than Hyphomicrobium methylovorum's SHMT (pH 6.0-7.5), In order to improve the catalytic efficiency of the wild type, the site-directed mutagenesis based on sequences alignment and bioinformatics prediction, was used and the catalytic efficiency of the mutant I249L was found to be 2.78-fold higher than that of the wild-type, with the replacement of isoleucine by leucine at the 249 position.

CONCLUSIONS

This research provides useful information about the interesting site, and the application of DOP-PCR in cloning a novel glyA gene.

Structure-based conversion of the coenzyme requirement of a short-chain dehydrogenase/reductase involved in bacterial alginate metabolism

The alginate-assimilating bacterium, Sphingomonas sp. strain A1, degrades the polysaccharides to monosaccharides through four alginate lyase reactions. The resultant monosaccharide, which is nonenzymatically converted to 4-deoxy-L-erythro-5-hexoseulose uronate (DEH), is further metabolized to 2-keto-3-deoxy-D-gluconate by NADPH-dependent reductase A1-R in the short-chain dehydrogenase/reductase (SDR) family. A1-R-deficient cells produced another DEH reductase, designated A1-R', with a preference for NADH. Here, we show the identification of a novel NADH-dependent DEH reductase A1-R' in strain A1, structural determination of A1-R' by x-ray crystallography, and structure-based conversion of a coenzyme requirement in SDR enzymes, A1-R and A1-R'. A1-R' was purified from strain A1 cells and enzymatically characterized. Except for the coenzyme requirement, there was no significant difference in enzyme characteristics between A1-R and A1-R'. Crystal structures of A1-R' and A1-R'·NAD(+) complex were determined at 1.8 and 2.7 Å resolutions, respectively. Because of a 64% sequence identity, overall structures of A1-R' and A1-R were similar, although a difference in the coenzyme-binding site (particularly the nucleoside ribose 2' region) was observed. Distinct from A1-R, A1-R' included a negatively charged, shallower binding site. These differences were caused by amino acid residues on the two loops around the site. The A1-R' mutant with the two A1-R-typed loops maintained potent enzyme activity with specificity for NADPH rather than NADH, demonstrating that the two loops determine the coenzyme requirement, and loop exchange is a promising method for conversion of coenzyme requirement in the SDR family.

Molecular determinants of the cofactor specificity of ribitol dehydrogenase, a short-chain dehydrogenase/reductase

Ribitol dehydrogenase from Zymomonas mobilis (ZmRDH) catalyzes the conversion of ribitol to d-ribulose and concomitantly reduces NAD(P)(+) to NAD(P)H. A systematic approach involving an initial sequence alignment-based residue screening, followed by a homology model-based screening and site-directed mutagenesis of the screened residues, was used to study the molecular determinants of the cofactor specificity of ZmRDH. A homologous conserved amino acid, Ser156, in the substrate-binding pocket of the wild-type ZmRDH was identified as an important residue affecting the cofactor specificity of ZmRDH. Further insights into the function of the Ser156 residue were obtained by substituting it with other hydrophobic nonpolar or polar amino acids. Substituting Ser156 with the negatively charged amino acids (Asp and Glu) altered the cofactor specificity of ZmRDH toward NAD(+) (S156D, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 10.9, where K(m)(,NAD) is the K(m) for NAD(+) and K(m)(,NADP) is the K(m) for NADP(+)). In contrast, the mutants containing positively charged amino acids (His, Lys, or Arg) at position 156 showed a higher efficiency with NADP(+) as the cofactor (S156H, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 0.11). These data, in addition to those of molecular dynamics and isothermal titration calorimetry studies, suggest that the cofactor specificity of ZmRDH can be modulated by manipulating the amino acid residue at position 156.

Re-engineering the discrimination between the oxidized coenzymes NAD+ and NADP+ in clostridial glutamate dehydrogenase and a thorough reappraisal of the coenzyme specificity of the wild-type enzyme

Clostridial glutamate dehydrogenase mutants, designed to accommodate the 2'-phosphate of disfavoured NADPH, showed the expected large specificity shifts with NAD(P)H. Puzzlingly, similar assays with oxidized cofactors initially revealed little improvement with NADP(+) , although rates with NAD(+) were markedly diminished. This article reveals that the enzyme's discrimination in favour of NAD(+) and against NADP(+) had been greatly underestimated and has indeed been abated by a factor of > 16,000 by the mutagenesis. Initially, stopped-flow studies of the wild-type enzyme showed a burst increase of A(340) with NADP(+) but not NAD(+), with amplitude depending on the concentration of the coenzyme, rather than enzyme. Amplitude also varied with the commercial source of the NADP(+). FPLC, HPLC and mass spectrometry identified NAD(+) contamination ranging from 0.04 to 0.37% in different commercial samples. It is now clear that apparent rates of NADP(+) utilization mainly reflected the reduction of contaminating NAD(+), creating an entirely false view of the initial coenzyme specificity and also of the effects of mutagenesis. Purification of the NADP(+) eliminated the burst. With freshly purified NADP(+), the NAD(+) : NADP(+) activity ratio under standard conditions, previously estimated as 300 : 1, is 11,000. The catalytic efficiency ratio is even higher at 80,000. Retested with pure cofactor, mutants showed marked specificity shifts in the expected direction, for example, 16 200 fold change in catalytic efficiency ratio for the mutant F238S/P262S, confirming that the key structural determinants of specificity have been successfully identified. Of wider significance, these results underline that, without purification, even the best commercial coenzyme preparations are inadequate for such studies.

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.

Computational design of Candida boidinii xylose reductase for altered cofactor specificity

In this study we introduce a computationally-driven enzyme redesign workflow for altering cofactor specificity from NADPH to NADH. By compiling and comparing data from previous studies involving cofactor switching mutations, we show that their effect cannot be explained as straightforward changes in volume, hydrophobicity, charge, or BLOSUM62 scores of the residues populating the cofactor binding site. Instead, we find that the use of a detailed cofactor binding energy approximation is needed to adequately capture the relative affinity towards different cofactors. The implicit solvation models Generalized Born with molecular volume integration and Generalized Born with simple switching were integrated in the iterative protein redesign and optimization (IPRO) framework to drive the redesign of Candida boidinii xylose reductase (CbXR) to function using the non-native cofactor NADH. We identified 10 variants, out of the 8,000 possible combinations of mutations, that improve the computationally assessed binding affinity for NADH by introducing mutations in the CbXR binding pocket. Experimental testing revealed that seven out of ten possessed significant xylose reductase activity utilizing NADH, with the best experimental design (CbXR-GGD) being 27-fold more active on NADH. The NADPH-dependent activity for eight out of ten predicted designs was either completely abolished or significantly diminished by at least 90%, yielding a greater than 10(4)-fold change in specificity to NADH (CbXR-REG). The remaining two variants (CbXR-RTT and CBXR-EQR) had dual cofactor specificity for both nicotinamide cofactors.

Structure-guided alteration of coenzyme specificity of formate dehydrogenase by saturation mutagenesis to enable efficient utilization of NADP+

Formate dehydrogenase from Candida boidinii (CboFDH) catalyses the oxidation of formate anion to carbon dioxide with concomitant reduction of NAD(+) to NADH. CboFDH is highly specific to NAD(+) and virtually fails to catalyze the reaction with NADP(+). Based on structural information for CboFDH, the loop region between beta-sheet 7 and alpha-helix 10 in the dinucleotide-binding fold was predicted as a principal determinant of coenzyme specificity. Sequence alignment with other formate dehydrogenases revealed two residues (Asp195 and Tyr196) that could account for the observed coenzyme specificity. Positions 195 and 196 were subjected to two rounds of site-saturation mutagenesis and screening and enabled the identification of a double mutant Asp195Gln/Tyr196His, which showed a more than 2 x 10(7)-fold improvement in overall catalytic efficiency with NADP(+) and a more than 900-fold decrease in the efficiency with NAD(+) as cofactors. The results demonstrate that the combined polar interactions and steric factors comprise the main structural determinants responsible for coenzyme specificity. The double mutant Asp195Gln/Tyr196His was tested for practical applicability in a cofactor recycling system composed of cytochrome P450 monooxygenase from Bacillus subtilis, (CYP102A2), NADP(+), formic acid and omega-(p-nitrophenyl)dodecanoic acid (12-pNCA). Using a 1250-fold excess of 12-pNCA over NADP(+) the first order rate constant was determined to be equal to k(obs) = 0.059 +/- 0.004 min(-1).

Localization and nucleotide specificity of Blastocystis succinyl-CoA synthetase

The anaerobic lifestyle of the intestinal parasite Blastocystis raises questions about the biochemistry and function of its mitochondria-like organelles. We have characterized the Blastocystis succinyl-CoA synthetase (SCS), a tricarboxylic acid cycle enzyme that conserves energy by substrate-level phosphorylation. We show that SCS localizes to the enigmatic Blastocystis organelles, indicating that these organelles might play a similar role in energy metabolism as classic mitochondria. Although analysis of residues inside the nucleotide-binding site suggests that Blastocystis SCS is GTP-specific, we demonstrate that it is ATP-specific. Homology modelling, followed by flexible docking and molecular dynamics simulations, indicates that while both ATP and GTP fit into the Blastocystis SCS active site, GTP is destabilized by electrostatic dipole interactions with Lys 42 and Lys 110, the side-chains of which lie outside the nucleotide-binding cavity. It has been proposed that residues in direct contact with the substrate determine nucleotide specificity in SCS. However, our results indicate that, in Blastocystis, an electrostatic gatekeeper controls which ligands can enter the binding site.

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.

Direct demonstration of an adaptive constraint

The role of constraint in adaptive evolution is an open question. Directed evolution of an engineered beta-isopropylmalate dehydrogenase (IMDH), with coenzyme specificity switched from nicotinamide adenine dinucleotide (NAD) to nicotinamide adenine dinucleotide phosphate (NADP), always produces mutants with lower affinities for NADP. This result is the correlated response to selection for relief from inhibition by NADPH (the reduced form of NADP) expected of an adaptive landscape subject to three enzymatic constraints: an upper limit to the rate of maximum turnover (kcat), a correlation in NADP and NADPH affinities, and a trade-off between NAD and NADP usage. Two additional constraints, high intracellular NADPH abundance and the cost of compensatory protein synthesis, have ensured the conserved use of NAD by IMDH throughout evolution. Our results show that selective mechanisms and evolutionary constraints are to be understood in terms of underlying adaptive landscapes.

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.

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.

Site-saturation mutagenesis is more efficient than DNA shuffling for the directed evolution of beta-fucosidase from beta-galactosidase

Protein engineers use a variety of mutagenic strategies to adapt enzymes to novel substrates. Directed evolution techniques (random mutagenesis and high-throughput screening) offer a systematic approach to the management of protein complexity. This sub-discipline was galvanized by the invention of DNA shuffling, a procedure that randomly recombines point mutations in vitro. In one influential study, Escherichia coli beta-galactosidase (BGAL) variants with enhanced beta-fucosidase activity (tenfold increase in k(cat)/K(M) in reactions with the novel para-nitrophenyl-beta-d-fucopyranoside substrate; 39-fold decrease in reactivity with the "native"para-nitrophenyl-beta-d-galactopyranoside substrate) were evolved in seven rounds of DNA shuffling and screening. Here, we show that a single round of site-saturation mutagenesis and screening enabled the identification of beta-fucosidases that are significantly more active (180-fold increase in k(cat)/K(M) in reactions with the novel substrate) and specific (700,000-fold inversion of specificity) than the best variants in the previous study. Site-saturation mutagenesis thus proved faster, less resource-intensive and more effective than DNA shuffling for this particular evolutionary pathway.

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.

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(+).

Catalytic mechanism and application of formate dehydrogenase

NAD+-dependent formate dehydrogenase (FDH) is an abundant enzyme that plays an important role in energy supply of methylotrophic microorganisms and in response to stress in plants. FDH belongs to the superfamily of D-specific 2-hydroxy acid dehydrogenases. FDH is widely accepted as a model enzyme to study the mechanism of hydride ion transfer in the active center of dehydrogenases because the reaction catalyzed by the enzyme is devoid of proton transfer steps and implies a substrate with relatively simple structure. FDH is also widely used in enzymatic syntheses of optically active compounds as a versatile biocatalyst for NAD(P)H regeneration consumed in the main reaction. This review covers the late developments in cloning genes of FDH from various sources, studies of its catalytic mechanism and physiological role, and its application for new chiral syntheses.

Identifying determinants of NADPH specificity in Baeyer-Villiger monooxygenases

The Baeyer-Villiger monooxygenase (BVMO), 4-hydroxyacetophenone monooxygenase (HAPMO), uses NADPH and O(2) to oxidize a variety of aromatic ketones and sulfides. The FAD-containing enzyme has a 700-fold preference for NADPH over NADH. Sequence alignment with other BVMOs, which are all known to be selective for NADPH, revealed three conserved basic residues, which could account for the observed coenzyme specificity. The corresponding residues in HAPMO (Arg339, Lys439 and Arg440) were mutated and the properties of the purified mutant enzymes were studied. For Arg440 no involvement in coenzyme recognition could be shown as mutant R440A was totally inactive. Although this mutant could still be fully reduced by NADPH, no oxygenation occurred, indicating that this residue is crucial for completing the catalytic cycle of HAPMO. Characterization of several Arg339 and Lys439 mutants revealed that these residues are indeed both involved in coenzyme recognition. Mutant R339A showed a largely decreased affinity for NADPH, as judged from kinetic analysis and binding experiments. Replacing Arg339 also resulted in a decreased catalytic efficiency with NADH. Mutant K439A displayed a 100-fold decrease in catalytic efficiency with NADPH, mainly caused by an increased K(m). However, the efficiency with NADH increased fourfold. Saturation mutagenesis at position 439 showed that the presence of an asparagine or a phenylalanine improves the catalytic efficiency with NADH by a factor of 6 to 7. All Lys439 mutants displayed a lower affinity for AADP(+), confirming a role of the lysine in recognizing the 2'-phosphate of NADPH. The results obtained could be extrapolated to the sequence-related cyclohexanone monooxygenase. Replacing Lys326 in this BVMO, which is analogous to Lys439 in HAPMO, again changed the coenzyme specificity towards NADH. These results indicate that the strict NADPH dependency of this class of monooxygenases is based upon recognition of the coenzyme by several basic residues.

A comparison of directed evolution approaches using the beta-glucuronidase model system

Protein engineers can alter the properties of enzymes by directing their evolution in vitro. Many methods to generate molecular diversity and to identify improved clones have been developed, but experimental evolution remains as much an art as a science. We previously used DNA shuffling (sexual recombination) and a histochemical screen to direct the evolution of Escherichia coli beta-glucuronidase (GUS) variants with improved beta-galactosidase (BGAL) activity. Here, we employ the same model evolutionary system to test the efficiencies of several other techniques: recursive random mutagenesis (asexual), combinatorial cassette mutagenesis (high-frequency recombination) and a versatile high-throughput microplate screen. GUS variants with altered specificity evolved in each trial, but different combinations of mutagenesis and screening techniques effected the fixation of different beneficial mutations. The new microplate screen identified a broader set of mutations than the previously employed X-gal colony screen. Recursive random mutagenesis produced essentially asexual populations, within which beneficial mutations drove each other into extinction (clonal interference); DNA shuffling and combinatorial cassette mutagenesis led instead to the accumulation of beneficial mutations within a single allele. These results explain why recombinational approaches generally increase the efficiency of laboratory evolution.

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.

Quantitative chimeric analysis of six specificity determinants that differentiate Escherichia coli aspartate from tyrosine aminotransferase

The six mutations, referred to as the Hex mutations, that together have been shown to convert Escherichia coli aspartate aminotransferase (AATase) specificity to be substantially like that of E. coli tyrosine aminotransferase (TATase) are dissected into two groups, (T109S/N297S) and (V39L/K41Y/T47I/N69L). The letters on the left and right of the numbers designate AATase and TATase residues, respectively. The T109S/N297S pair has been investigated previously. The latter group, the "Grease" set, is now placed in the AATase framework, and the retroGrease set (L39V/Y41K/I47T/L69N) is substituted into TATase. The Grease mutations in the AATase framework were found primarily to lower K(M)s for both aromatic and dicarboxylic substrates. In contrast, retroGrease TATase exhibits lowered k(cat)s for both substrates. The six retroHex mutations, combining retroGrease and S109T/S297N, were found to invert the substrate specificity of TATase, creating an enzyme with a nearly ninefold preference (k(cat)/K(M)) for aspartate over phenylalanine. The retroHex mutations perturb the electrostatic environment of the pyridoxal phosphate cofactor, as evidenced by a spectrophotometric titration of the internal aldimine, which uniquely shows two pK(a)s, 6.1 and 9.1. RetroHex was also found to have impaired dimer stability, with a K(D) for dimer dissociation of 350 nM compared with the wild type K(D) of 4 nM. Context dependence and additivity analyses demonstrate the importance of interactions of the Grease residues with the surrounding protein framework in both the AATase and TATase contexts, and with residues 109 and 297 in particular. Context dependence and cooperativity are particularly evident in the effects of mutations on k(cat)/K(M)(Asp). Effects on k(cat)/K(M)(Phe) are more nearly additive and context independent.

Engineering of coenzyme specificity of formate dehydrogenase from Saccharomyces cerevisiae

A eukaryotic formate dehydrogenase (EC 1.2.1.2, FDH) with its substrate specificity changed from NAD(+) to NADP(+) has been constructed by introducing two single-point mutations, Asp(196)-->Ala (D196A) and Tyr(197)-->Arg (Y197R). The mutagenesis was based on the results of homology modelling of a NAD(+)-specific FDH from Saccharomyces cerevisiae (SceFDH) using the Pseudomonas sp.101 FDH (PseFDH) crystal structure as a template. The resulting model structure suggested that Asp(196) and Tyr(197) mediate the absolute coenzyme specificity of SceFDH for NAD(+).

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.

Changing the target base specificity of the EcoRV DNA methyltransferase by rational de novo protein-design

The EcoRV DNA-(adenine-N(6))-methyltransferase (M.EcoRV) specifically modifies the first adenine residue within GATATC sequences. During catalysis, the enzyme flips its target base out of the DNA helix and binds it into a target base binding pocket which is formed in part by Lys16 and Tyr196. A cytosine residue is accepted by wild-type M.EcoRV as a substrate at a 31-fold reduced efficiency with respect to the k(cat)/K(M) values if it is located in a CT mismatch substrate (GCTATC/GATATC). Cytosine residues positioned in a CG base pair (GCTATC/GATAGC) are modified at much more reduced rates, because flipping out the target base is much more difficult in this case. We intended to change the target base specificity of M.EcoRV from adenine-N(6) to cytosine-N(4). To this end we generated, purified and characterized 15 variants of the enzyme, containing single, double and triple amino acid exchanges following different design approaches. One concept was to reduce the size of the target base binding pocket by site-directed mutagenesis. The K16R variant showed an altered specificity, with a 22-fold preference for cytosine as the target base in a mismatch substrate. This corresponds to a 680-fold change in specificity, which was accompanied by only a small loss in catalytic activity with the cytosine substrate. The K16R/Y196W variant no longer methylated adenine residues at all and its activity towards cytosine was reduced only 17-fold. Therefore, we have changed the target base specificity of M.EcoRV from adenine to cytosine by rational protein design. Because there are no natural paragons for the variants described here, a change of the target base specificity of a DNA interacting enzyme was possible by rational de novo design of its active site.

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.

Crystal structures of mutants of Thermus thermophilus IPMDH adapted to low temperatures

Random mutagenesis on thermophilic 3-isopropylmalate dehydrogenases (IPMDH; EC 1.1.1.85) produced mutant enzymes which adapt to low temperatures. These mutants had higher activity at lower temperatures than the wild-type enzyme without losing high thermostability. Here we report three structures of the mutants of Thermus thermophilus IPMDH determined by X-ray diffraction which was adapted to a low-temperature environment. Two of them have unstable coenzyme binding states and the other one probably has a stable substrate binding state. The present research suggests that the adaptation is correlated with the binding of either coenzyme or the substrate.

In vitro evolution of beta-glucuronidase into a beta-galactosidase proceeds through non-specific intermediates

The Escherichia coli beta-glucuronidase (GUS) was evolved in vitro to catalyze the hydrolysis of a beta-galactoside substrate 500 times more efficiently (k(cat)/K(m)) than the wild-type, with a 52 million-fold inversion in specificity. The amino acid substitutions that recurred among 32 clones isolated in three rounds of DNA shuffling and screening were mapped to the active site. The functional consequences of these mutations were investigated by introducing them individually or in combination into otherwise wild-type gusA genes. The kinetic behavior of the purified mutant proteins in reactions with a series of substrate analogues show that four mutations account for the changes in substrate specificity, and that they are synergistic. An evolutionary intermediate, unlike the wild-type and evolved forms, exhibits broadened specificity for substrates dissimilar to either glucuronides or galactosides. These results are consistent with the "patchwork" hypothesis, which postulates that modern enzymes diverged from ancestors with broad specificity.

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.

A highly specific monomeric isocitrate dehydrogenase from Corynebacterium glutamicum

The monomeric isocitrate dehydrogenase (IDH) of Corynebacterium glutamicum is compared to the topologically distinct dimeric IDH of Escherichia coli. Both IDHs have evolved to efficiently catalyze identical reactions with similar pH optimum as well as striking specificity toward NADP and isocitrate. However, the monomeric IDH is 10-fold more active (calculated as kcat/Km.isocitrate/Km.NADP) and 7-fold more NADP-specific than the dimeric enzyme, favoring NADP over NAD by a factor of 50,000. Such an extraordinary coenzyme specificity is not rivaled by any other characterized dehydrogenases. In addition, the monomeric enzyme is 10-fold more specific for isocitrate. The spectacular substrate specificity may be predominantly attributed to the isocitrate-assisted stabilization of catalytic complex during hydride transfer. No significant overall sequence identity is found between the monomeric and dimeric enzymes. However, structure-based alignment leads to the identification of three regions in the monomeric enzyme that match closely the three motifs located in the central region of dimeric IDHs and the homologous isopropylmalate dehydrogenases. The role of Lys253 as catalytic residue has been demonstrated by site-directed mutagenesis. Our results suggest that monomeric and dimeric forms of IDHs are functionally and structurally homologous.