Structural consequences of sequence patterns in the fingerprint region of the nucleotide binding fold. Implications for nucleotide specificity

In vitro synthesis of oryzamutaic acid H using recombinant L-lysine ε-dehydrogenase from Agrobacterium tumefaciens

Oryzamutaic acids, possessing a nitrogen-containing heterocyclic skeleton, have been isolated and identified from a rice mutant. Although oryzamutaic acids are expected to be functional ingredients, their functionality is difficult to evaluate, because of their wide variety and presence in trace amounts. Furthermore, how oryzamutaic acid is synthesized in vivo is unclear. Therefore, we developed a simple enzymatic synthesis method for these compounds in vitro. We focused on L-lysine ε-dehydrogenase (LysDH) from Agrobacterium tumefaciens, which synthesizes α-aminoadipate-δ-semialdehyde-a precursor of oryzamutaic acids. LysDH was cloned and expressed in Escherichia coli. Analysis of activity revealed that LysDH catalyzed the synthesis of oryzamutaic acid H at neutral pH in vitro. We synthesized 1.6 mg oryzamutaic acid H from 100 mg L-lysine. The synthesized oryzamutaic acid H exhibited UVA absorption, stability of temperature, and stability at a wide pH range. To our knowledge, this study is the first to report the enzymatic synthesis of oryzamutaic acid H in vitro and provides a basis for understanding the mechanisms of oryzamutaic acid synthesis in vivo.

Targeting Enterococcus faecalis HMG-CoA reductase with a non-statin inhibitor

HMG-CoA reductase (HMGR), a rate-limiting enzyme of the mevalonate pathway in Gram-positive pathogenic bacteria, is an attractive target for development of novel antibiotics. In this study, we report the crystal structures of HMGR from Enterococcus faecalis (efHMGR) in the apo and liganded forms, highlighting several unique features of this enzyme. Statins, which inhibit the human enzyme with nanomolar affinity, perform poorly against the bacterial HMGR homologs. We also report a potent competitive inhibitor (Chembridge2 ID 7828315 or compound 315) of the efHMGR enzyme identified by a high-throughput, in-vitro screening. The X-ray crystal structure of efHMGR in complex with 315 was determined to 1.27 Å resolution revealing that the inhibitor occupies the mevalonate-binding site and interacts with several key active site residues conserved among bacterial homologs. Importantly, 315 does not inhibit the human HMGR. Our identification of a selective, non-statin inhibitor of bacterial HMG-CoA reductases will be instrumental in lead optimization and development of novel antibacterial drug candidates.

Bioinformatics-Structural Approach to the Search for New D-Amino Acid Oxidases

D-amino acid oxidase (DAAO, EC 1.2.1.2) plays an important role in the functioning of prokaryotes as well as of lower (yeast and fungi) and higher eukaryotes (mammals). DAAO genes have not yet been found in archaean genomes. D-amino acid oxidase is increasingly used in various fields, which requires the development of new variants of the enzyme with specific properties. However, even within one related group (bacteria, yeasts and fungi, mammals), DAAOs show very low homology between amino acid sequences. In particular, this fact is clearly observed in the case of DAAO from bacteria. The high variability in the primary structures of DAAO severely limits the search for new enzymes in known genomes. As a result, many (if not most) DAAO genes remain either unannotated or incorrectly annotated. We propose an approach that uses bioinformatic methods in combination with general 3D structure and active center structure analysis to confirm that the gene found encodes D-amino acid oxidase and to predict the possible type of its substrate specificity. Using a homology search, we obtained a set of candidate sequences, modelled the tertiary structure of the selected enzymes, and compared them with experimental and model structures of known DAAOs. The effectiveness of the proposed approach for discrimination of DAAOs and glycine oxidases is shown. Using this approach, new DAAO genes were found in the genomes of six strains of extremophilic bacteria, and for the first time in the world, one gene was identified in the genome of halophilic archaea. Preliminary experiments confirmed the predicted specificity of DAAO from Natronosporangium hydrolyticum ACPA39 with D-Leu and D-Phe.

Structure-function characterization of an aldo-keto reductase involved in detoxification of the mycotoxin, deoxynivalenol

Deoxynivalenol (DON) is a mycotoxin, produced by filamentous fungi such as Fusarium graminearum, that causes significant yield losses of cereal grain crops worldwide. One of the most promising methods to detoxify this mycotoxin involves its enzymatic epimerization to 3-epi-DON. DepB plays a critical role in this process by reducing 3-keto-DON, an intermediate in the epimerization process, to 3-epi-DON. DepB_Rleg from Rhizobium leguminosarum is a member of the new aldo-keto reductase family, AKR18, and it has the unusual ability to utilize both NADH and NADPH as coenzymes, albeit with a 40-fold higher catalytic efficiency with NADPH compared to NADH. Structural analysis of DepB_Rleg revealed the putative roles of Lys-217, Arg-290, and Gln-294 in NADPH specificity. Replacement of these residues by site-specific mutagenesis to negatively charged amino acids compromised NADPH binding with minimal effects on NADH binding. The substrate-binding site of DepB_Rleg is larger than its closest structural homolog, AKR6A2, likely contributing to its ability to utilize a wide range of aldehydes and ketones, including the mycotoxin, patulin, as substrates. The structure of DepB_Rleg also suggests that 3-keto-DON can adopt two binding modes to facilitate 4-pro-R hydride transfer to either the re- or si-face of the C3 ketone providing a possible explanation for the enzyme's ability to convert 3-keto-DON to 3-epi-DON and DON in diastereomeric ratios of 67.2% and 32.8% respectively.

Agl24 is an ancient archaeal homolog of the eukaryotic N-glycan chitobiose synthesis enzymes

Protein N-glycosylation is a post-translational modification found in organisms of all domains of life. The crenarchaeal N-glycosylation begins with the synthesis of a lipid-linked chitobiose core structure, identical to that in Eukaryotes, although the enzyme catalyzing this reaction remains unknown. Here, we report the identification of a thermostable archaeal β-1,4-N-acetylglucosaminyltransferase, named archaeal glycosylation enzyme 24 (Agl24), responsible for the synthesis of the N-glycan chitobiose core. Biochemical characterization confirmed its function as an inverting β-D-GlcNAc-(1→4)-α-D-GlcNAc-diphosphodolichol glycosyltransferase. Substitution of a conserved histidine residue, found also in the eukaryotic and bacterial homologs, demonstrated its functional importance for Agl24. Furthermore, bioinformatics and structural modeling revealed similarities of Agl24 to the eukaryotic Alg14/13 and a distant relation to the bacterial MurG, which are catalyzing the same or a similar reaction, respectively. Phylogenetic analysis of Alg14/13 homologs indicates that they are ancient in Eukaryotes, either as a lateral transfer or inherited through eukaryogenesis.

Crystallographic analysis of TarI and TarJ, a cytidylyltransferase and reductase pair for CDP-ribitol synthesis in Staphylococcus aureus wall teichoic acid biogenesis

The cell wall of many pathogenic Gram-positive bacteria contains ribitol-phosphate wall teichoic acid (WTA), a polymer that is linked to virulence and regulation of essential physiological processes including cell division. CDP-ribitol, the activated precursor for ribitol-phosphate polymerization, is synthesized by a cytidylyltransferase and reductase pair known as TarI and TarJ, respectively. In this study, we present crystal structures of Staphylococcus aureus TarI and TarJ in their apo forms and in complex with substrates and products. The TarI structures illustrate the mechanism of CDP-ribitol synthesis from CTP and ribitol-phosphate and reveal structural changes required for substrate binding and catalysis. Insights into the upstream step of ribulose-phosphate reduction to ribitol-phosphate is provided by the structures of TarJ. Furthermore, we propose a general topology of the enzymes in a heterotetrameric form built using restraints from crosslinking mass spectrometry analysis. Together, our data present molecular details of CDP-ribitol production that may aid in the design of inhibitors against WTA biosynthesis.

