Reversal of coenzyme specificity and improvement of catalytic efficiency of Pichia stipitis xylose reductase by rational site-directed mutagenesis

Altered Cofactor Preference of Thermostable StDAPDH by a Single Mutation at K159

D-amino acid production from 2-keto acid by reductive amination is an attractive pathway because of its high yield and environmental safety. StDAPDH, a meso-diaminopimelate dehydrogenase (meso-DAPDH) from Symbiobacterium thermophilum, was the first meso-DAPDH to show amination of 2-keto acids. Furthermore, StDAPDH shows excellent thermostability compared to other meso-DAPDHs. However, the cofactor of StDAPDH is NADP(H), which is less common than NAD(H) in industrial applications. Therefore, cofactor engineering for StDAPDH is needed. In this study, the highly conserved cofactor binding sites around the adenosine moiety of NADPH were targeted to determine cofactor specificity. Lysine residues within a loop were found to be critical for the cofactor specificity of StDAPDH. Replacement of lysine with arginine resulted in the activity of pyruvic acid with NADH as the cofactor. The affinity of K159R to pyruvic acid was equal with NADH or NADPH as the cofactor, regardless of the mutation. Molecular dynamics simulations revealed that the large steric hindrance of arginine and the interaction of the salt bridge between NADH and arginine may have restricted the free movement of NADH, which prompted the formation of a stable active conformation of mutant K159R. These results provide further understanding of the catalytic mechanism of StDAPDH and guidance for the cofactor engineering of StDAPDH.

Industrially produced pullulanases with thermostability: Discovery, engineering, and heterologous expression

Pullulanases (EC 3.2.1.41) are well-known starch-debranching enzymes widely used to hydrolyze α-1,6-glucosidic linkages in starch, pullulan, amylopectin, and other oligosaccharides, with application potentials in food, brewing, and pharmaceutical industries. Although extensive studies are done to discover and express pullulanases, only few are available with desirable characteristics for industrial applications. This raises the challenge to mine new enzyme sources, engineer proteins based on sequence/structure, and regulate expressions. We review here the identification of extremophilic and mesophilic microbes as sources of industrial pullulanases with desirable characteristics, including acid-resistance, thermostability, and psychrotrophism. We present current advances in site-directed mutagenesis and sequence/structure-guided protein engineering of pullulanases. In addition, we discuss heterologous expression of pullulanases in prokaryotic and eukaryotic microbial systems, and address the effectiveness of the expression elements and their regulation of enzyme production. Finally, we indicate future research needs to develop desired industrial pullulanases.

Kinetics and Predicted Structure of a Novel Xylose Reductase from Chaetomium thermophilum

While in search of an enzyme for the conversion of xylose to xylitol at elevated temperatures, a xylose reductase (XR) gene was identified in the genome of the thermophilic fungus Chaetomium thermophilum. The gene was heterologously expressed in Escherichia coli as a His6-tagged fusion protein and characterized for function and structure. The enzyme exhibits dual cofactor specificity for NADPH and NADH and prefers D-xylose over other pentoses and investigated hexoses. A homology model based on a XR from Candida tenuis was generated and the architecture of the cofactor binding site was investigated in detail. Despite the outstanding thermophilicity of its host the enzyme is, however, not thermostable.

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 General Tool for Engineering the NAD/NADP Cofactor Preference of Oxidoreductases

The ability to control enzymatic nicotinamide cofactor utilization is critical for engineering efficient metabolic pathways. However, the complex interactions that determine cofactor-binding preference render this engineering particularly challenging. Physics-based models have been insufficiently accurate and blind directed evolution methods too inefficient to be widely adopted. Building on a comprehensive survey of previous studies and our own prior engineering successes, we present a structure-guided, semirational strategy for reversing enzymatic nicotinamide cofactor specificity. This heuristic-based approach leverages the diversity and sensitivity of catalytically productive cofactor binding geometries to limit the problem to an experimentally tractable scale. We demonstrate the efficacy of this strategy by inverting the cofactor specificity of four structurally diverse NADP-dependent enzymes: glyoxylate reductase, cinnamyl alcohol dehydrogenase, xylose reductase, and iron-containing alcohol dehydrogenase. The analytical components of this approach have been fully automated and are available in the form of an easy-to-use web tool: Cofactor Specificity Reversal-Structural Analysis and Library Design (CSR-SALAD).

