Multiparticle computer simulation of protein interactions in the photosynthetic membrane

The basic principles of the design of direct multiparticle models and the results of multiparticle computer simulation of electron transfer by mobile protein carriers in the photosynthetic membrane of a chloroplast thylakoid are presented. The reactions of complex formation of the plastocyanin with cytochrome f and the pigment-protein complex of photosystem I, as well as of ferredoxin with FNR and photosystem I are considered. The regulatory role of diffusion and electrostatic interactions as well as the effect of the shape of the reaction volume and ionic strength on the rate of electron transport are discussed.     

Modeling of primary photosynthetic processes using the kinetic Monte Carlo method

Processes that occur in the ensemble of photosynthetic electron transport systems have been modeled using a kinetic Monte Carlo method. The size of a simulated ensemble (3–5 million elementary photosynthetic chains) corresponds to the number of photosynthetic reaction centers in a plant cell. The method enables one to modify the structure of a model system according to different concepts of the organization of processes in a photosynthetic membrane. Using this model, the experimental kinetics of the chlorophyll fluorescence induction associated with the Photosystem II and the redox transformations of a photoactive pigment of the Photosystem I have been successfully reproduced. The model was verified by comparing the calculated fluorescence induction curves to experimental curves, obtained in the presence of various photosynthesis inhibitors and under temperature inactivation of the Photosystem II donor side.     

Multiparticle computer simulation of plastocyanin diffusion and interaction with cytochrome <Emphasis Type="Italic">f</Emphasis> in the electrostatic field of the thylakoid membrane

A multiparticle computer model of plastocyanin-cytochrome f complex formation in the thylakoid lumen has been designed, which takes into account the electrostatic interactions of proteins and membrane. The Poisson-Boltzmann formalism was used to determine the electrostatic potentials of the electron carrier proteins and the thylakoid membrane at different ionic strengths. The membrane electrostatic field was shown to influence plastocyanin diffusion and interaction with cytochrome f. The rate constants for plastocyanin-cytochrome f complex formation were calculated as a function of ionic strength and membrane surface charge.     

Modelling of pattern formation and oscillations in pH and transmembrane potential near the cell membrane of Chara corallina

A mathematical model of potencial-dependent proton transfer across the membrane of Chara corallina cells is considered. To construct the model, partial differential equations describing the system dynamics in time and in space were used. The variables of the model are the proton concentration and membrane potential. The model describes the experimentally observed inhomogeneous distribution of transmembrane potential and pH along the membrane and oscillations of the potential and pH in time. A mechanism of the distribution of pH and membrane potential along the Chara corallina cell is suggested.     

Modeling of hysteresis in pH pattern formation along the cell membrane of algae Chara corallina

It is known that illumination of the algae Chara corallina results in the formation along the membrane of regions with inhomogeneous distribution of pH. It was shown that, in a particular range of illumination intensities, two states with different pH distribution are realized at one and the same value of light intensity: an entirely homogeneous state and completely formed structures (pattern). The transition from the homogeneous state to the pattern formation takes place at one value of light intensity, and the back transition, at another light intensity, i.e., the hysteresis is observed. This phenomenon was studied by mathematical modeling. The mechanism of hysteresis is discussed.     

The role of electrostatic interactions in the formation of ferredoxin–ferredoxin NADP<Superscript>+</Superscript> reductase and ferredoxin–hydrogenase complexes

A competitive Brownian model for the interaction of ferredoxin, ferredoxin NADP⁺ reductase and hydrogenase has been built. In the model, molecules of three types of proteins are placed into a cubic reaction volume, where they move under Brownian and electrostatic forces created by neighboring molecules and the solution. It has been shown that the rate of ferredoxin binding with ferredoxin NADP⁺ reductase does not change at the pH range from 5.0 to 9.0. Thus, it may be suggested that regulation of ferredoxin NADP⁺ reductase activity is mediated by other processes. On the other hand, the rate of ferredoxin binding with hydrogenase in the model depends greatly on pH: if the pH value increases from 6.0 to 8.0 the rate increases by factor of three. The increase of the pH value in the stroma under illumination results in an increase of the rate of its interaction with ferredoxin, but decreases the level of protons that are the substrate for the reaction catalyzed by the protein. Thus, the rate of hydrogen production in the chloroplast stroma is low at low pH due to the reception of a small number of electrons by hydrogenase. When the pH increases, the number of electrons that are received by the enzyme from ferredoxin also increases; thus, the rate of hydrogen production increases as well.     

Modeling chlorophyll a fluorescence transient: Relation to photosynthesis

To honor Academician Alexander Abramovitch Krasnovsky, we present here an educational review on the relation of chlorophyll a fluorescence transient to various processes in photosynthesis. The initial event in oxygenic photosynthesis is light absorption by chlorophylls (Chls), carotenoids, and, in some cases, phycobilins; these pigments form the antenna. Most of the energy is transferred to reaction centers where it is used for charge separation. The small part of energy that is not used in photochemistry is dissipated as heat or re-emitted as fluorescence. When a photosynthetic sample is transferred from dark to light, Chl a fluorescence (ChlF) intensity shows characteristic changes in time called fluorescence transient, the OJIPSMT transient, where O (the origin) is for the first measured minimum fluorescence level; J and I for intermediate inflections; P for peak; S for semi-steady state level; M for maximum; and T for terminal steady state level. This transient is a real signature of photosynthesis, since diverse events can be related to it, such as: changes in redox states of components of the linear electron transport flow, involvement of alternative electron routes, the build-up of a transmembrane pH gradient and membrane potential, activation of different nonphotochemical quenching processes, activation of the Calvin-Benson cycle, and other processes. In this review, we present our views on how different segments of the OJIPSMT transient are influenced by various photosynthetic processes, and discuss a number of studies involving mathematical modeling and simulation of the ChlF transient. A special emphasis is given to the slower PSMT phase, for which many studies have been recently published, but they are less known than on the faster OJIP phase.     

A model of photosystem II for the analysis of fast fluorescence rise in plant leaves

The polyphasic patterns of fluorescence induction rise in pea leaves in vivo and after the treatment with ionophores have been studied using a Plant Efficiency Analyzer. To analyze in detail photosystem II (PS II) electron transfer processes, an extended PS II model was applied, which included the sums of exponential functions to specify explicitly the light-driven formation of the transmembrane electric potential (ΔΨ(t)) as well as pH in the lumen (pH_L(t)) and stroma (pH_S(t)). PS II model parameters and numerical coefficients in ΔΨ(t), pH_L(t), and pH_S(t) were evaluated to fit fluorescence induction data for different experimental conditions: leaf in vivo or after ionophore treatment at low or high light intensity. The model imitated changes in the pattern of fluorescence induction rise due to the elimination of transmembrane potential in the presence of ionophores, when ΔΨ = 0 and pH_L(t), pH_S(t) changed to small extent relative to control values in vivo, with maximum ΔΨ(t) ∼ 90 mV and ΔΨ(t) ∼ 40 mV for the stationary state at ΔpH ≅ 1.8. As the light intensity was increased from 300 to 1200 μmol m⁻² s⁻¹, the heat dissipation rate constants increased threefold for nonradiative recombination of P680⁺Phe⁻ and by ∼30% for P680⁺Q_A⁻. The parameters ΔΨ, pH_S and pH_L were analyzed as factors of PS II redox state populations and fluorescence yield. The kinetic mechanism of fluorescence quenching is discussed, which is related with light-induced lumen acidification, when ⁺Q_A⁻ and P680⁺ recombination probability increases to regulate the Q_A reduction.