Purification and Characterization of (2R,3R)-2,3-Butanediol Dehydrogenase of the Human Pathogen Neisseria gonorrhoeae FA1090 Produced in Escherichia coli

2,3-Butanediol dehydrogenase (BDH), also known as acetoin/diacetyl reductase, is a pivotal enzyme for the formation of 2,3-butanediol (2,3-BD), a chiral compound with potential roles in the virulence of certain pathogens. Here, a NAD(H)-dependent (2R,3R)-BDH from Neisseria gonorrhoeae FA1090 (NgBDH), the causative agent of gonorrhoea, was functionally characterized. Sequence analysis indicated that it belongs to zinc-containing medium-chain dehydrogenase/reductase family. The recombinant NgBDH migrated as a single band with a size of around 45 kDa on SDS-PAGE and could be confirmed by Western blotting and mass spectrometry. For the oxidation of either (2R,3R)-2,3-BD or meso-2,3-BD, the enzyme exhibited a broad pH optimum between pH 9.5 to 11.5. For the reduction of (3R/3S)-acetoin, the pH optimum was around 6.5. The enzyme could catalyze the stereospecific oxidation of (2R,3R)-2,3-BD (K_m = 0.16 mM, k_cat/K_m = 673 s^-1 · mM^-1) and meso-BD (K_m = 0.72 mM, k_cat/K_m = 165 s^-1 · mM^-1). Moreover, it could also reduce (3R/3S)-acetoin with a K_m of 0.14 mM and a k_cat/K_m of 885 s^-1 · mM^-1. The results presented here contribute to understand the 2,3-BD metabolism in N. gonorrhoeae and pave the way for studying the influence of 2,3-BD metabolism on the virulence of this pathogen in the future.

Structural basis for the catalytic mechanism and α-ketoglutarate cooperativity of glutamate dehydrogenase

Glutamate dehydrogenase (GDH) is a key enzyme connecting carbon and nitrogen metabolism in all living organisms. Despite extensive studies on GDHs from both prokaryotic and eukaryotic organisms in the last 40 years, the structural basis of the catalytic features of this enzyme remains incomplete. This study reports the structural basis of the GDH catalytic mechanism and allosteric behavior. We determined the first high-resolution crystal structures of glutamate dehydrogenase from the fungus Aspergillus niger (AnGDH), a unique NADP⁺-dependent allosteric enzyme that is forward-inhibited by the formation of mixed disulfide. We determined the structures of the active enzyme in its apo form and in binary/ternary complexes with bound substrate (α-ketoglutarate), inhibitor (isophthalate), coenzyme (NADPH), or two reaction intermediates (α-iminoglutarate and 2-amino-2-hydroxyglutarate). The structure of the forward-inhibited enzyme (fiAnGDH) was also determined. The hexameric AnGDH had three open subunits at one side and three partially closed protomers at the other, a configuration not previously reported. The AnGDH hexamers having subunits with different conformations indicated that its α-ketoglutarate-dependent homotropic cooperativity follows the Monod-Wyman-Changeux (MWC) model. Moreover, the position of the water attached to Asp-154 and Gly-153 defined the previously unresolved ammonium ion-binding pocket, and the binding site for the 2'-phosphate group of the coenzyme was also better defined by our structural data. Additional structural and mutagenesis experiments identified the residues essential for coenzyme recognition. This study reveals the structural features responsible for positioning α-ketoglutarate, NADPH, ammonium ion, and the reaction intermediates in the GDH active site.

Unique coenzyme binding mode of hyperthermophilic archaeal sn-glycerol-1-phosphate dehydrogenase from Pyrobaculum calidifontis

A gene encoding an sn-glycerol-1-phosphate dehydrogenase (G1PDH) was identified in the hyperthermophilic archaeon Pyrobaculum calidifontis. The gene was overexpressed in Escherichia coli, and its product was purified and characterized. In contrast to conventional G1PDHs, the expressed enzyme showed strong preference for NADH: the reaction rate (V_max ) with NADPH was only 2.4% of that with NADH. The crystal structure of the enzyme was determined at a resolution of 2.45 Å. The asymmetric unit consisted of one homohexamer. Refinement of the structure and HPLC analysis showed the presence of the bound cofactor NADPH in subunits D, E, and F, even though it was not added in the crystallization procedure. The phosphate group at C2' of the adenine ribose of NADPH is tightly held through the five biased hydrogen bonds with Ser40 and Thr42. In comparison with the known G1PDH structure, the NADPH molecule was observed to be pushed away from the normal coenzyme binding site. Interestingly, the S40A/T42A double mutant enzyme acquired much higher reactivity than the wild-type enzyme with NADPH, which suggests that the biased interactions around the C2'-phosphate group make NADPH binding insufficient for catalysis. Our results provide a unique structural basis for coenzyme preference in NAD(P)-dependent dehydrogenases. Proteins 2016; 84:1786-1796. © 2016 Wiley Periodicals, Inc.

Crystal Structures of a Hyperthermophilic Archaeal Homoserine Dehydrogenase Suggest a Novel Cofactor Binding Mode for Oxidoreductases

NAD(P)-dependent dehydrogenases differ according to their coenzyme preference: some prefer NAD, others NADP, and still others exhibit dual cofactor specificity. The structure of a newly identified archaeal homoserine dehydrogenase showed this enzyme to have a strong preference for NADP. However, NADP did not act as a cofactor with this enzyme, but as a strong inhibitor of NAD-dependent homoserine oxidation. Structural analysis and site-directed mutagenesis showed that the large number of interactions between the cofactor and the enzyme are responsible for the lack of reactivity of the enzyme towards NADP. This observation suggests this enzyme exhibits a new variation on cofactor binding to a dehydrogenase: very strong NADP binding that acts as an obstacle to NAD(P)-dependent dehydrogenase catalytic activity.

Enantioselective oxidation of galactitol 1-phosphate by galactitol-1-phosphate 5-dehydrogenase from Escherichia coli

Galactitol-1-phosphate 5-dehydrogenase (GPDH) is a polyol dehydrogenase that belongs to the medium-chain dehydrogenase/reductase (MDR) superfamily. It catalyses the Zn(2+)- and NAD(+)-dependent stereoselective dehydrogenation of L-galactitol 1-phosphate to D-tagatose 6-phosphate. Here, three crystal structures of GPDH from Escherichia coli are reported: that of the open state of GPDH with Zn(2+) in the catalytic site and those of the closed state in complex with the polyols Tris and glycerol, respectively. The closed state of GPDH reveals no bound cofactor, which is at variance with the conformational transition of the prototypical mammalian liver alcohol dehydrogenase. The main intersubunit-contacting interface within the GPDH homodimer presents a large internal cavity that probably facilitates the relative movement between the subunits. The substrate analogue glycerol bound within the active site partially mimics the catalytically relevant backbone of galactitol 1-phosphate. The glycerol binding mode reveals, for the first time in the polyol dehydrogenases, a pentacoordinated zinc ion in complex with a polyol and also a strong hydrogen bond between the primary hydroxyl group and the conserved Glu144, an interaction originally proposed more than thirty years ago that supports a catalytic role for this acidic residue.

