Reaction schemes, escape times and geminate recombinations in particle-based spatial simulations of biochemical reactions

Enzyme kinetics simulation at the scale of individual particles

Enzyme-catalyzed reactions involve two distinct timescales: a short timescale on which enzymes bind to substrate molecules to produce bound complexes and a comparatively long timescale on which the molecules of the complex are transformed into products. The uptake of the substrate in these reactions is the rate at which the product is made on the long timescale. Models often only consider the uptake to reduce the number of chemical species that need to be modeled and to avoid explicitly treating multiple timescales. Typically, the uptake rates cannot be described by mass action kinetics and are traditionally derived by applying singular perturbation theory to the system's governing differential equations. This analysis ignores short timescales by assuming that a pseudo-equilibrium between the enzyme and the enzyme-bound complex is maintained at all times. This assumption cannot be incorporated into current particle-based simulations of reaction-diffusion systems because they utilize proximity-based conditions to govern the instances of reactions that cannot maintain this pseudo-equilibrium for infinitely fast reactions. Instead, these methods must directly simulate the dynamics on the short timescale to accurately model the system. Due to the disparate timescales, such simulations require excessive amounts of computational time before the behavior on the long timescale can be observed. To resolve this problem, we use singular perturbation theory to develop a proximity-based reaction condition that enables us to ignore all fast reactions and directly reproduce non-mass action kinetics at long timescales. To demonstrate our approach, we implement simulations of a specific third order reaction with kinetics reminiscent of the prototypical Michaelis-Menten system.

Diffusion-limited attachment of large spherical particles to flexible membrane-immobilized receptors

Relatively large (~100 nm) spherical particles, e.g., virions, vesicles, or metal nanoparticles, often interact with short (<10 nm) flexible receptors immobilized in a lipid membrane or on other biologically relevant surfaces. The attachment kinetics of such particles may be limited globally by their diffusion toward a membrane or locally by diffusion around receptors. The detachment kinetics, also, can be limited by diffusion. Focusing on local diffusion limitations and using suitable approximations, we present expressions for the corresponding rate constants and identify their dependence on particle size and receptor length. We also illustrate features likely to be observed in such kinetics for particles (e.g., vesicles) with a substantial size distribution. The results obtained are generic and can be used to interpret a variety of situations. For example, we estimate upper values of virion attachment rate constants and clarify the likely effect of vesicle size distribution on previously observed non-exponential kinetics of vesicle detachment.

Cellular and molecular structure as a unifying framework for whole-cell modeling

Whole-cell modeling has the potential to play a major role in revolutionizing our understanding of cellular biology over the next few decades. A computational model of the entire cell would allow cellular biologists to integrate data from many disparate sources in a single consistent framework. Such a comprehensive model would be useful both for hypothesis testing and in the discovery of new behaviors that emerge from complex biological networks. Cellular and molecular structure can and should be a key organizing principle in a whole-cell model, connecting models across time and length scales in a multiscale approach. Here I present a summary of recent research centered around using molecular and cellular structure to model the behavior of cells.

Spatial simulations in systems biology: from molecules to cells

Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.