Evaluation by mutagenesis of the importance of 3 arginines in alpha, beta, and gamma subunits of human NAD-dependent isocitrate dehydrogenase

Two Different Isocitrate Dehydrogenases from Pseudomonas aeruginosa: Enzymology and Coenzyme-Evolutionary Implications

Pseudomonas aeruginosa PAO1, as an experimental model for Gram-negative bacteria, harbors two NADP+-dependent isocitrate dehydrogenases (NADP-IDHs) that were evolved from its ancient counterpart NAD-IDHs. For a better understanding of PaIDH1 and PaIDH2, we cloned the genes, overexpressed them in Escherichia coli and purified them to homogeneity. PaIDH1 displayed higher affinity to NADP+ and isocitrate, with lower Km values when compared to PaIDH2. Moreover, PaIDH1 possessed higher temperature tolerance (50 °C) and wider pH range tolerance (7.2-8.5) and could be phosphorylated. After treatment with the bifunctional PaIDH kinase/phosphatase (PaIDH K/P), PaIDH1 lost 80% of its enzymatic activity in one hour due to the phosphorylation of Ser115. Small-molecule compounds like glyoxylic acid and oxaloacetate can effectively inhibit the activity of PaIDHs. The mutant PaIDH1-D346I347A353K393 exhibited enhanced affinity for NAD+ while it lost activity towards NADP+, and the Km value (7770.67 μM) of the mutant PaIDH2-L589 I600 for NADP+ was higher than that observed for NAD+ (5824.33 μM), indicating a shift in coenzyme specificity from NADP+ to NAD+ for both PaIDHs. The experiments demonstrated that the mutation did not alter the oligomeric state of either protein. This study provides a foundation for the elucidation of the evolution and function of two NADP-IDHs in the pathogenic bacterium P. aeruginosa.

Molecular insights into the catalysis and regulation of mammalian NAD-dependent isocitrate dehydrogenases

Eukaryotic NAD-dependent isocitrate dehydrogenases (NAD-IDHs) are mitochondria-localized enzymes which catalyze the oxidative decarboxylation of isocitrate to α-ketoglutarate using NAD as a cofactor. In mammals, NAD-IDHs (or IDH3) consist of three types of subunits (α, β, and γ), and exist as (α2βγ)2 heterooctamer. Mammalian NAD-IDHs are regulated allosterically and/or competitively by a diversity of metabolites including citrate, ADP, ATP, NADH, and NADPH, which are associated with cellular metabolite flux, energy demands, and redox status. Proper assembly of the component subunits is essential for the catalysis and regulation of the enzymes. Recently, crystal structures of human IDH3 have been solved in apo form and in complex with various ligands, revealing the molecular mechanisms for the assembly, catalysis, and regulation of the enzyme.

Targeting the TCA cycle can ameliorate widespread axonal energy deficiency in neuroinflammatory lesions

Inflammation in the central nervous system can impair the function of neuronal mitochondria and contributes to axon degeneration in the common neuroinflammatory disease multiple sclerosis (MS). Here we combine cell-type-specific mitochondrial proteomics with in vivo biosensor imaging to dissect how inflammation alters the molecular composition and functional capacity of neuronal mitochondria. We show that neuroinflammatory lesions in the mouse spinal cord cause widespread and persisting axonal ATP deficiency, which precedes mitochondrial oxidation and calcium overload. This axonal energy deficiency is associated with impaired electron transport chain function, but also an upstream imbalance of tricarboxylic acid (TCA) cycle enzymes, with several, including key rate-limiting, enzymes being depleted in neuronal mitochondria in experimental models and in MS lesions. Notably, viral overexpression of individual TCA enzymes can ameliorate the axonal energy deficits in neuroinflammatory lesions, suggesting that TCA cycle dysfunction in MS may be amendable to therapy.

Structures of a constitutively active mutant of human IDH3 reveal new insights into the mechanisms of allosteric activation and the catalytic reaction

Human NAD-dependent isocitrate dehydrogenase or IDH3 (HsIDH3) catalyzes the decarboxylation of isocitrate into α-ketoglutarate in the tricarboxylic acid cycle. It consists of three types of subunits (α, β, and γ) and exists and functions as the (αβαγ)₂ heterooctamer. HsIDH3 is regulated allosterically and/or competitively by numerous metabolites including CIT, ADP, ATP, and NADH. Our previous studies have revealed the molecular basis for the activity and regulation of the αβ and αγ heterodimers. However, the molecular mechanism for the allosteric activation of the HsIDH3 holoenzyme remains elusive. In this work, we report the crystal structures of the αβ and αγ heterodimers and the (αβαγ)₂ heterooctamer containing an α-Q139A mutation in the clasp domain, which renders all the heterodimers and the heterooctamer constitutively active in the absence of activators. Our structural analysis shows that the α-Q139A mutation alters the hydrogen-bonding network at the heterodimer-heterodimer interface in a manner similar to that in the activator-bound αγ heterodimer. This alteration not only stabilizes the active sites of both α_Q139Aβ and α_Q139Aγ heterodimers in active conformations but also induces conformational changes of the pseudo-allosteric site of the α_Q139Aβ heterodimer enabling it to bind activators. In addition, the α_Q139A^ICT+Ca+NADβ^NAD structure presents the first pseudo-Michaelis complex of HsIDH3, which allows us to identify the key residues involved in the binding of cofactor, substrate, and metal ion. Our structural and biochemical data together reveal new insights into the molecular mechanisms for allosteric regulation and the catalytic reaction of HsIDH3.

