Davis 2017, Philos Trans.

edit Edit

The Davis 2017 model is a mechanistic, dynamic model of the photosynthetic electron transport chain and the thylakoid proton motive force (PMF), built to explain how the PMF can itself limit photosynthesis by driving photosystem II (PSII) photodamage. Its central feature is that it resolves the PMF into its two thermodynamically distinct components: the transmembrane electric field (Δψ) and the lumen pH gradient (ΔpH). The model follows electron flow from PSII through the plastoquinone (PQ) pool, the cytochrome b6f complex, plastocyanin, photosystem I (PSI), ferredoxin and FNR to NADPH and the Calvin-Benson-Bassham (CBB) cycle, while proton accumulation in the lumen and charge separation across the membrane build the PMF that powers ATP synthase. Counter-ion fluxes — a KEA3 K⁺/H⁺ antiporter and a K⁺ channel — partition the PMF between its Δψ and ΔpH components, and non-photochemical quenching is represented through PsbS protonation and the xanthophyll cycle (violaxanthin ↔ zeaxanthin).

The key result of the model is that it is the Δψ component, rather than lumen acidification, that drives elevated PSII charge recombination, producing singlet oxygen and the subsequent photodamage to PSII. This makes the model a useful tool for reasoning about photoprotection under fluctuating light, where rapid changes in the balance between Δψ and ΔpH can transiently expose PSII to damaging conditions. By including explicit ion transport, it also connects to the same questions of PMF partitioning later explored in the Li 2021 model, which builds directly on this work. As such it serves as a good mechanistic base within GreenSloth for studying how the thylakoid membrane manages energy and protects itself against its own light-harvesting machinery.

Scheme

Davis 2017 scheme

Analysis

Model definition

Variables
SymbolIDInitial value
QA_red0
PQH_20
pH_lumen7
Dpsi0.1
K_lu0.04
PC_ox0
Zx0
singO20
P700_ox0
Fd_red0
NADPH_st0
LEF0
ATP_made0
Parameters
SymbolIDValue
PPFD100
k_recomb0.33
phi_triplet0.45
phi_1O21
sigma0_II1
c_b6f1
pKa_reg6.5
Em_PC_pH70.37
Em_PQH2_pH70.11
Vmax_b6f500
pKa_PsbS6.4
NPQ_max5
pH_stroma7.8
PSI_antenna_size1
k_QA1000
Keq_QA200
k_PCtoP700500
k_FdtoNADP1000
K_st0.04
k_KEA30
P_K6000
lumen_protons_per_turnover0.000014
n4.666666666666667
DeltaGATP42
Vmax_ATPsynth1000
b_H0.03
volt_per_charge0.033
k_EZ0.03
nh_VDE4
pKa_VDE5.8
Vmax_VDE1
QA_total1
PQ_tot6
P700_total1
PC_tot2
Fd_tot1
NADP_tot1
Xanthophyll_tot1
k_CBB3000
Derived quantities
SymbolIDEquation
QA
P700_red
PQ
PC_red
Fd_ox
NADP_st
Vx
DeltaGATP_V
PsbSP
NPQ
PhiPSII
H_lumen
H_stroma
pmf
delta_pH
delta_pH_inVolts
ATP_synthase_driving_force
k_b6f
Reactions
SymbolIDRateStoichiometry
vPSII_recomb
vPSII_ChSep
v_PSII
v_PQ
v_b6f
PSI_ChSep
v_PSI_PCoxid
v_FNR
vATPsynthase
v_CBB
v_KEA3
v_VKC
v_Epox
v_Deepox

Curation

Curator's note

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

Figures
Fig3
Page figure
Fig4
Page figure
Fig5
Page figure