Substitutions at the cofactor phosphate-binding site of a clostridial alcohol dehydrogenase lead to unexpected changes in substrate specificity

Changing the cofactor specificity of an enzyme from nicotinamide adenine dinucleotide 2′-phosphate (NADPH) to the more abundant NADH is a common strategy for increasing overall enzyme efficiency in microbial metabolic engineering. The aim of this study was to switch the cofactor specificity of the primary–secondary alcohol dehydrogenase from Clostridium autoethanogenum, a bacterium with considerable promise for the bio-manufacturing of fuels and other petrochemicals, from strictly NADPH-dependent to NADH-dependent. We used insights from a homology model to build a site-saturation library focussed on residue S199, the position deemed most likely to disrupt binding of the 2′-phosphate of NADPH. Although the CaADH(S199X) library did not yield any NADH-dependent enzymes, it did reveal that substitutions at the cofactor phosphate-binding site can cause unanticipated changes in the substrate specificity of the enzyme. Using consensus-guided site-directed mutagenesis, we were able to create an enzyme that was stringently NADH-dependent, albeit with a concomitant reduction in activity. This study highlights the role that distal residues play in substrate specificity and the complexity of enzyme–cofactor interactions.

New insights into the evolutionary history of plant sorbitol dehydrogenase

BACKGROUND

Sorbitol dehydrogenase (SDH, EC 1.1.1.14) is the key enzyme involved in sorbitol metabolism in higher plants. SDH genes in some Rosaceae species could be divided into two groups. L-idonate-5-dehydrogenase (LIDH, EC 1.1.1.264) is involved in tartaric acid (TA) synthesis in Vitis vinifera and is highly homologous to plant SDHs. Despite efforts to understand the biological functions of plant SDH, the evolutionary history of plant SDH genes and their phylogenetic relationship with the V. vinifera LIDH gene have not been characterized.

RESULTS

A total of 92 SDH genes were identified from 42 angiosperm species. SDH genes have been highly duplicated within the Rosaceae family while monocot, Brassicaceae and most Asterid species exhibit singleton SDH genes. Core Eudicot SDHs have diverged into two phylogenetic lineages, now classified as SDH Class I and SDH Class II. V. vinifera LIDH was identified as a Class II SDH. Tandem duplication played a dominant role in the expansion of plant SDH family and Class II SDH genes were positioned in tandem with Class I SDH genes in several plant genomes. Protein modelling analyses of V. vinifera SDHs revealed 19 putative active site residues, three of which exhibited amino acid substitutions between Class I and Class II SDHs and were influenced by positive natural selection in the SDH Class II lineage. Gene expression analyses also demonstrated a clear transcriptional divergence between Class I and Class II SDH genes in V. vinifera and Citrus sinensis (orange).

CONCLUSIONS

Phylogenetic, natural selection and synteny analyses provided strong support for the emergence of SDH Class II by positive natural selection after tandem duplication in the common ancestor of core Eudicot plants. The substitutions of three putative active site residues might be responsible for the unique enzyme activity of V. vinifera LIDH, which belongs to SDH Class II and represents a novel function of SDH in V. vinifera that may be true also of other Class II SDHs. Gene expression analyses also supported the divergence of SDH Class II at the expression level. This study will facilitate future research into understanding the biological functions of plant SDHs.

Residues that influence coenzyme preference in the aldehyde dehydrogenases

To find out the residues that influence the coenzyme preference of aldehyde dehydrogenases (ALDHs), we reviewed, analyzed and correlated data from their known crystal structures and amino-acid sequences with their published kinetic parameters for NAD(P)(+). We found that the conformation of the Rossmann-fold loops participating in binding the adenosine ribose is very conserved among ALDHs, so that coenzyme specificity is mainly determined by the nature of the residue at position 195 (human ALDH2 numbering). Enzymes with glutamate or proline at 195 prefer NAD(+) because the side-chains of these residues electrostatically and/or sterically repel the 2'-phosphate group of NADP(+). But contrary to the conformational rigidity of proline, the conformational flexibility of glutamate may allow NADP(+)-binding in some enzymes by moving the carboxyl group away from the 2'-phosphate group, which is possible if a small neutral residue is located at position 224, and favored if the residue at position 53 interacts with Glu195 in a NADP(+)-compatible conformation. Of the residues found at position 195, only glutamate interacts with the NAD(+)-adenosine ribose; glutamine and histidine cannot since their side-chain points are opposite to the ribose, probably because the absence of the electrostatic attraction by the conserved nearby Lys192, or its electrostatic repulsion, respectively. The shorter side-chains of other residues-aspartate, serine, threonine, alanine, valine, leucine, or isoleucine-are distant from the ribose but leave room for binding the 2'-phosphate group. Generally, enzymes having a residue different from Glu bind NAD(+) with less affinity, but they can also bind NADP(+) even sometimes with higher affinity than NAD(+), as do enzymes containing Thr/Ser/Gln195. Coenzyme preference is a variable feature within many ALDH families, consistent with being mainly dependent on a single residue that apparently has no other structural or functional roles, and therefore can easily be changed through evolution and selected in response to physiological needs.

The shikimate dehydrogenase family: functional diversity within a conserved structural and mechanistic framework

Shikimate dehydrogenase (SDH) catalyzes the NADPH-dependent reduction of 3-deydroshikimate to shikimate, an essential reaction in the biosynthesis of the aromatic amino acids and a large number of other secondary metabolites in plants and microbes. The indispensible nature of this enzyme makes it a potential target for herbicides and antimicrobials. SDH is the archetypal member of a large protein family, which contains at least four additional functional classes with diverse metabolic roles. The different members of the SDH family share a highly similar three-dimensional structure and utilize a conserved catalytic mechanism, but exhibit distinct substrate preferences, making the family a particularly interesting system for studying modes of substrate recognition used by enzymes. Here, we review our current understanding of the biochemical and structural properties of each of the five previously identified SDH family functional classes.

Uncovering rare NADH-preferring ketol-acid reductoisomerases

All members of the ketol-acid reductoisomerase (KARI) enzyme family characterized to date have been shown to prefer the nicotinamide adenine dinucleotide phosphate hydride (NADPH) cofactor to nicotinamide adenine dinucleotide hydride (NADH). However, KARIs with the reversed cofactor preference are desirable for industrial applications, including anaerobic fermentation to produce branched-chain amino acids. By applying insights gained from structural and engineering studies of this enzyme family to a comprehensive multiple sequence alignment of KARIs, we identified putative NADH-utilizing KARIs and characterized eight whose catalytic efficiencies using NADH were equal to or greater than NADPH. These are the first naturally NADH-preferring KARIs reported and demonstrate that this property has evolved independently multiple times, using strategies unlike those used previously in the laboratory to engineer a KARI cofactor switch.

Improving the NADH-cofactor specificity of the highly active AdhZ3 and AdhZ2 from Escherichia coli K-12

Biocatalysis is a promising tool for the sustainable production of chemicals. When cofactor depending enzymatic reactions are involved the applicability of the right cofactor is a central issue. One important example in this regard is the production of alcohols by nicotinamide cofactor (NAD(P)(+)) depending alcohol dehydrogenases. AdhZ3 from Escherichia coli, which is important for the production of alcohols from biomass, has a preference for NADPH as cofactor. We used a structure guided site-specific random approach, to change the cofactor preference towards NADH and to deduce more general rules for redesigning the cofactor specificity. Transfer of a triplet motif from NADH preferring horse liver ADH to AdhZ3 showed an insufficient switch in the preference towards NADH. A combinatorial site saturation mutagenesis altering three residues at once was applied. Library screening with two different cofactor concentrations (0.1 and 0.3mM) resulted in nine improved variants with AdhZ3-LND having the highest vmax and AdhZ3-CND having the lowest K(m). Asparagine was the most frequent amino acid found in eight of nine triplet motifs. To verify the triplet-motif, two variants of E. coli AdhZ2 DIN and LND were designed and confirmed for improved activity with NADH.

