Comparative fluxome and metabolome analysis for overproduction of succinate in Escherichia coli

An aerobic succinate-producing Escherichia coli mutant was compared to its wild-type by quantitatively analyzing both the metabolome and fluxome, during glucose-limited steady-state and succinate excess dynamic conditions, in order to identify targets for further strain engineering towards more efficient succinate production. The mutant had four functional mutations under the conditions investigated: increased expression of a succinate exporter (DcuC), deletion of a succinate importer (Dct), deletion of succinate dehydrogenase (SUCDH) and expression of a PEP carboxylase (PPC) with increased capacity due to a point mutation. The steady-state and dynamic patterns of the intracellular metabolite levels and fluxes in response to changes were used to locate the quantitative differences in the physiology/metabolism of the mutant strain. Unexpectedly the mutant had a higher energy efficiency, indicated by a much lower rate of oxygen consumption, under glucose-limited conditions, caused by the deletion of the transcription factors IclR and ArcA. Furthermore the mutant had a much lower uptake capacity for succinate (26-fold) and oxygen (17-fold under succinate excess) compared to the wild-type strain. The mutant strain produced 7.9 mmol.CmolX(-1).h(-1) succinate during chemostat cultivation, showing that the choice of the applied genetic modifications was a successful strategy. Furthermore, the applied genetic modifications resulted in multiple large changes in metabolite levels (FBP, pyruvate, 6PG, NAD(+) /NADH ratio, α-ketogluarate) corresponding to large changes in fluxes. Compared to the wild-type a considerable flux shift occurred from the tricarboxylic acid (TCA) cycle to the oxidative part of the pentose phosphate pathway, including an inversion of the pyruvate kinase flux. The mutant responded very differently to excess of succinate, with a remarkable possible reversal of the TCA cycle. The mutant and the wild-type both showed homeostatic behaviour with respect to the energy charge. In contrast, large changes in redox ratios (NAD(+) /NADH) occurred in the wild-type, while the mutant showed even larger changes. This large redox change can be associated to the reversal of flux directions. The observed large flexibility in the central metabolism following genetic (deletions) and environmental (substrate excess) perturbations of the mutant, indicates that introducing a more efficient succinate exporter could result in an even higher succinate production rate.

Tracking the sterol biosynthesis pathway of the diatom Phaeodactylum tricornutum

Diatoms are unicellular photosynthetic microalgae that play a major role in global primary production and aquatic biogeochemical cycling. Endosymbiotic events and recurrent gene transfers uniquely shaped the genome of diatoms, which contains features from several domains of life. The biosynthesis pathways of sterols, essential compounds in all eukaryotic cells, and many of the enzymes involved are evolutionarily conserved in eukaryotes. Although well characterized in most eukaryotes, the pathway leading to sterol biosynthesis in diatoms has remained hitherto unidentified. Through the DiatomCyc database we reconstructed the mevalonate and sterol biosynthetic pathways of the model diatom Phaeodactylum tricornutum in silico. We experimentally verified the predicted pathways using enzyme inhibitor, gene silencing and heterologous gene expression approaches. Our analysis revealed a peculiar, chimeric organization of the diatom sterol biosynthesis pathway, which possesses features of both plant and fungal pathways. Strikingly, it lacks a conventional squalene epoxidase and utilizes an extended oxidosqualene cyclase and a multifunctional isopentenyl diphosphate isomerase/squalene synthase enzyme. The reconstruction of the P. tricornutum sterol pathway underscores the metabolic plasticity of diatoms and offers important insights for the engineering of diatoms for sustainable production of biofuels and high-value chemicals.

Detection of flavonoids in microalgae from different evolutionary lineages

Flavonoids are important secondary plant metabolites believed to be present mainly in land plants. As phenolics were detected previously in microalgae using photometric assays, we wanted to investigate the nature of these phenolics and verify whether flavonoids are present. Therefore, in this study, we used state-of-the-art ultra-high performance liquid chromatography-two-dimensional mass spectrometry (UHPLC-MS/MS) technology to investigate whether microalgae also contain flavonoids. For this, representative microalgal biomass samples from divergent evolutionary lineages (Cyanobacteria, Rhodophyta, Chlorophyta, Haptophyta, Ochrophyta) were screened for a set of carefully selected precursors, intermediates, and end products of the flavonoid biosynthesis pathways. Our data unequivocally showed that microalgae contain a wide range of flavonoids and thus must possess the enzyme pool required for their biosynthesis. Further, some of the microalgae displayed an intricate flavonoid pattern that is compatible with the established basic flavonoid pathway as observed in higher plants. This implies that the flavonoid biosynthesis pathway arose much earlier in evolution compared to what is generally accepted.

CathaCyc, a metabolic pathway database built from Catharanthus roseus RNA-Seq data