Structure and allosteric regulation of human NAD-dependent isocitrate dehydrogenase

Human NAD-dependent isocitrate dehydrogenase or HsIDH3 catalyzes the decarboxylation of isocitrate into α-ketoglutarate in the TCA cycle. HsIDH3 exists and functions as a heterooctamer composed of the αβ and αγ heterodimers, and is regulated allosterically and/or competitively by numerous metabolites including CIT, ADP, ATP, and NADH. In this work, we report the crystal structure of HsIDH3 containing a β mutant in apo form. In the HsIDH3 structure, the αβ and αγ heterodimers form the α₂βγ heterotetramer via their clasp domains, and two α₂βγ heterotetramers form the (α₂βγ)₂ heterooctamer through insertion of the N-terminus of the γ subunit of one heterotetramer into the back cleft of the β subunit of the other heterotetramer. The functional roles of the key residues at the allosteric site, the pseudo allosteric site, the heterodimer and heterodimer-heterodimer interfaces, and the N-terminal of the γ subunit are validated by mutagenesis and kinetic studies. Our structural and biochemical data together demonstrate that the allosteric site plays an important role but the pseudo allosteric site plays no role in the allosteric activation of the enzyme; the activation signal from the allosteric site is transmitted to the active sites of both αβ and αγ heterodimers via the clasp domains; and the N-terminal of the γ subunit plays a critical role in the formation of the heterooctamer to ensure the optimal activity of the enzyme. These findings reveal the molecular mechanism of the assembly and allosteric regulation of HsIDH3.

Comparative Proteomic Analysis of Rana chensinensis Oviduct

As one of most important traditional Chinese medicine resources, the oviduct of female Rana chensinensis (Chinese brown frog) was widely used in the treatment of asthenia after sickness or delivery, deficiency in vigor, palpitation, and insomnia. Unlike other vertebrates, the oviduct of Rana chensinensis oviduct significantly expands during prehibernation, in contrast to the breeding period. To explain this phenomenon at the molecular level, the protein expression profiles of Rana chensinensis oviduct during the breeding period and prehibernation were observed using isobaric tags for relative and absolute quantitation (iTRAQ) technique. Then, all identified proteins were used to obtain gene ontology (GO) annotation. Ultimately, KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis was performed to predict the pathway on differentially expressed proteins (DEPs). A total of 4479 proteins were identified, and 312 of them presented different expression profiling between prehibernation and breeding period. Compared with prehibernation group, 86 proteins were upregulated, and 226 proteins were downregulated in breeding period. After KEGG enrichment analysis, 163 DEPs were involved in 6 pathways, which were lysosome, RNA transport, glycosaminoglycan degradation, extracellular matrix (ECM)–receptor interaction, metabolic pathways and focal adhesion. This is the first report on the protein profiling of Rana chensinensis oviduct during the breeding period and prehibernation. Results show that this distinctive physiological phenomenon of Rana chensinensis oviduct was mainly involved in ECM–receptor interaction, metabolic pathways, and focal adhesion.

Insights into the inhibitory mechanisms of NADH on the αγ heterodimer of human NAD-dependent isocitrate dehydrogenase

Human NAD-dependent isocitrate dehydrogenase (NAD-IDH) catalyzes the oxidative decarboxylation of isocitrate in the citric acid cycle. In the α₂βγ heterotetramer of NAD-IDH, the γ subunit plays the regulatory role and the β subunit the structural role. Previous biochemical data have shown that mammalian NAD-IDHs can be inhibited by NADH; however, the molecular mechanism is unclear. In this work, we show that the αβ, αγ and α₂βγ enzymes of human NAD-IDH can be inhibited by NADH, and further determine the crystal structure of the αγ heterodimer bound with an Mg²⁺ and an NADH at the active site and an NADH at the allosteric site, which resembles that of the inactive α^Mgγ heterodimer. The NADH at the active site occupies the binding site for NAD⁺ and prevents the binding of the cofactor. The NADH at the allosteric site occupies the binding sites for ADP and citrate and blocks the binding of the activators. The biochemical data confirm that the NADH binding competes with the binding of NAD⁺ and the binding of citrate and ADP, and the two effects together contribute to the NADH inhibition on the activity. These findings provide insights into the inhibitory mechanisms of the αγ heterodimer by NADH.

