Saadat 2021, Front. Plant Sci.

edit Edit

The Saadat 2021 model builds upon previous models, particularly the Matuszynska2019 model, by incorporating and modifying various reactions and aspects of photosynthesis. Overall, the model can be divided into three modules: the ascorbate-glutathione cycle, the Calvin-Benson-Bassham (CBB) cycle and thioredoxin reductase-regulated reactions, and the photosynthetic electron transport chain (PETC).

The model is primarily used to investigate the electron flows around PSI and their relevance to photosynthetic efficiency. Several different analyses have been conducted to validate the model in both steady-state and dynamic environment conditions. The most interesting is the direct comparison of a knockout mutant of the protein PGR5. This protein is known to catalyse the reduction of plastoquinone by ferredoxin. The results of this comparison align with experimental values, which are, however, not presented in the publication but are referenced. Additionally, it is noted that the results should not be interpreted as accurate quantitative data, but rather as a proof of concept for the model.

Overall, the model has one advantage over other photosynthesis models, as it also highlights the importance of other electron flows, not just the PETC. Additionally, not only are the authors open about the model's flaws, but they are also insistent on making their code and analyses available on GitHub. Therefore, this model serves as a good stepping stone for more complex models that aim to incorporate aspects of photosynthesis, which are often simplified in other models.

Scheme

Saadat 2021 scheme

Analysis

Model definition

Variables
SymbolIDInitial value
_3PGA0.9167729479368978
BPGA0.0003814495319659031
GAP0.00580821050261484
DHAP0.1277806166216142
FBP0.005269452472931973
F6P0.2874944558066638
G6P0.6612372482712676
G1P0.03835176039761378
SBP0.011101373736607443
S7P0.1494578301900007
E4P0.00668295494870102
X5P0.020988553174809618
R5P0.035155825913785584
RUBP0.11293260727162346
RU5P0.014062330254191594
atp1.4612747767895344
fd_ox3.715702384326767
protons_lumen0.002086128887296243
lhc0.7805901436176024
nadph0.5578718406315588
pc_ox1.8083642974980014
pq_ox10.251099271612473
psbs_de0.9667381262477079
vx0.9629870646993118
MDA2.0353396709300447e-7
H2O21.2034405327140102e-7
DHA1.0296456279861962e-11
GSSG4.99986167652437e-12
tr_ox0.9334426859846461
E_inactive3.6023635680406634
Parameters
SymbolIDValue
PPFD100
CO2__dissolved_0.2
O2_lumen8
pH7.9
protons0.000012589254117941661
bH100
F96.485
E_0_PC0.38
E_0_P7000.48
E_0_FA-0.55
E_0_Fd-0.43
E_0_NADP-0.113
convf0.032
R0.0083
T298
Carotenoids_tot1
Fd_5
PC_tot4
PSBS_tot1
LHC_tot1
gamma00.1
gamma10.25
gamma20.6
gamma30.15
kZSat0.12
E_0_QA-0.14
E_0_PQ0.354
PQ_tot17.5
staticAntII0.1
staticAntI0.37
Thioredoxin_tot1
E_total6
NADP_0.8
A_P2.55
Pi_tot17.05
kf_ferredoxin_thioredoxin_reductase1
kf_tr_activation1
kf_tr_inactivation0.1
ASC_tot_10
Glutathion_tot10
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_fd_red1.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
kMehler1
E0_ferredoxin_reductase1
kcat_ferredoxin_reductase250000
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__dissolved_0.0107
ki_rubisco_carboxylase_3PGA0.04
ki_rubisco_carboxylase_FBP0.04
ki_rubisco_carboxylase_SBP0.075
ki_rubisco_carboxylase_pi0.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_pi12
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_pi12
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_pi4
ki_phosphoribulokinase_42.5
ki_phosphoribulokinase_50.4
kre_g6pi800000000
keq_g6pi2.3
kre_phosphoglucomutase800000000
keq_phosphoglucomutase0.058
pi_ext0.5
km_ex_pga0.25
km_ex_gap0.075
km_ex_dhap0.077
km_N_translocator_pi_ext0.74
km_N_translocator_pi0.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
kf_mda_reductase_1500
E0_mda_reductase_20.002
kcat_mda_reductase_2300
km_mda_reductase_2_nadph0.023
km_mda_reductase_2_MDA0.0014
kf110000
kr1220
kf210000
kr24000
kf32510
kf410000
kr44000
kf52510
XT0.07
E0_glutathion_reductase0.0014
kcat_glutathion_reductase595
km_glutathion_reductase_nadph0.003
km_glutathion_reductase_GSSG0.2
km_dehydroascorbate_reductase_DHA0.07
km_dehydroascorbate_reductase_GSH2.5
K0.5
E0_dehydroascorbate_reductase0.0017
kcat_dehydroascorbate_reductase142
kf_ex_atp0.2
kf_ex_nadph0.2
Derived quantities
SymbolIDEquation
RT
dG_pH
pH_lumen
zx
fd_red
pc_red
psbs_pr
lhc_prot
Q
keq_pq_red
pq_red
PSII_cross_section
tr_red
E_active
nadp
adp
pi
ascorbate
GSH
keq_atp_synthase
keq_b6f
keq_fnr
vmax_fnr
keq_PCP700
keq_ferredoxin_reductase
vmax_ferredoxin_reductase
E0_rubisco_active
vmax_rubisco_carboxylase
E0_fbpase_active
vmax_fbpase
E0_SBPase_active
vmax_SBPase
E0_phosphoribulokinase_active
vmax_phosphoribulokinase
vmax_ex_pga
N_translocator
E0_ex_g1p_active
vmax_ex_g1p
vmax_mda_reductase_2
vmax_glutathion_reductase
vmax_dehydroascorbate_reductase
B0
B1
B2
B3
A0
A1
A2
Reactions
SymbolIDRateStoichiometry
ferredoxin_thioredoxin_reductase
tr_activation
tr_inactivation
atp_synthase
b6f
lhc_protonation
lhc_deprotonation
cyclic_electron_flow
violaxanthin_deepoxidase
zeaxanthin_epoxidase
fnr
ndh
PSII
PSI
mehler
ferredoxin_reductase
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
mda_reductase_1
mda_reductase_2
ascorbate_peroxidase
glutathion_reductase
dehydroascorbate_reductase
ex_atp
ex_nadph
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
Fig3
Page figure
Fig4
Page figure
Fig5
Page figure
Fig6
Page figure