Poolman 2000, Exp. Botany

edit Edit

The Poolman 2000 model describes the Calvin-Benson-Bassham (CBB) cycle, covering the fixation of CO2 by rubisco through to the regeneration of RuBP and the export of triose phosphates, glycerate-3-phosphate and hexose phosphate from the chloroplast stroma.

This implementation is a variant of the classical Pettersson & Ryde-Petersson CBB cycle model. That original model was formulated as a differential-algebraic system (DAE): several near-equilibrium reactions (e.g. triose phosphate isomerase, the aldolase and transketolase steps, the pentose phosphate isomerase/epimerase, and phosphoglucomutase/PGI) were treated as instantaneous equilibria and enforced via algebraic constraints rather than rate equations. Here, each of those equilibrium constraints is instead replaced by a reversible mass-action rate law with a very fast rate constant (kre_*), driving the corresponding pool ratios to the same equilibrium constant (keq_*) dynamically rather than algebraically. This turns the DAE into a plain ODE system, at the cost of introducing stiffness from the disparate time scales of the fast quasi-equilibrium reactions and the slower, kinetically-limited steps (rubisco carboxylase, FBPase, SBPase, phosphoribulokinase, and the stromal export reactions).

This CBB cycle formulation is, in turn, the basis for the Calvin cycle in several later, more comprehensive photosynthesis models, including Matuszynska 2019 (10.1111/ppl.12962) and Saadat 2021 (10.3389/fpls.2021.750580).

Analysis

Model definition

Variables
SymbolIDInitial value
_3PGA0.6387788347932627
BPGA0.0013570885908749779
GAP0.011259431827358068
DHAP0.24770748227012374
FBP0.01980222074817044
F6P1.093666906864421
G6P2.5154338857582377
G1P0.14589516537322303
SBP0.09132688566151095
S7P0.23281380022778891
E4P0.02836065066520614
X5P0.03647242425941113
R5P0.06109130988031577
RUBP0.2672164362349537
RU5P0.0244365238237522
ATP0.43633201706180874
Parameters
SymbolIDValue
CO2_dissolved0.2
NADPH0.21
protons0.000012589254117941661
A_star_P0.5
NADP_star0.5
Pi_tot15
E0_rubisco_carboxylase1
kcat_rubisco_carboxylase2.72
km_rubisco_carboxylase_RUBP0.02
km_rubisco_carboxylase_CO2_dissolved0.0107
ki_rubisco_carboxylase_3PGA0.04
ki_rubisco_carboxylase_FBP0.04
ki_rubisco_carboxylase_SBP0.075
ki_rubisco_carboxylase_Orthophosphate0.9
ki_rubisco_carboxylase_NADPH0.07
kre_phosphoglycerate_kinase800000000
keq_phosphoglycerate_kinase0.00031
kre_gadph800000000
keq_gadph16000000
kre_triose_phosphate_isomerase800000000
keq_triose_phosphate_isomerase22
kre_aldolase_dhap_gap800000000
keq_aldolase_dhap_gap7.1
kre_aldolase_dhap_e4p800000000
keq_aldolase_dhap_e4p13
E0_fbpase1
kcat_fbpase1.6
km_fbpase_s0.03
ki_fbpase_F6P0.7
ki_fbpase_Orthophosphate12
kre_transketolase_gap_f6p800000000
keq_transketolase_gap_f6p0.084
kre_transketolase_gap_s7p800000000
keq_transketolase_gap_s7p0.85
E0_SBPase1
kcat_SBPase0.32
km_SBPase_s0.013
ki_SBPase_Orthophosphate12
kre_ribose_phosphate_isomerase800000000
keq_ribose_phosphate_isomerase0.4
kre_ribulose_phosphate_epimerase800000000
keq_ribulose_phosphate_epimerase0.67
E0_phosphoribulokinase1
kcat_phosphoribulokinase7.9992
km_phosphoribulokinase_RU5P0.05
km_phosphoribulokinase_ATP0.05
ki_phosphoribulokinase_3PGA2
ki_phosphoribulokinase_RUBP0.7
ki_phosphoribulokinase_Orthophosphate4
ki_phosphoribulokinase_42.5
ki_phosphoribulokinase_50.4
kre_g6pi800000000
keq_g6pi2.3
kre_phosphoglucomutase800000000
keq_phosphoglucomutase0.058
Orthophosphate_external0.5
km_ex_pga0.25
km_ex_gap0.075
km_ex_dhap0.077
km_N_translocator_Orthophosphate_external0.74
km_N_translocator_Orthophosphate0.63
kcat_N_translocator2
E0_N_translocator1
km_ex_g1p_G1P0.08
km_ex_g1p_ATP0.08
ki_ex_g1p10
ki_ex_g1p_3PGA0.1
ki_ex_g1p_F6P0.02
ki_ex_g1p_FBP0.02
E0_ex_g1p1
kcat_ex_g1p0.32
km_atp_synthase_ADP0.014
km_atp_synthase_Orthophosphate0.3
kcat_atp_synthase2.8
E0_atp_synthase1
Derived quantities
SymbolIDEquation
ADP
NADP
Orthophosphate
vmax_rubisco_carboxylase
vmax_fbpase
vmax_SBPase
vmax_phosphoribulokinase
vmax_ex_pga
N_translocator
vmax_ex_g1p
vmax_atp_synthase
Reactions
SymbolIDRateStoichiometry
rubisco_carboxylase
phosphoglycerate_kinase
gadph
triose_phosphate_isomerase
aldolase_dhap_gap
aldolase_dhap_e4p
fbpase
transketolase_gap_f6p
transketolase_gap_s7p
SBPase
ribose_phosphate_isomerase
ribulose_phosphate_epimerase
phosphoribulokinase
g6pi
phosphoglucomutase
ex_pga
ex_gap
ex_dhap
ex_g1p
atp_synthase

Curation

Curator's note

This model was validated by reproducing the following figures of the original publication.

Figures
Fig3
Page figure