The β and γ subunits play distinct functional roles in the α2βγ heterotetramer of human NAD-dependent isocitrate dehydrogenase

Human NAD-dependent isocitrate dehydrogenase existing as the α₂βγ heterotetramer, catalyzes the decarboxylation of isocitrate into α-ketoglutarate in the Krebs cycle, and is allosterically regulated by citrate, ADP and ATP. To explore the functional roles of the regulatory β and γ subunits, we systematically characterized the enzymatic properties of the holoenzyme and the composing αβ and αγ heterodimers in the absence and presence of regulators. The biochemical and mutagenesis data show that αβ and αγ alone have considerable basal activity but the full activity of α₂βγ requires the assembly and cooperative function of both heterodimers. α₂βγ and αγ can be activated by citrate or/and ADP, whereas αβ cannot. The binding of citrate or/and ADP decreases the S_{0.5,isocitrate} and thus enhances the catalytic efficiencies of the enzymes, and the two activators can act independently or synergistically. Moreover, ATP can activate α₂βγ and αγ at low concentration and inhibit the enzymes at high concentration, but has only inhibitory effect on αβ. Furthermore, the allosteric activation of α₂βγ is through the γ subunit not the β subunit. These results demonstrate that the γ subunit plays regulatory role to activate the holoenzyme, and the β subunit the structural role to facilitate the assembly of the holoenzyme.

Molecular mechanism of the allosteric regulation of the αγ heterodimer of human NAD-dependent isocitrate dehydrogenase

Human NAD-dependent isocitrate dehydrogenase catalyzes the decarboxylation of isocitrate (ICT) into α-ketoglutarate in the Krebs cycle. It exists as the α2βγ heterotetramer composed of the αβ and αγ heterodimers. Previously, we have demonstrated biochemically that the α2βγ heterotetramer and αγ heterodimer can be allosterically activated by citrate (CIT) and ADP. In this work, we report the crystal structures of the αγ heterodimer with the γ subunit bound without or with different activators. Structural analyses show that CIT, ADP and Mg2+ bind adjacent to each other at the allosteric site. The CIT binding induces conformational changes at the allosteric site, which are transmitted to the active site through the heterodimer interface, leading to stabilization of the ICT binding at the active site and thus activation of the enzyme. The ADP binding induces no further conformational changes but enhances the CIT binding through Mg2+-mediated interactions, yielding a synergistic activation effect. ICT can also bind to the CIT-binding subsite, which induces similar conformational changes but exhibits a weaker activation effect. The functional roles of the key residues are verified by mutagenesis, kinetic and structural studies. Our structural and functional data together reveal the molecular mechanism of the allosteric regulation of the αγ heterodimer.

Genome wide identification of aberrant alternative splicing events in myotonic dystrophy type 2

Myotonic dystrophy type 2 (DM2) is a genetic, autosomal dominant disease due to expansion of tetraplet (CCTG) repetitions in the first intron of the ZNF9/CNBP gene. DM2 is a multisystemic disorder affecting the skeletal muscle, the heart, the eye and the endocrine system. According to the proposed pathological mechanism, the expanded tetraplets have an RNA toxic effect, disrupting the splicing of many mRNAs. Thus, the identification of aberrantly spliced transcripts is instrumental for our understanding of the molecular mechanisms underpinning the disease. The aim of this study was the identification of new aberrant alternative splicing events in DM2 patients. By genome wide analysis of 10 DM2 patients and 10 controls (CTR), we identified 273 alternative spliced exons in 218 genes. While many aberrant splicing events were already identified in the past, most were new. A subset of these events was validated by qPCR assays in 19 DM2 and 15 CTR subjects. To gain insight into the molecular pathways involving the identified aberrantly spliced genes, we performed a bioinformatics analysis with Ingenuity system. This analysis indicated a deregulation of development, cell survival, metabolism, calcium signaling and contractility. In conclusion, our genome wide analysis provided a database of aberrant splicing events in the skeletal muscle of DM2 patients. The affected genes are involved in numerous pathways and networks important for muscle physio-pathology, suggesting that the identified variants may contribute to DM2 pathogenesis.

The role of NAD(+)-dependent isocitrate dehydrogenase 3 subunit α in AFB1 induced liver lesion

Aflatoxin B1 (AFB1) is a potent hepatocarcinogen that causes carcinogenesis in many animal species. In previous study, we found that isocitrate dehydrogenase 3α subunit (IDH3α) was upregulated in AFB1-induced carcinogenesis process. In this study, the sequences of IDH3α from various species were compared and the protein expression levels in different organs were examined, and the results showed that IDH3α was a widely distributed protein and shared highly conserved sequence in various species. In the same time, IDH3α was demonstrated to accumulate in a dose-dependent manner induced by AFB1 in cells, and was also up-regulated in the process of AFB1-induced liver lesion. Similar results were observed when H2O2 was used to replace AFB1. Over-expression of IDH3α increased the phosphorylation level of Akt (Protein kinase B) and neutralized the cellular toxicity induced by AFB1 or H2O2 and apoptosis induced by AFB1, while the reduced expression of IDH3α by siRNA decreased the phosphorylation, indicating that IDH3α played important roles in oxidative stress-induced PI3K/Akt pathway. Overall, the results suggested that AFB1 treatment could increase the expression of IDH3α, and the activated PI3K/Akt pathway by IDH3α eventually neutralized the apoptosis induced by AFB1.

