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

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

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

Metabolic engineering of Pichia pastoris for production of isopentanol (3-Methyl-1-butanol)

In recent years, the increasingly serious and clear effects of climate change have increased interest in renewable fuels and platform chemicals. Microbial platforms that can produce these compounds in an economically efficient way have emerged as an attractive alternative to the traditional production approaches. Here, we engineered the industrially-relevant yeast Pichia pastoris to produce the platform chemical 3-methyl-1-butanol (3M1B, isopentanol) directly from the renewable carbon source glucose. Specifically, we overexpressed the endogenous valine and leucine biosynthetic pathways to increase the production of the key pathway intermediate, 2-ketoisocaproate (2-KIC). Overexpression of the artificial keto-acid degradation pathway converted 2-KIC into 3M1B. Down-regulation of the side-product ethanol production using the CRISPR/Cas9 system led to a strain that is able to produce 3M1B at a titer of 191.0 ± 9.6 mg/L, the highest titer reported to date in a non-conventional yeast. We envision that our yeast system will pave the way for an efficient production system for this important class of platform compounds.

A Novel Type II NAD+-Specific Isocitrate Dehydrogenase from the Marine Bacterium Congregibacter litoralis KT71

In most living organisms, isocitrate dehydrogenases (IDHs) convert isocitrate into ɑ-ketoglutarate (ɑ-KG). Phylogenetic analyses divide the IDH protein family into two subgroups: types I and II. Based on cofactor usage, IDHs are either NAD⁺-specific (NAD-IDH) or NADP⁺-specific (NADP-IDH); NADP-IDH evolved from NAD-IDH. Type I IDHs include NAD-IDHs and NADP-IDHs; however, no type II NAD-IDHs have been reported to date. This study reports a novel type II NAD-IDH from the marine bacterium Congregibacter litoralis KT71 (ClIDH, GenBank accession no. EAQ96042). His-tagged recombinant ClIDH was produced in Escherichia coli and purified; the recombinant enzyme was NAD⁺-specific and showed no detectable activity with NADP⁺. The K_m values of the enzyme for NAD⁺ were 262.6±7.4 μM or 309.1±11.2 μM with Mg²⁺ or Mn²⁺ as the divalent cation, respectively. The coenzyme specificity of a ClIDH Asp487Arg/Leu488His mutant was altered, and the preference of the mutant for NADP⁺ was approximately 24-fold higher than that for NAD⁺, suggesting that ClIDH is an NAD⁺-specific ancestral enzyme in the type II IDH subgroup. Gel filtration and analytical ultracentrifugation analyses revealed the homohexameric structure of ClIDH, which is the first IDH hexamer discovered thus far. A 163-amino acid segment of CIIDH is essential to maintain its polymerization structure and activity, as a truncated version lacking this region forms a non-functional monomer. ClIDH was dependent on divalent cations, the most effective being Mn²⁺. The maximal activity of purified recombinant ClIDH was achieved at 35°C and pH 7.5, and a heat inactivation experiment showed that a 20-min incubation at 33°C caused a 50% loss of ClIDH activity. The discovery of a NAD⁺-specific, type II IDH fills a gap in the current classification of IDHs, and sheds light on the evolution of type II IDHs.

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

Background

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

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

Results

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

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

Conclusion

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

WRKY gene family evolution in Arabidopsis thaliana

The Arabidopsis thaliana WRKY proteins are characterized by a sequence of 60 amino acids including WRKY domain. It is well established that these proteins are involved in the regulation of various physiological programs unique to plants including pathogen defense, senescence and response to environmental stresses, which attracts attention of the scientific community as to how this family might have evolved. We tried to satisfy this curiosity and analyze reasons for duplications of these gene sequences leading to their diversified gene actions. The WRKY sequences available in Arabidopsis thaliana were used to evaluate selection pressure following duplication events. A phylogenetic tree was constructed and the WRKY family was divided into five sub-families. After that, tests were conducted to decide whether positive or purified selection played key role in these events. Our results suggest that purifying selection played major role during the evolution of this family. Some amino acid changes were also detected in specific branches of phylogeny suggesting that relaxed constraints might also have contributed to functional divergence among sub-families. Sites relaxed from purifying selection were identified and mapped onto the structural and functional regions of the WRKY1 protein. These analyses will enhance our understanding of the precise role played by natural selection to create functional diversity in WRKY family.

Evolutionary selection pressure of forkhead domain and functional divergence

Forkhead-box (Fox) genes encode a family of transcription factors defined by a "winged helix" DNA-binding domain which have been identified in many metazoans, and play important roles in diverse biological processes. Here we aimed to extend previous evolutionary selection analysis to fungi, using available sequences from E. cuniculi (Ec), Eremothecium gossypii (Eg), Saccharomyces cerevisiae (Sc), etc. The phylogeny of 335 Fox protein sequences was reconstructed, revealing the existence of 26 orthologous groups that were well supported by gene phylogeny which arose following a series of gene duplication events. Gene conversion events may also play important roles in the evolution of Fox genes. The nonsynonymous to synonymous substitution ratios (dN/dS) for orthologous groups suggested that after gene duplication and/or speciation of forkhead clusters, rapid differentiation and the negative selection have occurred, prompting the formation of distinct Fox subclasses and new functions. SDPpred was used to produce a set of the alignment positions (specificity determining positions) which is involved in conferring differential functional specificity. These findings explained the functional divergence of Fox gene family.