Reactions, diffusion, and volume exclusion in a conserved system of interacting particles

From random walks on networks to nonlinear diffusion

Mathematical models of motility are often based on random-walk descriptions of discrete individuals that can move according to certain rules. It is usually the case that large masses concentrated in small regions of space have a great impact on the collective movement of the group. For this reason, many models in mathematical biology have incorporated crowding effects and managed to understand their implications. Here, we build on a previously developed framework for random walks on networks to show that in the continuum limit, the underlying stochastic process can be identified with a diffusion partial differential equation. The diffusion coefficient of the emerging equation is, in general, density dependent, and can be directly related to the transition probabilities of the random walk. Moreover, the relaxation time of the stochastic process is directly linked to the diffusion coefficient and also to the network structure, as it usually happens in the case of linear diffusion. As a specific example, we study the equivalent of a porous-medium-type equation on networks, which shows similar properties to its continuum equivalent. For this equation, self-similar solutions on a lattice and on homogeneous trees can be found, showing finite speed of propagation in contrast to commonly used linear diffusion equations. These findings also provide insights into reaction-diffusion systems with general diffusion operators, which have appeared recently in some applications.

One-dimensional model for chemotaxis with hard-core interactions

In this paper we consider a biased velocity jump process with excluded-volume interactions for chemotaxis, where we account for the size of each particle. Starting with a system of N individual hard rod particles in one dimension, we derive a nonlinear kinetic model using two different approaches. The first approach is a systematic derivation for small occupied fraction of particles based on the method of matched asymptotic expansions. The second approach, based on a compression method that exploits the single-file motion of hard core particles, does not have the limitation of a small occupied fraction but requires constant tumbling rates. We validate our nonlinear model with numerical simulations, comparing its solutions with the corresponding noninteracting linear model as well as stochastic simulations of the underlying particle system.

Particle-Based Simulation Reveals Macromolecular Crowding Effects on the Michaelis-Menten Mechanism

Many computational models for analyzing and predicting cell physiology rely on in vitro data collected in dilute and controlled buffer solutions. However, this can mislead models because up to 40% of the intracellular volume—depending on the organism, the physiology, and the cellular compartment—is occupied by a dense mixture of proteins, lipids, polysaccharides, RNA, and DNA. These intracellular macromolecules interfere with the interactions of enzymes and their reactants and thus affect the kinetics of biochemical reactions, making in vivo reactions considerably more complex than the in vitro data indicates. In this work, we present a new, to our knowledge, type of kinetics that captures and quantifies the effect of volume exclusion and other spatial phenomena on the kinetics of elementary reactions. We further developed a framework that allows for the efficient parameterization of these kinetics using particle simulations. Our formulation, entitled generalized elementary kinetics, can be used to analyze and predict the effect of intracellular crowding on enzymatic reactions and was herein applied to investigate the influence of crowding on phosphoglycerate mutase in Escherichia coli, which exhibits prototypical reversible Michaelis-Menten kinetics. Current research indicates that many enzymes are reaction limited and not diffusion limited, and our results suggest that the influence of fractal diffusion is minimal for these reaction-limited enzymes. Instead, increased association rates and decreased dissociation rates lead to a strong decrease in the effective maximal velocities V_max and the effective Michaelis-Menten constants K_M under physiologically relevant volume occupancies. Finally, the effects of crowding were explored in the context of a linear pathway, with the finding that crowding can have a redistributing effect on the effective flux responses in the case of twofold enzyme overexpression. We suggest that this framework, in combination with detailed kinetics models, will improve our understanding of enzyme reaction networks under nonideal conditions.