Ronny Straube

Modeling the electron transport chain of purple non‐sulfur bacteria

Purple non‐sulfur bacteria (Rhodospirillaceae) have been extensively employed for studying principles of photosynthetic and respiratory electron transport phosphorylation and for investigating the regulation of gene expression in response to redox signals. Here, we use mathematical modeling to evaluate the steady‐state behavior of the electron transport chain (ETC) in these bacteria under different environmental conditions. Elementary‐modes analysis of a stoichiometric ETC model reveals nine operational modes. Most of them represent well‐known functional states, however, two modes constitute reverse electron flow under respiratory conditions, which has been barely considered so far. We further present and analyze a kinetic model of the ETC in which rate laws of electron transfer steps are based on redox potential differences. Our model reproduces well‐known phenomena of respiratory and photosynthetic operation of the ETC and also provides non‐intuitive predictions. As one key result, model simulations demonstrate a stronger reduction of ubiquinone when switching from high‐light to low‐light conditions. This result is parameter insensitive and supports the hypothesis that the redox state of ubiquinone is a suitable signal for controlling photosynthetic gene expression.

AN ASYMPTOTIC ANALYSIS OF THE MEAN FIRST PASSAGE TIME FOR NARROW ESCAPE PROBLEMS: PART II: THE SPHERE ∗

The mean first passage time (MFPT) is calculated for a Brownian particle in a spherical domain in

Diffusion on a Sphere with Localized Traps : Mean First Passage Time, Eigenvalue Asymptotics, and Fekete Points

\mathbb{R}^3

Spatial Desynchronization of Glycolytic Waves as Revealed by Karhunen−Loeve Analysis

that contains N small nonoverlapping absorbing windows, or traps, on its boundary. For the unit sphere, the method of matched asymptotic expansions is used to derive an explicit three-term asymptotic expansion for the MFPT for the case of N small locally circular absorbing windows. The third term in this expansion, not previously calculated, depends explicitly on the spatial configuration of the absorbing windows on the boundary of the sphere. The three-term asymptotic expansion for the average MFPT is shown to be in very close agreement with full numerical results. The average MFPT is shown to be minimized for trap configurations that minimize a certain discrete variational problem. This variational problem is closely related to the well-known optimization problem of determining the minimum energy configuration for N repelling point charges on the unit sphere. Numerical results, based on globa...

Reciprocal Enzyme Regulation as a Source of Bistability in Covalent Modification Cycles

A common scenario in cellular signal transduction is that a diffusing surface-bound molecule must arrive at a localized signaling region on the cell membrane before the signaling cascade can be completed. The question then arises of how quickly such signaling molecules can arrive at newly formed signaling regions. Here, we attack this problem by calculating asymptotic results for the mean first passage time for a diffusing particle confined to the surface of a sphere, in the presence of N partially absorbing traps of small radii. The rate at which the small diffusing molecule becomes captured by one of the traps is determined by asymptotically calculating the principal eigenvalue for the Laplace operator on the sphere with small localized traps. The asymptotic analysis relies on the method of matched asymptotic expansions, together with detailed properties of the Greens function for the Laplacian and the Helmholtz operators on the surface of the unit sphere. The asymptotic results compare favorably with ...

Sensitivity and Robustness in Covalent Modification Cycles with a Bifunctional Converter Enzyme

The dynamics of glycolytic waves in a yeast extract have been investigated in an open spatial reactor. At low protein contents in the extract, we find a transition from inwardly moving target patterns at the beginning of the experiment to outwardly moving spiral- or circular-shaped waves at later stages. These two phases are separated by a transition phase of more complex spatiotemporal dynamics. We have analyzed the pattern dynamics in these three intervals at different spatial scales by means of a Karhunen-Loeve (KL) decomposition. During the initial phase of the experiment, the observed patterns are sufficiently described by the two dominant KL modes independently of the spatial scale. However, during the last stage of the experiment, at least 6 KL modes are needed to account for the observed patterns at spatial scales larger than 3 mm, while for smaller scales, 2 KL modes are still sufficient. This indicates that in the course of the experiment, the local glycolytic oscillators become desynchronized at spatial scales larger than 3 mm. Possible reasons for the desynchronization of the glycolytic waves are discussed.

Reciprocal Regulation as a Source of Ultrasensitivity in Two-Component Systems with a Bifunctional Sensor Kinase

