Morales 2018, Plant Cell Environ.

edit Edit

The Morales 2018 model — "Dynamic modelling of limitations on improving leaf CO₂ assimilation under fluctuating irradiance" — extends the classical steady-state Farquhar-von Caemmerer-Berry framework into a fully dynamic description of C₃ photosynthesis, built to identify which regulatory processes limit carbon gain when irradiance fluctuates rather than staying constant. On top of the biochemical CO₂-fixation reactions (RuBP and PGA pools, Rubisco carboxylation/oxygenation, TPU and electron-transport limitations), it resolves the kinetics of the mechanisms that actually lag behind a change in light: energy-dependent non-photochemical quenching (qE) and its relaxation, photoinhibition and PSII repair, chloroplast movement (leaf absorptance), Rubisco activation state, and stomatal opening/closure. A four-compartment CO₂/H₂O diffusion pathway — ambient, boundary layer, substomatal and mesophyll/chloroplast — couples these regulatory dynamics back to the leaf gas-exchange measurements the model is validated against.

The model was parameterised and tested against combined gas-exchange and chlorophyll-fluorescence time courses recorded on Arabidopsis thaliana (Col-0) and a panel of mutants specifically impaired in Rubisco activation, qE or sucrose synthesis, which let the authors isolate the individual contribution of each kinetic limitation to the loss of carbon assimilation under fluctuating light. It predicted that the rate of Calvin-cycle enzyme activation and the rate of stomatal opening are the most limiting processes on their own (up to 17% each), with a combined improvement of up to 32% if all kinetic limitations were removed simultaneously — a result later used to guide breeding and engineering targets for photosynthetic efficiency under field-realistic, fluctuating light. It is included in GreenSloth as a detailed, experimentally validated dynamic gas-exchange model that links photoprotection, Rubisco regulation and stomatal behaviour within a single framework.

Analysis

Model definition

Variables
SymbolIDInitial value
PGA0.00005
RuBP0.00005
fRB0.25
fP0
fZ0
alphar1
PSIId0
fR0
PR0
Cc0.00038
Ccyt0.00038
Ci0.00038
Ca0.00038
H2OS0.02
gsw0.09
sumA0
Parameters
SymbolIDValue
Ib100
Ig0
Ir900
Ta298.15
Tl298.15
H2OR0.015
CO2R0.00038
sigma20.5
Jmax250.00013928
DHaJmax36210000
DHdJmax215900000
DsJmax690000
kD0455000000
kDinh5000000000
kf56000000
kp2654000000
fcyc0.1
fpseudo0.1
theta0.7
gamma10.2
gamma20.6
gamma30.2
PhiqEmax0.2
KiqEp0.0187
KdqEp0.0239
KiqEz0.00187
KdqEz0.00239
Kinh00.1
fprot0.1
Krep250.000192
DHaKrep160800000
DHdKrep233230000
DsKrep780000
alphar_alpha256550
DHaAlphar67320000
DHaKalpha90500000
DsKalpha1080000
DHdKalpha328000000
Iac0.0000016
alpharac0.05
alpharav0.25
thetaalphar0.36
Kialpha250.00149
Kdalpha250.00186
fR00.04
alphafR2500
thetafR0.96
KiR0.00628
KdR0.0075
Vrmax0.00011865
KmPGA0.000005
RB0.0000159
Kc254.16
DHaKc41820000
Ko251.26
DHaKo55150000
Kmc250.0002617
DHaKmc49430000
Kmo250.1985
DHaKmo29080000
KaRCA0.0102
ac0.27
bc14000
KdRB0.00068
Krca0.0863
fRBmin0.48
O20.21
RCA0.11737
DHdRCA290200000
ToRCA300.4
DHaRCA30000000
KmRuBP0.02
Vch0.00001
KiPGA0.84
TPU250.00000747
DHaTPU57500000
DHdTPU246700000
DsTPU790000
DHaGc70200000
DHdGc94000000
DsGc320000
DHaGw70200000
DHdGw94000000
DsGw320000
Rm259.9e-7
DHaRm56200000
kPR0.024
Vref0.000155
Scm7.1
Sm9.8
falphaSc0.93
gcm250.39
gw250.75
D0740000
fI00.39
gswm0.48
alphafI767
thetafI0.88
Kgsi0.00114
Kgsd0.00114
gbw9.2
volume_chamber0.00008
leaf_surface0.0002
Flow0.0005
alphab0.92
alphag0.72
alphared0.83
alphabp0.66
alphagp0.6
alpharp0.8
Derived quantities
SymbolIDEquation
PAR
PARa
PARaP
Jmax
Krep
TPU
Kc
Ko
Kmc
Kmo
Rm
gcm
gw
alphar_alpha
Kialpha
Kdalpha
PhiqE
Rp
Sc
gc
Kmapp_RuBP
fRCA
fRBmax
PARaP2
fRuBP
phi
Vc
VrTPU
VrE
A
gm
Fm_d
Fm_a
Fmp_d
Fo_d
Fo_a
Fop_d
fRBss_nr
qI
alpharss
Kinh
PhiIId
Fo
fRss
Fm
kD
Fmp_a
PhiIIop
Fmp
PhiIIo
J2pp
J2pm
Fop_a
qPp
NPQ
qPm
qM
qPno_qD
qE
fqEss
PhiqEss
PhiIIoss
J2qE
VrJ
J2
fRBss_r
Vr
qP
PhiII
reg_limit
fRBss
gbc
gsc
Mv
ea
Photo
Trmmol
es_leaf
VPDleaf
fvpd
gtw
transpiration
fI
gss
Cond
Reactions
SymbolIDRateStoichiometry
dPGA_dt
dRuBP_dt
dPR_dt
dfP_dt
dfZ_dt
dalphar_dt
dPSIId_dt
dfR_dt
dfRB_dt
dCc_dt
dCcyt_dt
dCi_dt
dCa_dt
dH2OS_dt
dgsw_dt

Curation

Curator's note

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

Figures
Fig6
Page figure