Cofactory: sequence-based prediction of cofactor specificity of Rossmann folds

Obtaining optimal cofactor balance to drive production is a challenge in metabolically engineered microbial production strains. To facilitate identification of heterologous enzymes with desirable altered cofactor requirements from native content, we have developed Cofactory, a method for prediction of enzyme cofactor specificity using only primary amino acid sequence information. The algorithm identifies potential cofactor binding Rossmann folds and predicts the specificity for the cofactors FAD(H2), NAD(H), and NADP(H). The Rossmann fold sequence search is carried out using hidden Markov models whereas artificial neural networks are used for specificity prediction. Training was carried out using experimental data from protein-cofactor structure complexes. The overall performance was benchmarked against an independent evaluation set obtaining Matthews correlation coefficients of 0.94, 0.79, and 0.65 for FAD(H2), NAD(H), and NADP(H), respectively. The Cofactory method is made publicly available at http://www.cbs.dtu.dk/services/Cofactory.

Molecular characterization of an NADPH-dependent acetoin reductase/2,3-butanediol dehydrogenase from Clostridium beijerinckii NCIMB 8052

Acetoin reductase is an important enzyme for the fermentative production of 2,3-butanediol, a chemical compound with a very broad industrial use. Here, we report on the discovery and characterization of an acetoin reductase from Clostridium beijerinckii NCIMB 8052. An in silico screen of the C. beijerinckii genome revealed eight potential acetoin reductases. One of them (CBEI_1464) showed substantial acetoin reductase activity after expression in Escherichia coli. The purified enzyme (C. beijerinckii acetoin reductase [Cb-ACR]) was found to exist predominantly as a homodimer. In addition to acetoin (or 2,3-butanediol), other secondary alcohols and corresponding ketones were converted as well, provided that another electronegative group was attached to the adjacent C-3 carbon. Optimal activity was at pH 6.5 (reduction) and 9.5 (oxidation) and around 68°C. Cb-ACR accepts both NADH and NADPH as electron donors; however, unlike closely related enzymes, NADPH is preferred (Km, 32 μM). Cb-ACR was compared to characterized close homologs, all belonging to the "threonine dehydrogenase and related Zn-dependent dehydrogenases" (COG1063). Metal analysis confirmed the presence of 2 Zn(2+) atoms. To gain insight into the substrate and cofactor specificity, a structural model was constructed. The catalytic zinc atom is likely coordinated by Cys37, His70, and Glu71, while the structural zinc site is probably composed of Cys100, Cys103, Cys106, and Cys114. Residues determining NADP specificity were predicted as well. The physiological role of Cb-ACR in C. beijerinckii is discussed.

Structure of NADP(+)-dependent glutamate dehydrogenase from Escherichia coli--reflections on the basis of coenzyme specificity in the family of glutamate dehydrogenases

Glutamate dehydrogenases (GDHs; EC 1.4.1.2-4) catalyse the oxidative deamination of L-glutamate to α-ketoglutarate, using NAD(+) and/or NADP(+) as a cofactor. Subunits of homo-hexameric bacterial enzymes comprise a substrate-binding domain I followed by a nucleotide-binding domain II. The reaction occurs in a catalytic cleft between the two domains. Although conserved residues in the nucleotide-binding domains of various dehydrogenases have been linked to cofactor preferences, the structural basis for specificity in the GDH family remains poorly understood. Here, the refined crystal structure of Escherichia coli GDH in the absence of reactants is described at 2.5-Å resolution. Modelling of NADP(+) in domain II reveals the potential contribution of positively charged residues from a neighbouring α-helical hairpin to phosphate recognition. In addition, a serine that follows the P7 aspartate is presumed to form a hydrogen bond with the 2'-phosphate. Mutagenesis and kinetic analysis confirms the importance of these residues in NADP(+) recognition. Surprisingly, one of the positively charged residues is conserved in all sequences of NAD(+)-dependent enzymes, but the conformations adopted by the corresponding regions in proteins whose structure has been solved preclude their contribution to the coordination of the 2'-ribose phosphate of NADP(+). These studies clarify the sequence-structure relationships in bacterial GDHs, revealing that identical residues may specify different coenzyme preferences, depending on the structural context. Primary sequence alone is therefore not a reliable guide for predicting coenzyme specificity. We also consider how it is possible for a single sequence to accommodate both coenzymes in the dual-specificity GDHs of animals.

Glutamate dehydrogenases: the why and how of coenzyme specificity

NAD(+) and NADP(+), chemically similar and with almost identical standard oxidation-reduction potentials, nevertheless have distinct roles, NAD(+) serving catabolism and ATP generation whereas NADPH is the biosynthetic reductant. Separating these roles requires strict specificity for one or the other coenzyme for most dehydrogenases. In many organisms this holds also for glutamate dehydrogenases (GDH), NAD(+)-dependent for glutamate oxidation, NADP(+)-dependent for fixing ammonia. In higher animals, however, GDH has dual specificity. It has been suggested that GDH in mitochondria reacts only with NADP(H), the NAD(+) reaction being an in vitro artefact. However, contrary evidence suggests mitochondrial GDH not only reacts with NAD(+) but maintains equilibrium using the same pool as accessed by β-hydroxybutyrate dehydrogenase. Another complication is the presence of an energy-linked dehydrogenase driving NADP(+) reduction by NADH, maintaining the coenzyme pools at different oxidation-reduction potentials. Its coexistence with GDH makes possible a futile cycle, control of which is not yet properly explained. Structural studies show NAD(+)-dependent, NADP(+)-dependent and dual-specificity GDHs are closely related and a few site-directed mutations can reverse specificity. Specificity for NAD(+) or for NADP(+) has probably emerged repeatedly during evolution, using different structural solutions on different occasions. In various GDHs the P7 position in the coenzyme-binding domain plays a key role. However, whereas in other dehydrogenases an acidic P7 residue usually hydrogen bonds to the 2'- and 3'-hydroxyls, dictating NAD(+) specificity, among GDHs, depending on detailed conformation of surrounding residues, an acidic P7 may permit binding of NAD(+) only, NADP(+) only, or in higher animals both.

Novel birch pollen specific immunotherapy formulation based on contiguous overlapping peptides

Background

Synthetic contiguous overlapping peptides (COPs) may represent an alternative to allergen extracts or recombinant allergens for allergen specific immunotherapy. In combination, COPs encompass the entire allergen sequence, providing all potential T cell epitopes, while preventing IgE conformational epitopes of the native allergen.

Methods

Individual COPs were derived from the sequence of Bet v 1, the major allergen of birch pollen, and its known crystal structure, and designed to avoid IgE binding. Three sets of COPs were tested in vitro in competition ELISA and basophil degranulation assays. Their in vivo reactivity was determined by intraperitoneal challenge in rBet v 1 sensitized mice as well as by skin prick tests in volunteers with allergic rhinoconjunctivitis to birch pollen.

Results