Covalent modification cycles (CMCs) are the building blocks of many regulatory networks in biological systems. Under proper kinetic conditions such mono-cyclic enzyme systems can show a higher sensitivity to effectors than enzymes subject to direct allosteric regulation. Using methods from reaction network theory it has been argued that CMCs can potentially exhibit multiple steady states if the converter enzymes are regulated in a reciprocal manner, but the underlying mechanism as well as the kinetic requirements for the emergence of such a behavior remained unclear. Here, we reinvestigate CMCs with reciprocal regulation of the converter enzymes for two common regulatory mechanisms: allosteric regulation and covalent modification. To analyze the steady state behavior of the corresponding mass-action equations, we derive reduced models by means of a quasi-steady state approximation (QSSA). We also derive reduced models using the total QSSA which often better reproduces the transient dynamics of enzyme-catalyzed reaction systems. Through a steady state analysis of the reduced models we show that the occurrence of bistability can be associated with the presence of a double negative feedback loop. We also derive constraints for the model parameters which might help to evaluate the potential significance of the mechanisms described here for the generation of bistability in natural systems. In particular, our results support the view of a possible bistable response in the metabolic PFK1/F1,6BPase cycle as observed experimentally in rat liver extracts, and it suggests an alternative view on the origin of bistability in the Cdk1-Wee1-Cdc25 system.

Modeling the Light- and Redox-Dependent Interaction of PpsR/AppA in Rhodobacter sphaeroides

Regulation by covalent modification is a common mechanism to transmit signals in biological systems. The modifying reactions are catalyzed either by two distinct converter enzymes or by a single bifunctional enzyme (which may employ either one or two catalytic sites for its opposing activities). The reason for this diversification is unclear, but contemporary theoretical models predict that systems with distinct converter enzymes can exhibit enhanced sensitivity to input signals whereas bifunctional enzymes with two catalytic sites are believed to generate robustness against variations in systems components. However, experiments indicate that bifunctional enzymes can also exhibit enhanced sensitivity due to the zero-order effect, raising the question whether both phenomena could be understood within a common mechanistic model. Here, I argue that this is, indeed, the case. Specifically, I show that bifunctional enzymes with two catalytic sites can exhibit both ultrasensitivity and concentration robustness, depending on the kinetic operating regime of the enzymes opposing activities. The model predictions are discussed in the context of experimental observations of ultrasensitivity and concentration robustness in the uridylylation cycle of the PII protein, and in the phosphorylation cycle of the isocitrate dehydrogenase, respectively.

Inward Rotating Spiral Waves in Glycolysis

Two-component signal transduction systems, where the phosphorylation state of a regulator protein is modulated by a sensor kinase, are common in bacteria and other microbes. In many of these systems, the sensor kinase is bifunctional catalyzing both, the phosphorylation and the dephosphorylation of the regulator protein in response to input signals. Previous studies have shown that systems with a bifunctional enzyme can adjust the phosphorylation level of the regulator protein independently of the total protein concentrations – a property known as concentration robustness. Here, I argue that two-component systems with a bifunctional enzyme may also exhibit ultrasensitivity if the input signal reciprocally affects multiple activities of the sensor kinase. To this end, I consider the case where an allosteric effector inhibits autophosphorylation and, concomitantly, activates the enzymes phosphatase activity, as observed experimentally in the PhoQ/PhoP and NRII/NRI systems. A theoretical analysis reveals two operating regimes under steady state conditions depending on the effector affinity: If the affinity is low the system produces a graded response with respect to input signals and exhibits stimulus-dependent concentration robustness – consistent with previous experiments. In contrast, a high-affinity effector may generate ultrasensitivity by a similar mechanism as phosphorylation-dephosphorylation cycles with distinct converter enzymes. The occurrence of ultrasensitivity requires saturation of the sensor kinases phosphatase activity, but is restricted to low effector concentrations, which suggests that this mode of operation might be employed for the detection and amplification of low abundant input signals. Interestingly, the same mechanism also applies to covalent modification cycles with a bifunctional converter enzyme, which suggests that reciprocal regulation, as a mechanism to generate ultrasensitivity, is not restricted to two-component systems, but may apply more generally to bifunctional enzyme systems.

An asymptotic analysis of intracellular signaling gradients arising from multiple small compartments

Facultative photosynthetic bacteria switch their energy generation mechanism from respiration to photosynthesis depending on oxygen tension and light. Part of this transition is mediated by the aerobic transcriptional repressor PpsR. In Rhodobacter sphaeroides, the repressive action of PpsR is antagonized by the redox- and blue-light-sensitive flavoprotein AppA which results in a unique phenotype: the repression of photosynthesis genes at intermediate oxygen levels and high light intensity, which is believed to reduce the risk of photooxidative stress. To analyze the underlying mechanism we developed a simple mathematical model based on the AppA-dependent reduction of a disulfide bond in PpsR and the light-sensitive complex formation between the reduced forms of AppA and PpsR. A steady-state analysis shows that high light repression can indeed occur at intermediate oxygen levels if PpsR is reduced on a faster timescale than AppA and if the electron transfer from AppA to PpsR is effectively irreversible. The model further predicts that if AppA copy numbers exceed those of PpsR by at least a factor of two, the transition from aerobic to anaerobic growth mode can occur via a bistable regime. We provide necessary conditions for the emergence of bistability and discuss possible experimental verifications.