Computational simulation of the reactive oxygen species and redox network in the regulation of chloroplast metabolism

Cells contain a thiol redox regulatory network to coordinate metabolic and developmental activities with exogenous and endogenous cues. This network controls the redox state and activity of many target proteins. Electrons are fed into the network from metabolism and reach the target proteins via redox transmitters such as thioredoxin (TRX) and NADPH-dependent thioredoxin reductases (NTR). Electrons are drained from the network by reactive oxygen species (ROS) through thiol peroxidases, e.g., peroxiredoxins (PRX). Mathematical modeling promises access to quantitative understanding of the network function and was implemented for the photosynthesizing chloroplast by using published kinetic parameters combined with fitting to known biochemical data. Two networks were assembled, namely the ferredoxin (FDX), FDX-dependent TRX reductase (FTR), TRX, fructose-1,6-bisphosphatase pathway with 2-cysteine PRX/ROS as oxidant, and separately the FDX, FDX-dependent NADP reductase (FNR), NADPH, NTRC-pathway for 2-CysPRX reduction. Combining both modules allowed drawing several important conclusions of network performance. The resting H2O2 concentration was estimated to be about 30 nM in the chloroplast stroma. The electron flow to metabolism exceeds that into thiol regulation of FBPase more than 7000-fold under physiological conditions. The electron flow from NTRC to 2-CysPRX is about 5.46-times more efficient than that from TRX-f1 to 2-CysPRX. Under severe stress (30 μM H2O2) the ratio of electron flow to the thiol network relative to metabolism sinks to 1:251 whereas the ratio of electron flow from NTRC to 2-CysPRX and TRX-f1 to 2-CysPRX rises up to 1:80. Thus, the simulation provides clues on experimentally inaccessible parameters and describes the functional state of the chloroplast thiol regulatory network.Authors summaryThe state of the thiol redox regulatory network is a fundamental feature of all cells and determines metabolic and developmental processes. However, only some parameters are quantifiable in experiments. This paper establishes partial mathematical models which enable simulation of electron flows through the regulatory system. This in turn allows for estimating rates and states of components of the network and to tentatively address previously unknown parameters such as the resting hydrogen peroxide levels or the expenditure of reductive power for regulation relative to metabolism. The establishment of such models for simulating the performance and dynamics of the redox regulatory network is of significance not only for photosynthesis but also, e.g., in bacterial and animal cells exposed to environmental stress or pathological disorders.