Metabolic Futile Cycles and Their Functions: A Systems Analysis of Energy and Control

It has long been hypothesized that futile cycles in cellular metabolism are involved in the regulation of biochemical pathways. Following the work of Newsholme and Crabtree, we develop a quantitative theory for this idea based on open-system thermodynamics and metabolic control analysis. It is shown that the {\it stoichiometric sensitivity} of an intermediary metabolite concentration with respect to changes in steady-state flux is governed by the effective equilibrium constant of the intermediate formation, and the equilibrium can be regulated by a futile cycle. The direction of the shift in the effective equilibrium constant depends on the direction of operation of the futile cycle. High stoichiometric sensitivity corresponds to ultrasensitivity of an intermediate concentration to net flow through a pathway; low stoichiometric sensitivity corresponds to super-robustness of concentration with respect to changes in flux. Both cases potentially play important roles in metabolic regulation. Futile cycles actively shift the effective equilibrium by expending energy; the magnitude of changes in effective equilibria and sensitivities is a function of the amount of energy used by a futile cycle. This proposed mechanism for control by futile cycles works remarkably similarly to kinetic proofreading in biosynthesis. The sensitivity of the system is also intimately related to the rate of concentration fluctuations of intermediate metabolites. The possibly different roles of the two major mechanisms for cellular biochemical regulation, namely reversible chemical modifications via futile cycles and shifting equilibrium by macromolecular binding, are discussed.Comment: 11 pages, 5 figure

Systems biology of lactic acid bacteria: a critical review

Understanding the properties of a system as emerging from the interaction of well described parts is the most important goal of Systems Biology. Although in the practice of Lactic Acid Bacteria (LAB) physiology we most often think of the parts as the proteins and metabolites, a wider interpretation of what a part is can be useful. For example, different strains or species can be the parts of a community, or we could study only the chemical reactions as the parts of metabolism (and forgetting about the enzymes that catalyze them), as is done in flux balance analysis. As long as we have some understanding of the properties of these parts, we can investigate whether their interaction leads to novel or unanticipated behaviour of the system that they constitute

Optimality principles in the regulation of metabolic networks.

One of the challenging tasks in systems biology is to understand how molecular networks give rise to emergent functionality and whether universal design principles apply to molecular networks. To achieve this, the biophysical, evolutionary and physiological constraints that act on those networks need to be identified in addition to the characterisation of the molecular components and interactions. Then, the cellular “task” of the network—its function—should be identified. A network contributes to organismal fitness through its function. The premise is that the same functions are often implemented in different organisms by the same type of network; hence, the concept of design principles. In biology, due to the strong forces of selective pressure and natural selection, network functions can often be understood as the outcome of fitness optimisation. The hypothesis of fitness optimisation to understand the design of a network has proven to be a powerful strategy. Here, we outline the use of several optimisation principles applied to biological networks, with an emphasis on metabolic regulatory networks. We discuss the different objective functions and constraints that are considered and the kind of understanding that they provide

Computing the shortest elementary flux modes in genome-scale metabolic networks

This article is available open access through the publisher’s website through the link below. Copyright @ The Author 2009.Motivation: Elementary flux modes (EFMs) represent a key concept to analyze metabolic networks from a pathway-oriented perspective. In spite of considerable work in this field, the computation of the full set of elementary flux modes in large-scale metabolic networks still constitutes a challenging issue due to its underlying combinatorial complexity. Results: In this article, we illustrate that the full set of EFMs can be enumerated in increasing order of number of reactions via integer linear programming. In this light, we present a novel procedure to efficiently determine the K-shortest EFMs in large-scale metabolic networks. Our method was applied to find the K-shortest EFMs that produce lysine in the genome-scale metabolic networks of Escherichia coli and Corynebacterium glutamicum. A detailed analysis of the biological significance of the K-shortest EFMs was conducted, finding that glucose catabolism, ammonium assimilation, lysine anabolism and cofactor balancing were correctly predicted. The work presented here represents an important step forward in the analysis and computation of EFMs for large-scale metabolic networks, where traditional methods fail for networks of even moderate size. Contact: [email protected] Supplementary information: Supplementary data are available at Bioinformatics online (http://bioinformatics.oxfordjournals.org/cgi/content/full/btp564/DC1).Fundação Calouste Gulbenkian, Fundação para a Ciência e a Tecnologia (FCT) and Siemens SA Portugal

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

Model evaluation for glycolytic oscillations in yeast biotransformations of xenobiotics

Anaerobic glycolysis in yeast perturbed by the reduction of xenobiotic ketones is studied numerically in two models which possess the same topology but different levels of complexity. By comparing both models' predictions for concentrations and fluxes as well as steady or oscillatory temporal behavior we answer the question what phenomena require what kind of minimum model abstraction. While mean concentrations and fluxes are predicted in agreement by both models we observe different domains of oscillatory behavior in parameter space. Generic properties of the glycolytic response to ketones are discussed

Виправте!

BACKGROUND: The aim of the current work was to perform a clinical trial simulation (CTS) analysis to optimize a drug-drug interaction (DDI) study of vincristine in children who also received azole antifungals, taking into account challenges of conducting clinical trials in this population, and, to provide a motivating example of the application of CTS in the design of pediatric oncology clinical trials. PROCEDURE: A pharmacokinetic (PK) model for vincristine in children was used to simulate concentration-time profiles. A continuous model for body surface area versus age was defined based on pediatric growth curves. Informative sampling time windows were derived using D-optimal design. The CTS framework was used to different magnitudes of clearance inhibition (10%, 25%, or 40%), sample size (30-500), the impact of missing samples or sampling occasions, and the age distribution, on the power to detect a significant inhibition effect, and in addition, the relative estimation error (REE) of the interaction effect. RESULTS: A minimum group specific sample size of 38 patients with a total sample size of 150 patients was required to detect a clearance inhibition effect of 40% with 80% power, while in the case of a lower effect of clearance inhibition, a substantially larger sample size was required. However, for the majority of re-estimated drug effects, the inhibition effect could be estimated precisely (REE < 25%) in even smaller sample sizes and with lower effect sizes. CONCLUSION: This work demonstrated the utility of CTS for the evaluation of PK clinical trial designs in the pediatric oncology population