Yokota 1985, Planta

edit Edit

The Yokota 1985 model is a kinetic simulation of the peroxisomal/glycolate branch of photorespiration in C₃ plants and the alga Euglena gracilis. It grew out of the authors' measurement that glyoxylate reacts non-enzymatically with hydrogen peroxide as a second-order reaction (rate constant 2.27 l mol⁻¹ s⁻¹ at pH 8.0, 25 °C), a reaction that had been proposed as a non-enzymatic, decarboxylating source of photorespiratory CO₂. Using this rate constant together with reported Michaelis-Menten parameters for glycolate oxidase, glyoxylate:glutamate/alanine aminotransferase, glycine decarboxylase, serine:glyoxylate aminotransferase, hydroxypyruvate (glycerate) reductase and catalase, the authors built a reaction-kinetic model of the pathway glycolate → glyoxylate → glycine → serine → hydroxypyruvate, with H₂O₂ produced alongside glyoxylate and removed by catalase.

The GreenSloth implementation mirrors this reaction network directly, with each enzyme step given as a Michaelis-Menten (or ping-pong bi-bi, for serine:glyoxylate aminotransferase) rate law and glycolate supplied at a constant rate from phosphoglycolate phosphatase. Simulating the full system let the authors show that essentially all of the H₂O₂ generated during glycolate oxidation is removed by catalase rather than by reaction with glyoxylate, meaning the non-enzymatic glyoxylate–H₂O₂ reaction cannot be a significant CO₂ source under physiological conditions — so photorespiratory CO₂ must instead originate from the glycine decarboxylase step, not from glyoxylate decarboxylation. In Euglena, which has comparatively low catalase activity, the simulation showed the same non-enzymatic reaction becomes fast enough to matter, explaining species-specific differences in glycolate metabolism. Together with Hahn 1987 and Zhu 2009, this model rounds out GreenSloth's photorespiration models, focusing specifically on the downstream glycolate pathway rather than the Rubisco oxygenation step itself.

Analysis

Model definition

Variables
SymbolIDInitial value
glycolate0.09
glyoxylate0.7964601770483386
glycine8.999999999424611
serine2.5385608670239126
hydroxypyruvate0.009782608695111009
H2O20.010880542843616855
Parameters
SymbolIDValue
kf_phosphoglycolate_phosphatase60
E0_glycolate_oxidase1
kcat_glycolate_oxidase100
km_glycolate_oxidase_s0.06
E0_glycine_transaminase1
kcat_glycine_transaminase143
km_glycine_transaminase_s3
E0_glycine_decarboxylase0.5
kcat_glycine_decarboxylase100
km_glycine_decarboxylase_s6
E0_serine_glyoxylate_transaminase1
kcat_serine_glyoxylate_transaminase159
km_serine_glyoxylate_transaminase_glyoxylate0.15
km_serine_glyoxylate_transaminase_serine2.72
E0_glycerate_dehydrogenase1
kcat_glycerate_dehydrogenase398
km_glycerate_dehydrogenase_s0.12
E0_catalase1
kcat_catalase760500
km_catalase_s137.9
Derived quantities
SymbolIDEquation
vmax_glycolate_oxidase
vmax_glycine_transaminase
vmax_glycine_decarboxylase
vmax_serine_glyoxylate_transaminase
vmax_glycerate_dehydrogenase
vmax_catalase
Reactions
SymbolIDRateStoichiometry
phosphoglycolate_phosphatase
glycolate_oxidase
glycine_transaminase
glycine_decarboxylase
serine_glyoxylate_transaminase
glycerate_dehydrogenase
catalase

Curation

Curator's note

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

Figures
Fig6a
Page figure