The combination, named AllerT, of three COPs selected for undetectable IgE binding in competition assays and for the absence of basophil activation in vitro was unable to induce anaphylaxis in sensitized mice in contrast to rBet v 1. In addition no positive reactivity to AllerT was observed in skin prick tests in human volunteers allergic to birch pollen. In contrast, a second set of COPs, AllerT4-T5 displayed some residual IgE binding in competition ELISA and a weak subliminal reactivity to skin prick testing.

Conclusions

The hypoallergenicity of contiguous overlapping peptides was confirmed by low, if any, IgE binding activity in vitro, by the absence of basophil activation and the absence of in vivo induction of allergic reactions in mouse and human.

Trial registration

ClinicalTrials.gov NCT01719133

Homology cloning, sequence characterization, and expression analysis of cDNA encoding electron transfer flavoprotein beta polypeptide in mud crab (Scylla paramamosain)

Electron transfer flavoproteins (ETFs) are αβ-heterodimers found in eukaryotic mitochondria and bacteria. Herein we report a full-length complementary DNA of a mud crab (Scylla paramamosain) ETF β subunit (Scpa-ETFB) isolated with a homology cloning strategy. The complete complementary DNA of the Scpa-ETFB contains a 17-nt 5'-untranslated region, a 765-nt open reading frame encoding 254 amino acids, and a 248-nt 3'-untranslated region. The high identity of Scpa-ETFB with ETFB in other organisms indicated that Scpa-ETFB is a new member of the ETFB family. Although the conserved motif associated with flavin adenine dinucleotide binding is absent in Scpa-ETFB, the signature sequences of the ETF superfamily were identified. Using reverse transcriptase polymerase chain reaction, we detected the messenger RNA transcript of Scpa-ETFB in high levels in the tissues of the hepatopancreas, ovary, heart, and muscle. Phylogenetic analysis showed that Scpa-ETFB is most closely related to the ETFB genes of Caligus rogercresseyi and Lepeophtheirus salmonis. These results provided basic information for elucidating the molecular mechanism of energy production in the mud crab.

Crystal structure of NAD+-dependent Peptoniphilus asaccharolyticus glutamate dehydrogenase reveals determinants of cofactor specificity

Glutamate dehydrogenases (EC 1.4.1.2-4) catalyse the oxidative deamination of l-glutamate to α-ketoglutarate using NAD(P) as a cofactor. The bacterial enzymes are hexamers and each polypeptide consists of an N-terminal substrate-binding (Domain I) followed by a C-terminal cofactor-binding segment (Domain II). The reaction takes place at the junction of the two domains, which move as rigid bodies and are presumed to narrow the cleft during catalysis. Distinct signature sequences in the nucleotide-binding domain have been linked to NAD(+) vs. NADP(+) specificity, but they are not unambiguous predictors of cofactor preferences. Here, we have determined the crystal structure of NAD(+)-specific Peptoniphilus asaccharolyticus glutamate dehydrogenase in the apo state. The poor quality of native crystals was resolved by derivatization with selenomethionine, and the structure was solved by single-wavelength anomalous diffraction methods. The structure reveals an open catalytic cleft in the absence of substrate and cofactor. Modeling of NAD(+) in Domain II suggests that a hydrophobic pocket and polar residues contribute to nucleotide specificity. Mutagenesis and isothermal titration calorimetry studies of a critical glutamate at the P7 position of the core fingerprint confirms its role in NAD(+) binding. Finally, the cofactor binding site is compared with bacterial and mammalian enzymes to understand how the amino acid sequences and three-dimensional structures may distinguish between NAD(+) vs. NADP(+) recognition.

An Examination by Site-Directed Mutagenesis of Putative Key Residues in the Determination of Coenzyme Specificity in Clostridial NAD-Dependent Glutamate Dehydrogenase

Sequence and structure comparisons of various glutamate dehydrogenases (GDH) and other nicotinamide nucleotide-dependent dehydrogenases have potentially implicated certain residues in coenzyme binding and discrimination. We have mutated key residues in Clostridium symbiosum NAD⁺-specific GDH to investigate their contribution to specificity and to enhance acceptance of NADPH. Comparisons with E. coli NADPH-dependent GDH prompted design of mutants F238S, P262S, and F238S/P262S, which were purified and assessed at pH 6.0, 7.0, and 8.0. They showed markedly increased catalytic efficiency with NADPH, especially at pH 8.0 (∼170-fold for P262S and F238S/P262S with relatively small changes for NADH). A positive charge introduced through the D263K mutation also greatly increased catalytic efficiency with NADPH (over 100-fold at pH 8) and slightly decreased activity with NADH. At position 242, “P6” of the “core fingerprint,” where NAD⁺- and NADP⁺-dependent enzymes normally have Gly or Ala, respectively, clostridial GDH already has Ala. Replacement with Gly produced negligible shift in coenzyme specificity.

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.

A hypoxia-induced positive feedback loop promotes hypoxia-inducible factor 1alpha stability through miR-210 suppression of glycerol-3-phosphate dehydrogenase 1-like

Oxygen-dependent regulation of the transcription factor HIF-1α relies on a family of prolyl hydroxylases (PHDs) that hydroxylate hypoxia-inducible factor 1α (HIF-1α) protein at two prolines during normal oxygen conditions, resulting in degradation by the proteasome. During low-oxygen conditions, these prolines are no longer hydroxylated and HIF-1α degradation is blocked. Hypoxia-induced miRNA-210 (miR-210) is a direct transcriptional target of HIF-1α, but its complete role and targets during hypoxia are not well understood. Here, we identify the enzyme glycerol-3-phosphate dehydrogenase 1-like (GPD1L) as a novel regulator of HIF-1α stability and a direct target of miR-210. Expression of miR-210 results in stabilization of HIF-1α due to decreased levels of GPD1L resulting in an increase in HIF-1α target genes. Altering GPD1L levels by overexpression or knockdown results in a decrease or increase in HIF-1α stability, respectively. GPD1L-mediated decreases in HIF-1α stability can be reversed by pharmacological inhibition of the proteasome or PHD activity. When rescued from degradation by proteasome inhibition, elevated amounts of GPD1L cause hyperhydroxylation of HIF-1α, suggesting increases in PHD activity. Importantly, expression of GPD1L attenuates the hypoxic response, preventing complete HIF-1α induction. We propose a model in which hypoxia-induced miR-210 represses GPD1L, contributing to suppression of PHD activity, and increases of HIF-1α protein levels.

Combinatorial analysis for sequence and spatial motif discovery in short sequence fragments

Motifs are overrepresented sequence or spatial patterns appearing in proteins. They often play important roles in maintaining protein stability and in facilitating protein function. When motifs are located in short sequence fragments, as in transmembrane domains that are only 6-20 residues in length, and when there is only very limited data, it is difficult to identify motifs. In this study, we introduce combinatorial models based on permutation for assessing statistically significant sequence and spatial patterns in short sequences. We show that our method can uncover previously unknown sequence and spatial motifs in beta-barrel membrane proteins and that our method outperforms existing methods in detecting statistically significant motifs in this data set. Last, we discuss implications of motif analysis for problems involving short sequences in other families of proteins.

Molecular modeling studies of L-arabinitol 4-dehydrogenase of Hypocrea jecorina: its binding interactions with substrate and cofactor

