Three-state mathematical model for the assessment of DCMU-treated photosystem II heterogeneity.

Photosystem II (PSII) is one of the main pigment-protein complexes of photosynthesis which is highly sensitive to unfavorable environmental factors. The heterogeneity of PSII properties is essential for the resistance of autotrophic organisms to stress factors. Assessment of the PSII heterogeneity may be used in environmental monitoring for on-line detection of contamination of the environment. We propose an approach to assess PSII oxygen-evolving complex and light-harvesting antenna heterogeneity that is based on mathematical modeling of the shape of chlorophyll a fluorescence rise of 3-(3,4-dichlorophenyl)-1,1-dimethylurea-treated samples. The hierarchy of characteristic times of the processes considered in the model makes it possible to reduce the model to a system of three ordinary differential equations. The analytic solution of the reduced three-state model is expressed as a sum of two exponential functions, and it exactly reproduces the solution of the complete system within the time range from microseconds to hundreds of milliseconds. The combination of several such models for reaction centers with different properties made it possible to use it as an instrument to study PSII heterogeneity. PSII heterogeneity was studied for Chlamydomonas at different intensities of actinic light, for Scenedesmus under short-term heating, and for Chlorella grown in nitrate-enriched and nitrate-depleted media.

Mathematical Simulation of Electron Transport in the Primary Photosynthetic Processes.

Summarized results of investigation of regulation of electron transport and associated processes in the photosynthetic membrane using methods of mathematical and computer modeling carried out at the Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, are presented in this review. Detailed kinetic models of processes in the thylakoid membrane were developed using the apparatus of differential equations. Fitting of the model curves to the data of spectral measurements allowed us to estimate the values of parameters that were not determined directly in experiments. The probabilistic method of agent-based Monte Carlo modeling provides ample opportunities for studying dynamics of heterogeneous systems based on the rules for the behavior of individual elements of the system. Algorithms for simplified representation of Big Data make it possible to monitor changes in the photosynthetic apparatus in the course of culture growth in a photobioreactor and for the purpose of environmental monitoring. Brownian and molecular models describe movement and interaction of individual electron carrier proteins and make it possible to study electrostatic, hydrophobic, and other interactions leading to regulation of conformational changes in the reaction complexes. Direct multiparticle models explicitly simulate Brownian diffusion of the mobile protein carriers and their electrostatic interactions with multienzyme complexes both in solution and in heterogeneous interior of a biomembrane. The combined use of methods of kinetic and Brownian multiparticle and molecular modeling makes it possible to study the mechanisms of regulation of an integral system of electron transport processes in plants and algae at molecular and subcellular levels.

'Ecological photobiology' session of the Russian Photobiology Society 9th Congress (Shepsi, September 12-19, 2021).

'Ecological photobiology' session of the Russian Photobiology Society 9th Congress was devoted to a wide range of problems related to the assessment of the environmental state by photobiological methods and included oral presentations and a poster session. A short survey of these presentations is given.

Molecular, Brownian, kinetic and stochastic models of the processes in photosynthetic membrane of green plants and microalgae.

The paper presents the results of recent work at the Department of Biophysics of the Biological Faculty, Lomonosov Moscow State University on the kinetic and multiparticle modeling of processes in the photosynthetic membrane. The detailed kinetic models and the rule-based kinetic Monte Carlo models allow to reproduce the fluorescence induction curves and redox transformations of the photoactive pigment P700 in the time range from 100 ns to dozens of seconds and make it possible to reveal the role of individual carriers in their formation for different types of photosynthetic organisms under different illumination regimes, in the presence of inhibitors, under stress conditions. The fitting of the model curves to the experimental data quantifies the reaction rate constants that cannot be directly measured experimentally, including the non-radiative thermal relaxation reactions. We use the direct multiparticle models to explicitly describe the interactions of mobile photosynthetic carrier proteins with multienzyme complexes both in solution and in the biomembrane interior. An analysis of these models reveals the role of diffusion and electrostatic factors in the regulation of electron transport, the influence of ionic strength and pH of the cellular environment on the rate of electron transport reactions between carrier proteins. To describe the conformational intramolecular processes of formation of the final complex, in which the actual electron transfer occurs, we use the methods of molecular dynamics. The results obtained using kinetic and molecular models supplement our knowledge of the mechanisms of organization of the photosynthetic electron transport processes at the cellular and molecular levels.