Cloning, differential expression, and association analysis with fat traits of porcine IDH3γ gene

Mitochondrial NAD⁺-dependent isocitrate dehydrogenase (IDH3) catalyzes the allosterically regulated rate-limiting step of the tricarboxylic acid cycle activated. In pigs, very little is known about this gene. Here, we cloned 1,346 bp full-length cDNA and 8,778 bp genomic sequence of porcine γ subunit of IDH3 (IDH3γ). IDH3γ contains 12 exons separated by 11 introns. Real-time PCR revealed that IDH3γ mRNA were upregulated in backfat of Large White compared with Meishan and F1 hybrids, and most abundant in small intestine via tissue distribution profile. A microsatellite ("GT" repeats) in second intron was found. The selected pigs were genotyped at this microsatellite. The IDH3γ genotypes showed a significant effect on backfat thickness at thorax-waist (P < 0.05), backfat thickness at sixth to seventh thorax (P < 0.01), and average backfat thickness (P < 0.05). This site seemed to be significantly dominant in action (P < 0.05 for backfat thickness at sixth to seventh thorax, backfat thickness at thorax-waist, and average backfat thickness), and allele B was associated with increase of thickness values of these traits. This locus is possibly considered as a marker for adipose deposition traits.

Ligand binding and structural changes associated with allostery in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an octameric enzyme composed of four each of regulatory IDH1 and catalytic IDH2 subunits that share 42% sequence identity. IDH2 contains catalytic isocitrate/Mg2+ and NAD+ binding sites whereas IDH1 contains homologous binding sites, respectively, for cooperative binding of isocitrate and for allosteric binding of AMP. Ligand binding is highly ordered in vitro, and IDH exhibits the unusual property of half-site binding for all ligands. The structures of IDH solved in the absence or presence of ligands have shown: (a) a heterodimer to be the basic structural/functional unit of the enzyme, (b) the organization of heterodimers to form tetramer and octamer structures, (c) structural differences that may underlie cooperative and allosteric regulatory mechanisms, and (d) the possibility for formation of a disulfide bond that could reduce catalytic activity. In vivo analyses of mutant enzymes have elucidated the physiological importance of catalytic activity and allosteric regulation of this tricarboxylic acid cycle enzyme. Other studies have established the importance of a disulfide bond in regulation of IDH activity in vivo, as well as contributions of this bond to the property of half-site ligand binding exhibited by the wild-type enzyme.

Basis for half-site ligand binding in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase is an allosterically regulated octameric enzyme composed of four heterodimers of a catalytic IDH2 subunit and a regulatory IDH1 subunit. Despite structural predictions that the enzyme would contain eight isocitrate binding sites, four NAD(+) binding sites, and four AMP binding sites, only half of the sites for each ligand can be measured in binding assays. On the basis of a potential interaction between side chains of Cys-150 residues in IDH2 subunits in each tetramer of the enzyme, ligand binding assays of wild-type (IDH1/IDH2) and IDH1/IDH2(C150S) octameric enzymes were conducted in the presence of dithiothreitol. These assays demonstrated the presence of eight isocitrate and four AMP binding sites for the wild-type enzyme in the presence of dithiothreitol and for the IDH1/IDH2(C150S) enzyme in the absence or presence of this reagent, suggesting that interactions between sulfhydryl side chains of IDH2 Cys-150 residues limit access to these sites. However, only two NAD(+) sites could be measured for either enzyme. A tetrameric form of IDH (an IDH1(G15D)/IDH2 mutant enzyme) demonstrated half-site binding for isocitrate (two sites) in the absence of dithiothreitol and full-site binding (four sites) in the presence of dithiothreitol. Only one NAD(+) site could be measured for the tetramer under both conditions. In the context of the structure of the enzyme, these results suggest that an observed asymmetry between heterotetramers in the holoenzyme contributes to interactions between IDH2 Cys-150 residues and to half-site binding of isocitrate, but that a form of negative cooperativity may limit access to apparently equivalent NAD(+) binding sites.

A 304 bp insertion/deletion mutation in promoter region induces the increase of porcine IDH3β gene expression