L-arabinitol 4-dehydrogenase (LAD1; EC 1.1.1.12) is an enzyme in the L-arabinose catabolic pathway of fungi that catalyzes the conversion of L-arabinitol into L-xylulose. The primary objective of this work is to identify the catalytic and coenzyme binding domains of LAD1 from Hypocrea jecorina in order to provide better insight into the possible catalytic events in these domains. The 3D structure of NAD(+)-dependent LAD1 was developed based on the crystal structure of human sorbitol dehydrogenase as a template. A series of molecular mechanics and dynamics operations were performed to find the most stable binding interaction for the enzyme and its ligands. Using the verified model, a docking study was performed with the substrate L-arabinitol, Zn(2+) and NAD(+). This study found a catalytic Zn(2+) binding domain (Cys66, His91, Glu92 and Glu176) and a cofactor NAD(+) binding domain (Gly202, ILeu204, Gly205, Cys273, Arg229 and Val298) with strong hydrogen bonding contacts with the substrate and cofactor. The binding pockets of the enzyme for l-arabinitol, NAD(+), and Zn(2+) have been explicitly defined. The results from this study should guide future mutagenesis studies and provide useful clues for engineering enzymes to improve the utilization of polyols for rare sugar production.

First crystal structure of L-lysine 6-dehydrogenase as an NAD-dependent amine dehydrogenase

A gene encoding an L-lysine dehydrogenase was identified in the hyperthermophilic archaeon Pyrococcus horikoshii. The gene was overexpressed in Escherichia coli, and its product was purified and characterized. The expressed enzyme is the most thermostable L-lysine dehydrogenase yet described, with a half-life of 180 min at 100 degrees C. The product of the enzyme's catalytic activity is Delta(1)-piperideine-6-carboxylate, which makes this enzyme an L-lysine 6-dehydrogenase (EC 1.4.1.18) that catalyzes the reductive deamination of the epsilon- amino group and a type of NAD-dependent amine dehydrogenase. The three-dimensional structure of the enzyme was determined using the mercury-based multiple-wavelength anomalous dispersion method at a resolution of 2.44 A in the presence of NAD and sulfate ion. The asymmetric unit consisted of two subunits, and a crystallographic 2-fold axis generated the functional dimer. Each monomer consisted of a Rossmann fold domain and a C-terminal catalytic domain, and the fold of the catalytic domain showed similarity to that of saccharopine reductase. Notably, the structures of subunits A and B differed significantly. In subunit A, the active site contained a sulfate ion that was not seen in subunit B. Consequently, subunit A adopted a closed conformation, whereas subunit B adopted an open one. In each subunit, one NAD molecule was bound to the active site in an anti-conformation, indicating that the enzyme makes use of pro-R-specific hydride transfer between the two hydrides at C-4 of NADH (type A specificity). This is the first description of the three-dimensional structure of l-lysine 6-dehydrogenase as an NAD-dependent amine dehydrogenase.

Refolding, characterization and crystal structure of (S)-malate dehydrogenase from the hyperthermophilic archaeon Aeropyrum pernix

Tartrate oxidation activity was found in the crude extract of an aerobic hyperthermophilic archaeon Aeropyrum pernix, and the enzyme was identified as (S)-malate dehydrogenase (MDH), which, when produced in Escherichia coli, was mainly obtained as an inactive inclusion body. The inclusion body was dissolved in 6 M guanidine-HCl and gradually refolded to the active enzyme through dilution of the denaturant. The purified recombinant enzyme consisted of four identical subunits with a molecular mass of about 110 kDa. NADP was preferred as a coenzyme over NAD for (S)-malate oxidation and, unlike MDHs from other sources, this enzyme readily catalyzed the oxidation of (2S,3S)-tartrate and (2S,3R)-tartrate. The tartrate oxidation activity was also observed in MDHs from the hyperthermophilic archaea Methanocaldococcus jannaschii and Archaeoglobus fulgidus, suggesting these hyperthermophilic MDHs loosely bind their substrates. The refolded A. pernix MDH was also crystallized, and the structure was determined at a resolution of 2.9 A. Its overall structure was similar to those of the M. jannaschii, Chloroflexus aurantiacus, Chlorobium vibrioforme and Cryptosporidium parvum [lactate dehydrogenase-like] MDHs with root-mean-square-deviation values between 1.4 and 2.1 A. Consistent with earlier reports, Ala at position 53 was responsible for coenzyme specificity, and the next residue, Arg, was important for NADP binding. Structural comparison revealed that the hyperthermostability of the A. pernix MDH is likely attributable to its smaller cavity volume and larger numbers of ion pairs and ion-pair networks, but the molecular strategy for thermostability may be specific for each enzyme.

Three-dimensional structure and enzymatic function of proapoptotic human p53-inducible quinone oxidoreductase PIG3

Tumor suppressor p53 regulates the expression of p53-induced genes (PIG) that trigger apoptosis. PIG3 or TP53I3 is the only known member of the medium chain dehydrogenase/reductase superfamily induced by p53 and is used as a proapoptotic marker. Although the participation of PIG3 in the apoptotic pathway is proven, the protein and its mechanism of action were never characterized. We analyzed human PIG3 enzymatic function and found NADPH-dependent reductase activity with ortho-quinones, which is consistent with the classification of PIG3 in the quinone oxidoreductase family. However, the activity is much lower than that of zeta-crystallin, a better known quinone oxidoreductase. In addition, we report the crystallographic structure of PIG3, which allowed the identification of substrate- and cofactor-binding sites, with residues fully conserved from bacteria to human. Tyr-59 in zeta-crystallin (Tyr-51 in PIG3) was suggested to participate in the catalysis of quinone reduction. However, kinetics of Tyr/Phe and Tyr/Ala mutants of both enzymes demonstrated that the active site Tyr is not catalytic but may participate in substrate binding, consistent with a mechanism based on propinquity effects. It has been proposed that PIG3 contribution to apoptosis would be through oxidative stress generation. We found that in vitro activity and in vivo overexpression of PIG3 accumulate reactive oxygen species. Accordingly, an inactive PIG3 mutant (S151V) did not produce reactive oxygen species in cells, indicating that enzymatically active protein is necessary for this function. This supports that PIG3 action is through oxidative stress produced by its enzymatic activity and provides essential knowledge for eventual control of apoptosis.

Predicting small ligand binding sites in proteins using backbone structure

Motivation: Specific non-covalent binding of metal ions and ligands, such as nucleotides and cofactors, is essential for the function of many proteins. Computational methods are useful for predicting the location of such binding sites when experimental information is lacking. Methods that use structural information, when available, are particularly promising since they can potentially identify non-contiguous binding motifs that cannot be found using only the amino acid sequence. Furthermore, a prediction method that can utilize low-resolution models is advantageous because high-resolution structures are available for only a relatively small fraction of proteins.

Results: SitePredict is a machine learning-based method for predicting binding sites in protein structures for specific metal ions or small molecules. The method uses Random Forest classifiers trained on diverse residue-based site properties including spatial clustering of residue types and evolutionary conservation. SitePredict was tested by cross-validation on a set of known binding sites for six different metal ions and five different small molecules in a non-redundant set of protein–ligand complex structures. The prediction performance was good for all ligands considered, as reflected by AUC values of at least 0.8. Furthermore, a more realistic test on unbound structures showed only a slight decrease in the accuracy. The properties that contribute the most to the prediction accuracy of each ligand were also examined. Finally, examples of predicted binding sites in homology models and uncharacterized proteins are discussed.

Availability: Binding site prediction results for all PDB protein structures and human protein homology models are available at http://sitepredict.org/.

Contact:bordner.andrew@mayo.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Structure and mechanism of GumK, a membrane-associated glucuronosyltransferase

