Matuszynska 2019, Physiol. Plant.

edit Edit

The Matuszynska 2019 model is a kinetic model of photosynthesis that couples the photosynthetic electron transport chain (PETC) to the carbon-fixing reactions, developed to explore — from a theoretical perspective — how a plant balances the supply of energy (ATP) and reducing power (NADPH) against their demand. It resolves the full electron-transport chain — photosystem II (PSII), the plastoquinone (PQ) pool, the cytochrome b6f complex, plastocyanin, photosystem I (PSI), FNR and ATP synthase — together with the regulatory layers of non-photochemical quenching (PsbS protonation and the xanthophyll cycle) and light-harvesting complex state transitions, and connects this machinery to the Calvin-Benson-Bassham (CBB) cycle that consumes ATP and NADPH. By tracking lumen and stroma proton concentrations and the resulting pH, it links the proton motive force directly to ATP synthesis and to the activation of photoprotective quenching.

The model builds on the earlier Matuszyńska light-acclimation models (see the Matuszynska 2016 NPQ model and the fuller PhD-thesis version) by adding explicit carbon metabolism, which lets it address the mismatch between the fixed ATP:NADPH stoichiometry produced by linear electron flow and the different ratio required by the CBB cycle — a central question in understanding why mechanisms such as cyclic electron flow and quenching exist. It went on to serve as the foundation for the Saadat 2021 model, which extends it with the ascorbate-glutathione cycle and additional electron flows around PSI. Because it captures the coupling between light reactions and carbon fixation in a single, well-documented framework, this model is included in GreenSloth as a versatile base for studying the energetic regulation of photosynthesis.

Scheme

Matuszynska 2019 scheme

Analysis

Model definition

Variables
SymbolIDInitial value
_3PGA0.9928653922138561
BPGA0.0005297732935310749
GAP0.0062663539939955834
DHAP0.13785977143668732
FBP0.006133532145409954
F6P0.31271973359685457
G6P0.719255387166192
G1P0.041716812452951633
SBP0.013123745088361893
S7P0.15890073845176905
E4P0.007322797350442026
X5P0.022478763225333428
R5P0.037651927659696716
RUBP0.13184790283048484
RU5P0.015060770937455408
ATP1.612922506604933
Ferredoxine_oxidised3.8624032084329674
protons_lumen0.002208423037307405
Light_minus_harvesting_complex0.80137477470646
NADPH0.491395685599137
Plastocyanine_oxidised1.885391998090184
Plastoquinone_oxidised10.991562708096392
PsbS_de_minus_protonated0.9610220887579118
Violaxanthin0.9514408605906095
Parameters
SymbolIDValue
protons0.000012589254117941661
pH7.9
CO2_dissolved0.2
O2_dissolved_lumen8
PPFD100
bH100
F96.485
E0_PC0.38
E0_P7000.48
E0_FA-0.55
E0_Fd-0.43
E0_NADP-0.113
convf0.032
R0.0083
T298
Carotenoids_tot1
Fd_star5
PC_tot4
PSBS_tot1
LHC_tot1
gamma00.1
gamma10.25
gamma20.6
gamma30.15
kZSat0.12
E0_QA-0.14
E0_PQ0.354
PQ_tot17.5
staticAntII0.1
staticAntI0.37
NADP_star0.8
A_star_P2.55
Pi_tot17.05
kf_atp_synthase20
HPR4.666666666666667
Pi_mol0.01
DeltaG0_ATP30.6
kcat_b6f2.5
kh_lhc_protonation3
kf_lhc_protonation0.0096
ksat_lhc_protonation5.8
kf_lhc_deprotonation0.0096
kf_cyclic_electron_flow1
kf_violaxanthin_deepoxidase0.0024
kh_violaxanthin_deepoxidase5
ksat_violaxanthin_deepoxidase5.8
kf_zeaxanthin_epoxidase0.00024
km_fnr_Ferredoxine_reduced1.56
km_fnr_NADP0.22
E0_fnr3
kcat_fnr500
kf_ndh0.002
PSII_total2.5
PSI_total2.5
kH0500000000
kPQred250
kPCox2500
kFdred250000
k25000000000
kH5000000000
kF625000000
kf_proton_leak10
kPTOX0.01
kStt70.0035
km_lhc_state_transition_120.2
n_ST2
kPph10.0013
E0_rubisco1
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
E0_ex_g1p1
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
kcat_ex_g1p0.32
Derived quantities
SymbolIDEquation
RT
dG_pH
pH_lumen
Zeaxanthin
Ferredoxine_reduced
Plastocyanine_reduced
PsbS_protonated
Light_minus_harvesting_complex_protonated
Q
keq_Plastoquinone_reduced
Plastoquinone_reduced
PSII_cross_section
NADP
ADP
Orthophosphate
keq_atp_synthase
keq_b6f
keq_fnr
vmax_fnr
keq_PCP700
keq_ferredoxin_reductase
A1
vmax_rubisco_carboxylase
vmax_fbpase
vmax_SBPase
vmax_phosphoribulokinase
vmax_ex_pga
N_translocator
vmax_ex_g1p
B0
B1
B2
B3
PQ_ox_div_tot
Fd_ox_div_tot
PC_ox_div_tot
NADPH_div_tot
ATP_div_tot
Fluo
Reactions
SymbolIDRateStoichiometry
atp_synthase
b6f
lhc_protonation
lhc_deprotonation
cyclic_electron_flow
violaxanthin_deepoxidase
zeaxanthin_epoxidase
fnr
ndh
PSII
PSI
proton_leak
PTOX
lhc_state_transition_12
lhc_state_transition_21
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
Changes

The original model calculates the pH as log10(proton_concentration). To avoid numerical instabilities, we clamped the pH to the range of 1-14 using log10(clamp(0.1, proton_concentration, 1e-14))..

Curation

Curator's note

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

Figures
Fig2
Page figure