The medicinal plant Madagascar periwinkle (Catharanthus roseus) synthesizes numerous terpenoid indole alkaloids (TIAs), such as the anticancer drugs vinblastine and vincristine. The TIA pathway operates in a complex metabolic network that steers plant growth and survival. Pathway databases and metabolic networks reconstructed from 'omics' sequence data can help to discover missing enzymes, study metabolic pathway evolution and, ultimately, engineer metabolic pathways. To date, such databases have mainly been built for model plant species with sequenced genomes. Although genome sequence data are not available for most medicinal plant species, next-generation sequencing is now extensively employed to create comprehensive medicinal plant transcriptome sequence resources. Here we report on the construction of CathaCyc, a detailed metabolic pathway database, from C. roseus RNA-Seq data sets. CathaCyc (version 1.0) contains 390 pathways with 1,347 assigned enzymes and spans primary and secondary metabolism. Curation of the pathways linked with the synthesis of TIAs and triterpenoids, their primary metabolic precursors, and their elicitors, the jasmonate hormones, demonstrated that RNA-Seq resources are suitable for the construction of pathway databases. CathaCyc is accessible online (http://www.cathacyc.org) and offers a range of tools for the visualization and analysis of metabolic networks and 'omics' data. Overlay with expression data from publicly available RNA-Seq resources demonstrated that two well-characterized C. roseus terpenoid pathways, those of TIAs and triterpenoids, are subject to distinct regulation by both developmental and environmental cues. We anticipate that databases such as CathaCyc will become key to the study and exploitation of the metabolism of medicinal plants.

The metabolic blueprint of Phaeodactylum tricornutum reveals a eukaryotic Entner-Doudoroff glycolytic pathway

Diatoms are one of the most successful groups of unicellular eukaryotic algae. Successive endosymbiotic events contributed to their flexible metabolism, making them competitive in variable aquatic habitats. Although the recently sequenced genomes of the model diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana have provided the first insights into their metabolic organization, the current knowledge on diatom biochemistry remains fragmentary. By means of a genome-wide approach, we developed DiatomCyc, a detailed pathway/genome database of P. tricornutum. DiatomCyc contains 286 pathways with 1719 metabolic reactions and 1613 assigned enzymes, spanning both the central and parts of the secondary metabolism of P. tricornutum. Central metabolic pathways, such as those of carbohydrates, amino acids and fatty acids, were covered. Furthermore, our understanding of the carbohydrate model in P. tricornutum was extended. In particular we highlight the discovery of a functional Entner-Doudoroff pathway, an ancient alternative for the glycolytic Embden-Meyerhof-Parnas pathway, and a putative phosphoketolase pathway, both uncommon in eukaryotes. DiatomCyc is accessible online (http://www.diatomcyc.org), and offers a range of software tools for the visualization and analysis of metabolic networks and 'omics' data. We anticipate that DiatomCyc will be key to gaining further understanding of diatom metabolism and, ultimately, will feed metabolic engineering strategies for the industrial valorization of diatoms.

Genome-scale metabolic models: reconstruction and analysis

Metabolism can be defined as the complete set of chemical reactions that occur in living organisms in order to maintain life. Enzymes are the main players in this process as they are responsible for catalyzing the chemical reactions. The enzyme-reaction relationships can be used for the reconstruction of a network of reactions, which leads to a metabolic model of metabolism. A genome-scale metabolic network of chemical reactions that take place inside a living organism is primarily reconstructed from the information that is present in its genome and the literature and involves steps such as functional annotation of the genome, identification of the associated reactions and determination of their stoichiometry, assignment of localization, determination of the biomass composition, estimation of energy requirements, and definition of model constraints. This information can be integrated into a stoichiometric model of metabolism that can be used for detailed analysis of the metabolic potential of the organism using constraint-based modeling approaches and hence is valuable in understanding its metabolic capabilities.

Expression of phosphofructokinase in Neisseria meningitidis

Neisseria meningitidis serogroup B is a pathogen that can infect diverse sites within the human host. According to the N. meningitidis genomic information and experimental observations, glucose can be completely catabolized through the Entner-Doudoroff pathway and the pentose phosphate pathway. The Embden-Meyerhof-Parnas pathway is not functional, because the gene for phosphofructokinase (PFK) is not present. The phylogenetic distribution of PFK indicates that in most obligate aerobic organisms, PFK is lacking. We conclude that this is because of the limited contribution of PFK to the energy supply in aerobically grown organisms in comparison with the energy generated through oxidative phosphorylation. Under anaerobic or microaerobic conditions, the available energy is limiting, and PFK provides an advantage, which explains the presence of PFK in many (facultatively) anaerobic organisms. In accordance with this, in silico flux balance analysis predicted an increase of biomass yield as a result of PFK expression. However, analysis of a genetically engineered N. meningitidis strain that expressed a heterologous PFK showed that the yield of biomass on substrate decreased in comparison with a pfkA-deficient control strain, which was associated mainly with an increase in CO(2) production, whereas production of by-products was similar in the two strains. This might explain why the pfkA gene has not been obtained by horizontal gene transfer, since it is initially unfavourable for biomass yield. No large effects related to heterologous expression of pfkA were observed in the transcriptome. Although our results suggest that introduction of PFK does not contribute to a more efficient strain in terms of biomass yield, achievement of a robust, optimal metabolic network that enables a higher growth rate or a higher biomass yield might be possible after adaptive evolution of the strain, which remains to be investigated.

Modeling Neisseria meningitidis B metabolism at different specific growth rates

Neisseria meningitidis is a human pathogen that can infect diverse sites within the human host. The major diseases caused by N. meningitidis are responsible for death and disability, especially in young infants. At the Netherlands Vaccine Institute (NVI) a vaccine against serogroup B organisms is currently being developed. This study describes the influence of the growth rate of N. meningitidis on its macro-molecular composition and its metabolic activity and was determined in chemostat cultures. In the applied range of growth rates, no significant changes in RNA content and protein content with growth rate were observed in N. meningitidis. The DNA content in N. meningitidis was somewhat higher at the highest applied growth rate. The phospholipid and lipopolysaccharide content in N. meningitidis changed with growth rate but no specific trends were observed. The cellular fatty acid composition and the amino acid composition did not change significantly with growth rate. Additionally, it was found that the PorA content in outer membrane vesicles was significantly lower at the highest growth rate. The metabolic fluxes at various growth rates were calculated using flux balance analysis. Errors in fluxes were calculated using Monte Carlo Simulation and the reliability of the calculated flux distribution could be indicated, which has not been reported for this type of analysis. The yield of biomass on substrate (Y(x/s)) and the maintenance coefficient (m(s)) were determined as 0.44 (+/-0.04) g g(-1) and 0.04 (+/-0.02) g g(-1) h(-1), respectively. The growth associated energy requirement (Y(x/ATP)) and the non-growth associated ATP requirement for maintenance (m(ATP)) were estimated as 0.13 (+/-0.04) mol mol(-1) and 0.43 (+/-0.14) mol mol(-1) h(-1), respectively. It was found that the split ratio between the Entner-Doudoroff and the pentose phosphate pathway, the sole glucose utilizing pathways in N. meningitidis, had a minor effect on ATP formation rate but a major effect on the fluxes going through for instance the citric-acid cycle. For this reason, we presented flux ranges for underdetermined parts of metabolic network rather than presenting single flux values, which is more commonly done in literature.

Scale-up for bulk production of vaccine against meningococcal disease

At the Netherlands Vaccine Institute (NVI) a vaccine against Neisseria meningitidis serogroup B organisms based on different porA subtypes contained in outer membrane vesicles (OMVs) is in advanced stage of development and will be evaluated in clinical trial studies in the near future. In order to meet the expected demand for product, the current biopharmaceutical production process is being scaled-up. This study describes the scale-up approach for the upstream process and the resulting bioreactor design and operation strategy leading towards a feasible solution for bulk production of a vaccine against meningococcal disease. The technically realized 1.2 m(3) bioreactor, equipped with a turbine impeller for gas dispersion, was complemented with an upward pumping impeller and a rotary plate foam breaker to contain foam inside the bioreactor. Aeration and ventilation in the culture broth were controlled by increasing the stirrer speed and gas flow rate simultaneously at increasing oxygen demand. The scale-up was successful and comparable growth curves and nutrient consumption profiles were reached on 0.06 and 1.2 m(3).

Identification of Formaldehyde-Induced Modifications in Proteins: Reactions with Insulin

Formaldehyde is frequently used to inactivate, stabilize, or immobilize proteins. The treatment results in a large variety of chemical modifications in proteins, such as the formation of methylol groups, Schiff bases, and methylene bridges. The purpose of the present study was to identify the stable formaldehyde-induced modifications in a small protein. Therefore, insulin was treated with excess formaldehyde (CH2O) or deuterated formaldehyde (CD2O). In a separate experiment, insulin was modified by formaldehyde (CH2O vs CD2O) and glycine. The mixture of CH2O-treated and CD2O-treated insulin was digested by the proteinase Glu-C. The peptide fragments obtained were analyzed by liquid chromatography-mass spectrometry (LC-MS). Seven intramolecular cross-links were identified in formaldehyde-treated insulin. Furthermore, eight out of the sixteen potentially reactive sites of the insulin molecule were modified by incubation with formaldehyde and glycine. Both the location and the chemical nature of the modifications could be assigned based on the mass increase of potential adducts as elucidated in our previous study (B. Metz et al. (2004) J. Biol. Chem. 279, 6235-6243). To confirm the assigned structures, LC-MS measurements with collision-induced dissociation (LC-MS/MS) were performed on insulin fragments. The results of the LC-MS/MS analyses agreed excellently with the assignments. The study showed that arginine, tyrosine, and lysine residues were very reactive. However, eight theoretically reactive residues did not show detectable modifications, probably because of their low intrinsic reactivity, inaccessibility, or both. The asparagine, glutamine, and histidine residues were not converted in insulin. The N-termini of insulin were partly converted to the expected imidazolidinone adducts, indicating that the protein conformation affects the accessibility and reactivity of these residues. In conclusion, this study shows that, based on our current insights in the chemistry of the reactions between proteins and formaldehyde, we are able to elucidate the location and nature of formaldehyde-induced modifications in a small protein. The approach followed in this study may be generally applicable to larger formaldehyde-treated proteins, such as toxoids used in vaccines.