title: Simulating metabolism with statistical thermodynamics.
abstract: New methods are needed for large scale modeling of metabolism that predict metabolite levels and characterize the thermodynamics of individual reactions and pathways. Current approaches use either kinetic simulations, which are difficult to extend to large networks of reactions because of the need for rate constants, or flux-based methods, which have a large number of feasible solutions because they are unconstrained by the law of mass action. This report presents an alternative modeling approach based on statistical thermodynamics. The principles of this approach are demonstrated using a simple set of coupled reactions, and then the system is characterized with respect to the changes in energy, entropy, free energy, and entropy production. Finally, the physical and biochemical insights that this approach can provide for metabolism are demonstrated by application to the tricarboxylic acid (TCA) cycle of Escherichia coli. The reaction and pathway thermodynamics are evaluated and predictions are made regarding changes in concentration of TCA cycle intermediates due to 10- and 100-fold changes in the ratio of NAD+:NADH concentrations. Finally, the assumptions and caveats regarding the use of statistical thermodynamics to model non-equilibrium reactions are discussed.
authors.affiliation: Computational Biology and Bioinformatics Group, Fundamental and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington, United States of America.
authors.forename: William R
authors.initials: WR
authors.lastname: Cannon
mesh_terms: D002952:Citric Acid Cycle; D003198:Computer Simulation; D019277:Entropy; D004926:Escherichia coli / Q000378:metabolism; D016013:Likelihood Functions; D008660:Metabolism; D011336:Probability; D013223:Statistics as Topic; D013816:Thermodynamics
references.citation: PLoS Comput Biol. 2014 Feb 20;10(2):e1003483
references.pmid: 24586134
references.citation: Eur J Biochem. 2003 Feb;270(3):415-21
references.pmid: 12542691
references.citation: Proc Natl Acad Sci U S A. 2004 Sep 7;101(36):13168-73
references.pmid: 15340153
references.citation: Bioinformatics. 2004 Nov 1;20(16):2874-7
references.pmid: 15145806
references.citation: Proc Natl Acad Sci U S A. 1922 Jun;8(6):151-4
references.pmid: 16576643
references.citation: PLoS Comput Biol. 2012;8(6):e1002562
references.pmid: 22737065
references.citation: Phys Rev E Stat Nonlin Soft Matter Phys. 2010 May;81(5 Pt 1):051133
references.pmid: 20866211
references.citation: Nucleic Acids Res. 2012 Jan;40(Database issue):D770-5
references.pmid: 22064852
references.citation: Proc Natl Acad Sci U S A. 2000 Jan 4;97(1):32-7
references.pmid: 10618366
references.citation: BMC Syst Biol. 2007 Jun 01;1:23
references.pmid: 17543097
references.citation: Nat Chem Biol. 2009 Aug;5(8):593-9
references.pmid: 19561621
references.citation: Mol Syst Biol. 2006;2:2006.0034
references.pmid: 16788595
references.citation: PLoS One. 2007 Jan 03;2(1):e144
references.pmid: 17206279
references.citation: J Chem Theory Comput. 2010 Mar 4;6(4):1109-1117
references.pmid: 20419072
references.citation: J R Soc Interface. 2014 Feb 19;11(94):20131017
references.pmid: 24554573
references.citation: Database (Oxford). 2011 Apr 11;2011:bar005
references.pmid: 21482578
references.citation: Biophys J. 2007 Mar 1;92(5):1792-805
references.pmid: 17172310
references.citation: Nucleic Acids Res. 2013 Jan;41(Database issue):D764-72
references.pmid: 23203881
references.citation: Proc Natl Acad Sci U S A. 1922 Jun;8(6):147-51
references.pmid: 16576642
references.citation: Biophys J. 2012 Apr 18;102(8):1703-11
references.pmid: 22768925
references.citation: Protein Eng. 2000 Jun;13(6):413-21
references.pmid: 10877852