Obese and lean pig breeds show obvious differences in adipose metabolism/fat deposition; however, the molecular mechanism underlying phenotype variation remains unknown. In order to understand it, we analyzed the differences of gene expression in backfat between Meishan (a typical Chinese indigenous obese breed) and Large White (a lean Western breed) pigs. Here, we cloned porcine β subunit of IDH3 (IDH3B) and 2447 bp 5'-flanking sequence of this gene, and determined the genomic structure. Porcine IDH3B contains three isoforms, IDH3B ( 1 ), IDH3B ( 2 ) and IDH3B ( 3 ). Real-time RT-PCR revealed that these three isoforms were prevalently up-regulated in backfat of western commercial pigs, Large White, Landrace and Duroc, compared with Chinese indigenous breeds, Meishan and Tongcheng pigs. A 304 bp insertion/deletion variant was found in the 5'-flanking region. Dual-luciferase reporter assays showed that in vitro the promoter of IDH3B gene with the insertion had higher luciferase activity as compared with the wild type. Three genotypes AA, AB and BB, due to this insertion, were detected, and the frequency of allele A was dominant in western commercial pigs, whereas allele B predominated in Chinese indigenous breeds. IDH3B mRNA expression in Meishan pigs was more abundant with genotype AA than with genotype AB or BB, as in Large White pigs. In addition, the polymorphism was detected in 317 pigs of a Large White × Meishan F2 resource population. Association analysis showed that pigs with genotype AA possessed higher backfat thickness at buttocks than those with genotype AB (P < 0.05) or BB. These data suggested that the 304 bp insertion mutation in promoter region increased the expression of porcine IDH3β transcripts and this mutation might be a candidate marker for marker assistant selection in swine breeding.

The tricarboxylic acid cycle in L₃ Teladorsagia circumcincta: metabolism of acetyl CoA to succinyl CoA

Nematodes, like other species, derive much of the energy for cellular processes from mitochondrial pathways including the TCA cycle. Previously, we have shown L₃ Teladorsagia circumcincta consume oxygen and so may utilise a full TCA cycle for aerobic energy metabolism. We have assessed the relative activity levels and substrate affinities of citrate synthase, aconitase, isocitrate dehydrogenase (both NAD+ and NADP+ specific) and α-ketoglutarate dehydrogenase in homogenates of L₃ T. circumcincta. All of these enzymes were present in homogenates. Compared with citrate synthase, low levels of enzyme activity and low catalytic efficiency was observed for NAD+ isocitrate dehydrogenase and especially α-ketoglutarate dehydrogenase. Therefore, it is likely that the activity of these to enzymes regulate overall metabolite flow through the TCA cycle, especially when [NAD+] limits enzyme activity. Of the enzymes tested, only citrate synthase had substrate affinities which were markedly different from values obtained from mammalian species. Overall, the results are consistent with the suggestion that a full TCA cycle exists withinL₃ T. circumcincta. While there may subtle variations in enzyme properties, particularly for citrate synthase, the control points for the TCA cycle inL₃ T. circumcincta are probably similar to those in the tissues of their host species.

Construction and analyses of tetrameric forms of yeast NAD+-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an octameric enzyme composed of four heterodimers of regulatory IDH1 and catalytic IDH2 subunits. The crystal structure suggested that the interactions between tetramers in the octamer are restricted to defined regions in IDH1 subunits from each tetramer. Using truncation and mutagenesis, we constructed three tetrameric forms of IDH. Truncation of five residues from the amino terminus of IDH1 did not alter the octameric form of the enzyme, but this truncation with an IDH1 G15D or IDH1 D168K residue substitution produced tetrameric enzymes as assessed by sedimentation velocity ultracentrifugation. The IDH1 G15D substitution in the absence of any truncation of IDH1 was subsequently found to be sufficient for production of a tetrameric enzyme. The tetrameric forms of IDH exhibited ∼50% reductions in V(max) and in cooperativity with respect to isocitrate relative to those of the wild-type enzyme, but they retained the property of allosteric activation by AMP. The truncated (-5)IDH1/IDH2 and tetrameric enzymes were much more sensitive than the wild-type enzyme to inhibition by the oxidant diamide and concomitant formation of a disulfide bond between IDH2 Cys-150 residues. Binding of ligands reduced the sensitivity of the wild-type enzyme to diamide but had no effect on inhibition of the truncated or tetrameric enzymes. These results suggest that the octameric structure of IDH has in part evolved for regulation of disulfide bond formation and activity by ensuring the proximity of the amino terminus of an IDH1 subunit of one tetramer to the IDH2 Cys-150 residues in the other tetramer.

Each conserved active site tyr in the three subunits of human isocitrate dehydrogenase has a different function

The human NAD-dependent isocitrate dehydrogenase (IDH) is a heterotetrameric mitochondrial enzyme with 2alpha:1beta:1gamma subunit ratio. The three subunits share 40-52% identity in amino acid sequence and each includes a tyrosine in a comparable position: alphaY126, betaY137, and gammaY135. To study the role of the corresponding tyrosines of each of the subunits of human NAD-IDH, the tyrosines were mutated (one subunit at a time) to Ser, Phe, or Glu. Enzymes were expressed with one mutant and two wild-type subunits. The results of characterization of the mutant enzymes suggest that betaY137 is involved in NAD binding and allosteric activation by ADP. The alphaY126 is required for catalytic activity and likely acts as a general acid in the reaction. The gammaY135 is also required for catalytic activity and may be involved in proper folding of the enzyme. The corresponding tyrosines in the three dissimilar subunits of NAD-IDH thus have distinctive functions.

