Stable, Precise, and Reproducible Patterning of Bicoid and Hunchback Molecules in the Early Drosophila Embryo

Precise patterning of morphogen molecules and their accurate reading out are of key importance in embryonic development. Recent experiments have visualized distributions of proteins in developing embryos and shown that the gradient of concentration of Bicoid morphogen in Drosophila embryos is established rapidly after fertilization and remains stable through syncytial mitoses. This stable Bicoid gradient is read out in a precise way to distribute Hunchback with small fluctuations in each embryo and in a reproducible way, with small embryo-to-embryo fluctuation. The mechanisms of such stable, precise, and reproducible patterning through noisy cellular processes, however, still remain mysterious. To address these issues, here we develop the one- and three-dimensional stochastic models of the early Drosophila embryo. The simulated results show that the fluctuation in expression of the hunchback gene is dominated by the random arrival of Bicoid at the hunchback enhancer. Slow diffusion of Hunchback protein, however, averages out this intense fluctuation, leading to the precise patterning of distribution of Hunchback without loss of sharpness of the boundary of its distribution. The coordinated rates of diffusion and transport of input Bicoid and output Hunchback play decisive roles in suppressing fluctuations arising from the dynamical structure change in embryos and those arising from the random diffusion of molecules, and give rise to the stable, precise, and reproducible patterning of Bicoid and Hunchback distributions.Author Summary: For developing embryos, the precise, position-specific regulation of molecular processes is of fatal importance. As the mechanism of such regulation, widely accepted has been the notion of the intraembryonic distribution of regulatory molecules called “morphogens”. One of the best-studied morphogens is Bicoid in the early developmental stage of the Drosophila embryo. Synthesized around the anterior pole of the embryo, Bicoid forms an exponential gradient of concentration to initiate expression of a target gene, hunchback, in nuclei at the periphery of the embryo. This invariably forms a concentration boundary of the product protein Hunchback at around 49% embryo length. Remarkably, the embryo-embryo variability in the boundary position is less than 5%. Reactions in embryos, however, should be intrinsically noisy because the number of molecules involved is small, and those reactions are governed by randomly diffusing molecules. The mechanisms to generate the invariable Hunchback distribution by filtering the intense noise remain mysterious, and here we construct models to shed light on this problem. Stochastic simulations show that the slow diffusion of Hunchback averages out the intense noise, so that the coordinated rates of diffusion and transport of input Bicoid and output Hunchback play decisive roles in suppressing fluctuations.