Protons in small spaces: Discrete simulations of vesicle acidification

Abstract

The lumenal pH of an organelle is one of its defining characteristics and central to its biological function. Experiments have elucidated many of the key pH regulatory elements and how they vary from compartment-to-compartment, and continuum mathematical models have played an important role in understanding how these elements (proton pumps, counter-ion fluxes, membrane potential, buffering capacity, etc.) work together to achieve specific pH setpoints. While continuum models have proven successful in describing ion regulation at the cellular length scale, it is unknown if they are valid at the subcellular level where volumes are small, ion numbers may fluctuate wildly, and biochemical heterogeneity is large. Here, we create a discrete, stochastic (DS) model of vesicular acidification to answer this question. We used this simplified model to analyze pH measurements of isolated vesicles containing single proton pumps and compared these results to solutions from a continuum, ordinary differential equations (ODE)-based model. Both models predict similar parameter estimates for the mean proton pumping rate, membrane permeability, etc., but, as expected, the ODE model fails to report on the fluctuations in the system. The stochastic model predicts that pH fluctuations decrease during acidification, but noise analysis of single-vesicle data confirms our finding that the experimental noise is dominated by the fluorescent dye, and it reveals no insight into the true noise in the proton fluctuations. Finally, we again use the reduced DS model explore the acidification of large, lysosome-like vesicles to determine how stochastic elements, such as variations in proton-pump copy number and cycling between on and off states, impact the pH setpoint and fluctuations around this setpoint. PLOS Computational Biology | https://doi.

Stimulus-Response signaling dynamics characterize macrophage polarization states

“…To describe differences between experimental conditions, summary statistics are easy to compute and interpret, however they are insufficient. For example, taking the timepoint by timepoint mean of the single-cell trajectories can obscure asynchronous oscillatory dynamics observed at the single-cell level 66 . Further, average behavior descriptions mask how the distributions of responses actually overlap.…”

Pump-Rest-Leak-Repeat: regulation of the mammalian-brain V-ATPase via ultra-slow mode-switching

“…1b (t = 5 min). The plateau reached after ~15 min reflects a dynamic equilibrium between protons being actively pumped into the lumen and protons passively permeating through the membrane down the concentration gradient 19,25 . Addition of the V-ATPase specific inhibitor bafilomycin (Fig.…”

Collective Behavior of Urease pH Clocks in Nano- and Microvesicles Controlled by Fast Ammonia Transport

“…Hence, some insight can be gained from simulations with a simple ODE model of the vesicles and the external solution (SI 3.4). The reaction can be modeled by taking into account stochastic effects; however, in other work it was determined that population-level behavior was retained in the ODE models. , The rate of change of the concentration of a species A i in a vesicle was determined by the reaction and mass transfer rate as follows: − where f ( A i ) contains the enzyme reaction and solution equilibria terms (SI 3.1–3.2) taken from earlier work, ,, A 0 is the concentration of the species in the outer solution, P i is the permeability coefficient of species i, and r is the radius of the vesicle. In the outer solution, the rate of change of the concentration of each species, A 0 , was given by the reaction rate ( g ( A 0 )) and the mass transfer rate (for identical vesicles) as follows: where ϕ = NV j / V 0 = the vesicle volume fraction, which takes into account dilution as a result of the volume change from vesicle to solution.…”