Stochastic fluctuations in reaction-diffusion processes often have substantial effect on spatial and temporal dynamics of signal transductions in complex biological systems. One popular approach for simulating these processes is to divide the system into small spatial compartments assuming that molecules react only within the same compartment and jump between adjacent compartments driven by the diffusion. While the approach is convenient in terms of its implementation, its computational cost may become prohibitive when diffusive jumps occur significantly more frequently than reactions, as in the case of rapid diffusion. Here, we present a hybrid continuous-discrete method in which diffusion is simulated using continuous approximation while reactions are based on the Gillespie algorithm. Specifically, the diffusive jumps are approximated as continuous Gaussian random vectors with time-dependent means and covariances, allowing use of a large time step, even for rapid diffusion. By considering the correlation among diffusive jumps, the approximation is accurate for the second moment of the diffusion process. In addition, a criterion is obtained for identifying the region in which such diffusion approximation is required to enable adaptive calculations for better accuracy. Applications to a linear diffusion system and two nonlinear systems of morphogens demonstrate the effectiveness and benefits of the new hybrid method.

Morphogen gradients induce sharply defined domains of gene expression in a concentration-dependent manner, yet how cells interpret these signals in the face of spatial and temporal noise remains unclear. Using fluorescence lifetime imaging microscopy (FLIM) and phasor analysis to measure endogenous retinoic acid (RA) directly in vivo, we have investigated the amplitude of noise in RA signaling, and how modulation of this noise affects patterning of hindbrain segments (rhombomeres) in the zebrafish embryo. We demonstrate that RA forms a noisy gradient during critical stages of hindbrain patterning and that cells use distinct intracellular binding proteins to attenuate noise in RA levels. Increasing noise disrupts sharpening of rhombomere boundaries and proper patterning of the hindbrain. These findings reveal novel cellular mechanisms of noise regulation, which are likely to play important roles in other aspects of physiology and disease.

Positional information derived from local morphogen concentration plays an important role in patterning. A key question is how morphogen diffusion and gene expression regulation shape positional information into an appropriate profile with suitably low noise. We address this question using a model system--the C. elegans germline--whose regulatory network has been well characterized genetically but whose spatiotemporal dynamics are poorly understood. We show that diffusion within the germline syncytium is a critical control of stem cell differentiation and that semi-permeable diffusion barriers present at key locations make it possible--in combination with a feedback loop in the germline regulatory network--for mitotic zone size to be robust against spatial noise in Notch signaling. Spatial averaging within compartments defined by diffusion barriers is an advantageous patterning strategy, which attenuates noise while still allowing for sharp transitions between compartments. This strategy could apply to other organs.