Simulating the Interplay between the Uptake of Inorganic Phosphate and the Cell Phosphate Metabolism under Phosphorus Feast and Famine Conditions in Chlorella vulgaris.

Using a mathematical simulation approach, we studied the dynamics of the green microalga Chlorella vulgaris phosphate metabolism response to shortage and subsequent replenishing of inorganic phosphate in the medium. A three-pool interaction model was used to describe the phosphate uptake from the medium, its incorporation into the cell organic compounds, its storage in the form of polyphosphates, and culture growth. The model comprises a system of ordinary differential equations. The distribution of phosphorous between cell pools was examined for three different stages of the experiment: growth in phosphate-rich medium, incubation in phosphate-free medium, and phosphate addition to the phosphorus-starving culture. Mathematical modeling offers two possible scenarios for the appearance of the peak of polyphosphates (PolyP). The first scenario explains the accumulation of PolyP by activation of the processes of its synthesis, and the decline in PolyP is due to its redistribution between dividing cells during growth. The second scenario includes a hysteretic mechanism for the regulation of PolyP hydrolysis, depending on the intracellular content of inorganic phosphate. The new model of the dynamics of P pools in the cell allows one to better understand the phenomena taking place during P starvation and re-feeding of the P-starved microalgal cultures with inorganic phosphate such as transient PolyP accumulation. Biotechnological implications of the observed dynamics of the polyphosphate pool of the microalgal cell are considered. An approach enhancing the microalgae-based wastewater treatment method based on these scenarios is proposed.

Call for contributions to the Special Issue on the 9th Congress of the Russian Photobiological Society held in Shepsi, Krasnodar region, Russia, on September 12-19, 2021.

This Commentary describes a call for submissions for the upcoming Special Issue focused on the science presented at the 9th Congress of The Russian Photobiological Society entitled "Actual Problems in Photobiology."

Comparative analysis of plastocyanin-cytochrome f complex formation in higher plants, green algae and cyanobacteria.

Mechanisms of the complex formation between plastocyanin and cytochrome f in higher plants (Spinacia oleracea and Brassica rapa), green microalgae Chlamydomonas reinhardtii and two species of cyanobacteria (Phormidium laminosum and Nostoc sp.) were investigated using combined Brownian and molecular dynamics simulations and hierarchical cluster analysis. In higher plants and green algae, electrostatic interactions force plastocyanin molecule close to the heme of cytochrome f. In the subsequent rotation of plastocyanin molecule around the point of electrostatic contact in the vicinity of cytochrome f, copper (Cu) atom approaches cytochrome heme forming a stable configuration where cytochrome f molecule behaves as a rather rigid body without conformational changes. In Nostoc plastocyanin molecule approaches cytochrome f in a different orientation (head-on) where the stabilization of the plastocyanin-cytochrome f complex is accompanied by the conformational changes of the G188E189D190 loop that stabilizes the whole complex. In cyanobacterium P. laminosum, electrostatic preorientation of the approaching molecules was not detected, thus indicating that random motions rather than long-range electrostatic interactions are responsible for the proper mutual orientation. We demonstrated that despite the structural similarity of the investigated electron transport proteins in different photosynthetic organisms, the complexity of molecular mechanisms of the complex formation increases in the following sequence: non-heterocystous cyanobacteria - heterocystous cyanobacteria - green algae - flowering plants.

Multiparticle Brownian dynamics simulation of experimental kinetics of cytochrome bf oxidation and photosystem I reduction by plastocyanin.

A model of electron transport from cytochrome f to photosystem I mediated by plastocyanin was designed on the basis of the multiparticle Brownian dynamics method. The model combines events which occur over a wide time range, including protein diffusion along the thylakoid membrane, long-distance interactions between proteins, formation of a multiprotein complex, electron transfer within a complex and complex dissociation. Results of the modeling were compared with the experimental kinetics measured in chloroplast thylakoids. Computer simulation demonstrated that the complex interior of the photosynthetic membrane, electrostatic interactions and Brownian diffusion provide physical conditions for the directed electron flow along the photosynthetic electron transport chain.