Xanthomonas campestris GumK (beta-1,2-glucuronosyltransferase) is a 44-kDa membrane-associated protein that is involved in the biosynthesis of xanthan, an exopolysaccharide crucial for this bacterium's phytopathogenicity. Xanthan also has many important industrial applications. The GumK enzyme is the founding member of the glycosyltransferase family 70 of carbohydrate-active enzymes, which is composed of bacterial glycosyltransferases involved in exopolysaccharide synthesis. No x-ray structures have been reported for this family. To better understand the mechanism of action of the bacterial glycosyltransferases in this family, the x-ray crystal structure of apo-GumK was solved at 1.9 angstroms resolution. The enzyme has two well defined Rossmann domains with a catalytic cleft between them, which is a typical feature of the glycosyltransferase B superfamily. Additionally, the crystal structure of GumK complexed with UDP was solved at 2.28 angstroms resolution. We identified a number of catalytically important residues, including Asp157, which serves as the general base in the transfer reaction. Residues Met231, Met273, Glu272, Tyr292, Met306, Lys307, and Gln310 interact with UDP, and mutation of these residues affected protein activity both in vitro and in vivo. The biological and structural data reported here shed light on the molecular basis for donor and acceptor selectivity in this glycosyltransferase family. These results also provide a rationale to obtain new polysaccharides by varying residues in the conserved alpha/beta/alpha structural motif of GumK.

Combinatorial model for sequence and spatial motif discovery in short sequence fragments: examples from beta-barrel membrane proteins

Motifs are over-represented sequence or spatial patterns appearing in proteins. They often play important roles in maintaining protein stability and in facilitating protein functions. When motifs are located in short sequence fragments, as in transmembrane domains that are only 10-20 residues in length, and when there is only very limited data, it is difficult to identify motifs. In this study, we develop combinatorial models for assessing statistically significant sequence and spatial patterns. We show our method can uncover previously unknown sequence and spatial motifs in beta-barrel membrane proteins.

Role of the amino acid invariants in the active site of MurG as evaluated by site-directed mutagenesis

To evaluate their role in the active site of the MurG enzyme from Escherichia coli, 13 residues conserved in the sequences of 73 MurG orthologues were submitted to site-directed mutagenesis. All these residues lay within, or close to, the active site of MurG as defined by its tridimensional structure [Ha et al., Prot. Sci. 9 (2000) 1045-1052, and Hu et al., Proc. Natl. Acad. Sci. USA 100 (2003) 845-849]. Thirteen mutants proteins, in which residues T15, H18, Y105, H124, E125, N127, N134, S191, N198, R260, E268, Q288 or N291 have been replaced by alanine, were obtained as the C-terminal His-tagged forms. The effects of the mutations on the activity were checked: (i) by functional complementation of an E. coli murG mutant strain by the mutated genes; and (ii) by the determination of the steady-state kinetic parameters of the purified proteins. Most mutations resulted in an important loss of activity and, in the case of N134A, in the production of a highly unstable protein. The results correlated with the assigned or putative functions of the residues based on the tridimensional structure.

Probing the determinants of coenzyme specificity in Peptostreptococcus asaccharolyticus glutamate dehydrogenase by site-directed mutagenesis

Glutamate dehydrogenase (EC 1.4.1.2-4) from Peptostreptococcus asaccharolyticus has a strong preference for NADH over NADPH as a coenzyme, over 1000-fold in terms of kcat/Km values. Sequence alignments across the wider family of NAD(P)-dependent dehydrogenases might suggest that this preference is mainly due to a negatively charged glutamate at position 243 (E243) in the adenine ribose-binding pocket. We have examined the possibility of altering coenzyme specificity of the Peptostreptococcus enzyme, and, more specifically, the role of residue 243 and neighbouring residues in coenzyme binding, by introducing a range of point mutations. Glutamate dehydrogenases are unusual among dehydrogenases in that NADPH-specific forms usually have aspartate at this position. However, replacement of E243 with aspartate led to only a nine-fold relaxation of the strong discrimination against NADPH. By contrast, replacement with a more positively charged lysine or arginine, as found in NADPH-dependent members of other dehydrogenase families, allows a more than 1000-fold shift toward NADPH, resulting in enzymes equally efficient with NADH or NADPH. Smaller shifts in the same direction were also observed in enzymes where a neighboring tryptophan, W244, was replaced by a smaller alanine (approximately six-fold) or Asp245 was changed to lysine (32-fold). Coenzyme binding studies confirm that the mutations result in the expected major changes in relative affinities for NADH and NADPH, and pH studies indicate that improved affinity for the extra phosphate of NADPH is the predominant reason for the increased catalytic efficiency with this coenzyme. The marked difference between the results of replacing E243 with aspartate and with positive residues implies that the mode of NADPH binding in naturally occurring NADPH-dependent glutamate dehydrogenases differs from that adopted in E243K or E243D and in other dehydrogenases.

Crystal structure of archaeal highly thermostable L-aspartate dehydrogenase/NAD/citrate ternary complex

The crystal structure of the highly thermostable L-aspartate dehydrogenase (L-aspDH; EC 1.4.1.21) from the hyperthermophilic archaeon Archaeoglobus fulgidus was determined in the presence of NAD and a substrate analog, citrate. The dimeric structure of A. fulgidus L-aspDH was refined at a resolution of 1.9 A with a crystallographic R-factor of 21.7% (R(free) = 22.6%). The structure indicates that each subunit consists of two domains separated by a deep cleft containing an active site. Structural comparison of the A. fulgidus L-aspDH/NAD/citrate ternary complex and the Thermotoga maritima L-aspDH/NAD binary complex showed that A. fulgidus L-aspDH assumes a closed conformation and that a large movement of the two loops takes place during substrate binding. Like T. maritima L-aspDH, the A. fulgidus enzyme is highly thermostable. But whereas a large number of inter- and intrasubunit ion pairs are responsible for the stability of A. fulgidus L-aspDH, a large number of inter- and intrasubunit aromatic pairs stabilize the T. maritima enzyme. Thus stabilization of these two L-aspDHs appears to be achieved in different ways. This is the first detailed description of substrate and coenzyme binding to L-aspDH and of the molecular basis of the high thermostability of a hyperthermophilic L-aspDH.

Crystal structure and biochemical properties of the D-arabinose dehydrogenase from Sulfolobus solfataricus

Sulfolobus solfataricus metabolizes the five-carbon sugar d-arabinose to 2-oxoglutarate by an inducible pathway consisting of dehydrogenases and dehydratases. Here we report the crystal structure and biochemical properties of the first enzyme of this pathway: the d-arabinose dehydrogenase. The AraDH structure was solved to a resolution of 1.80 A by single-wavelength anomalous diffraction and phased using the two endogenous zinc ions per subunit. The structure revealed a catalytic and cofactor binding domain, typically present in mesophilic and thermophilic alcohol dehydrogenases. Cofactor modeling showed the presence of a phosphate binding pocket sequence motif (SRS-X2-H), which is likely to be responsible for the enzyme's preference for NADP+. The homo-tetrameric enzyme is specific for d-arabinose, l-fucose, l-galactose and d-ribose, which could be explained by the hydrogen bonding patterns of the C3 and C4 hydroxyl groups observed in substrate docking simulations. The enzyme optimally converts sugars at pH 8.2 and 91 degrees C, and displays a half-life of 42 and 26 min at 85 and 90 degrees C, respectively, indicating that the enzyme is thermostable at physiological operating temperatures of 80 degrees C. The structure represents the first crystal structure of an NADP+-dependent member of the medium-chain dehydrogenase/reductase (MDR) superfamily from Archaea.