Bioethanol production by recycled Scheffersomyces stipitis in sequential batch fermentations with high cell density using xylose and glucose mixture

Here, it is shown three-step investigative procedures aiming to improve pentose-rich fermentations performance, involving a simple system for elevated mass production by Scheffersomyces stipitis (I), cellular recycle batch fermentations (CRBFs) at high cell density using two temperature strategies (fixed at 30°C; decreasing from 30 to 26°C) (II), and a short-term adaptation action seeking to acclimatize the microorganism in xylose rich-media (III). Cellular propagation provided 0.52gdrycellweightgRS(-1), resulting in an expressive value of 45.9gdrycellweightL(-1). The yeast robustness in CRBF was proven by effective ethanol production, reaching high xylose consumption (81%) and EtOH productivity (1.53gL(-1)h(-1)). Regarding the short-term adaptation, S. stipitis strengthened its robustness, as shown by a 6-fold increase in xylose reductase (XR) activity. The short fermentation time (20h for each batch) and the fermentation kinetics for ethanol production from xylose are quite promising.

In vitro Engineering of Novel Bioactivity in the Natural Enzymes

Enzymes catalyze various biochemical functions with high efficiency and specificity. In vitro design of the enzyme leads to novel bioactivity in this natural biomolecule that give answers of some vital questions like crucial residues in binding with substrate, molecular evolution, cofactor specificity etc. Enzyme engineering technology involves directed evolution, rational designing, semi-rational designing, and structure-based designing using chemical modifications. Similarly, combined computational and in vitro evolution approaches together help in artificial designing of novel bioactivity in the natural enzyme. DNA shuffling, error prone PCR and staggered extension process are used to artificially redesign active site of enzyme, which can alter its efficiency and specificity. Modifications of the enzyme can lead to the discovery of new path of molecular evolution, designing of efficient enzymes, locating active sites and crucial residues, shift in substrate, and cofactor specificity. The methods and thermodynamics of in vitro designing of the enzyme are also discussed. Similarly, engineered thermophilic and psychrophilic enzymes attain substrate specificity and activity of mesophilic enzymes that may also be beneficial for industry and therapeutics.

Improving ethanol and xylitol fermentation at elevated temperature through substitution of xylose reductase in Kluyveromyces marxianus

Thermo-tolerant yeast Kluyveromyces marxianus is able to utilize a wide range of substrates, including xylose; however, the xylose fermentation ability is weak because of the redox imbalance under oxygen-limited conditions. Alleviating the intracellular redox imbalance through engineering the coenzyme specificity of NADPH-preferring xylose reductase (XR) and improving the expression of XR should promote xylose consumption and fermentation. In this study, the native xylose reductase gene (Kmxyl1) of the K. marxianus strain was substituted with XR or its mutant genes from Pichia stipitis (Scheffersomyces stipitis). The ability of the resultant recombinant strains to assimilate xylose to produce xylitol and ethanol at elevated temperature was greatly improved. The strain YZB014 expressing mutant PsXR N272D, which has a higher activity with both NADPH and NADH as the coenzyme, achieved the best results, and produced 3.55 g l(-1) ethanol and 11.32 g l(-1) xylitol-an increase of 12.24- and 2.70-fold in product at 42 °C, respectively. A 3.94-fold increase of xylose consumption was observed compared with the K. marxianus YHJ010 harboring KmXyl1. However, the strain YZB015 expressing a mutant PsXR K21A/N272D, with which co-enzyme preference was completely reversed from NADPH to NADH, failed to ferment due to the low expression. So in order to improve xylose consumption and fermentation in K. marxianus, both higher activity and co-enzyme specificity change are necessary.

