The effect of reactions on the formation and readout of the gradient of Bicoid

During early development, the establishment of gradients of transcriptional factors determines the patterning of cell fates. The case of Bicoid (Bcd) in Drosophila melanogaster embryos is well documented and studied. There are still controversies as to whether SDD models in which Bcd is Synthesized at one end, then Diffuses and is Degraded can explain the gradient formation within the timescale observed experimentally. The Bcd gradient is observed in embryos that express a Bicoid-eGFP fusion protein (Bcd-GFP) which cannot differentiate if Bcd is freely diffusing or bound to immobile sites. In this work we analyze an SDID model that includes the Interaction of Bcd with binding sites. We simulate numerically the resulting full reaction-diffusion system in a cylindrical domain using previously determined biophysical parameters and a simplified version of the Bcd source. In this way we obtain solutions that depend on the spatial location approximately as observed experimentally and that reach such dependence at a time that is also compatible with the experimental observations. Analyzing the differences between the free and bound Bcd distributions we observe that the latter spans over a longer lengthscale. We conclude that deriving the lengthscale from the distribution of Bcd-GFP can lead to an overestimation of the gradient lengthscale and of the Hill coefficient that relates the concentrations of Bcd and of the protein, Hunchback, whose production is regulated by Bcd.