Narrow escape in composite domains forming heterogeneous networks

Abstract: Cellular networks are often composed of thin tubules connecting much larger node compartments. These structures serve for active or diffusion transport of proteins. Examples are glial networks in the brain, the endoplasmic reticulum in cells or dendritic spines located on dendrites. In this latter case, a large ball forming the head is connected by a narrow passage. In all cases, how the transport of molecules, ions or proteins is regulated determines the time scale of chemical reactions or signal transduction. In the present study, based on modeling diffusion in three dimensions, we compute the mean time for a Brownian particle to reach a narrow target inside such a composite network made of tubules connected to spherical nodes. We derive asymptotic formulas by solving a mixed Neumann-Dirichlet boundary value problem with small Dirichlet part. We first consider the case of a network domain organized in a 2-D lattice structure that consists of spherical ball compartments connected via narrow cylindrical passages. When there is a single target we derive a matrix equation for each Mean First Passage Time (MFPT) averaged over each spherical compartment. We then consider a composite domain consisting of a spherical head-like domain connected to a large cylinder via a narrow cylindrical neck. For Brownian particles starting within the narrow neck, we derive formulas for the MFPT to reach a target on the spherical head. When diffusing particles can be absorbed upon hitting additional absorbing boundaries of the large cylinder, we compute the probability and conditional MFPT to reach a target. We compare these formulas with numerical solutions of the mixed boundary value problem and with Brownian simulations. To conclude, the present analysis reveals that the mean arrival time, driven by diffusion in heterogeneous networks, is controlled by the target and narrow passage sizes.