The stoichiometry and energetics of oxygenic phototrophic growth.

The values of gross metabolic flows in cells are essentially interconnected due to conservation laws of chemical elements and interrelations of biochemical coupling. Therefore, the overall stoichiometry of cellular metabolism, such as the biomass quantum yield, the ratio between linear and circular flows via the electron transport chain, etc., can be calculated using balances of metabolic flows in the network branching points and coupling ratios related to ATP formation and expenditures. This work has studied the energetic stoichiometry of photosynthetic cells by considering the transfer of reductivity in the course of biochemical reactions. This approach yielded rigorous mathematical expressions for biomass quantum yield and other integral bioenergetic indices of cellular growth as functions of ATP balance parameters. The effect of cellular substance turnover has been taken into account. The obtained theoretical estimation of biomass quantum yield is rather close to experimental data which confirms the predictive capacity of this approach.

Investigation of the maximum quantum yield of PS II in Haematococcus pluvialis cell cultures during growth: effects of chemical or high-intensity light treatment.

In this study, we investigated the increase in photosynthetic quantum yield that occurs in advance of increased microalgal growth. Haematococcus pluvialis was cultivated under normal conditions; the number of cells, the maximum quantum yield of photosystem II (F(v)/F(m)), and optical density were measured. We observed an increase in F(v)/F(m) approximately 72h prior to the cell growth phase. To confirm the relationship between photosynthetic yield and growth, samples were treated with several chemicals under high-intensity light illumination and control conditions to inhibit photosystem II and induce a decrease in the quantum photosynthetic yield. The samples were exposed to high-intensity light at an irradiance of 400µmol photonsm(-2)s(-1) for varied amount of time and were treated with chemicals such as 3-(3,4-dichlorophenyl)-1,1-dimethylurea, nigericin sodium salt and valinomycin. We observed that both the photooxidation of photosystem II reaction centers and the formation of transmembrane electrochemical gradients led to an initial decrease in fluorescence yield after the onset of high-intensity light illumination. We also observed that treatment of high-intensity light illuminated cells with antibiotics after adaptation to moderate light intensities caused a difference in photosynthetic activity. In conclusion, the maximum quantum yield of photosystem II is obtained prior to the cell growth phase and can therefore be used as a prediction parameter for cell growth.

Computer simulation of interaction of photosystem 1 with plastocyanin and ferredoxin.

We designed 3D multiparticle computer models to simulate diffusion and interactions of spinach plastocyanin and ferredoxin with plant photosystem 1 in a solution. Using these models we studied kinetic characteristics of plastocyanin-photosystem 1 and ferredoxin-photosystem 1 complex formation at a variety of ionic strength values. The computer multiparticle models demonstrate non-monotonic dependences of complex formation rates on the ionic strength as the result of long-range electrostatic interactions. Our calculations show that the decrease in the association second-order rate constant at low values of the ionic strength is caused by the protein pairs spending more time in "wrong" orientations which do not satisfy the docking conditions and so do not form the final complex capable of the electron transfer.

New direct dynamic models of protein interactions coupled to photosynthetic electron transport reactions.

This review covers the methods of computer simulation of protein interactions taking part in photosynthetic electron transport reactions. A direct multiparticle simulation method that simulates reactions describing interactions of ensembles of molecules in the heterogeneous interior of a cell is developed. In the models, protein molecules move according to the laws of Brownian dynamics, mutually orient themselves in the electrical field, and form complexes in the 3D scene. The method allows us to visualize the processes of molecule interactions and to calculate the rate constants for protein complex formation reactions in the solution and in the photosynthetic membrane. Three-dimensional multiparticle computer models for simulating the complex formation kinetics for plastocyanin with photosystem I and cytochrome bf complex, and ferredoxin with photosystem I and ferredoxin:NADP+-reductase are considered. Effects of ionic strength are featured for wild type and mutant proteins. The computer multiparticle models describe nonmonotonic dependences of complex formation rates on the ionic strength as the result of long-range electrostatic interactions.