Dynamics and control of the central carbon metabolism in hepatoma cells

BACKGROUND

The liver plays a major role in metabolism and performs a number of vital functions in the body. Therefore, the determination of hepatic metabolite dynamics and the analysis of the control of the respective biochemical pathways are of great pharmacological and medical importance. Extra- and intracellular time-series data from stimulus-response experiments are gaining in importance in the identification of in vivo metabolite dynamics, while dynamic network models are excellent tools for analyzing complex metabolic control patterns. This is the first study that has been undertaken on the data-driven identification of a dynamic liver central carbon metabolism model and its application in the analysis of the distribution of metabolic control in hepatoma cells.

RESULTS

Dynamic metabolite data were collected from HepG2 cells after they had been deprived of extracellular glucose. The concentration of 25 extra- and intracellular intermediates was quantified using HPLC, LC-MS-MS, and GC-MS. The in silico metabolite dynamics were in accordance with the experimental data. The central carbon metabolism of hepatomas was further analyzed with a particular focus on the control of metabolite concentrations and metabolic fluxes. It was observed that the enzyme glucose-6-phosphate dehydrogenase exerted substantial negative control over the glycolytic flux, whereas oxidative phosphorylation had a significant positive control. The control over the rate of NADPH consumption was found to be shared between the NADPH-demand itself (0.65) and the NADPH supply (0.38).

CONCLUSIONS

Based on time-series data, a dynamic central carbon metabolism model was developed for the investigation of new and complex metabolic control patterns in hepatoma cells. The control patterns found support the hypotheses that the glucose-6-phosphate dehydrogenase and the Warburg effect are promising targets for tumor treatment. The systems-oriented identification of metabolite dynamics is a first step towards the genome-based assessment of potential risks posed by nutrients and drugs.

Long-term prediction of fish growth under varying ambient temperature using a multiscale dynamic model

Background

Feed composition has a large impact on the growth of animals, particularly marine fish. We have developed a quantitative dynamic model that can predict the growth and body composition of marine fish for a given feed composition over a timespan of several months. The model takes into consideration the effects of environmental factors, particularly temperature, on growth, and it incorporates detailed kinetics describing the main metabolic processes (protein, lipid, and central metabolism) known to play major roles in growth and body composition.

Results

For validation, we compared our model's predictions with the results of several experimental studies. We showed that the model gives reliable predictions of growth, nutrient utilization (including amino acid retention), and body composition over a timespan of several months, longer than most of the previously developed predictive models.

Conclusion

We demonstrate that, despite the difficulties involved, multiscale models in biology can yield reasonable and useful results. The model predictions are reliable over several timescales and in the presence of strong temperature fluctuations, which are crucial factors for modeling marine organism growth. The model provides important improvements over existing models.

Disulfide bond formation in yeast NAD+-specific isocitrate dehydrogenase

The tricarboxylic acid cycle NAD+-specific isocitrate dehydrogenase (IDH) of Saccharomyces cerevisiae is an octameric enzyme composed of four heterodimers of regulatory IDH1 and catalytic IDH2 subunits. Recent structural analyses revealed the close proximity of Cys-150 residues from IDH2 in adjacent heterodimers, and features of the structure for the ligand-free enzyme suggested that formation of a disulfide bond between these residues might stabilize an inactive form of the enzyme. We constructed two mutant forms of IDH, one containing a C150S substitution in IDH2 and the other containing C56S/C242S substitutions in IDH2 leaving Cys-150 as the sole cysteine residue. Treatment of the affinity-purified enzymes with diamide resulted in the formation of disulfide bonds and in decreased activities for the wild-type and C56S/C242S enzymes. Both effects were reversible by the addition of dithiothreitol. Diamide had no effect on the C150S mutant enzyme, suggesting that Cys-150 is essential for the formation of a disulfide bond that inhibits IDH activity. Diamide-induced formation of the Cys-150 disulfide bond was also observed in vivo for yeast transformants expressing the wild-type or C56S/C242S enzymes but not for a transformant expressing the C150S enzyme. Finally, natural formation of the Cys-150 disulfide bond with a concomitant decrease in cellular IDH activity was observed during the stationary phase for the parental strain and for transformants expressing wild-type or C56S/C242S enzymes but not for a transformant expressing the C150S enzyme. A reduction in viability for the latter strain suggests that a decrease in IDH activity is important for metabolic changes in stationary phase cells.

Insights from retinitis pigmentosa into the roles of isocitrate dehydrogenases in the Krebs cycle

