Amplified biochemical oscillations in cellular systems

We describe a mechanism for pronounced biochemical oscillations, relevant to microscopic systems, such as the intracellular environment. This mechanism operates for reaction schemes which, when modeled using deterministic rate equations, fail to exhibit oscillations for any values of rate constants. The mechanism relies on amplification of the underlying stochasticity of reaction kinetics within a narrow window of frequencies. This amplification allows fluctuations to beat the central limit theorem, having a dominant effect even though the number of molecules in the system is relatively large. The mechanism is quantitatively studied within simple models of self-regulatory gene expression, and glycolytic oscillations.Comment: 35 pages, 6 figure

Characterization of growth and metabolism of the haloalkaliphile Natronomonas pharaonis

Natronomonas pharaonis is an archaeon adapted to two extreme conditions: high salt concentration and alkaline pH. It has become one of the model organisms for the study of extremophilic life. Here, we present a genome-scale, manually curated metabolic reconstruction for the microorganism. The reconstruction itself represents a knowledge base of the haloalkaliphile's metabolism and, as such, would greatly assist further investigations on archaeal pathways. In addition, we experimentally determined several parameters relevant to growth, including a characterization of the biomass composition and a quantification of carbon and oxygen consumption. Using the metabolic reconstruction and the experimental data, we formulated a constraints-based model which we used to analyze the behavior of the archaeon when grown on a single carbon source. Results of the analysis include the finding that Natronomonas pharaonis, when grown aerobically on acetate, uses a carbon to oxygen consumption ratio that is theoretically near-optimal with respect to growth and energy production. This supports the hypothesis that, under simple conditions, the microorganism optimizes its metabolism with respect to the two objectives. We also found that the archaeon has a very low carbon efficiency of only about 35%. This inefficiency is probably due to a very low P/O ratio as well as to the other difficulties posed by its extreme environment

Systems analysis of bioenergetics and growth of the extreme halophile Halobacterium salinarum

Halobacterium salinarum is a bioenergetically flexible, halophilic microorganism that can generate energy by respiration, photosynthesis, and the fermentation of arginine. In a previous study, using a genome-scale metabolic model, we have shown that the archaeon unexpectedly degrades essential amino acids under aerobic conditions, a behavior that can lead to the termination of growth earlier than necessary. Here, we further integratively investigate energy generation, nutrient utilization, and biomass production using an extended methodology that accounts for dynamically changing transport patterns, including those that arise from interactions among the supplied metabolites. Moreover, we widen the scope of our analysis to include phototrophic conditions to explore the interplay between different bioenergetic modes. Surprisingly, we found that cells also degrade essential amino acids even during phototropy, when energy should already be abundant. We also found that under both conditions considerable amounts of nutrients that were taken up were neither incorporated into the biomass nor used as respiratory substrates, implying the considerable production and accumulation of several metabolites in the medium. Some of these are likely the products of forms of overflow metabolism. In addition, our results also show that arginine fermentation, contrary to what is typically assumed, occurs simultaneously with respiration and photosynthesis and can contribute energy in levels that are comparable to the primary bioenergetic modes, if not more. These findings portray a picture that the organism takes an approach toward growth that favors the here and now, even at the cost of longer-term concerns. We believe that the seemingly "greedy" behavior exhibited actually consists of adaptations by the organism to its natural environments, where nutrients are not only irregularly available but may altogether be absent for extended periods that may span several years. Such a setting probably predisposed the cells to grow as much as possible when the conditions become favorable

A Data Integration and Visualization Resource for the Metabolic Network of Synechocystis sp. PCC 6803

Data integration is a central activity in systems biology. The integration of genomic, transcript, protein, metabolite, flux, and computational data yields unprecedented information about the system level functioning of organisms. Often, data integration is done purely computationally, leaving the user with little insight besides statistical information. In this article, we present a visualization tool for the metabolic network of Synechocystis PCC6803, an important model cyanobacterium for sustainable biofuel production. We illustrate how this metabolic map can be used to integrate experimental and computational data for Synechocystis systems biology and metabolic engineering studies. Additionally, we discuss how this map, and the software infrastructure that we supply with it, can be used in the development of other organism-specific metabolic network visualizations. Besides a Python console package VoNDA (http://vonda.sf.net), we provide a working demonstration of the interactive metabolic map and the associated Synechocystis genome-scale stoichiometric model, as well as various ready-to-visualize microarray data sets, at http://f-a-m-e.org/synechocystis/

Understanding start-up problems in yeast glycolysis

Yeast glycolysis has been the focus of research for decades, yet a number of dynamical aspects of yeast glycolysis remain poorly understood at present. If nutrients are scarce, yeast will provide its catabolic and energetic needs with other pathways, but the enzymes catalysing upper glycolytic fluxes are still expressed. We conjecture that this overexpression facilitates the rapid transition to glycolysis in case of a sudden increase in nutrient concentration. However, if starved yeast is presented with abundant glucose, it can enter into an imbalanced state where glycolytic intermediates keep accumulating, leading to arrested growth and cell death. The bistability between regularly functioning and imbalanced phenotypes has been shown to depend on redox balance. We shed new light on these phenomena with a mathematical analysis of an ordinary differential equation model, including NADH to account for the redox balance. In order to gain qualitative insight, most of the analysis is parameter-free, i.e., without assigning a numerical value to any of the parameters. The model has a subtle bifurcation at the switch between an inviable equilibrium state and stable flux through glycolysis. This switch occurs if the ratio between the flux through upper glycolysis and ATP consumption rate of the cell exceeds a fixed threshold. If the enzymes of upper glycolysis would be barely expressed, our model predicts that there will be no glycolytic flux, even if external glucose would be at growth-permissable levels. The existence of the imbalanced state can be found for certain parameter conditions independent of the mentioned bifurcation. The parameter-free analysis proved too complex to directly gain insight into the imbalanced states, but the starting point of a branch of imbalanced states can be shown to exist in detail. Moreover, the analysis offers the key ingredients necessary for successful numerical continuation, which highlight the existence of this bistability and the influence of the redox balance

Optimal flux spaces of genome-scale stoichiometric models are determined by a few subnetworks

The metabolism of organisms can be studied with comprehensive stoichiometric models of their metabolic networks. Flux balance analysis (FBA) calculates optimal metabolic performance of stoichiometric models. However, detailed biological interpretation of FBA is limited because, in general, a huge number of flux patterns give rise to the same optimal performance. The complete description of the resulting optimal solution spaces was thus far a computationally intractable problem. Here we present CoPE-FBA: Comprehensive Polyhedra Enumeration Flux Balance Analysis, a computational method that solves this problem. CoPE-FBA indicates that the thousands to millions of optimal flux patterns result from a combinatorial explosion of flux patterns in just a few metabolic sub-networks. The entire optimal solution space can now be compactly described in terms of the topology of these sub-networks. CoPE-FBA simplifies the biological interpretation of stoichiometric models of metabolism, and provides a profound understanding of metabolic flexibility in optimal states