Computational design of glutamate dehydrogenase in Bacillus subtilis natto

Bacillus subtilis natto is widely used in industry to produce natto, a traditional and popular Japanese soybean food. However, during its secondary fermentation, high amounts of ammonia are released to give a negative influence on the flavor of natto. Glutamate dehydrogenase (GDH) is a key enzyme for the ammonia produced and released, because it catalyzes the oxidative deamination of glutamate to alpha-ketoglutarate using NAD(+) or NADP(+) as co-factor during carbon and nitrogen metabolism processes. To solve this problem, we employed multiple computational methods model and re-design GDH from Bacillus subtilis natto. Firstly, a structure model of GDH with cofactor NADP(+) was constructed by threading and ab initio modeling. Then the substrate glutamate were flexibly docked into the structure model to form the substrate-binding mode. According to the structural analysis of the substrate-binding mode, Lys80, Lys116, Arg196, Thr200, and Ser351 in the active site were found could form a significant hydrogen bonding network with the substrate, which was thought to play a crucial role in the substrate recognition and position. Thus, these residues were then mutated into other amino acids, and the substrate binding affinities for each mutant were calculated. Finally, three single mutants (K80A, K116Q, and S351A) were found to have significant decrease in the substrate binding affinities, which was further supported by our biochemical experiments.

Scaffold-based pan-agonist design for the PPARα, PPARβ and PPARγ receptors