Here we describe two families with retinitis pigmentosa, a hereditary neurodegeneration of rod and cone photoreceptors in the retina. Affected family members were homozygous for loss-of-function mutations in IDH3B, encoding the beta-subunit of NAD-specific isocitrate dehydrogenase (NAD-IDH, or IDH3), which is believed to catalyze the oxidation of isocitrate to alpha-ketoglutarate in the citric acid cycle. Cells from affected individuals had a substantial reduction of NAD-IDH activity, with about a 300-fold increase in the K(m) for NAD. NADP-specific isocitrate dehydrogenase (NADP-IDH, or IDH2), an enzyme that catalyzes the same reaction, was normal in affected individuals, and they had no health problems associated with the enzyme deficiency except for retinitis pigmentosa. These findings support the hypothesis that mitochondrial NADP-IDH, rather than NAD-IDH, serves as the main catalyst for this reaction in the citric acid cycle outside the retina, and that the retina has a particular requirement for NAD-IDH.

Detailed kinetics and regulation of mammalian NAD-linked isocitrate dehydrogenase

A mathematical model is presented to describe the catalytic mechanism of mammalian NAD-linked isocitrate dehydrogenase (NAD-IDH), a highly regulated enzyme in the tricarboxylic acid cycle, a crucial pathway in energy metabolism and biosynthesis. The mechanism accounts for allosteric regulation by magnesium-bound isocitrate and EGTA and calcium-bound ATP and ADP. The developed model is used to analyze kinetic data for the cardiac enzyme and to estimate kinetic parameter values. Since the kinetic mechanism is expressed in terms of chemical species (rather than biochemical reactants), the model explicitly accounts for the effects of biochemical state (ionic strength, pH, temperature, and metal cation concentration) on the kinetics. Because the substrate isocitrate competes with allosteric activators (ATP and ADP) and an inhibitor (EGTA) for metal ion cofactors (Ca(2+) and Mg(2+)), the observed kinetic relationships between reactants, activator and inhibitor concentrations, and catalytic flux are complex. Our analysis reveals that under physiological conditions, the ADP/ATP ratio plays a more significant role than Ca(2+) concentration in regulating the enzyme's activity. In addition, the enzyme is highly sensitive to Mg(2+) concentration in the physiological range, pointing to a potential regulatory role of [Mg(2+)] in mitochondrial energy metabolism.

Allosteric motions in structures of yeast NAD+-specific isocitrate dehydrogenase

Mitochondrial NAD(+)-specific isocitrate dehydrogenases (IDHs) are key regulators of flux through biosynthetic and oxidative pathways in response to cellular energy levels. Here we present the first structures of a eukaryotic member of this enzyme family, the allosteric, hetero-octameric, NAD(+)-specific IDH from yeast in three forms: 1) without ligands, 2) with bound analog citrate, and 3) with bound citrate + AMP. The structures reveal the molecular basis for ligand binding to homologous but distinct regulatory and catalytic sites positioned at the interfaces between IDH1 and IDH2 subunits and define pathways of communication between heterodimers and heterotetramers in the hetero-octamer. Disulfide bonds observed at the heterotetrameric interfaces in the unliganded IDH hetero-octamer are reduced in the ligand-bound forms, suggesting a redox regulatory mechanism that may be analogous to the "on-off" regulation of non-allosteric bacterial IDHs via phosphorylation. The results strongly suggest that eukaryotic IDH enzymes are exquisitely tuned to ensure that allosteric activation occurs only when concentrations of isocitrate are elevated.

Role of alpha-Asp181, beta-Asp192, and gamma-Asp190 in the distinctive subunits of human NAD-specific isocitrate dehydrogenase

Human NAD-dependent isocitrate dehydrogenase (IDH) is allosterically activated by ADP by lowering the Km for isocitrate. The enzyme has three subunit types with distinguishable sequences present in the approximate ratio 2alpha:1beta:1gamma and, per tetramer, binds 2 mol of each ligand. To evaluate whether the subunits also have distinct functions, we replaced equivalent aspartates, one subunit at a time, by asparagines; each expressed, purified enzyme was composed of one mutant and two wild-type subunits. The aspartates were chosen because beta-Asp192 and gamma-Asp190 had previously been affinity labeled by a reactive ADP analogue and alpha-Asp181 is equivalent based on sequence alignments. The alpha-D181N IDH mutant exhibits a 2000-fold decrease in Vmax, with increases of 15-fold in the Kms for Mn(II) and NAD and a much smaller change in the Km for isocitrate. In contrast, the Vmax values of the beta-D192N and gamma-D190N IDHs are only reduced 4-5-fold as compared to wild-type enzyme. The Km for NAD of the beta-D192N enzyme is 9 times that of the normal enzyme with little or no effect on the affinity for Mn(II) or isocitrate, while the Kms for coenzyme and for Mn(II) of the gamma-D190N enzyme are 19 and 72 times, respectively, that of the normal enzyme with a much smaller effect on the Km for isocitrate. Finally, all three mutant enzymes fail to respond to ADP by lowering the Km for isocitrate, although they do bind ADP. Thus, these aspartates are close to but not in the ADP site and are required for communication between the ADP and isocitrate sites. These results demonstrate that alpha-Asp181 is the only one of these aspartates essential for catalysis. Beta-Asp192 is a determinant of the enzyme's affinity for NAD, as is gamma-Asp190, while gamma-Asp190 also influences the enzyme's affinity for metal ion. We conclude that the NAD and ADP sites are shared between alpha- and beta- and alpha- and gamma-subunits, and the Mn(II) site is shared between alpha- and gamma-subunits, while the alpha-subunit is essential for catalysis. Although alpha-Asp181, beta-Asp192, and gamma-Asp190 may have derived from a common progenitor, these aspartates of the three subunits have evolved distinct functions.

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