A computational analysis of the three isoforms of glutamate dehydrogenase reveals structural features of the isoform EC 1.4.1.4 supporting a key role in ammonium assimilation by plants

BACKGROUND

There are three isoforms of glutamate dehydrogenase. The isoform EC 1.4.1.4 (GDH4) catalyses glutamate synthesis from 2-oxoglutarate and ammonium, using NAD(P)H. Ammonium assimilation is critical for plant growth. Although GDH4 from animals and prokaryotes are well characterized, there are few data concerning plant GDH4, even from those whose genomes are well annotated.

RESULTS

A large set of the three GDH isoforms was built resulting in 116 non-redundant full polypeptide sequences. A computational analysis was made to gain more information concerning the structure-function relationship of GDH4 from plants (Eukaryota, Viridiplantae). The tested plant GDH4 sequences were the two ones known to date, those of Chlorella sorokiniana. This analysis revealed several structural features specific of plant GDH4: (i) the lack of a structure called "antenna"; (ii) the NAD(P)-binding motif GAGNVA; and (iii) a second putative coenzyme-binding motif GVLTGKG together with four residues involved in the binding of the reduced form of NADP.

CONCLUSION

A number of structural features specific of plant GDH4 have been found. The results reinforce the probable key role of GDH4 in ammonium assimilation by plants.

REVIEWERS

This article was reviewed by Tina Bakolitsa (nominated by Eugene Koonin), Martin Jambon (nominated by Laura Landweber), Sandor Pangor and Franck Eisenhaber.

Molecular modeling and site-directed mutagenesis of plant chloroplast monogalactosyldiacylglycerol synthase reveal critical residues for activity

Monogalactosyldiacylglycerol (MGDG), the major lipid of plant and algal plastids, is synthesized by MGD (or MGDG synthase), a dimeric and membrane-bound glycosyltransferase of the plastid envelope that catalyzes the transfer of a galactosyl group from a UDP-galactose donor onto a diacylglycerol acceptor. Although this enzyme is essential for biogenesis, and therefore an interesting target for herbicide design, no structural information is available. MGD monomers share sequence similarity with MURG, a bacterial glycosyltransferase catalyzing the transfer of N-acetyl-glucosamine on Lipid 1. Using the x-ray structure of Escherichia coli MURG as a template, we computed a model for the fold of Spinacia oleracea MGD. This structural prediction was supported by site-directed mutagenesis analyses. The predicted monomer architecture is a double Rossmann fold. The binding site for UDP-galactose was predicted in the cleft separating the two Rossmann folds. Two short segments of MGD (beta2-alpha2 and beta6-beta7 loops) have no counterparts in MURG, and their structure could not be determined. Combining the obtained model with phylogenetic and biochemical information, we collected evidence supporting the beta2-alpha2 loop in the N-domain as likely to be involved in diacylglycerol binding. Additionally, the monotopic insertion of MGD in one membrane leaflet of the plastid envelope occurs very likely at the level of hydrophobic amino acids of the N-terminal domain.

Reactive amino acid residues involved in glutamate-binding of human glutamate dehydrogenase isozymes

In the present study, the cassette mutagenesis at several putative positions (K94, G96, K118, K130, or D172) was performed to examine the residues involved in the glutamate-binding of the human glutamate dehydrogenase isozymes (hGDH1 and hGDH2). None of the mutations tested affected the expression or stability of the proteins. There was dramatic reduction in the catalytic efficiency in mutant proteins at K94, G96, K118, or K130 site, but not at D172 site. The K(M) values for glutamate were 4-10-fold greater for the mutants at K94, G96, or K118 site than for the wild-type hGDH1 and hGDH2, whereas no differences in the K(M) values for NAD(+) were detected between the mutant and wild-type enzymes. For K130Y mutant, the K(M) value for glutamate increased 1.6-fold, whereas the catalytic efficiency (k(cat)/K(M)) showed only 2-3% of the wild-type. Therefore, the decreased catalytic efficiency of the K130 mutant mainly results from the reduced k(cat) value, suggesting a possibility that the K130Y residue may be involved in the catalysis rather than in the glutamate-binding. The D172Y mutant did not show any changes in k(cat) value and K(M) values for glutamate and NAD(+), indicating that D172Y is not directly involved in catalysis and substrates binding of the hGDH isozymes. For sensitivity to ADP activation, only the D172Y mutant showed a reduced sensitivity to ADP activation. The reduction of ADP activation in D172Y mutant was more profoundly observed in hGDH2 than in hGDH1. There were no differences in their sensitivities to GTP inhibition between the wild-type and mutant GDHs at all positions tested. Our results suggest that K94, G96, and K118 residues play an important role, although at different degrees, in the binding of glutamate to hGDH isozymes.

The first crystal structure of hyperthermostable NAD-dependent glutamate dehydrogenase from Pyrobaculum islandicum

The extremely thermostable NAD-dependent glutamate dehydrogenase (NAD-GluDH) from Pyrobaculum islandicum, a member of the Crenarchaeota, was crystallized, and its 3D structure has been determined by X-ray diffraction methods. The homohexameric structure of Pb. islandicum glutamate dehydrogenase (Pis-GluDH) was solved and refined at a resolution of 2.9A with a crystallographic R-factor of 19.9% (Rfree 26.0%). The structure indicates that each subunit consists of two domains separated by a deep cleft containing an active site. The secondary structural elements and catalytically important residues of the enzyme were highly conserved among the NAD(P)-dependent GluDHs from other sources. A structural comparison of Pis-GluDH with other NAD(P)-dependent GluDHs suggests that a significant difference in the alpha8-loop-alpha9 region of this enzyme is associated with its coenzyme specificity. From the analysis of the 3D structure, hydrophobic interactions between intersubunits were found to be important features for the enzyme oligomerization. It has been reported that Pis-GluDH is highly thermostable, like the GluDH of the hyperthermophilic archaeum Pyrococcus furiosus, and the increase in the intersubunit ion pair networks is responsible for the extreme thermostability of the Pc. furiosus enzyme. However, the number of intersubunit ion pairs in the Pis-GluDH molecules is much smaller than those of the Pc. furiosus GluDH. The number of hydrophobic interactions at the intersubunit interfaces were increased and responsible for the extremely high thermostability. This indicates that the major molecular strategy for high thermostability of the GluDHs may be different for each hyperthermophile.

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

Methanoarchaeal sulfolactate dehydrogenase: prototype of a new family of NADH-dependent enzymes

The crystal structure of the sulfolactate dehydrogenase from the hyperthermophilic and methanogenic archaeon Methanocaldococcus jannaschii was solved at 2.5 A resolution (PDB id. 1RFM). The asymmetric unit contains a tetramer of tight dimers. This structure, complexed with NADH, does not contain a cofactor-binding domain with 'Rossmann-fold' topology. Instead, the tertiary and quaternary structures indicate a novel fold. The NADH is bound in an extended conformation in each active site, in a manner that explains the pro-S specificity. Cofactor binding involves residues belonging to both subunits within the tight dimers, which are therefore the smallest enzymatically active units. The protein was found to be a homodimer in solution by size-exclusion chromatography, analytical ultracentrifugation and small-angle neutron scattering. Various compounds were tested as putative substrates. The results indicate the existence of a substrate discrimination mechanism, which involves electrostatic interactions. Based on sequence homology and phylogenetic analyses, several other enzymes were classified as belonging to this novel family of homologous (S)-2-hydroxyacid dehydrogenases.