Photons, Photosynthesis, and High-Performance Computing: Challenges, Progress, and Promise of Modeling Metabolism in Green Algae

The complexity associated with biological metabolism considered at a kinetic level presents a challenge to quantitative modeling. In particular, the relatively sparse knowledge of parameters for enzymes with known kinetic responses is problematic. The possible space of these parameters is of high-dimension, and sampling of such a space typifies a combinatorial explosion of possible dynamic states. However, with sufficient quantitative transcriptomics, proteomics, and metabolomics data at hand, these challenges could be met by high-performance software with sampling, fitting, and optimization capabilities. With this in mind, we present the High-Performance Systems Biology Toolkit HiPer SBTK, an evolving software package to simulate, fit, and optimize metabolite concentrations and fluxes within the space of rate and binding parameters associated with detailed enzyme kinetic models. We present our chosen modeling paradigm for the formulation of metabolic pathway models, the means to address the challenge of representing such models in a precise and persistent fashion using the standardized Systems Biology Markup Language, and our second-generation model of H2-associated Chlamydomonas metabolism. Processing of such models for hierarchically parallelized simulation and optimization, job specification by the user through a GUI interface, software capabilities and initial scaling data, and the mapping of the computation to biological questions is also discussed. Moreover, we present near-term future software and model development goals.

[1]  Xin-Guang Zhu,et al.  Optimizing the Distribution of Resources between Enzymes of Carbon Metabolism Can Dramatically Increase Photosynthetic Rate: A Numerical Simulation Using an Evolutionary Algorithm1[W][OA] , 2007, Plant Physiology.

[2]  K. Vinnakota,et al.  Computer Modeling of Mitochondrial Tricarboxylic Acid Cycle, Oxidative Phosphorylation, Metabolite Transport, and Electrophysiology* , 2007, Journal of Biological Chemistry.

[3]  S. Schnell,et al.  A systematic investigation of the rate laws valid in intracellular environments. , 2006, Biophysical chemistry.

[4]  Michael Hucka,et al.  LibSBML: an API Library for SBML , 2008, Bioinform..

[5]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[6]  Michael Darsow,et al.  ChEBI: a database and ontology for chemical entities of biological interest , 2007, Nucleic Acids Res..

[7]  M. Ghirardi,et al.  Hydrogen Photoproduction Is Attenuated by Disruption of an Isoamylase Gene in Chlamydomonas reinhardtii , 2004, The Plant Cell Online.

[8]  H. Qian,et al.  Thermodynamic constraints for biochemical networks. , 2004, Journal of theoretical biology.

[9]  Daniel A. Beard,et al.  A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation , 2005, PLoS Comput. Biol..

[10]  H. Qian,et al.  Energy balance for analysis of complex metabolic networks. , 2002, Biophysical journal.

[11]  Peter Graf,et al.  Addressing unknown constants and metabolic network behaviors through petascale computing: understanding H2 production in green algae , 2007 .

[12]  Arthur R. Grossman,et al.  Anaerobic Acclimation in Chlamydomonas reinhardtii , 2007, Journal of Biological Chemistry.

[13]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[14]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[15]  B. Palsson,et al.  Genome-scale Reconstruction of Metabolic Network in Bacillus subtilis Based on High-throughput Phenotyping and Gene Essentiality Data* , 2007, Journal of Biological Chemistry.

[16]  Stefanie Widder,et al.  The SBML ODE Solver Library: a native API for symbolic and fast numerical analysis of reaction networks , 2006, Bioinform..

[17]  Barbara M. Bakker,et al.  Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. , 2000, European journal of biochemistry.

[18]  B. Rumberg,et al.  Kinetic modeling of the photosynthetic electron transport chain. , 2001, Bioelectrochemistry.

[19]  Bernhard O. Palsson,et al.  Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets , 2007 .

[20]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..