Communication: Clusters of absorbing disks on a reflecting wall: competition for diffusing particles

Permeability and diffusion resistance of porous membranes: Analytical theory and its numerical test

This study is devoted to the transport of neutral solutes through porous flat membranes, driven by the solute concentration difference in the reservoirs separated by the membrane. Transport occurs through membrane channels, which are assumed to be non-overlapping, identical, straight cylindrical pores connecting the reservoirs. The key quantities characterizing transport are membrane permeability and its diffusion resistance. Such transport problems arising in very different contexts, ranging from plant physiology and cell biology to chemical engineering, have been studied for more than a century. Nevertheless, an expression giving the permeability for a membrane of arbitrary thickness at arbitrary surface densities of the channel openings is still unknown. Here, we fill in the gap and derive such an expression. Since this expression is approximate, we compare its predictions with the permeability obtained from Brownian dynamics simulations and find good agreement between the two.

Saturating Receiver and Receptor Competition in Synaptic DMC: Deterministic and Statistical Signal Models

Synaptic communication is based on a biological Molecular Communication (MC) system which may serve as a blueprint for the design of synthetic MC systems. However, the physical modeling of synaptic MC is complicated by the possible saturation of the molecular receiver caused by the competition of neurotransmitters (NTs) for postsynaptic receptors. Receiver saturation renders the system behavior nonlinear in the number of released NTs and is commonly neglected in existing analytical models. Furthermore, due to the ligands' competition for receptors (and vice versa), the individual binding events at the molecular receiver are in general not statistically independent and the commonly used binomial model for the statistics of the received signal does not apply. Hence, in this work, we propose a novel deterministic model for receptor saturation in terms of a state-space description based on an eigenfunction expansion of Fick's diffusion equation. The presented solution is numerically stable and computationally efficient. Employing the proposed deterministic model, we show that saturation at the molecular receiver effectively reduces the peak-value of the expected received signal and accelerates the clearance of NTs as compared to the case when receptor occupancy is neglected. We further derive a statistical model for the received signal in terms of the hypergeometric distribution which accounts for the competition of NTs for receptors and the competition of receptors for NTs. The proposed statistical model reveals how the signal statistics are shaped by the number of released NTs, the number of receptors, and the binding kinetics of the receptors, respectively, in the presence of competition. In particular, we show that the impact of these parameters on the signal variance is qualitatively different depending on the relative numbers of NTs and receptors. Finally, the accuracy of the proposed deterministic and statistical models is verified by particle-based computer simulations.

Mean first passage time for a particle diffusing on a disk with two absorbing traps at the boundary

The problem of survival of a Brownian particle diffusing on a disk with a reflective boundary that has two absorbing arcs is treated analytically. The framework of boundary homogenization is applied to calculate the effective trapping rate of the disk boundary, and this enables estimation of the mean first passage time. The method of conformal mapping is applied to transform the original system to a simpler geometrical configuration (a flat reflective boundary with a periodic configuration of identical absorbing strips) for which the analytical solution is known. The expression for the mean first passage time is simplified for some limiting cases (small arc or small gap). The derived analytical expressions compare favorably with the results of Brownian particle simulations and other analytical results from the literature.

Trapping of diffusing particles by short absorbing spikes periodically protruding from reflecting base

We study trapping of diffusing particles by a periodic non-uniform boundary formed by absorbing spikes protruding from a reflecting flat base. It is argued that such a boundary can be replaced by a flat uniform partially absorbing boundary with a properly chosen effective trapping rate. Assuming that the spikes are short compared to the inter-spike distance, we propose an approximate expression which gives the trapping rate in terms of geometric parameters of the boundary and the particle diffusivity. To validate this result, we compare some theoretical predictions based on the expression for the effective trapping rate with corresponding quantities obtained from Brownian dynamics simulations.

Effect of surface curvature on diffusion-limited reactions on a curved surface

To investigate how the curvature of a reactive surface can affect reaction kinetics, we use a simple model in which a diffusion-limited bimolecular reaction occurs on a curved surface that is hollowed inward, flat, or extended outward while keeping the reactive area on the surface constant. By numerically solving the diffusion equation for this model using the finite element method, we find that the rate constant is a non-linear function of the surface curvature and that there is an optimal curvature providing the maximum value of the rate constant, which indicates that a spherical reactant whose entire surface is reactive (a uniformly reactive sphere) is not the most reactive species for a given reactive surface area. We discuss how this result arises from the interplay between two opposing effects: the exposedness of the reactive area to its partner reactants, which causes the rate constant to increase as the curvature increases, and the competition occurring on the reactive surface, which decreases the rate constant. This study helps us to understand the role of curvature in surface reactions and allows us to rationally design reactants that provide a high reaction rate.

Trapping by clusters of channels, receptors, and transporters: quantitative description

Various membrane functional units such as receptors, transporters, and channels, whose action necessarily involves capturing diffusing molecules, are often organized into multimeric complexes forming clusters on the cell and organelle membranes. These functional units themselves are usually oligomers of several integral proteins, which have their own symmetry. Depending on the symmetry, they form clusters on different packing lattices. Moreover, local membrane inhomogeneities, e.g., the so-called membrane domains, rafts, stalks, etc., lead to different patterns even within the structures on the same packing lattice. Units in the cluster compete for diffusing molecules and screen each other. Here we propose a general approach that allows one to quantify the screening effects. The approach is used to derive simple approximate formulas giving the trapping rates of diffusing molecules by clusters of absorbers on lattices of different packing symmetries. The obtained results describe smooth variation of the trapping rate from the sum of the rates of individual absorbers forming the cluster to the effective collective rate. The latter shows how the trapping efficiency of an individual absorber decreases as the number of absorbers in the cluster increases and/or the inter-absorber distance decreases. Numerical tests demonstrate good agreement between the rates predicted by the theory and obtained from Brownian dynamics simulations for clusters of different shapes and sizes.

Trapping of diffusing particles by clusters of absorbing disks on a reflecting wall with disk centers on sites of a square lattice

A simple approximate formula is derived for the rate constant that describes steady-state flux of diffusing particles through a cluster of perfectly absorbing disks on the otherwise reflecting flat wall, assuming that the disk centers occupy neighboring sites of a square lattice. A distinctive feature of trapping by a disk cluster is that disks located at the cluster periphery shield the disks in the center of the cluster. This competition of the disks for diffusing particles makes it impossible to find an exact analytical solution for the rate constant in the general case. To derive the approximate formula, we use a recently suggested approach [A. M. Berezhkovskii, L. Dagdug, V. A. Lizunov, J. Zimmerberg, and S. M. Bezrukov, J. Chem. Phys. 136, 211102 (2012)], which is based on the replacement of the disk cluster by an effective uniform partially absorbing spot. The formula shows how the rate constant depends on the size and shape of the cluster. To check the accuracy of the formula, we compare its predictions with the values of the rate constant obtained from Brownian dynamics simulations. The comparison made for 18 clusters of various shapes and sizes shows good agreement between the theoretical predictions and numerical results.