As important members of nuclear receptor superfamily, Peroxisome proliferator-activated receptors (PPAR) play essential roles in regulating cellular differentiation, development, metabolism, and tumorigenesis of higher organisms. The PPAR receptors have 3 identified subtypes: PPARα, PPARβ and PPARγ, all of which have been treated as attractive targets for developing drugs to treat type 2 diabetes. Due to the undesirable side-effects, many PPAR agonists including PPARα/γ and PPARβ/γ dual agonists are stopped by US FDA in the clinical trials. An alternative strategy is to design novel pan-agonist that can simultaneously activate PPARα, PPARβ and PPARγ. Under such an idea, in the current study we adopted the core hopping algorithm and glide docking procedure to generate 7 novel compounds based on a typical PPAR pan-agonist LY465608. It was observed by the docking procedures and molecular dynamics simulations that the compounds generated by the core hopping and glide docking not only possessed the similar functions as the original LY465608 compound to activate PPARα, PPARβ and PPARγ receptors, but also had more favorable conformation for binding to the PPAR receptors. The additional absorption, distribution, metabolism and excretion (ADME) predictions showed that the 7 compounds (especially Cpd#1) hold high potential to be novel lead compounds for the PPAR pan-agonist. Our findings can provide a new strategy or useful insights for designing the effective pan-agonists against the type 2 diabetes.

Computational studies on the substrate interactions of influenza A virus PB2 subunit

Influenza virus, which spreads around the world in seasonal epidemics and leads to large numbers of deaths every year, has several ribonucleoproteins in the central core of the viral particle. These viral ribonucleoproteins can specifically bind the conserved 3' and 5' caps of the viral RNAs with responsibility for replication and transcription of the viral RNA in the nucleus of infected cells. A fundamental question of most importance is that how the cap-binding proteins in the influenza virus discriminates between capped RNAs and non-capped ones. To get an answer, we performed molecular dynamics simulations and free energy calculations on the influenza A virus PB2 subunit, an important component of the RNP complexes, with a cap analog m7GTP. Our calculations showed that some key residues in the active site, such as Arg355, His357, Glu361 as well as Gln406, could offer significant hydrogen bonding and hydrophobic interactions with the guanine ring of the cap analog m7GTP to form an aromatic sandwich mechanism for the cap recognition and positioning in the active site. Subsequently, we applied this idea to a virtual screening procedure and identified 5 potential candidates that might be inhibitors against the PB2 subunit. Interestingly, 2 candidates Cpd1 and Cpd2 have been already reported to have inhibitory activities to the influenza virus cap-binding proteins. Further calculation also showed that they had comparatively higher binding affinities to the PB2 subunit than that of m7GTP. We believed that our findings could give an atomic insight into the deeper understanding of the cap recognition and binding mechanism, providing useful information for searching or designing novel drugs against influenza viruses.

Autoinhibitory mechanism for the mutation-induced impaired FGF9 signaling

Fibroblast growth factor 9 (FGF9), an important member of the fibroblast growth factor (FGF) family, can bind with high affinity to FGFR3 in a heparin-dependent approach. In humans, the deletions and mutations resulting in dysfunction of the FGF9 signaling can cause human skeletal dysplasia and cancers. A mutation (S99N) in this protein has been identified to be associated with significantly impaired FGF signaling considered as a potential cause of synostoses syndrome. However, the detailed mechanism for this observation still remains unknown. In this study, we used molecular dynamics simulations and free energy calculations to study the interactions of FGF9(WT/S99N), FGFR3c, and heparin, with an aim of providing atomic sights into the detailed mechanism for the impaired FGF signaling caused by the S99N mutation. We found that the S99N mutation has a well-ordered C-terminal structure, which can reduce its homodimerization ability so as to break the monomer-dimer equilibrium in the FGF signaling, which is considered as a key factor to regulate extracellular matrix affinity and tissue diffusion in the FGF signaling pathway. The FGF9(WT) monomer can preferentially form a homodimer owing to its comparatively favorable binding free energy. In contrast, the FGF9(S99N) monomer is preferred to bind with the FGFR3c receptor to form an inactive complex, leading to impair FGF signaling. To support our computational findings, we also performed biochemical experiments, which confirm the computational results mentioned above. The impaired FGF signaling is believed to be a potential cause of human synostoses syndrome, implicating an important role for FGF9 in normal joint development.

Exploration of conformational transition in the aryl-binding site of human FXa using molecular dynamics simulations

Human coagulation Factor X (FX), a member of the vitamin K-dependent serine protease family, is a crucial component of the human coagulation cascade. Activated FX (FXa) participates in forming the prothrombinase complex on activated platelets to convert prothrombin to thrombin in coagulation reactions. In the current study, 30-ns MD simulations were performed on both the open and closed states of human FXa. Root mean squares (RMS) fluctuations showed that structural fluctuations concentrated on the loop regions of FXa, and the presence of a ligand in the closed system resulted in larger fluctuations of the gating residues. The open system had a gating distance from 9.23 to 11.33 Å, i.e., significantly larger than that of the closed system (4.69-6.35 Å), which allows diversified substrates of variable size to enter. Although the solvent accessible surface areas (SASA) of FXa remained the same in both systems, the open system generally had a larger total SASA or hydrophobic SASA (or both) for residues surrounding the S4 pocket. Additionally, more hydrogen bonds were formed in the closed state than in the open state of FXa, which is believed to play a significant role in maintaining the closed confirmation of the aryl-binding site. Based on the results of MD simulations, we propose that an induced-fit mechanism governs the functioning of human coagulation FX, which helps provide a better understanding of the interactions between FXa and its substrate, and the mechanism of the conformational changes involved in human coagulation.

Identification of a xylose reductase gene in the xylose metabolic pathway of Kluyveromyces marxianus NBRC1777

Kluyveromyces marxianus is thermotolerant yeast that is able to utilize a wider range of substrates and has greater thermal tolerance than most other yeast species. K. marxianus can assimilate xylose, but its ability to produce ethanol from xylose in oxygen-limited environments is poor. In the present study, the K. marxianus xylose reductase (KmXR) gene (Kmxyl1) was cloned and the recombinant enzyme was characterized to clarify the factors that limit xylose fermentation in K. marxianus NBRC1777. KmXR is a key enzyme in the xylose metabolism of K. marxianus, which was verified by disruption of the Kmxyl1 gene. The Km of the recombinant KmXR for NADPH is 65.67 μM and KmXR activity is 1.295 U/mg, which is lower than those of most reported yeast XRs, and the enzyme has no activity with coenzyme NADH. This result demonstrates that the XR from K. marxianus is highly coenzyme specific; combined with the extremely low XDH activity of K. marxianus with NADP+, the limitation of xylose fermentation is due to a redox imbalance under anaerobic conditions and low KmXR activity.

Insights from modeling the 3D structure of New Delhi metallo-β-lactamse and its binding interactions with antibiotic drugs

New Delhi metallo-beta-lactamase (NDM-1) is an enzyme that makes bacteria resistant to a broad range of beta-lactam antibiotic drugs. This is because it can inactivate most beta-lactam antibiotic drugs by hydrolyzing them. For in-depth understanding of the hydrolysis mechanism, the three-dimensional structure of NDM-1 was developed. With such a structural frame, two enzyme-ligand complexes were derived by respectively docking Imipenem and Meropenem (two typical beta-lactam antibiotic drugs) to the NDM-1 receptor. It was revealed from the NDM-1/Imipenem complex that the antibiotic drug was hydrolyzed while sitting in a binding pocket of NDM-1 formed by nine residues. And for the case of NDM-1/Meropenem complex, the antibiotic drug was hydrolyzed in a binding pocket formed by twelve residues. All these constituent residues of the two binding pockets were explicitly defined and graphically labeled. It is anticipated that the findings reported here may provide useful insights for developing new antibiotic drugs to overcome the resistance problem.

Investigating the coenzyme specificity of phenylacetone monooxygenase from Thermobifida fusca

Type I Baeyer–Villiger monooxygenases (BVMOs) strongly prefer NADPH over NADH as an electron donor. In order to elucidate the molecular basis for this coenzyme specificity, we have performed a site-directed mutagenesis study on phenylacetone monooxygenase (PAMO) from Thermobifida fusca. Using sequence alignments of type I BVMOs and crystal structures of PAMO and cyclohexanone monooxygenase in complex with NADP⁺, we identified four residues that could interact with the 2′-phosphate moiety of NADPH in PAMO. The mutagenesis study revealed that the conserved R217 is essential for binding the adenine moiety of the nicotinamide coenzyme while it also contributes to the recognition of the 2′-phosphate moiety of NADPH. The substitution of T218 did not have a strong effect on the coenzyme specificity. The H220N and H220Q mutants exhibited a ~3-fold improvement in the catalytic efficiency with NADH while the catalytic efficiency with NADPH was hardly affected. Mutating K336 did not increase the activity of PAMO with NADH, but it had a significant and beneficial effect on the enantioselectivity of Baeyer–Villiger oxidations and sulfoxidations. In conclusion, our results indicate that the function of NADPH in catalysis cannot be easily replaced by NADH. This finding is in line with the complex catalytic mechanism and the vital role of the coenzyme in BVMOs.

Structure and engineering of L-arabinitol 4-dehydrogenase from Neurospora crassa

L-arabinitol 4-dehydrogenase (LAD) catalyzes the conversion of l-arabinitol into l-xylulose with concomitant NAD(+) reduction. It is an essential enzyme in the development of recombinant organisms that convert l-arabinose into fuels and chemicals using the fungal l-arabinose catabolic pathway. Here we report the crystal structure of LAD from the filamentous fungus Neurospora crassa at 2.6 A resolution. In addition, we created a number of site-directed variants of N. crassa LAD that are capable of utilizing NADP(+) as cofactor, yielding the first example of LAD with an almost completely switched cofactor specificity. This work represents the first structural data on any LAD and provides a molecular basis for understanding the existing literature on the substrate specificity and cofactor specificity of this enzyme. The engineered LAD mutants with altered cofactor specificity should be useful for applications in industrial biotechnology.