TY - JOUR A1 - Straube, Arthur A1 - Winkelmann, Stefanie A1 - Schütte, Christof A1 - Höfling, Felix T1 - Stochastic pH oscillations in a model of the urea–urease reaction confined to lipid vesicles JF - J. Phys. Chem. Lett. N2 - The urea-urease clock reaction is a pH switch from acid to basic that can turn into a pH oscillator if it occurs inside a suitable open reactor. We numerically study the confinement of the reaction to lipid vesicles, which permit the exchange with an external reservoir by differential transport, enabling the recovery of the pH level and yielding a constant supply of urea molecules. For microscopically small vesicles, the discreteness of the number of molecules requires a stochastic treatment of the reaction dynamics. Our analysis shows that intrinsic noise induces a significant statistical variation of the oscillation period, which increases as the vesicles become smaller. The mean period, however, is found to be remarkably robust for vesicle sizes down to approximately 200 nm, but the periodicity of the rhythm is gradually destroyed for smaller vesicles. The observed oscillations are explained as a canard-like limit cycle that differs from the wide class of conventional feedback oscillators. Y1 - 2021 U6 - https://doi.org/10.1021/acs.jpclett.1c03016 VL - 12 SP - 9888 EP - 9893 ER - TY - JOUR A1 - Ernst, Ariane A1 - Unger, Nathalie A1 - Schütte, Christof A1 - Walter, Alexander A1 - Winkelmann, Stefanie T1 - Rate-limiting recovery processes in neurotransmission under sustained stimulation JF - Mathematical Biosciences N2 - At chemical synapses, an arriving electric signal induces the fusion of vesicles with the presynaptic membrane, thereby releasing neurotransmitters into the synaptic cleft. After a fusion event, both the release site and the vesicle undergo a recovery process before becoming available for reuse again. Of central interest is the question which of the two restoration steps acts as the limiting factor during neurotrans-mission under high-frequency sustained stimulation. In order to investigate this question, we introduce a novel non-linear reaction network which involves explicit recovery steps for both the vesicles and the release sites, and includes the induced time-dependent output current. The associated reaction dynamics are formulated by means of ordinary differential equations (ODEs), as well as via the associated stochastic jump process. While the stochastic jump model describes a single release site, the average over many release sites is close to the ODE solution and shares its periodic structure. The reason for this can be traced back to the insight that recovery dynamics of vesicles and release sites are statistically almost independent. A sensitivity analysis on the recovery rates based on the ODE formulation reveals that neither the vesicle nor the release site recovery step can be identified as the essential rate-limiting step but that the rate- limiting feature changes over the course of stimulation. Under sustained stimulation the dynamics given by the ODEs exhibit transient dynamics leading from an initial depression of the postsynaptic response to an asymptotic periodic orbit, while the individual trajectories of the stochastic jump model lack the oscillatory behavior an asymptotic periodicity of the ODE-solution. Y1 - 2023 U6 - https://doi.org/10.1016/j.mbs.2023.109023 VL - 362 ER - TY - JOUR A1 - Straube, Arthur A1 - Winkelmann, Stefanie A1 - Höfling, Felix T1 - Accurate reduced models for the pH oscillations in the urea-urease reaction confined to giant lipid vesicles JF - The Journal of Physical Chemistry B N2 - This theoretical study concerns a pH oscillator based on the urea-urease reaction confined to giant lipid vesicles. Under suitable conditions, differential transport of urea and hydrogen ion across the unilamellar vesicle membrane periodically resets the pH clock that switches the system from acid to basic, resulting in self-sustained oscillations. We analyse the structure of the phase flow and of the limit cycle, which controls the dynamics for giant vesicles and dominates the pronouncedly stochastic oscillations in small vesicles of submicrometer size. To this end, we derive reduced models, which are amenable to analytic treatments that are complemented by numerical solutions, and obtain the period and amplitude of the oscillations as well as the parameter domain, where oscillatory behavior persists. We show that the accuracy of these predictions is highly sensitive to the employed reduction scheme. In particular, we suggest an accurate two-variable model and show its equivalence to a three-variable model that admits an interpretation in terms of a chemical reaction network. The faithful modeling of a single pH oscillator appears crucial for rationalizing experiments and understanding communication of vesicles and synchronization of rhythms. Y1 - 2023 U6 - https://doi.org/10.1021/acs.jpcb.2c09092 VL - 127 IS - 13 SP - 2955 EP - 2967 ER - TY - GEN A1 - Straube, Arthur A1 - Winkelmann, Stefanie A1 - Höfling, Felix T1 - Accurate reduced models for the pH oscillations in the urea-urease reaction confined to giant lipid vesicles N2 - Our theoretical study concerns an urea-urease-based pH oscillator confined to giant lipid vesicles. Under suitable conditions, differential transport of urea and hydrogen ion across the unilamellar vesicle membrane periodically resets the pH clock that switches the system from acid to basic, resulting in self-sustained oscillations. We analyse the structure of the limit cycle, which controls the dynamics for giant vesicles and dominates the strongly stochastic oscillations in small vesicles of submicrometer size. To this end, we derive reduced models, amenable to analytic treatments, and show that the accuracy of predictions, including the period of oscillations, is highly sensitive to the choice of the reduction scheme. In particular, we suggest an accurate two-variable model and show its equivalence to a three-variable model that admits an interpretation in terms of a chemical reaction network. The accurate description of a single pH oscillator appears crucial for rationalizing experiments and understanding communication of vesicles and synchronization of rhythms. T3 - ZIB-Report - 22-21 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88179 SN - 1438-0064 ER - TY - JOUR A1 - del Razo, Mauricio A1 - Winkelmann, Stefanie A1 - Klein, Rupert A1 - Höfling, Felix T1 - Chemical diffusion master equation: formulations of reaction-diffusion processes on the molecular level JF - Journal of Mathematical Physics N2 - The chemical diffusion master equation (CDME) describes the probabilistic dynamics of reaction--diffusion systems at the molecular level [del Razo et al., Lett. Math. Phys. 112:49, 2022]; it can be considered the master equation for reaction--diffusion processes. The CDME consists of an infinite ordered family of Fokker--Planck equations, where each level of the ordered family corresponds to a certain number of particles and each particle represents a molecule. The equations at each level describe the spatial diffusion of the corresponding set of particles, and they are coupled to each other via reaction operators --linear operators representing chemical reactions. These operators change the number of particles in the system, and thus transport probability between different levels in the family. In this work, we present three approaches to formulate the CDME and show the relations between them. We further deduce the non-trivial combinatorial factors contained in the reaction operators, and we elucidate the relation to the original formulation of the CDME, which is based on creation and annihilation operators acting on many-particle probability density functions. Finally we discuss applications to multiscale simulations of biochemical systems among other future prospects. Y1 - 2023 U6 - https://doi.org/10.1063/5.0129620 VL - 64 IS - 1 ER - TY - JOUR A1 - Lücke, Marvin A1 - Heitzig, Jobst A1 - Koltai, Péter A1 - Molkethin, Nora A1 - Winkelmann, Stefanie T1 - Large population limits of Markov processes on random networks JF - Stochastic Processes and their Applications N2 - We consider time-continuous Markovian discrete-state dynamics on random networks of interacting agents and study the large population limit. The dynamics are projected onto low-dimensional collective variables given by the shares of each discrete state in the system, or in certain subsystems, and general conditions for the convergence of the collective variable dynamics to a mean-field ordinary differential equation are proved. We discuss the convergence to this mean-field limit for a continuous-time noisy version of the so-called "voter model" on Erdős-Rényi random graphs, on the stochastic block model, as well as on random regular graphs. Moreover, a heterogeneous population of agents is studied. For each of these types of interaction networks, we specify the convergence conditions in dependency on the corresponding model parameters. Y1 - 2023 U6 - https://doi.org/10.1016/j.spa.2023.09.007 VL - 166 ER - TY - JOUR A1 - Steudle, Gesine A1 - Winkelmann, Stefanie A1 - Fürst, Steffen A1 - Wolf, Sarah T1 - Understanding Memory Mechanisms in Socio-Technical Systems: the Case of an Agent-based Mobility Model JF - Advances in Complex Systems N2 - This paper explores memory mechanisms in complex socio-technical systems, using a mobility demand model as an example case. We simplified a large-scale agent-based mobility model into a Markov process and discover that the mobility decision process is non-Markovian. This is due to its dependence on the system’s history, including social structure and local infrastructure, which evolve based on prior mobility decisions. To make the process Markovian, we extend the state space by incorporating two history-dependent components. Although our model is a very much reduced version of the original one, it remains too complex for the application of usual analytic methods. Instead, we employ simulations to examine the functionalities of the two history-dependent components. We think that the structure of the analyzed stochastic process is exemplary for many socio-technical, -economic, -ecological systems. Additionally, it exhibits analogies with the framework of extended evolution, which has previously been used to study cultural evolution. Y1 - 2024 U6 - https://doi.org/10.1142/S0219525924400034 VL - 27 ER - TY - JOUR A1 - Boltz, Horst-Holger A1 - Sirbu, Alexei A1 - Stelzer, Nina A1 - de Lanerolle, Primal A1 - Winkelmann, Stefanie A1 - Annibale, Paolo T1 - The Impact of Membrane Protein Diffusion on GPCR Signaling JF - Cells N2 - Spatiotemporal signal shaping in G protein-coupled receptor (GPCR) signaling is now a well-established and accepted notion to explain how signaling specificity can be achieved by a superfamily sharing only a handful of downstream second messengers. Dozens of Gs-coupled GPCR signals ultimately converge on the production of cAMP, a ubiquitous second messenger. This idea is almost always framed in terms of local concentrations, the differences in which are maintained by means of spatial separation. However, given the dynamic nature of the reaction-diffusion processes at hand, the dynamics, in particular the local diffusional properties of the receptors and their cognate G proteins, are also important. By combining some first principle considerations, simulated data, and experimental data of the receptors diffusing on the membranes of living cells, we offer a short perspective on the modulatory role of local membrane diffusion in regulating GPCR-mediated cell signaling. Our analysis points to a diffusion-limited regime where the effective production rate of activated G protein scales linearly with the receptor–G protein complex’s relative diffusion rate and to an interesting role played by the membrane geometry in modulating the efficiency of coupling Y1 - 2022 U6 - https://doi.org/10.3390/cells11101660 VL - 11 IS - 10 SP - 1660 ER - TY - JOUR A1 - Engel, Maximilian A1 - Olicón-Méndez, Guillermo A1 - Wehlitz, Nathalie A1 - Winkelmann, Stefanie T1 - Synchronization and random attractors in reaction jump processes JF - Journal of Dynamics and Differential Equations N2 - This work explores a synchronization-like phenomenon induced by common noise for continuous-time Markov jump processes given by chemical reaction networks. Based on Gillespie’s stochastic simulation algorithm, a corresponding random dynamical system is formulated in a two-step procedure, at first for the states of the embedded discrete-time Markov chain and then for the augmented Markov chain including random jump times. We uncover a time-shifted synchronization in the sense that—after some initial waiting time—one trajectory exactly replicates another one with a certain time delay. Whether or not such a synchronization behavior occurs depends on the combination of the initial states. We prove this partial time-shifted synchronization for the special setting of a birth-death process by analyzing the corresponding two-point motion of the embedded Markov chain and determine the structure of the associated random attractor. In this context, we also provide general results on existence and form of random attractors for discrete-time, discrete-space random dynamical systems. Y1 - 2024 U6 - https://doi.org/10.1007/s10884-023-10345-4 ER - TY - JOUR A1 - Lücke, Marvin A1 - Winkelmann, Stefanie A1 - Koltai, Peter T1 - Accurate mean-field equation for voter model dynamics on scale-free networks JF - Physical Review E N2 - Understanding the emergent macroscopic behavior of dynamical systems on networks is a crucial but challenging task. One of the simplest and most effective methods to construct a reduced macroscopic model is given by mean-field theory. The resulting approximations perform well on dense and homogeneous networks but poorly on scale-free networks, which, however, are more realistic in many applications. In this paper, we introduce a modified version of the mean-field approximation for voter model dynamics on scale-free networks. The two main deviations from classical theory are that we use degree-weighted shares as coarse variables and that we introduce a correlation factor that can be interpreted as slowing down dynamics induced by interactions. We observe that the correlation factor is only a property of the network and not of the state or of parameters of the process. This approach achieves a significantly smaller approximation error than standard methods without increasing dimensionality. Y1 - 2026 U6 - https://doi.org/10.1103/vkpx-5cvt VL - 113 ER - TY - JOUR A1 - Wehlitz, Nathalie A1 - Pavliotis, Grigorios A1 - Schütte, Christof A1 - Winkelmann, Stefanie T1 - Data-driven Reduction of Transfer Operators for Particle Clustering Dynamics N2 - We develop an operator-based framework to coarse-grain interacting particle systems that exhibit clustering dynamics. Starting from the particle-based transfer operator, we first construct a sequence of reduced representations: the operator is projected onto concentrations and then further reduced by representing the concentration dynamics on a geometric low-dimensional manifold and an adapted finite-state discretization. The resulting coarse-grained transfer operator is finally estimated from dynamical simulation data by inferring the transition probabilities between the Markov states. Applied to systems with multichromatic and Morse interaction potentials, the reduced model reproduces key features of the clustering process, including transitions between cluster configurations and the emergence of metastable states. Spectral analysis and transition-path analysis of the estimated operator reveal implied time scales and dominant transition pathways, providing an interpretable and efficient description of particle-clustering dynamics. Y1 - 2026 ER -