Isocitrate and isopropylmalalte dehydrogenases are homologous enzymes important for the cell metabolism. They oxidize their substrates using NAD or NADP as cofactors. Thus, they have two specificities, towards the substrate and the cofactor, appearing in three combinations. Although many three-dimensional (3D) structures are resolved, identification of amino acids determining these specificities remains a challenge. We present computational identification and analysis of specificity-determining positions (SDPs). Besides many experimentally proven SDPs, we predict new SDPs, for example, four substrate-specific positions (103Leu, 105Thr, 337Ala, and 341Thr in IDH from E. coli) that contact the cofactor and may play a role in the recognition process.

Identification of Mn2+-binding aspartates from alpha, beta, and gamma subunits of human NAD-dependent isocitrate dehydrogenase

The human NAD-dependent isocitrate dehydrogenase (IDH), with three types of subunits present in the ratio of 2alpha:1beta:1gamma, requires a divalent metal ion to catalyze the oxidative decarboxylation of isocitrate. With the aim of identifying ligands of the enzyme-bound Mn(2+), we mutated aspartates on the alpha, beta, or gamma subunits. Mutagenesis target sites were based on crystal structures of metal-isocitrate complexes of Escherichia coli and pig mitochondrial NADP-IDH and sequence alignments. Aspartates replaced by asparagine or cysteine were 206, 230, and 234 of the alpha subunit and those corresponding to alpha-Asp-206: 217 of the beta subunit and 215 of the gamma subunit. Each expressed, purified mutant enzyme has two wild-type subunits and one subunit with a single mutation. Specific activities of WT, alpha-D206N, alpha-D230C, alpha-D234C, beta-D217N, and gamma-D215N enzymes are 22, 29, 1.4, 0.2, 7.3 and 3.7 micromol of NADH/min/mg, respectively, whereas alpha-D230N and alpha-D234N enzymes showed no activity. The K(m,Mn(2+)) for alpha-D230C and gamma-D215N are increased 32- and 100-fold, respectively, along with elevations in K(m,isocitrate). The K(m,NAD) of alpha-D230C is increased 16-fold, whereas that of beta-D217N is elevated 10-fold. For all the mutants K(m,isocitrate) is decreased by ADP, indicating that these aspartates are not needed for normal ADP activation. This study demonstrates that alpha-Asp-230 and alpha-Asp-234 are critical for catalytic activity, but alpha-Asp-206 is not needed; alpha-Asp-230 and gamma-Asp-215 may interact directly with the Mn(2+); and alpha-Asp-230 and beta-Asp-217 contribute to the affinity of the enzyme for NAD. These results suggest that the active sites of the human NAD-IDH are shared between alpha and gamma subunits and between alpha and beta subunits.

Potential use of sugar binding proteins in reactors for regeneration of CO2 fixation acceptor D-Ribulose-1,5-bisphosphate

Sugar binding proteins and binders of intermediate sugar metabolites derived from microbes are increasingly being used as reagents in new and expanding areas of biotechnology. The fixation of carbon dioxide at emission source has recently emerged as a technology with potentially significant implications for environmental biotechnology. Carbon dioxide is fixed onto a five carbon sugar D-ribulose-1,5-bisphosphate. We present a review of enzymatic and non-enzymatic binding proteins, for 3-phosphoglycerate (3PGA), 3-phosphoglyceraldehyde (3PGAL), dihydroxyacetone phosphate (DHAP), xylulose-5-phosphate (X5P) and ribulose-1,5-bisphosphate (RuBP) which could be potentially used in reactors regenerating RuBP from 3PGA. A series of reactors combined in a linear fashion has been previously shown to convert 3-PGA, (the product of fixed CO₂ on RuBP as starting material) into RuBP (Bhattacharya et al., 2004; Bhattacharya, 2001). This was the basis for designing reactors harboring enzyme complexes/mixtures instead of linear combination of single-enzyme reactors for conversion of 3PGA into RuBP. Specific sugars in such enzyme-complex harboring reactors requires removal at key steps and fed to different reactors necessitating reversible sugar binders. In this review we present an account of existing microbial sugar binding proteins and their potential utility in these operations.