Stable, precise, and reproducible patterning of bicoid and hunchback molecules in the early Drosophila embryo.  
Jump to Full Text  
MedLine Citation:

PMID: 19714200 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

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 embryotoembryo 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 threedimensional 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. 
Authors:

Yurie OkabeOho; Hiroki Murakami; Suguru Oho; Masaki Sasai 
Publication Detail:

Type: Journal Article; Research Support, NonU.S. Gov't Date: 20090828 
Journal Detail:

Title: PLoS computational biology Volume: 5 ISSN: 15537358 ISO Abbreviation: PLoS Comput. Biol. Publication Date: 2009 Aug 
Date Detail:

Created Date: 20090828 Completed Date: 20091116 Revised Date: 20130531 
Medline Journal Info:

Nlm Unique ID: 101238922 Medline TA: PLoS Comput Biol Country: United States 
Other Details:

Languages: eng Pagination: e1000486 Citation Subset: IM 
Affiliation:

Department of Computational Science and Engineering, Nagoya University, Nagoya, Japan. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Animals Body Patterning* Computational Biology / methods* Computer Simulation DNABinding Proteins / genetics, physiology* Databases, Protein Drosophila / embryology*, physiology* Drosophila Proteins / genetics, physiology* Gene Expression Profiling Gene Expression Regulation, Developmental Homeodomain Proteins / genetics, physiology* Mitosis Models, Biological Models, Genetic Models, Statistical Reproducibility of Results TransActivators / genetics, physiology* Transcription Factors / genetics, physiology* 
Chemical  
Reg. No./Substance:

0/DNABinding Proteins; 0/Drosophila Proteins; 0/Homeodomain Proteins; 0/TransActivators; 0/Transcription Factors; 0/bicoid protein, Drosophila; 0/hunchback protein, Drosophila 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): PLoS Comput Biol Journal ID (publisherid): plos Journal ID (pmc): ploscomp ISSN: 1553734X ISSN: 15537358 Publisher: Public Library of Science, San Francisco, USA 
Article Information Download PDF OkabeOho et al. This is an openaccess article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Received Day: 23 Month: 3 Year: 2009 Accepted Day: 23 Month: 7 Year: 2009 collection publication date: Month: 8 Year: 2009 Print publication date: Month: 8 Year: 2009 Electronic publication date: Day: 28 Month: 8 Year: 2009 Volume: 5 Issue: 8 Elocation ID: e1000486 ID: 2720536 PubMed Id: 19714200 Publisher Id: 09PLCBRA0298R3 DOI: 10.1371/journal.pcbi.1000486 
Stable, Precise, and Reproducible Patterning of Bicoid and Hunchback Molecules in the Early Drosophila Embryo Alternate Title:Stable, Precise, and Reproducible Patterning  
Yurie OkabeOho^{1}  
Hiroki Murakami^{1}  
Suguru Oho^{2}  
Masaki Sasai^{1}^{3}^{4}*  
Daniel A. Beardedit1 
Role: Editor 
1Department of Computational Science and Engineering, Nagoya University, Nagoya, Japan 

2Department of Environmental Engineering and Architecture, Nagoya University, Nagoya, Japan 

3Department of Applied Physics, Nagoya University, Nagoya, Japan 

4School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Korea 

Medical College of Wisconsin, United States of America 

Correspondence: * Email: sasai@nuap.nagoyau.ac.jp Contributed by footnote: Conceived and designed the experiments: YOO HM MS. Performed the experiments: YOO HM. Analyzed the data: YOO HM SO MS. Contributed reagents/materials/analysis tools: YOO HM SO. Wrote the paper: YOO HM MS. 
Pattern formation of multicellular organisms requires the accurate spatial regulation of gene expression. Such regulation has been attributed to the position dependent distribution of signaling molecules called morphogens ^{[1]},^{[2]}. One of the best studied morphogens is Bicoid (Bcd) in embryos of Drosophila melanogaster, which works as a transcription factor to regulate the expression of hunchback (hb) and other downstream genes ^{[3]}?^{[5]}. Bcd is synthesized from the maternal mRNA which is localized at around the anterior pole of embryo, and the concentration of Bcd exhibits a gradient with approximately exponential decay with distance from the anterior pole ^{[6]}. Although the decisive importance of the gradient in concentration of Bcd has been established ^{[7]},^{[8]}, the mechanism for forming the stable Bcd gradient and that for the reliable readout of the gradient remain largely unresolved.
Important advances have been made by visualizing spatial patterns of distributions of the input protein Bcd and the output protein Hunchback (Hb) with the antibody staining technique ^{[6]},^{[9]},^{[10]}. These visualized data, however, have been taken from fixed embryos, which provided only the snapshots of developing embryos. Further insights were gained by combining the method using the enhanced green fluorescent protein (eGFP) with the antibody staining technique: In two papers ^{[11]},^{[12]}, Gregor et al. reported the data of dynamically changing Bcd patterns in living embryos by observing fluorescence of a BcdeGFP fusion protein. These data provoked fundamental questions as to the formation and readout of the Bcd gradient, which presents challenges in computational biology.
In their first paper ^{[11]}, Gregor et al. showed that the Bcd distribution is approximately an exponential gradient in intranuclear concentration, which is established rather rapidly in less than 90 min after fertilization. As Bcd is a DNAbinding protein, Bcd is localized in nuclei. However, by observing responses of distributions of BcdeGFP to the photobleaching stimuli, it was shown that Bcd is not simply trapped in nuclei but is in dynamic equilibrium between influx and efflux with the cytoplasm and intranuclear degradation. The early embryo in the blastodermal stage undergoes a series of syncytial mitoses from nuclear cycles 10 to 14, during which the nuclei are not separated by cell membranes. It was observed that Bcd which is stored in nuclei is released into cytoplasm as nuclei lose their envelopes during each mitosis, and the released Bcd is concentrated into nuclei again as envelopes are reformed during the interphase period. These observations raised the question on how the stable Bcd gradient is established in spite of the repeated dynamical structural change of embryo due to mitoses of nuclei. Though the effects of dynamically changing nuclei have begun to be to be investigated theoretically ^{[13]}?^{[15]}, most of the hitherto developed theoretical models of the Bcd gradient have been based on the assumptions of diffusion and degradation of Bcd in a static continuum ^{[16]}?^{[20]}.
In the second paper ^{[12]}, the authors reported the data of fluctuations in concentration profiles of Bcd and Hb and showed that formation of those profiles is precise and reproducible: It was shown that fluctuations of profiles of Bcd and Hb in each embryo are about 10%, which are small enough to distinguish adjacent nuclei, and that fluctuations from embryo to embryo are also as small as 10%. These data imply that the system exerts precise control over concentration of Bcd, and that Hb is produced as a reliable response to the small concentration difference of Bcd. Response of the hb activity to Bcd, however, is the process to receive randomly arriving Bcd molecules at the hb enhancer on DNA. The random arrival of Bcd at the nanometer scale region of DNA should induce intense noise in the hb activity and Gregor et al. argued that this effect should be so large that we have to assume some unconfirmed mechanism working to reduce the noise of the readout process to the observed level of precision. There are further sources of intrinsic noise in the readout process arising from the bursting production of proteins and the nonadiabatic gene switching ^{[21]}?^{[28]}, so that the observed precise readout of the Bcd gradient implies that the noise arising from those sources is also suppressed to the level of 10% fluctuations.
In this paper, we construct two computational models, onedimensional and threedimensional models of Drosophila embryo, to discuss two major issues raised by Gregor et al.; how the stable patterning of the Bcd gradient is realized in the dynamically changing syncytial embryo; and how the Bcd gradient is read in a precise and reproducible way in spite of intense noise induced by the random diffusion of Bcd molecules. Our onedimensional model of embryo focuses on the latter issue on the noise induced by random diffusion, and effects of the dynamically changing embryonic structure are examined with the threedimensional model.
Importance of the noise induced by random diffusion has been stressed by Berg and Purcell in their early study of bacterial chemotaxis ^{[29]}. They have argued that certain physical limit of sensitivity should exist when the bacterial receptor receives the randomly arriving chemical signals. Bialek and colleagues ^{[30]}?^{[32]} extended this analysis in a theoretically more systematic way and showed that the physical limit arising from random diffusion of signaling molecules is not the special feature of bacterial chemotaxis, but is common to many biological processes of receiving chemical signals. Gregor et al. ^{[12]} applied Bialek's theory to the problem of Drosophila embryo and showed that the limit of accuracy in sensing the Bcd concentration is
where c is the intranuclear concentration of Bcd, and ?c is fluctuation in the Bcd concentration sensed by the hb enhancer. The hb enhancer is assumed to have the typical length a. D is the intranuclear diffusion constant of Bcd, and ? is the time length during which Bcd can bind to the enhancer. In Drosophila embryo, decision of whether hb is activated or not is made at the point of threshold concentration ; hb is activated in the nuclei of and is kept silent if . Since the observed values of D and are and ^{[12]}, and the size of the enhancer binding site should be as small as , Equation 1 implies at the decision point of the hb activation ^{[12]}. When we regard ? as the typical time length of one nuclear cycle, we have , which leads to . In Ref.12, ?c was estimated from the data of fluctuation of the Hb concentration by assuming that ?c is reflected to the hb gene activity and hence is reflected to the fluctuation in concentration of Hb. The value observed in this way was , which is much smaller than the expected value of . To explain the small value of , ? should be almost two hours, which is unacceptably long and near to the entire time available for development from fertilization up to cellularization. Thus, some mechanism lacking in this argument is needed to reduce the expected large noise of to the level of the observed precision . We may refer to this problem as the paradox of signal interpretation.Though Bialek and colleagues derived Equation 1 by solving the coupled equations of diffusion, degradation, and binding/unbinding processes, the same result can be obtained by the intuitive argument as illustrated in Figure 1: Since Bcd molecules diffuse over a distance during the allowed time ?, Bcd molecules interacting with the target of size a during the period ? are distributed in the volume . We call this volume ?interaction volume?. The number of molecules in the interaction volume is in average. If Bcd molecules diffuse as independent particles, fluctuation in n should obey the Poisson statistics to exhibit . Then, we have , which gives the lower limit of fluctuation, and in more generic cases with other additional sources of noise, we have the criterion of Equation 1. In the following, we explicitly use the idea of interaction volume to model the readout process in embryo.
In the next section, we first discuss the onedimensional stochastic model of embryo. To focus on the problem of the paradox of signal interpretation, we neglect the dynamical structural change of embryo and assume that nuclei stay still without exhibiting nuclear cycles. With this model we perform the stochastic simulation of reactions and diffusion of Bcd and Hb. Bcd molecules are produced at the anterior pole, diffuse, degrade, and bind/unbind to/from the hb enhancer, and Hb molecules are produced with the rate depending on the hb enhancer status, and diffuse and degrade. With this onedimensional model, fluctuations in distributions of Bcd and Hb can be examined in a wide range of parameters and we show that the selfaveraging process due to the Hb diffusion should resolve the paradox of signal interpretation.
Effects of the dynamically changing embryonic structure are investigated by constructing a more realistic threedimensional stochastic model of embryo. In the threedimensional model, an embryo is represented by a cylinder consisting of 43,750 hexagonal sites. Bcd molecules are synthesized at the most front slice of this cylinder. The inside structure of this cylinder is modified during simulation by following the rules governing nuclear mitoses in the model: At the initial instance, one site at the core of the cylinder is assumed to be a nucleus. This nucleus starts to divide, and at the nuclear cycle 8?10, the multiplicated nuclei move to the peripheral. At the nuclear cycle 14, about 6000 sites at the surface layer of the cylinder are nuclear sites. Within the nuclear site the same reactions are assumed as in the onedimensional model, and in other cytosolic sites degradation and diffusion of Bcd and Hb are considered. We show that patterning of the Bcd profile is stable and precise in spite of dynamic mitoses of nuclei in embryo. We further confirm that the same mechanism of selfaveraging of Hb molecules as discussed with the onedimensional model also works in the threedimensional model to resolve the paradox of signal interpretation.
A simplified onedimensional model of embryo is developed to focus on the problem of precise readout of the Bcd gradient. In this model, nuclei are assumed to stay static to simulate the stable distribution of Bcd during the interphase period of nuclear cycle 14. In this period, about 6000 nuclei are located in the vicinity of the surface of the embryo whose size varies around 500 ?m ^{[33]}. Each of these nuclei is not wrapped by the cell membrane but is exposed to the common cytoplasm in a syncytial embryo. These nuclei are arranged in about 100 rows along the anteriorposterior (AP) axis of embryo ^{[34]}. Here, such a distribution of nuclei is modeled by the array of nuclei placed along the onedimensional AP axis. As shown in Figure 2A, the onedimensional system of length L is composed of N sites each of which has the length of ?x?=?L/N with and N?=?100. Each site contains a nucleus, and numbers of Bcd and Hb molecules in the ith site are denoted by N_{Bcd}(i) and N_{Hb}(i), respectively. With this system, we performed the stochastic simulation of the following molecular processes: (i) binding/unbinding of Bcd and Hb to/from the hb enhancer; (ii) diffusion and degradation; and (iii) synthesis of Bcd and Hb. See Text S1 for more detailed explanation of the model.
(i) Binding/unbinding to/from the hb enhancer: Random binding of diffusive molecules to the hb enhancer is the source of noise which should lead to the signal interpretation paradox. Since the size of the hb enhancer is in nanometer scale and a nucleus extends to the micrometer scale, each nucleus should be represented by the mesh of 1000 or more refined lattice points to directly describe this diffusionbinding process in the model. Such refinement requires the high computational cost, which should impair the merit of the onedimensional model in survey of widely different parameter regions. Instead of refining the mesh size, we here use the concept of ?interaction volume? to describe the random arrival of molecules at the nanometer region. We assume that Bcd and Hb have the same diffusion constant inside the nucleus and define the interaction volume of size , where is the intranuclear diffusion constant of Bcd or Hb, a is the size of the hb enhancer, and ? is the time length of the interphase of nucleus. Use of different intranuclear diffusion constants for Bcd and Hb does not add significant effects in the present simplified model.
N_{Bcd}(i) is divided into the part outside the interaction volume , and the part inside the interaction volume , where is the number of Bcd molecules that remain in the interaction volume without binding to the hb enhancer, and is the number of Bcd molecules that are bound to the hb enhancer. N_{Hb}(i) is also divided in the same way, so that
We assume that only the molecules in the interaction volume can bind to the hb enhancer. Due to the intranuclear diffusion, molecules are stochastically transferred into/from the interaction volume. We write the volume of a nucleus as , so thatrepresents the relative size of the interaction volume. The limit of V_{r}?=?1 corresponds to the case that the inside of nucleus is so well stirred that the noise induced by the random arrival of Bcd at the hb enhancer is negligible. The smaller V_{r} is, the smaller is allowed to stay in the interaction volume and fluctuations of become larger, which enhances the noise in the binding reaction of Bcd. For V_{r}?1, the noise in the hb expression is dominated by this effect.The hb enhancer contains multiple binding sites of Bcd and Hb, which should allow the various different states of the hb enhancer depending on the number of bound Bcd and Hb. As suggested in experiment ^{[35]},^{[36]}, we assume that the hb enhancer has six sites to bind Bcd and two sites to bind Hb. Since there can be cooperative interactions among bound Bcd and Hb ^{[35]}?^{[38]}, the most populated states should be categorized into the following four states; the state with no bound Bcd or Hb (S(i)?=?0), the state with two bound Hb molecules and with no bound Bcd (S(i)?=?1), the state with six bound Bcd molecules and with no bound Hb (S(i)?=?2), and the state with six bound Bcd and two bound Hb molecules (S(i)?=?3). We consider that the cooperative binding process of Bcd takes place as
with the rate and that the cooperative binding process of Hb takes place aswith the rate . Due to the large fluctuation of and , binding processes of Equation 4 and Equation 5 should bear intense noise as argued in the paradox of signal interpretation. In this paper, embryos homozygous for the mutated hb allele ^{[36]},^{[39]},^{[40]}, which codes for Hb having noDNA binding affinity, are also simulated by assuming reactions of Equation 4. For such mutant embryos, we put h_{h}?=?0 and only consider S(i)?=?0 and S(i)?=?2 as possible states of the hb enhancer.Unbinding reactions are the reverse processes of Equations 4 and 5 taking place with the rate for the unbinding of Bcd and with the rate for the unbinding of Hb. The ratio or determines the affinity of Bcd or Hb to the hb enhancer site. We should also note that and can be used as measures of time scale of change in the hb enhancer state.
(ii) Diffusion and degradation: Due to the diffusive movement in cytoplasm, Bcd molecules are transferred between adjacent sites with the rate d_{b}N_{Bcd}(i), where d_{b} is proportional to the diffusion constant of Bcd, D_{b}, as d_{b}?=?D_{b}/?x^{2}. Bcd is degraded at each site with the rate or . Diffusion and degradation of Hb are treated in the same way as those of Bcd but with different rate constants: Hb molecules are transferred between adjacent sites with the rate d_{h}N_{Hb}(i), where d_{h} is proportional to the diffusion constant D_{h} as d_{h}?=?D_{h}/?x^{2}. Hb is degraded at each site with the rate or .
(iii) Synthesis: Bcd is synthesized from the maternal mRNA distributed at around the anterior pole, which is simulated by assuming that Bcd is generated only at the site i?=?1 with the rate J. Reactions involved in transcription and translation of hb are coarsegrained and described as a combined single reaction step, which is the approximation adopted in studies of gene expression ^{[26]}?^{[28]}. Then, the process to synthesize Hb around the ith nucleus is
where N_{burst} is the number of Hb molecules synthesized from a shortlived hb mRNA, representing the bursting production of Hb ^{[25]},^{[28]}. The bursting production with N_{burst}?1 can be a source of intrinsic noise of hb expression ^{[22]},^{[28]}. As illustrated in Figure 2B, we consider that Equation 6 is executed with the rate g(S(i)), where S(i)?=?0?3 is the state of the hb enhancer in the ith site. As suggested by the experimental data ^{[36]},^{[40]}, binding of Hb and binding of Bcd should cooperatively activate synthesis of Hb, and hence we assumeFor a mutant of loosing binding affinity of Hb with h_{h}?=?0, we only consider the case of S?=?0 and S?=?2, so thatThe scheme of BcdHb cooperation of Equation 7 implies that Hb activates expression of hb, which constitutes a positive feedback loop of Hb production. In the scheme of Equation 8, on the other hand, the positive feedback loop is absent. In this paper we compare two schemes of Equation 7 and Equation 8 to examine the effect of the positive feedback.In the ensemble of many embryos, the average system size is approximately , and L fluctuates from embryo to embryo with the relative amplitude of ^{[12]} where is the average over many embryos. Corresponding to these average and fluctuation, we repeat stochastic simulations 100 times to take the ensemble average of patterning of Bcd and Hb by randomly varying ?x at each simulation run with , and . To analyze the calculated results of different trajectories in a unified way, , , , , , , , and , obtained at the last step of each simulation run are normalized by the factor to be , , , , , and . We consider that the intranuclear diffusion constant is similar to the one for the shortlength diffusion in cytoplasm as D_{nuc}?=?0.3?1.0 ?m^{2}/s ^{[11]}. Since a should be a typical length in DNA and ? should be a typical duration of interphase, we assume as the length of 10 base pairs and , which leads to . Diameter of a nucleus at nuclear cycle 14 is about 6.5 ?m in average, so that . We use .
Other parameters used in the model are explained in Method section and are summarized in Table 1. We refer to the parameterization of Table 1 as the standard parameterization in the onedimensional model. Results of different parameter values from those of Table 1 are also examined extensively with the onedimensional model. Especially to be mentioned in Table 1 is the value of the diffusion constant D_{b} of Bcd in cytoplasm. In the standard parameterization, we assume D_{b}?=?10 ?m^{2}/s from the following two reasons: (i) The diffusion constant of inert molecules moving through cytoplasm was observed to be as large as 10 ?m^{2}/s ^{[10]}; (ii) If bcd mRNA, and hence synthesis of Bcd, are localized at the anterior pole of embryo, Bcd has to diffuse as fast as 10 ?m^{2}/s to form a suitable Bcd gradient in a realistic time scale. The smaller value is used, on the other hand, for the intranuclear diffusion constant as D_{nuc}<1.0 ?m^{2}/s. The use of a larger value of will result in , and the fluctuation of the hb expression will be largely reduced with such a large value of . D_{nuc}, however, is the diffusion constant in a short time scale in a much denser medium than cytoplasm, which should be nearer to the observed diffusion constant of 0.3?1.0 ?m^{2}/s in cortex ^{[11]},^{[12]}, so that the large value of seems not appropriate and instead we should adopt the value , which leads to a small as used in Table 1.
After finishing preparation of the present manuscript, we became aware of the recent experimental report ^{[41]} suggesting that bcd mRNA is not strictly localized at the anterior pole but forms a gradient along the AP axis. With such a gradient of bcd mRNA, the smaller D_{b} might be enough to explain the observed Bcd gradient. Intense interest has been focused on the diffusion constant of Bcd ^{[11]},^{[14]},^{[41]}, but consensus has not yet been achieved on either of the value of diffusion constant or the role of cytoplasmic diffusion of Bcd. We here adopt assumptions of the localized bcd mRNA and the large D_{b} as have been traditionally assumed in theoretical models ^{[10]}, ^{[13]}, ^{[16]}?^{[20]},^{[36]}. Further discussion on the value of diffusion constant will be given in Discussion section of the present paper.
We emphasize that embryo is not in a steady state but changes dynamically during development. We decompose problems, however, into two categories; dynamical problems and problems of accurate reading out. The dynamical problems will be considered in the more realistic threedimensional model, and with the simplified onedimensional model, we concentrate on the problems of accurate reading out by assuming the stationary nuclear configuration. This simplification in the onedimensional model allows the extensive survey of the input/output relation over the wide parameter range. With the onedimensional model, stochastic simulations were performed to reach the steady state and the data obtained at the last step of simulation were used for analyses. Then, the simulation was repeated 100 times to take the ensemble average. Dynamical quantities will be examined in the threedimensional model introduced later in this paper, which has the appropriate time scale and is suited to examine the temporal change of quantities. Results of the onedimensional model will be checked by the threedimensional model to confirm whether the decomposition of problems into dynamic and static ones correctly describes the important features of embryos.
In Figure 3A calculated with the onedimensional model with the standard parameterization is plotted as a function of position, . Here, results of 100 simulation runs are superposed and their average is compared with the observed data (Figure 5A of Ref.12). The panel shows that adjustment of a parameter J enables the model to consistently reproduce the observed Bcd profile in a semiquantitative way. Plotted in Figure 3B is for 100 simulation runs and their average . Plotted in Figure 3C is , which measures the relative gene activity at the ith site, and its average . and are plotted in Figure 3D and and are plotted in Figure 3E. We see that fluctuations in and are much more significant than in and as expected from the smallness of and . These large fluctuations are responsible for the large fluctuation of shown in Figure 3C. In spite of such large fluctuation in , fluctuation in seems to be moderate.
Fluctuation in can be quantitatively analyzed with the input/output (IO) relation by regarding Bcd as input and Hb as output. We define the mean input as , where i_{half} is the position satisfying and is the maximal value of in the system. The mean output is defined by . Thus defined mean IO relation is plotted in Figure 4A, which agrees well with the experimental data (Figure 4A of Ref.12). Fluctuation in the IO relation is shown in Figure 4B, where we plot as a function of with . We find that the amplitude of fluctuation shown in Figure 4B is about 10%, which is as small as was observed in experiment (Figure 4B in Ref.12). The fluctuation of output data shown in Figure 4B can be converted to the form of equivalent input data as in Figure 4C by using the relation , where , and . We can compare with of Equation 1 to see whether the output fluctuation is as large as that expected in Bialek's argument: As plotted in Figure 4C, the calculated fluctuation in output, , is around 10%, which is as small as was observed in experiment (Figure 4C of Ref.12) and is much smaller than the value expected in the argument of Equation 1. Thus, the present onedimensional model shows that either the random diffusion of Bcd and Hb or the noisy expression of hb does not bring about large fluctuations in the IO relation.
It should be noted that in the onedimensional model the IO relation is calculated by averaging quantities over the ensemble which is different from that in the experimental data: In experiment, many stripes along the AP axis of a single embryo can be used as an ensemble of data to measure the fluctuation of Bcd or Hb concentration in single embryo, and the small IO fluctuation observed in this ensemble is denoted by ?preciseness? in Ref.12. In the present onedimensional model, many runs of simulation are necessary to obtain the ensemble of enough size, and we here regard the smallness of fluctuations in this ensemble of many simulation runs as ?preciseness? in the readout process. The same ensemble as in the experimental measurement can be taken in the threedimensional model, and we will show later that the results of the threedimensional model are essentially same as those of the present onedimensional model, which should validate the use of this ensemble in the onedimensional model to analyze precision in the readout process.
Effects of the parameter variation on the precision can be examined by changing parameters one by one from values in the standard parameterization. In Figure 5, the IO relation is analyzed by changing parameters, V_{r}, f_{b}, f_{h}, N_{burst}, and D_{h}. Parameters f_{b} and f_{h} define the rate of change in the hb enhancer status. Instead of using the bare values of f_{b} and f_{h}, we here use the normalized forms of ?_{b}?=?f_{b}/k_{b} and ?_{h}?=?f_{h}/k_{h} to explain the results. From computational ^{[25]},^{[28]} and theoretical ^{[24]} investigations, we can expect that decrease in the value of ?_{b} or ?_{h} leads to increase in noise in the Hb synthesis. The large N_{burst} can also be a source of noise in the Hb synthesis ^{[22]}. Results of Figure 5 indicate, however, that changes in these parameters do not strongly affect the IO relation unless too small ?_{b} or ?_{h} is used. When V_{r} or D_{h} is varied, on the other hand, the clear change of the IO relation can be found: Smaller V_{r} gives the larger noise in the IO relation, showing that the random diffusion of Bcd and Hb are dominant sources to limit the precision of readout as was suggested by Gregor et al. ^{[12]}. D_{h}?=?0 corresponds to the case that the spatial heterogeneity of hb expression is directly reflected in the Hb distribution. In this case the fluctuation in the IO relation is large as was expected from the argument of Gregor et al. ^{[12]}: The random diffusion and random reception at the hb enhancer lead to the fluctuation larger than the observed data, which exceeds the required precision to distinguish neighboring nuclei.
These dependences on parameters are further examined in Figure 6. In Figure 6A, sensitive dependence of the fluctuation of Hb concentration on V_{r} and D_{h} is shown. When we assume the nondiffusive Hb with D_{h}?=?0, the calculated results show the small enough amplitude of fluctuation only when V_{r} is larger than 10^{?2}. With the realistic value of , small but finite D_{h} with D_{h} is necessary to explain the observed 10?20% fluctuation in the readout process. Effects of changing ?_{b} and ?_{h} were also examined (Figure S1): Small enough values of ?_{b} and ?_{h} bring about the large fluctuation, but in the parameter region of ?_{b}>1 and ?_{h}>1, the fluctuation amplitude is moderate and does not much depend on ?_{b} or ?_{h}. In Figures 6B and 6C, dependence of the profile of Hb distribution on V_{r} and D_{h} is examined. The number of Hb molecules become very small when V_{r}?10^{?3} and D_{h}>0.5 ?m^{2}/s, but does not sensitively depend on V_{r} or D_{h} when V_{r}>10^{?3} or D_{h}<0.5 ?m^{2}/s (Figure 6B). The angle of gradient of Hb distribution, on the other hand, sensitively depends on D_{h} (Figure 6C). From this dependence of the angle on D_{h}, we see that D_{h} should be less than 0.3 ?m^{2}/s to explain the observed experimental profile of Hb distribution ^{[12]}. Diffusion constant of D_{h}<0.3 ?m^{2}/s is consistent with the data that Hb does not diffuse over the long distance during a nuclear cycle ^{[36]}.
Sensitive dependence of the fluctuation of Hb distribution on V_{r} is the result expected from the argument of Gregor et al. ^{[12]} and Bialek et al. ^{[30]}?^{[32]}, but the suppression of fluctuation by nonzero D_{h} is rather unexpected. To understand the reason why the small but finite D_{h} suppresses the Hb fluctuation, we performed simulation by using Equation 8, the assumption of disability of Hb to bind to DNA. With this assumption, the positive feedback loop of the hb regulation is lost and the cooperativity at binding sites of the hb enhancer is reduced. By putting h_{h}?=?0 and prohibiting the state S?=?3, the number of produced Hb decreases to about 25% of the amount calculated with the standard parameterization as was observed in the experiment of Ref.40. The slope angle of the Hb concentration at the threshold position becomes slightly smaller due to the decrease of cooperativity in the hb activation, which is consistent with the observation of Ref.36. The features of fluctuation in the IO relation, however, do not show a significant difference from the results obtained with the standard parameterization. In Figure 7, sensitive dependence of fluctuation in the IO relation on V_{r} and D_{h} in the case of h_{h}?=?0 is shown. The fluctuation is dominated by the small V_{r} effect, but is substantially suppressed by the small but finite diffusion constant of Hb with 0<D_{h}<0.3 ?m^{2}/s. Thus, we conclude that not the feedback loop of the hb regulation but the selfaveraging due to the diffusion of Hb is sufficient to suppress the large fluctuations in the hb expression. With D_{h}?0.1 ?m^{2}/s, Hb molecules synthesized at around each nucleus are averaged over positions of several neighbor nuclei through diffusion. This simple mechanism is sufficient to suppress the fluctuation in the Hb distribution and hence resolves the paradox of signal interpretation.
Effects of the dynamically changing nuclei are examined by developing the threedimensional model. As shown in Figures 8A and 8B, an embryo is represented by a cylinder of length L extending along the AP axis. 15 embryos were sampled in Refs.11 and 12, and their average size was found to be about 490 ?m. In these 15 embryos, about 70 nuclei were distributed at the intersection of the midsagittal plane. In order to quantitatively compare the simulated results with these observed data, we use the system having the similar size: We consider the cylinder composed of 70 plates each of which has the thickness of ?x?=?L/70 with L?490 ?m. Each plate has an approximately circular shape of radius 90 ?m, and consists of 625 hexagonal sites whose side length is ??4 ?m. The system, therefore, consists of 70?625?=?43750 sites. Starting from t?=?0 at the time of oviposition, the dynamical change of distribution of Bcd and Hb is followed with the stochastic simulation until t?=?144 min at the interphase of nuclear cycle 14. This simulation was repeated 10 times by varying ? and ?x around averages and with deviations of and to simulate the observed fluctuation of 4.1% in size of embryo ^{[12]}. The system has three types of sites; cortical, core, and nuclear sites. Each nuclear site contains a single nucleus and other types of sites do not contain nucleus. Position of the site is designated by (i, j), where i is the number of slice from i?=?1 at the anterior slice to i?=?70 at the posterior slice. j?=?1?625 designates position of the site in the ith slice. The numbers of Bcd and Hb in each site at time t are denoted by and , where ??=?cortical, core, or nuclear specifies the type of site. Diffusion and reactions defined at sites depend on types of sites.
We should keep in mind that the shape of embryo is more like a prolate ellipsoid and that cylinder is just a 0th order approximation to it. Due to this disagreement in shape of the system, the simulated Bcd and Hb distributions at the anterior pole should deviate from the observed distributions. These deviations, however, should be localized at the anterior pole and are not important to examine the global features of the Bcd gradient and the boundary of the Hb distribution observed at around the middle of embryo. The present cylindrical model, therefore, provides a computationally efficient tool to study the global Bcd and Hb distributions in embryo. See Text S1 for more detailed explanation of the model.
Diffusion and reactions of Bcd and Hb defined in the threedimensional model are the straightforward extension of those considered in the onedimensional model. Bcd and Hb diffuse between adjacent lattice sites in the threedimensional space as in the same way as in the onedimensional model. Diffusion constants, however, depend on types of starting or arriving sites, which are explained in more details in Figure 8C and in Method section. Bcd and Hb are degraded with the rate of and . Here, and are assumed to be small in cortical or core sites but large in nuclear sites. At the core sites in the slice i?=?1, Bcd is generated with the rate J(m) with m representing the nuclear cycle. Since the observed production rate of Bcd increases as nuclear cycles proceed ^{[42]}, we assume J(m) to be an increasing function of m.
As in the onedimensional model, effects of random reception of Bcd or Hb molecules at the hb enhancer are simulated by introducing the nanometer scale spatial region of interaction volume. We use V_{r} to designate the ratio of interaction volume to the volume of a nucleus, so that the smallness of V_{r} is the source of fluctuations in the hb expression. and are decomposed into three parts as in the onedimensional model depending on whether the molecule is in or out of the interaction volume, and on whether the molecule is bound to the hb enhancer or not, so that
The hb enhancer state in the threedimensional model is S(i, j, t)?=?0?3. Hb is synthesized at around nuclei, i.e. in nuclear sites with the rate g(S(i, j, t)). Values of g(S(i, j, t)) can be determined in the same way as in Equation 7 or Equation 8. Further details on parameters are explained in Method section and are summarized in Table 2. We refer to the parameterization of Table 2 as the standard parameterization in the threedimensional model.Nuclei undergo nuclear cycles, increase in number, and change their locations in embryo. One nuclear cycle consists of the interphase period and the mitotic period. During the interphase period, nucleus is stable, which is described by a nuclear site in the model. During the mitotic period, nuclear membrane is dissolved and nuclear sites are turned into cortical or core sites in the model. After the mitosis, nuclear membrane is rebuilt to enter into the next nuclear cycle, and the number of nuclear sites is doubled. Since it was observed that somatic nuclei in Drosophila embryo undergo nearly synchronous nuclear divisions as a syncytium ^{[43]}, we treat these nuclear cycles as synchronized events in the model.
At t?=?0, one nuclear site is defined at around the center of embryo and nuclear cycle 1 is initiated. During nuclear cycle n?=?2 to 7, the number of nuclei increases as 2^{n}^{?1}. The increased nuclei are placed randomly at around the center of embryo. At the 8th nuclear cycle, about 90 nuclei of total 2^{7}?=?128 nuclei begin to move to the surface region of embryo and the remaining nuclei are left in the middle of embryo ^{[43]},^{[44]}. In the model, about 90 nuclear sites are redistributed randomly in the region somewhat apart from the center of the cylinder. The nuclei left in the center are yolk nuclei ^{[43]},^{[44]}, which are turned into the core sites in the model. About 90 nuclei departing from the center reach the surface region at nuclear cycle 10 and are doubled in number through nuclear cycle 14. With this rule of mitoses and movement, the surface of the cylinder in the model is occupied by about 5500?6200 nuclear sites at the interphase period of nuclear cycle 14 depending on the number of nuclei selected to move toward the cylinder surface at nuclear cycle 8. More detailed explanation of rules is given in Text S1. During the interphase period, volume of nuclei increases in Drosophila embryo ^{[11]} but here we neglect effects of the nuclear volume change and keep the size of nuclear site constant for the sake of simplicity.
Sites near the cylinder surface are defined as cortical sites if they are not nuclear sites. The cortical sites form a single layer of sites on the surface of cylinder until nuclear cycle 10, and the thickness of cortical layer increases as nuclear cycle proceeds from 11 to 14. See Text S1 for the further explanation on cortical sites. Sites other than cortical or nuclear sites are core sites. In this way the dynamical change in embryonic structure is represented by rules defined at each time step in the model. We perform stochastic simulations of diffusion and reactions of Bcd and Hb under these rules to see how such dynamically changing embryonic environment affects the stability, precision and reproducibility of distributions of Bcd and Hb.
With the threedimensional model, stochastic simulations of diffusion and reactions are performed. 10 simulation runs were repeated by allowing the fluctuation of embryonic size with . In order to compare the calculated results with the observed data, quantities are normalized by the factor as , , , , , and .
Shown in Figure 9 are snapshots of a simulation run with the threedimensional model, in which the Bcd concentration is represented by green shaded colors. As the mitoses are repeated from nuclear cycle 10 to 14, the number of nuclei on the surface of embryo increases and the cortical layer becomes thick. Bcd is concentrated in nuclei which can be recognized as the bright green dots and is also accumulated in the cortical layer to exhibit the green shaded region. These features are similar to the observed snapshots in distribution of BcdeGFP (Figure 6A of Ref.11).
Profiles of distribution of Bcd and Hb and the distribution of the hb enhancer state are examined in Figure 10. In the left panels of Figure 10 (Figures 10A, 10D, 10G), results of 10 simulated embryos are superposed, each of which is obtained by monitoring nuclear sites along a stripe of the cylinder surface at an instance in the interphase period of nuclear cycle 14. In the middle panels of Figure 10 (Figures 10B, 10E, 10H), the embryotoembryo fluctuations of the Bcd distribution, , the Hb distribution, , and the hb enhancer state, are plotted, where and is the ensemble average over 10 embryos. We can see from Figure 10 that the embryotoembryo fluctuation of Bcd distribution is as small as 10% as in the observed data (Figure 5B of Ref.12). , , and can be transformed to the positional fluctuations as , , and , which are plotted in the right panels of Figure 10 (Figures 10C, 10F, 10I). Positional fluctuation of Bcd, , is comparable with the observed data (Figure 5C of Ref.12) showing that the Bcd gradient is generated in a reproducible way with small embryotoembryo fluctuation at the interphase of nuclear cycle 14. Fluctuations of the Hb distribution (Figure 10E) and the hb enhancer state (Figure 10H) are larger than fluctuation of Bcd (Figure 10B); At the threshold position of x/L?0.48, and amount to 40?50%. The positional fluctuations, and , however, are less than 5% (Figures 10F and 10I) because of the sharp change in and at the threshold position. Thus, Bcd, Hb, and the hb enhancer state are reproducible with the small embryotoembryo fluctuation.
We calculated the IO relation by taking the average over many nuclei in a single embryo at the interphase of nuclear cycle 14. Lines of the IO relation obtained from each of 10 runs of simulation are superposed in Figure 11A. These 10 lines do not deviate much from each other, which is the feature consistent with the experimental data (Figure 4A of Ref.12). Fluctuations in the IO relation are quantified in Figures 11B and 11C, which show a good agreement with the observed data (Figures 4B and 4C of Ref.12). Small values of each IO fluctuation shown by each line of Figures 11B and 11C imply that the Bcd gradient is read in a precise way in each embryo, and the small deviation of multiple IO lines from each other implies that the Bcd gradient is read out in a reproducible way with the small embryotoembryo fluctuation. In this way, results of the threedimensional model show that both the generation of the Bcd gradient and the readout of the Bcd gradient are precise and reproducible at the nuclear cycle 14 in spite of the large dynamical structure change of embryo through nuclear cycle 1?14.
Dependence of the IO relation on parameters has also similar features to those discussed with the onedimensional model (Figure S2). The fluctuation of the IO relation, ?_{h} at nuclear cycle 14, does not sensitively depend on ?_{b}, ?_{h} or N_{burst} as far as ?_{b} or ?_{h} is not extremely small. The IO fluctuation is sensitive, on the other hand, to V_{r} and D_{h}. Thus, the qualitatively same results as the onedimensional model were obtained with the threedimensional model, which supports the view that the fluctuation in the IO relation is dominated by the small value of V_{r}, and the fluctuation is masked by the small but finite D_{h}. When the positive feedback of Hb is turned off by posing h_{h}?=?0, the amplitude of N_{h} decreases as in the onedimensional model, but features of preciseness and reproducibility are not altered much (Figure S3), which implies that the positive feedback mechanism is not necessary to suppress the fluctuations.
This mechanism to resolve the paradox of signal interpretation can be visually confirmed by taking snapshots of the sideview of the model embryo (Figure S4). The hb enhancer state, ?, shows large spatial heterogeneity at the threshold region with about half of nuclei being randomly activated to produce Hb molecules for both two cases of D_{h}?=?0.1 ?m^{2}/s (Figure S4A) and D_{h}?=?0 (Figure S4B). When D_{h}?=?0.1 ?m^{2} s^{?1}, thus produced Hb diffuses over several nuclei, which smoothes out this heterogeneity to produce a boundary of the Hb distribution (Figure S4C), which is in sharp contrary to the case of D_{h}?=?0 for which the Hb distribution directly reflects the heterogeneity of the gene activity (Figure S4D).
Though hb mRNA is not explicitly taken into account in the present model, we may be able to expect that the noise level of the hb mRNA distribution is in between the noise level of distribution of the hb gene activity and that of the Hb distribution since mRNA first localizes in nuclei in which the hb gene is active but gradually diffuses out of nuclei to represent both characteristics of the localized hb gene activity and the more delocalized diffusive hb products. This expectation is consistent with the observed larger fluctuation of the hb mRNA distribution than the Hb distribution (Figure S8 of Ref.36).
With the threedimensional model, the simulated temporal change of the Bcd distribution during nuclear cycle 10?14 can be compared with the experimental data. In Figure 12A, the simulated Bcd concentration at the anterior part is compared with the experimental data ^{[11]}. Plotted are numbers of Bcd molecules in nuclear and cortical sites at the outmost layer of the cylinder in the slice i?=?10 (x/L?=?0.14). The Bcd concentration in nuclei is large during the interphase and sharply drops to the smaller values during mitoses. During the mitotic period, Bcd is released to the cortical layer and the Bcd concentration in cortical sites increases. In the observed data, the Bcd concentration in nuclei sharply rises in the beginning of the interphase period and then gradually decreases during the interphase period. Gregor et al. showed that this decrease in the Bcd concentration should be attributed to increase in size of each nucleus during the interphase period: With increase of the nuclear volume, Bcd is diluted in nuclei, which is partially offset by increase of incoming flux to the nuclei. Such dynamical balance among incoming/outgoing fluxes and dilution of Bcd should result in the gradual decrease of the Bcd concentration during the interphase period ^{[11]}. Since the growth of the nuclear size is not taken into account in the present model, decrease of Bcd concentration during the interphase period is not reproduced in simulation but the Bcd concentration gradually increases due to the active transport of Bcd into nuclei. Apart from such a difference, the overall features of the simulated results agree well with the observed data (Figure 3D of Ref.11). The Bcd concentration is stably reproduced in the next nuclear cycle even after Bcd is lost during the mitotic period. This stability is quantified by comparing the peak values of the Bcd concentration in successive nuclear cycles. In Figure 12B, the peak values in successive nuclear cycles are shown to differ only about 10%, showing that the Bcd concentration profile is as stable as was observed in experiment (Figure 3E of Ref.11). Such stable Bcd profile is realized because Bcd molecules released to cortical sites during the mitotic period are efficiently concentrated again into nuclei in the next nuclear cycle. The rapid transport of Bcd from cortex to nuclei is the necessary condition to prevent Bcd from escaping into core sites to achieve this efficient concentration. In Figure 12C, the simulated accumulation of Bcd in the cortical layer is compared with the experimental data (Figure 6D of Ref.11). The factor Q(m) defined in Table 2, which represents the gradual decrease of diffusion constant from cortex to core sites, was determined to reproduce this experimental plotting. See Method section for further explanation of Q(m). In Figure S5, the snapshots of the Bcd distribution are shown to visualize how Bcd is released in the mitotic period and is attracted into nuclei again in the interphase period.
The Hb distribution is also stable in dynamically changing embryo. By calculating the distribution of the number of Hb molecules along a strip on the cylinder for t?=?126?143 min in the interphase period of nuclear cycle 14, we can see that the number of Hb molecules increases in a coherent manner along the AP axis as Bcd is being concentrated in nuclei and the Hb synthesis is promoted (Figure S6A). The distribution profile of Hb is established soon after the Hb synthesis starts, so that the threshold position x_{1/2}/L?=?i_{half}?x/L of the Hb distribution only slightly shifts during nuclear cycle 14 (Figure S6B), which is consistent with the small shift of x_{1/2}/L in the experimental data (Figure 2C of Ref.36). The slope of the Hb distribution at x_{1/2} is kept steep during nuclear cycle 14 (Figure S6C), which roughly agrees with the experimental data (Figure 2D of Ref.36).
The experimentally observed IO relation around x_{1/2} has been fitted by the Hill relation as
with the Hill coefficient n?=?3?5 ^{[12]},^{[35]},^{[36]},^{[38]}. Results of the same fitting of the simulated data are shown in Figure 13A during nuclear cycle 11?14. Through these nuclear cycles, n stays constant in interphase of different nuclear cycles in spite of the repeated mitoses between interphase periods. We should note that such stable Hb distribution should be realized in dynamically changing embryo only when the stable and precise patterning of Bcd distribution is established.Temporal change of fluctuation of the IO relation in nuclear sites near the threshold position x_{1/2} is shown in Figure 13B: As shown in Figure 11B, the fluctuation of the IO relation shows a peak at x_{1/2}. Plotted in Figure 13B is the time course of the peak height of the fluctuation of the IO relation. When D_{h}?=?0, the fluctuation sharply increases at the onset of interphase in each nuclear cycle because multiplicated nuclei are placed at positions away from Hb molecules inherited from the previous nuclear cycle. As the busting production takes place multiple times, the Hb molecules are accumulated around nuclei and the Hb concentration approaches the steady level due to the balance between the Hb production and degradation, but the remaining fluctuation is still larger than the fluctuation observed in experiment ^{[12]}. When D_{h} is finite with D_{h}?=?0.1?0.5 ?m^{2}/s, on the other hand, the fluctuation is much suppressed and reaches the experimentally observed level ^{[12]}. Here, the important feature is that the fluctuation for nonzero D_{h} decreases by taking a long time from 50?35% in nuclear cycle 11 to 15?10% in nuclear cycle 14. This longtime decrease is due to the slow diffusion of Hb. See also Video S1 to grasp the intuitive picture of how the slow diffusion of Hb averages the Hb distribution through multiple nuclear cycles to reach the homogeneous distribution in the surface layer of a slice of embryo. Though the larger D_{h} is more effective to reduce the fluctuation, too large D_{h} smears out the Hb distribution along the AP axis and flattens the slope of the profile at the boundary of hb activation as shown in Figure 6 for the onedimensional model. To keep the sharpness of the Hb boundary, D_{h} has to be smaller than 0.3 ?m^{2}/s. As D_{h} is small, the selfaveraging process of Hb should take time longer than the duration of a single nuclear cycle and multiple cycles are needed to reduce the fluctuation. In this way, the stable Bcd and Hb profiles lasting through multiple nuclear cycles is the basis to assure the effective reduction of fluctuation through multiple nuclear cycles.
One and threedimensional models of Drosophila embryo were developed and two major issues raised by the experimental visualization of embryonic development ^{[11]},^{[12]} were examined with stochastic simulations: One of the issues is on the mechanism to generate the profile of Bcd gradient. We showed that the stochastic processes of synthesis, diffusion and degradation of Bcd give rise to the Bcd distribution whose entire profile remains stable through multiple nuclear cycles. This stable profile is precise to exhibit only the small fluctuation within each embryo and is reproducible with small fluctuation among multiple embryos. The stable profile of Bcd is realized by the rapid transport of Bcd to nuclei. The other issue is on the readout of the Bcd gradient. Random diffusion and reception of Bcd molecules at the nanometer scale region of DNA induce the intense fluctuation in the hb gene activity. Our models showed that this fluctuation indeed dominates fluctuation in the hb expression, but the fluctuation in the Hb distribution is masked by the slow diffusion of Hb molecules: The slow diffusion of Hb through multiple nuclear cycles averages the Hb distribution without losing the sharp boundary of distribution at around the threshold position. In this way, selfaveraging due to the slow diffusion of the output protein resolves the paradox of signal interpretation and enables the precise and reproducible readout of the gradient. Since the selfaveraging process of Hb has to continue over multiple nuclear cycles to achieve the sufficient accuracy, the stable Bcd distribution over multiple nuclear cycles is a necessary condition for realizing the accurate Hb distribution. Thus, the coordinated diffusion of input and output molecules is the basis to generate the stable, precise, and reproducible patterning of both input and output molecules.
Our prediction of the Hb selfaveraging mechanism should be tested experimentally. Because mobility of the fused protein of Hb and the other protein domain should be affected by the characteristics of the added domain, diffusion constant of the Hbfusion protein could be systematically changed by modulating the mass, size, and surface charge of the added domain. One possible example is to add a small tag such as histidine residues to eGFP, which can reduce velocity of the HbeGFP fusion protein ^{[45]}. Although such modulation may change the affinity of Hb to the hb enhancer to reduce effects of the feedback regulation of the hb expression, we should note that this feedback regulation has little effect on fluctuation in the Hb distribution as shown in Figure 7 and Figure S3, so that the effects of the Hb selfaveraging mechanism on fluctuation can be tested by systematically modulating the fused Hb proteins.
Gregor et al. ^{[12]} also discussed the possible resolution of the paradox of signal interpretation through the selfaveraging diffusion of Hb molecules. In their proposed mechanism, however, the positive feedback due to the binding of Hb protein to the hb enhancer was the necessary condition: They assumed that the hb genes in N?50 neighboring nuclei are activated synchronously with the positive feedback interaction of the diffusing Hb and hb enhancers, thereby the effective size of the target DNA region is enlarged from a to aN, which leads to the modification of the criterion of Equation 1 to be . In the present simulations, however, we have not found an evidence of synchronized gene switching in multiple nuclei, so that this ?selfaveraging of input to the feedback loop? proposed in Ref.12 is not necessary to resolve the paradox of signal interpretation. The paradox was resolved even in models without the positive feedback, which clearly shows that the simpler assumption of ?selfaveraging of output? is sufficient to resolve the paradox. In Drosophila embryos, molecules other than Bcd or Hb may also affect the hb gene activity ^{[18]},^{[46]}. The present results showed that assumption of the other input molecules is not necessary but the Bcd input and Hb output are sufficient to explain the stable, precise, and reproducible profiles of Bcd and Hb molecules at least in the anterior to middle regions of embryo.
Several different measures have been used to evaluate the accuracy of the Hb distribution ^{[6]}, ^{[9]}?^{[12]},^{[17]},^{[36]},^{[47]},^{[48]}. In this paper, we analyzed ?preciseness? ^{[12]} of the Hb distribution by evaluating the smallness of fluctuation in each embryo, and ?reproducibility? ^{[6]},^{[10]},^{[12]} by evaluating the smallness of fluctuation among multiple embryos. Other than these measures, ?sharpness? ^{[9]},^{[36]}, and ?robustness? ^{[17]},^{[47]},^{[48]} have been used to evaluate the accuracy. Our results showed that the sharp boundary of the Hb distribution is realized when D_{h} is within an optimal range: D_{h} should be smaller than 0.3 ?m^{2}/s to keep the average slope of the Hb boundary steep and D_{h} should be nonzero to suppress fluctuation around the average. In this way, the small fluctuation, or preciseness is a necessary condition for the sharpness. The present simulation showed that input of multiple kinds of transcriptional factors to hb is not necessary for the sharp Hb distribution. An experimental observation which is consistent with this result is the expression pattern of the reporter gene embedded into the zygotic genome ^{[9]}: The expressed pattern of the reporter gene which contains binding sites only for Bcd had a sharp boundary in embryos, which indicates that the other transcription factors than Bcd are not necessary for sharpness. This result is also consistent with the theoretical investigation: By solving the coupled equations of diffusion and reactions of Bialek and others ^{[30]}?^{[32]}, we can see that the physical limit of Equation 1 is not changed when the hb enhancer receives multiple input molecules. Thus, the multiple input molecules do not by themselves resolve the paradox of signal interpretation but the mechanism of selfaveraging diffusion is necessary for preciseness and hence for sharpness of the Hb distribution. Interactions among multiple genes, however, should be important for the Hb distribution to be robust against variations of conditions. Comparison between the observed and simulated data has suggested that the feedback regulations among gap genes are necessary for the robust Hb distribution ^{[47]},^{[48]}. This ?canalization? mechanism, therefore, should be a necessary condition for robustness but our analyses suggest that the canalization mechanism by itself is not sufficient to assure preciseness and sharpness.
We emphasize that the system described by the threedimensional model is not in a steady state, so that the absolute values of concentrations at each position vary during nuclear cycles. Effects of nonsteadiness were also focused on by Bergmann et al. ^{[15]},^{[49]}. Since the number of nuclei is small in the earlier stage of development, the fluctuation allowed to distinguish adjacent nuclei is about 50% in nuclear cycle 9 and 30% in nuclear cycle 11 by assuming an even distribution of nuclei during all nuclear cycles. Bergmann et al ^{[49]} suggested that the fate of each nucleus is determined as early as in nuclear cycle 9. Such presteady state decoding should tolerate the large fluctuation because of the large spatial distance between nuclei. Our simulation data showed that the system is not in the steady state all through nuclear cycles as pointed out by Bergmann et al, but the fluctuation in the simulated IO relation was larger than 30% in nuclear cycle 11 when D_{h}?=?0. The presteady state decoding, therefore, does not resolve the paradox of signal interpretation by itself, but the slow diffusion of Hb is necessary for the precise decoding.
The mechanism of diffusion of Bcd remains as an unsolved problem. From the experimental data of the fluorescence relaxation of BcdeGFP, diffusion constant was estimated to be D_{b}<1 ?m^{2}/s ^{[11]}, but D_{b}?=?10?20 ?m^{2}/s was suggested by monitoring diffusion of inert molecules in embryo ^{[10]}. One possible explanation of this discrepancy is based on the difference of environment in these measurements: The fluorescence relaxation giving small D_{b} was measured at the surface layer of embryo. Diffusion should be slow in such measurement because the complex and dense cytoskeletal structures in cortical layer trap the diffusing molecules ^{[50]},^{[51]}. D_{b}?10?20 ?m^{2}/s, on the other hand, was a result of diffusion through the unstructured core of embryo. The fluorescence relaxation in less structured unfertilized embryo, however, also showed D_{b}<1 ?m^{2}/s, which casts doubts on such simple interpretation based on the difference of environment ^{[11]}. Another possible explanation is based on the notion of anomalous diffusion. Particles diffusing through medium of crowded obstacles should show the diffusion as for the short distance, and for the long distance with ?<?^{[52]}. If both and are fitted by and , the effective diffusion constant, , should become much smaller than , and hence it is possible to assume and . We can also expect that the crossover distance between these two behaviors and other parameters such as and depend on the densities of cytoskeletal obstacles and take different values in cortical layer and core of embryos. Also possible is the mechanism that the cytoplasmic flow ^{[14]},^{[53]} enhances transport of molecules and increases the effective value of D_{b} for the long spatiotemporal scale. Either type of the longlived coherent flow ^{[14]} or the randomized flow should accelerate the movement of Bcd in embryo. Assumptions on values of diffusion constant in our models have not yet been checked from the microscopic viewpoint of the diffusion process, and further theoretical and experimental investigations are necessary to define diffusion constants on a sound basis.
The recent experimental report suggested that bcd mRNA is not strictly localized at the anterior pole but forms a gradient by the mechanism of quasirandom active transport through a nonpolar microtublular network in cortex of the embryo ^{[41]}. If this is the case, Bcd can form the gradient by following its mRNA gradient with the relatively small diffusion constant, so that the assumption of the anomalous diffusion of Bcd may not be necessary. Test of this mechanism of dictation of the bcd mRNA gradient is possible with the present one and threedimensional models by changing the definition of sites which can synthesize Bcd and adjusting diffusion constants of Bcd. With such modifications, the mechanism of Bcd gradient formation will be altered in the model but we expect that the main results of the paper, such as the stable Bcd distribution through the rapid transportation of Bcd into nuclei and the precise and reproducible Hb distribution through resolution of the paradox, will not be changed because these results should be independent of whether the Bcd gradient is formed by following the bcd mRNA gradient or it is formed by the localized synthesis and diffusion.
Another important problem is the scaling of distributions of developmental factors in embryo ^{[54]}. Several dipteran species have different size of embryos varying up to five fold, in which the distribution of Bcd scales with the length of embryo ^{[10]},^{[17]},^{[34]},^{[55]},^{[56]}. A possible mechanism of this scaling is the domination of degradation of Bcd at nuclei to determine the Bcd profile ^{[10]},^{[11]}. Because the positional distribution of nuclei scales with the length of embryo, one might consider that this assumption leads to the scaling of Bcd distribution. We examined this hypothesis with our threedimensional model, but the results indicated that the Bcd distribution does not scale even when nuclei distribution scales (data not shown). To solve this problem, the assumption of the larger rate of trapping to the nuclei might be needed ^{[57]}, or the scaling of transport or diffusion mechanism of Bcd or bcd mRNA should be necessary. Examination of these hypotheses is left for future investigations.
In both one and threedimensional models, stochastic simulations were performed with the Gillespie algorithm ^{[58]}. In each step in simulation with the Gillespie algorithm, one molecular process has to be picked out efficiently from a large number of diffusion and reaction processes defined in the model. Such efficient picking out was realized by sorting a large number of processes with the optimized direct method ^{[59]}.
Most of the standard parameter values were determined by consulting the relevant experimental observation, and several parameters such as J, g(3), h_{b} and h_{h} in the onedimensional model or J(m) and Q(m) in the threedimensional model were calibrated to reproduce the observed averaged values of distributions of Bcd and Hb. The calculated results with this model, therefore, were derived by using the experimental data for the averaged values of distributions, but the experimental data for fluctuations of distributions were not used in the definition of the model or in determination of parameter values. The main objective of the model was not to show the consistency of the model to explain the observed average values of distributions of Bcd and Hb, but to show that the model provides consistent data on precision to read out the Bcd gradient.
Standard parameters used in the onedimensional model are summarized in Table 1. These standard parameter values were determined by following the guideline mentioned above. To determine the standard value of D_{b} in Table 1, however, there remains a fundamental problem: When inert molecules whose molecular weight is comparable with Bcd were injected into embryo, the diffusion constant of those molecules obtained from the measurement for an hour was D_{b}?10?20 ?m^{2}/s ^{[10]}. This value of D_{b} can well explain the observed exponential gradient of Bcd in embryo ^{[10]}. When the diffusion constant was measured by the relaxation of fluorescence strength in the BcdeGFP distribution at the cortical cytoplasm for a minute, on the other hand, the obtained value was D_{b}<1 ?m^{2}/s ^{[11]}, which is far different from the value observed from the long spatiotemporal measurement. A possible explanation of this discrepancy is that Bcd diffuses anomalously because Bcd is frequently trapped by structures in cytoplasm during its diffusion ^{[50]},^{[51]}. This anomalous diffusion may make D_{b} dependent on the distance l to be measured: For l>10 ?m, D_{b} should be as large as 10?20 ?m^{2}/s, but for l<1 ?m, D_{b} should be small enough to be less than 1 ?m^{2}/s. Another possible explanation is that the cytoplasmic flow in embryo ^{[14]},^{[53]} effectively enhances the diffusion over long distance. Though either of these hypotheses has not yet been confirmed, we here rely on such possibilities and use D_{b}?=?10 ?m^{2}/s as a standard parameter to model the cytosolic diffusion process.
Since s varied in every simulation run, V_{r} is also varied as and the other volume dependent parameters are , , and . and in Table 1 are related to each other by . We consider that , where is a distance inside the nucleus, so that we use in Table 1. Diffusion constant of Hb in cytoplasm is assumed to be small because the observed data suggests that Hb diffuses for much smaller distance than Bcd ^{[36]}, and hence we assume that D_{h} is as small as D_{nuc}. We use D_{h}?=?0.1 ?m^{2}/s as a standard parameter.
Bcd is stable in cytoplasm but could be degraded in nucleus during the interphase period, so that we assume that the lifetime of Bcd is similar to or longer than the duration of the interphase. Hb may also have the same order of lifetime. We use k_{b}?=?7.7?10^{?4} s^{?1}?=?1/21.6 min (halflife period of 15 min) and k_{h}?=?3.9?10^{?4} s^{?1} in Table 1. The rate of synthesis of Bcd at the anterior pole J was adjusted to fit the calculated value of the Bcd concentration to the observed data. Considering that the diameter of nucleus is about 6.5 ?m, V_{nuc} is around 144 ?m^{3}. The observed concentration of Bcd in nuclei shows a peak of 55 nM near the anterior pole and decays to 8 nM at around the middle point of the AP axis and becomes almost 0 at around the posterior pole ^{[12]}, which correspond to 4700 molecules/nucleus at the anterior part, 700 molecules/nucleus at the middle part and almost 0 at the posterior part. In Table 1, J was chosen to fit the calculated N_{Bcd}(i) to these values. For N_{burst}, we use a typical value N_{burst}?=?10 found in the singlemolecule observation of a gene switch ^{[60]}. The synthesis rate of Hb becomes very low when Bcd does not bind to the hb enhancer ^{[40]}, so that g(0) and g(1) have to be much smaller than g(2) or g(3). In Table 1 we assume g(0)/g(3)?=?g(1)/g(3)?=?10^{?5}. The experimental data showed that in the embryo which lacks ability of Hb to bind the hb enhancer, the concentration of Hb at the anterior region is reduced to ?20% of that in the wild type ^{[36]},^{[40]}, so that we assume g(2)/g(3)?=?0.2 in the standard parameterization. As g(3)N_{burst}/k_{b} should give an estimation of N_{Hb}(i) around the nucleus of the state S(i)?=?3, we determined g(3) by adjusting the calculated results of N_{Hb}(i) at the anterior part to have the same order of magnitude as N_{Bcd}(i?=?1).
The rate of unbinding Bcd or Hb from the hb enhancer can be of order of the duration of the interphase, so that we determine f_{b} and f_{h} to be 1/3 min?5?10^{?3} s^{?1}. The binding affinity of Bcd and Hb to the hb enhancer is determined by f_{b}/h_{b} and f_{h}/h_{h}, respectively. Parameters h_{b} and h_{h}, therefore, determine which nuclei should be S(i)?=?2 or 3, and hence determine the region where hb actively expresses. We adjusted the values of h_{b} and h_{h} to fit the location of the sharp increase in the calculated N_{Hb}(i) to the observed threshold position of x/L?48% ^{[12]}.
Standard parameter values are summarized in Table 2. Diffusion constants in the threedimensional model depend on types of starting or arriving sites, which are explained in Figure 2C. In cortical sites, for example, microtubules or other skeletal structures are developed, which should trap molecules and slower the diffusive movement of molecules, and hence diffusion constant between cortical sites is set to be smaller than that between core sites. Since the cytoskeltal structure in cortical sites develops as embryo develops, Bcd should become more easily trapped to such structure in the later stage of nuclear cycles, which should decrease the diffusion constant D_{Bcd, corticalcore} in Figure 8C. We represent this tendency by introducing a decreasing factor Q(m) as a function of the nuclear cycle m as defined in Table 2. Q(m) was determined so as to reproduce the observed increase of concentration of Bcd at the cortical layer (Figure 6D of Ref.11). This gradually decreasing diffusion constant prevented Bcd from escaping into core sites to accumulate in the cortical layer. Values of other parameter in Table 2 were chosen according to the same guiding principle as in Table 1, but we make several comments on some specific features of parameters in the threedimensional model: We assume that Bcd is actively transported into nuclear, so that the rate constant for transport of Bcd into the nuclear sites are larger than that for transporting Bcd to the outside of the nuclear sites in Table 2. Bcd and Hb are degraded with the rate of and in cortical or core sites. In the nuclear sites, the degradation rates are and at the outside of the interaction volume, and and at the inside of the interaction volume. Here, and are assumed to be small in cortical or core sites but large in nuclear sites. At the core sites in the slice i?=?1, Bcd is generated with the rate J(m) as in the onedimensional model. Since the observed production rate of Bcd increases as nuclear cycles proceed ^{[42]}, we assume that J(m) is an increasing function of m. The value of J(14) is adjusted so that the calculated results reproduce the observed concentration profile of Bcd at nuclear cycle 14.
The hb enhancer state in the threedimensional model is S(i, j, t)?=?0?3. Also assumed is that Hb is synthesized at around nuclei, i.e. in the nuclear site with the rate g(S(i, j, t)). Values of g(S(i, j, t)) can be determined in the same way as Equation 7 or Equation 8. In Table 2, g(3) was adjusted so that the calculated number of Hb molecules becomes of same order as that of Bcd molecules.
One and ThreeDimensional Models.
(0.28 MB PDF)
Click here for additional data file (pcbi.1000486.s001.pdf)
Figure S1
Noise in the Input/Output Relation in the OneDimensional Model with Varied Frequencies of the State Change of the hb Enhancer. ?_{h}(i)/max[<N_{h}(i)>] is plotted in the two dimensional plane of the relative rates of gene switching, ?_{b}?=?f_{b}/k_{b} and ?_{h}?=?f_{h}/k_{h}. V_{r}?=?10^{?2.5} is used. Other parameters are set to have standard values.
(0.12 MB PDF)
Click here for additional data file (pcbi.1000486.s002.pdf)
Figure S2
Input/Output Relations and Their Noise in the ThreeDimensional Model with Varied Parameters. Parameterization is changed by choosing one parameter from Table 2 to be varied from its standard value. For one parameterization, lines of 10 simulated runs are calculated as in each panel of Figure 11 and then these 10 lines are averaged to obtain one curve. Thus calculated averaged curves for different parameterizations are superposed. (Top row) Mean input/output relations. (Middle row) The standard deviation of input/output relations. (Bottom row) Translation of the data of Middle row into lines of equivalent input noise. (Left) The relative size of the interaction volume is varied as V_{r}?=?10^{?1} (navy), 10^{?1.5} (dark green), 10^{?2} (light green), 10^{?2.5} (blue, standard value), 10^{?2.9} (orange), 10^{?3} (magenta), and 10^{?3.5} (red). (Second left) The dissociation rate of Bcd from the hb enhancer is varied as ?_{b}?=?f_{b}/k_{b, nuclear}?=?10^{4} (dark green), 10^{3} (light green), 10^{2} (blue, standard value), 10 (orange), 1 (magenta), and 10^{?1} (red). (Middle) The dissociation rate of Hb from the hb enhancer is varied as ?_{h}?=?f_{h}/k_{h, nuclear}?=?10^{4} (dark green), 10^{3} (light green), 10^{2} (blue, standard value), 10 (orange), 1 (magenta), and 10^{?1} (red). (Second Right) The number of Hb molecules synthesized in a burst is varied as N_{burst}?=?100 (dark green), 50 (light green), 10 (blue, standard value), and 1 (red). The frequency of bursting, g(S), is also varied to keep g(S)N_{burst} constant. (Right) The diffusion constant of Hb is varied as D_{h}?=?5 ?m^{2}/s (navy), 1 ?m^{2}/s (dark green), 0.5 ?m^{2}/s (light green), 0.25 ?m^{2}/s (light blue), 0.1 ?m^{2}/s (blue, standard value), 0.05 ?m^{2}/s (orange), 0.01 ?m^{2}/s (magenta), and 0 (red).
(0.26 MB PDF)
Click here for additional data file (pcbi.1000486.s003.pdf)
Figure S3
Input/Output Relations and Their Noise in the ThreeDimensional Model with Varied Parameters for the Case of No Feedback Regulation of hb. Parameterization is changed by choosing one parameter from Table 2 to be varied from its standard value. For one parameterization, lines of 10 simulated runs are calculated and then these 10 lines are averaged to obtain one curve. Thus calculated averaged curves for different parameterizations are superposed. The same plot as in Figure S2 but with h_{h}?=?0.
(0.26 MB PDF)
Click here for additional data file (pcbi.1000486.s004.pdf)
Figure S4
SideView of Simulated Embryos. Snapshots of the hb expression pattern at t?=?95 min, which is in the interphase of nuclear cycle 12, are shown by highlighting nuclear sites of large ? with white dots. (A) D_{h}?=?0.1 ?m^{2}/s, and (B) D_{h}?=?0. Snapshots of the distribution pattern of Hb molecule are shown by green dots (C) D_{h}?=?0.1 ?m^{2}/s at the same instance as in A, and (D) D_{h}?=?0 at the same instance as in B.
(0.72 MB PDF)
Click here for additional data file (pcbi.1000486.s005.pdf)
Figure S5
Import and Export Dynamics of Nuclear Bcd. Snapshots of temporal changes of Bcd concentration at the slice i?=?5 of the threedimensional model. Bcd concentration is expressed by green shaded color. Nuclear cycle 13 starts from 103 min 00 sec and lasts until 119 min 00 sec. Bcd molecules start to accumulate in nuclei at the onset of the interphase (103 min 00 sec), and quickly flow out of nuclei at the onset of mitosis in nuclear cycle 14 (119 min 00 sec).
(0.17 MB PDF)
Click here for additional data file (pcbi.1000486.s006.pdf)
Figure S6
Development of Hb Profiles During Nuclear Cycle 14. (A) Numbers of Hb molecules averaged over 10 simulation runs, <N_{h, nuclear}(i,100,t)>, are shown at every one minute during t?=?126?143 min in nuclear cycle 14 as functions of x/L?=?i?x/L. The j?=?100th sites are on the surface of the cylinder. (B) Position of x_{1/2}/L?=?i_{half}?x/L is plotted for t?=?126?143 min in nuclear cycle 14 (blue line) Red dots are the experimental data of x_{1/2}/L (Figure 2C of Ref.36). (C) Angle of slope of <N_{h, nuclear}(i_{half},100,t)> is plotted at every one minute for t?=?126?143 min in nuclear cycle 14 with the standard parameterization (blue dots) and with the parametrization of h_{h}?=?0 (green dots). The lack of positive feedback with h_{h}?=?0 makes the angle slightly smaller. Red dots are the experimental data of angle at x_{1/2}/L (Figure 2D of Ref.36).
(0.17 MB PDF)
Click here for additional data file (pcbi.1000486.s007.pdf)
Video S1
Temporal Change of Distributions of Bcd and Hb in the ThreeDimensional Model. Concentrations of Bcd and Hb at the i?=?5th slice are shown in the green shaded color. The simulated 80minute process from nuclear cycle 8 to nuclear cycle 14 is shown in oneminute video run. (Left) Bcd concentration, and (Right) Hb concentration.
(0.89 MB WMV)
Click here for additional data file (pcbi.1000486.s008.wmv)
Notes
The authors have declared that no competing interests exist.
This work was supported by GrantinAid for Scientific Research in Priority Area from MEXT and GrantinAid for Scientific Research (A) from JSPS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Y.O. and H.M. thank Dr. Aude Porcher for fruitful discussions on experiments of Drosophila embryos. The authors thank Prof. Motonori Ota for discussions on the possible experimental verification of our results.
References
1.  Wolpert L. Year: 1969Positional information and the spatial pattern of cellular differentiation.J Theor Biol251474390734 
2.  Crick F. Year: 1970Diffusion in embryogenesis.Nature2254204225411117 
3.  Driever W,N?ssleinVolhard C. Year: 1988A gradient of Bicoid protein in Drosophila embryos.Cell5483933383244 
4.  Driever W,N?ssleinVolhard C. Year: 1988The Bicoid protein determines position in the Drosophila embryo.Cell54951043383245 
5.  Struhl G,Struhl K,Macdonald PM. Year: 1989The gradient morphogen Bicoid is a concentrationdependent transcriptional activator.Cell57125912732567637 
6.  Houchmandzadeh B,Wieschaus E,Leibler S. Year: 2002Establishment of developmental precision and proportions in the early Drosophila embryo.Nature41579880211845210 
7.  Ephrussi A,St Johnston D. Year: 2004Seeing is believing: The Bicoid morphogen gradient matures.Cell11614315214744427 
8.  Peel AD,Chipman AD,Akam M. Year: 2005Arthropod segmentation: beyond the Drosophila paradigm.Nat Rev Genet690591616341071 
9.  Crauk O,Dostatni N. Year: 2005Bicoid determines sharp and precise target gene expression in the Drosophila embryo.Curr Biol151888189816271865 
10.  Gregor T,Bialek W,de Ruyter van Steveninck RR,Tank DW,Wieschaus EF. Year: 2005Diffusion and scaling during early embryonic pattern formation.Proc Natl Acad Sci USA102184031840716352710 
11.  Gregor T,Wieschaus EF,McGregor AP,Bialek W,Tank DW. Year: 2007Stability and nuclear dynamics of the Bicoid morphogen gradient.Cell1301415217632061 
12.  Gregor T,Tank DW,Wieschaus EF,Bialek W. Year: 2007Probing the limits to positional information.Cell13015316417632062 
13.  Coopey M,Berezhkovskii AM,Kim Y,Boettiger AN,Shvartsman SY. Year: 2007Modeling the bicoid gradient: Diffusion and reversible nuclear trapping of a stable protein.Dev Biol31262363018001703 
14.  Hecht I,Rappel WJ,Levine H. Year: 2009Determining the scale of the Bicoid morphogen gradient.Proc Natl Acad Sci USA1066 1710171519118201 
15.  Bergmann S,Sandler O,Sberro H,Shnider S,Schejter E,et al. Year: 2007Presteadystate decoding of the Bicoid morphogen gradient.PLoS Biol5e4617298180 
16.  Reeves GT,Muratov CB,Sch?pbach T,Shvartsman SY. Year: 2006Quantitative models of developmental pattern formation.Dev Cell1128930016950121 
17.  Umulis D,O'Connor MB,Othmer HG. Year: 2008Robustness of embryonic spatial patterning in Drosophila melanogaster.Curr Topics in Dev Biol8165111 
18.  Howard M,ten Wolde PR. Year: 2005Finding the center reliably: Robust patterns of developmental gene expression.Phys Rev Lett9520810316384103 
19.  Tostevin F,ten Wolde PR,Howard M. Year: 2007Fundamental limits to position determination by concentration gradients.PLoS Comput Biol3e7817465676 
20.  Lepzelter D,Wang J. Year: 2008Exact probabilistic solution of spatialdependent stochastics and associated spatial potential landscape for the bicoid protein.Phys Rev E77041917 
21.  Raser JM,O'Shea EK. Year: 2005Noise in gene expression: origins, consequences, and control.Science3092010201316179466 
22.  K?rn M,Elston TC,Blake WJ,Collins JJ. Year: 2005Stochasticity in gene expression: From theories to phenotypes.Nat Rev Genet645146415883588 
23.  Paulsson J. Year: 2004Summing up the noise in gene networks.Nature42741541814749823 
24.  Okabe Y,Yagi Y,Sasai M. Year: 2007Effects of the DNA state fluctuation on singlecell dynamics of selfregulating gene.J Chem Phys12710510717867786 
25.  Yoda M,Ushikubo T,Inoue W,Sasai M. Year: 2007Roles of noise in single and coupled multiple genetic oscillators.J Chem Phys12611510117381233 
26.  Ushikubo T,Inoue W,Yoda M,Sasai M. Year: 2006Testing the transition state theory in stochastic dynamics of a genetic switch.Chem Phys Lett430139143 
27.  Sasai M,Wolynes PG. Year: 2003Stochastic gene expression as a many body problem.Proc Natl Acad Sci USA1002374237912606710 
28.  Walczak AM,Sasai M,Wolynes PG. Year: 2005Self consistent proteomic field theory of stochastic gene switches.Biophys J8882885015542546 
29.  Berg HC,Purcell EM. Year: 1977Physics of chemoreception.Biophys J20193219911982 
30.  Bialek W,Setayeshgar S. Year: 2005Physical limits to biochemical signaling.Proc Natl Acad Sci USA102100401004516006514 
31.  Tkacik G,Bialek W. Year: 2007Diffusion, dimensionality and noise in transcriptional regulation.arXiv0712.1852 [qbio.MN]. 
32.  Bialek W,Setayeshgar S. Year: 2008Cooperativity, sensitivity, and noise in biochemical signaling.Phys Rev Lett10025810118643705 
33.  Fowlkes CC,Hendriks CLL,Ker?nen SVE,Weber GH,R?bel O,et al. Year: 2008A quantitative spatiotemporal atlas of gene expression in the Drosophila blastoderm.Cell13336437418423206 
34.  Gregor T,McGregor AP,Wieschaus EF. Year: 2008Shape and function of the Bicoid morphogen gradient in dipteran species with different sized embryos.Develop Biol31635035818328473 
35.  Ma XG,Yuan D,Diepold K,Scarborough T,Ma J. Year: 1996The Drosophila morphogenetic protein Bicoid binds DNA cooperatively.Development122119512068620846 
36.  Lopes FJP,Vieira FMC,Holloway DM,Bisch PM,Spirov AV. Year: 2008Spatial bistability generates hunchback expression sharpness in the Drosophila embryo.PLoS Comput Biol4e100018418818726 
37.  Driever W,Thoma G,NussleinVolhard C. Year: 1989Determination of spatial domains of zygotic geneexpression in the Drosophila embryo by the affinity of bindingsites for the Bicoid morphogen.Nature3403633672502714 
38.  Burz DS,RiveraPomar R,J?ckle H,Hanes SD. Year: 1998Cooperative DNAbinding by Bicoid provides a mechanism for thresholddependent gene activation in the Drosophila embryo.EMBO J17599860099774343 
39.  H?lskamp M,Lukowitz W,Beermann A,Glaser G,Tautz D. Year: 1994Differential regulation of target genes by different alleles of the segmentation gene Hunchback in Drosophila.Genetics1381251348001780 
40.  SimpsonBrose M,Treisman J,Desplan C. Year: 1994Synergy between the Hunchback and Bicoid morphogens is required for anterior patterning in Drosophila.Cell788558658087852 
41.  Spirov A,Fahmy K,Schneider M,Frei E,Noll M,et al. Year: 2009Formation of the bicoid morphogen gradient: an mRNA gradient dictates the protein gradient.Development13660561419168676 
42.  Salles FJ,Lieberfarb ME,Wreden C,Gergen JP,Strickland S. Year: 1994Coordinate initialtion of Drosophila developement by regulated polyadenylation of maternal messenger RNAs.Science266199619997801127 
43.  Foe VE,Alberts BM. Year: 1983Studies of nuclear and cytoplasmic behaviour during the five mitotic cycles that precede gastrulation in Drosophila embryogenesis.J Cell Sci6131706411748 
44.  Foe VE,Odell GM,Edgar BA. Year: 1993Mitosis and morphogenesis in the Drosophila embryo: point and counterpoint. The development of Drosophila melanogaster,Bate M,MartinezArias ANYCold Spring Harbor Laboratory Press149300 
45.  Elowitz MB,Surette MG,Wolf PE,Stock JB,Leibler S. Year: 1999Protein mobility in the cytoplasm of Escherichia coli.J Bacteriol1811972039864330 
46.  McHale P,Rappel WJ,Levine H. Year: 2006Embryonic pattern scaling achieved by oppositely directed morphogen gradients.Phys Biol310712016829697 
47.  Manu,Surkova S,Spirov AV,Gursky VV,Janssens H,et al. Year: 2009Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation.PLoS Biol7e1000049 
48.  Manu,Surkova S,Spirov AV,Gursky VV,Janssens H,et al. Year: 2009Canalization of gene expression and domain shifts in the Drosophila blastoderm by dynamical attractors.PLoS Comput Biol5e100030319282965 
49.  Bergmann S,Tamari Z,Schejter E,Shilo BZ,Barkai N. Year: 2008Reexamining the stability of the Bicoid morphogen gradient.Cell132151718191212 
50.  Seksek O,Biwersi J,Verkman AS. Year: 1997Translational diffusion of macromoleculesized solutes in cytoplasm and nucleus.J Cell Biol1381311429214387 
51.  LippincottSchwartz J,Snapp E,Kenworthy A. Year: 2001Studying protein dynamics in living cells.Nat Rev Mol Cell Biol244445611389468 
52.  Bouchaud JP,Georges A. Year: 1990Anomalous diffusion in disordered media ? statistical mechanics, models, and physical applications.Phys Rep195127293 
53.  von Dassow G,Schubiger G. Year: 1994How an actin network might cause fountain streaming and nuclear migration in the syncytial Drosophila embryo.J Cell Biol127163716537798318 
54.  West GB,Brown JH,Enquist BJ. Year: 1997A general model for the origin of allometric scaling laws in biology.Science2761221269082983 
55.  Lott SE,Kreitman M,Palsson A,Alekseeva E,Ludwig MZ. Year: 2007Canalization of segmentation and its evolution in Drosophila.Proc Natl Acad Sci USA104109261093117569783 
56.  Sommer R,Tautz D. Year: 1991Segmentation gene expression in the housefly Musca domestica.Development1134194301685986 
57.  Umulis DM. Year: 2009Analysis of dynamic morphogen scale invariance.J R Soc Interface 2009.0015 
58.  Gillespie DT. Year: 1977Exact stochastic simulation of coupled chemical reactions.J Phys Chem8123402361 
59.  Gibson MA,Bruck J. Year: 2000Efficient exact stochastic simulation of chemical systems with many species and many channels.J Phys ChemA 10418761889 
60.  Cai L,Friedman N,Xie XS. Year: 2006Stochastic protein expression in individual cells at the single molecule level.Nature44035836216541077 
Figures
[Figure ID: pcbi1000486g001] 
doi: 10.1371/journal.pcbi.1000486.g001. Figure 1 Interaction volume. Effects of the random arrival of Bcd at the hb enhancer on DNA can be intuitively explained with the concept of interaction volume. Since molecules diffuse over a distance during the allowed time ?, Bcd molecules which can interact with the target of size a during the period ? are distributed in the volume of size , which is called the interaction volume. The interaction volume is small when a is in a nanometer lengthscale. With the average Bcd concentration c, the average number of Bcd molecules in the interaction volume is . If Bcd molecules diffuse as independent particles, fluctuation in n should obey the Poisson statistics to exhibit . Then, the Bcd concentration sensed at the target fluctuates to the extent at least . Smallness of leads to the large . 
[Figure ID: pcbi1000486g002] 
doi: 10.1371/journal.pcbi.1000486.g002. Figure 2 Onedimensional model of Drosophila embryo. (A) The onedimensional system of length L is composed of 100 sites each of which has the length of ?x?=?L/100. In each site diffusion and reactions of Bcd molecules (brown spheres) and Hb molecules (yellow spheres) are defined. Bcd diffuses from or into adjacent sites, degrades, and diffuses from or into the small region of interaction volume (dark pink, whose relative size is enlarged in the figure for explanation but its volume is about 1/1000 of nuclear volume). Hb is synthesized at each site and diffuses from or into adjacent sites, degrades, and diffuses from or into the small region of interaction volume. Bcd is synthesized only at the most anterior site. The rate of each process is written on each arrow. (B) The rate of Hb synthesis g(S) depends on the state of the hb enhancer, S?=?0?3; the state with no bound Bcd or Hb (S?=?0), the state with two bound Hb molecules and with no bound Bcd (S?=?1), the state with six bound Bcd molecules and with no bound Hb (S?=?2), and the state with six bound Bcd and two bound Hb molecules (S?=?3). g(S) satisfies Equation 7 or Equation 8. 
[Figure ID: pcbi1000486g003] 
doi: 10.1371/journal.pcbi.1000486.g003. Figure 3 Profiles of Bcd, Hb, and the hb state in the onedimensional model. All data (A?E) are plotted along positions on the AP axis, and all statistics are calculated from the ensemble of 100 simulation runs. (A) Bcd profiles in the ith site, , for 100 simulation runs are superposed and their average is shown with the yellow line. The observed data of the number of Bcd molecules in nuclei are plotted with red dots, which are averages taken over the ensemble of 15 embryos (Figure 5A of Ref.12). (B) Hb profiles in the ith site, , for 100 simulation runs are superposed and their average is shown with the yellow line. (C) Profiles of the hb enhencer state, , for 100 simulation runs are superposed and their average is shown with the yellow line. (D) Profiles of Bcd inside the interaction volume in the ith site, , for 100 simulation runs are superposed and their average is shown with the yellow line. (E) Profiles of Hb inside the interaction volume in the ith site, , for 100 simulation runs are superposed and their average is shown with the yellow line. 
[Figure ID: pcbi1000486g004] 
doi: 10.1371/journal.pcbi.1000486.g004. Figure 4 Input/output relation and its noise in the onedimensional model with the standard parameterization. (A) The mean output is plotted as a function of logarithm of the mean input, (blue line), where is the maximal value of and i_{half} is defined by . Red dots are the experimental data of the mean input/output relation averaged over nine embryos (Figure 4A of Ref.12). (B) Standard deviation of Hb levels for 100 simulation runs, , is plotted as a function of the input Bcd level (blue line). Red dots are the experimental data (Figure 4B of Ref.12). (C) Translation of the output noise of (B) into an equivalent input noise (blue line). Red dots are the experimental data (Figure 4C of Ref.12). 
[Figure ID: pcbi1000486g005] 
doi: 10.1371/journal.pcbi.1000486.g005. Figure 5 Input/output relations and their noise in the onedimensional model with varied parameters. The same plot as in Figure 4 but with different parameterizations. For each simulation, one parameter is chosen from Table 1 to be varied from its standard value. (Left) The relative size of the interaction volume is varied as V_{r}?=?10^{?1} (navy), 10^{?1.5} (dark green), 10^{?2} (light green), 10^{?2.5} (light blue), 10^{?2.9} (blue, standard value), 10^{?3} (magenta), and 10^{?3.5} (red). (Second Left) The dissociation rate of Bcd from the hb enhancer is varied as ?_{b}?=?f_{b}/k_{b}?=?10^{4} (dark green), 10^{3} (light green), 10^{2} (light blue), 10 (blue, standard value), 1 (magenta), and 10^{?1} (red). (Middle) The dissociation rate of Hb from the hb enhancer is varied as ?_{h}?=?f_{h}/k_{h}?=?10^{4} (dark green), 10^{3} (light green), 10^{2} (light blue), 10 (blue, standard value), 1 (magenta), and 10^{?1} (red). (Second Right) The number of Hb molecules synthesized in a burst is varied as N_{burst}?=?100 (dark green), 50 (light green), 10 (blue, standard value), and 1 (red). The frequency of bursting, g(S), is also varied to keep g(S)N_{burst} constant. (Right) The diffusion constant of Hb is varied as D_{h}?=?5 ?m^{2}/s (navy), 1 ?m^{2}/s (dark green), 0.5 ?m^{2}/s (light green), 0.25 ?m^{2}/s (light blue), 0.1 ?m^{2}/s (blue, standard value), 0.05 ?m^{2}/s (orange), 0.01 ?m^{2}/s (magenta), and 0 (red). 
[Figure ID: pcbi1000486g006] 
doi: 10.1371/journal.pcbi.1000486.g006. Figure 6 The input/output noise and the mean Hb profile in the onedimensional model with varied parameters. The noise in the Input/Output relation and features of the mean Hb profile are plotted in the two dimensional plane of the Hb diffusion constant D_{h} and the relative size of the interaction volume V_{r}. Other parameters are set to have standard values. (A) The maximal value of the standard deviation of Hb levels, evaluated from the ensemble of 100 simulation runs, is plotted in the contour map. (B) at the threshold position x?=?x_{1/2}?=?i_{half}?x. (C) Angle ? of the Hb profile at the threshold position x?=?x_{1/2}. Here, ? is calculated as , and is a function of x?=?i?x obtained by fitting it to the calculated data of . 
[Figure ID: pcbi1000486g007] 
doi: 10.1371/journal.pcbi.1000486.g007. Figure 7 Input/output relations and their noise in the onedimensional model with varied parameters for the case of no feedback regulation of hb. The same plot as in Figure 5 but Hb is assumed not to bind to the hb enhancer with h_{h}?=?0. 
[Figure ID: pcbi1000486g008] 
doi: 10.1371/journal.pcbi.1000486.g008. Figure 8 Threedimensional model of Drosophila embryo. (A) An embryo is modeled by a cylinder of radius R?90 ?m and length L?490 ?m extending along the anteriorposterior axis. The cylinder is composed of 70 plates, each of which has thickness of ?x?=?L/70 and consists of 625 hexagonal sites: The system consists of 70?625?=?43750 sites in total. These sites are classified into three types; nuclear sites (red), cortical sites (yellow), and core sites (white). Sites in the most anterior slice are core sites, from which Bcd is synthesized (blue). (B) Closeup of a part of the most anterior sites and sites in the next slice in the model. (C) Diffusion constants of molecule X can have different values D_{X, YZ}, depending on the type of starting site Y, and the type of arriving site Z. Types of sites are represented by the same colors as in A and B. We write X?=?b for Bcd and X?=?h for Hb molecules. D_{X, corenuclear} is denoted by D_{X, import} and D_{X, nuclearcore} is denoted by D_{X, export}. Values of these constants are listed in Table 2. 
[Figure ID: pcbi1000486g009] 
doi: 10.1371/journal.pcbi.1000486.g009. Figure 9 Snapshots of the Bcd distribution in the threedimensional model. Concentration of Bcd in each site is expressed by green shaded color. Slices at i?=?3, 7, 16, 32, and 60 are shown at nuclear cycle 10?14. 
[Figure ID: pcbi1000486g010] 
doi: 10.1371/journal.pcbi.1000486.g010. Figure 10 Reproducibility of profiles of Bcd, Hb, and the hb state in the threedimensional model. All data are plotted along positions on the AP axis. In simulation, all statistics are calculated from the ensemble of 10 simulation runs. The compared experimental statistics of Ref.12 were obtained from the ensemble of 15 embryos. (A) Bcd profiles in a stripe of the j?=?100th sites in slices i?=?1?70 at time t?=?144 min. Data of 10 simulation runs are superposed. The j?=?100th sites are on the surface of the cylinder. The average of these 10 curves, , is drawn with the yellow line. The observed data of averaged Bcd numbers in nuclei (Figure 5A of Ref.12) is plotted with red dots. (B) The standard deviation of is plotted in the normalized form as (blue line). The observed standard deviation (Figure 5B of Ref.12) is plotted with red dots. (C) Variability of Bcd profiles is translated into an effective positional fluctuation (blue dots). The observed positional fluctuation (Figure 5C of Ref.12) is plotted with red dots. (D) Hb profiles in a stripe of the j?=?100th site in slices i?=?1?70 at time t?=?144 min for 10 simulation runs are superposed. The average of these 10 curves, , is drawn with the yellow line. (E) The standard deviation of is plotted in the normalized form as . (F) Variability of Hb profiles is translated into an effective positional fluctuation . (G) Profiles of hb gene state, in a stripe of the j?=?100th site in slices i?=?1?70 at time t?=?144 min for 10 simulation runs are superposed. The average of these 10 curves, , is drawn with the yellow line. (H) The standard deviation of is plotted in the normalized form as . (I) Variability of profiles of hb state is translated into an effective positional fluctuation . Horizontal lines in middle and right panels are drawn for guides of eyes. 
[Figure ID: pcbi1000486g011] 
doi: 10.1371/journal.pcbi.1000486.g011. Figure 11 Input/output relations and their noise in the threedimensional model with the standard parameterization. (A) Mean input/output relations for 10 simulation runs are superposed (blue lines). A mean input/output relation is obtained by plotting in one simulation as a function of for all nuclei (about 6000 nuclei) in the system at t?=?144 min and by binning them into 50 bins. Red dots are the experimental data (Figure 4A of Ref.12). (B) For each bin from (A) the standard deviation of is calculated. Drawn are 10 lines of the standard deviation for 10 simulation runs. Red dots are experimental data of standard deviation of input/output relations for nine embryos (Figure 4B of Ref.12). (C) Translation of 10 lines of the output noise of (B) into lines of equivalent input noise (blue lines). Red dots are the experimental data (Figure 4C of Ref.12). 
[Figure ID: pcbi1000486g012] 
doi: 10.1371/journal.pcbi.1000486.g012. Figure 12 Temporal development of the Bcd distribution. (A) Nuclear and cytoplasmic development of the number of Bcd molecules from nuclear cycle 10 to 14. Numbers of Bcd molecules in nuclear sites at the i?=?10th slice (x/L?=?0.14), , are averaged over j in the outmost layer of the cylinder and over 10 simulation runs (blue line), and numbers of Bcd molecules in cortical sites at the i?=?10th slice (x/L?=?0.14), , are averaged over j in the outmost layer of the cylinder and over 10 simulation runs (magenta line). Experimental data of Bcd concentrations (Figure 3D of Ref.11) at nuclei (blue dots) and at cortex (red dots) are compared. (B) Scatter plot of peak numbers of Bcd molecules in successive nuclear cycles. The maximum value of during nuclear cycle n is recorded as for slice i?=?5, 10, 15?, 65, and then the recorded data are plotted at the coordinate ( , ) in the twodimensional plane for n?=?8?10. The j?=?100th sites are on the surface of the cylinder. Line with the slope one corresponds to the perfect stability across nuclear cycles and dotted lines correspond to 10% deviation. (C) Development of ratio of concentrations. Numbers of Bcd molecules averaged over the whole system and over the time duration of interphase of nuclear cycle n are denoted by , , and , for nuclear, cortical and core sites, respectively. Plotted are (blue line), and (green line). 10 lines for 10 simulation runs are superposed. The experimental data for (Bcd concentration in cortex)/(Bcd concentration in core) are plotted with red open circles and the experimental data for (Bcd concentration in nuclei)/(Bcd concentration in cortex) are plotted with red dots, both of which are taken from Figure 6D of Ref.11. 
[Figure ID: pcbi1000486g013] 
doi: 10.1371/journal.pcbi.1000486.g013. Figure 13 Development of the Hill coefficient and input/output noise. (A) Temporal change of the Hill coefficient n obtained by fitting the Hill relation to the input/output relation of Figure 11A. n is put to be 0 during mitoses. (B) Temporal change of the peak height of the average standard deviation of input/output relations. The standard deviation is obtained from ensemble of nuclei in the single embryo as lines in Figure 11B, and the average was taken for 10 data of the standard deviation of 10 simulation runs. The peak height is put to be 0 during mitoses. D_{h}?=?0.5 ?m^{2}/s (blue), 0.1 ?m^{2}/s (green, standard value), and 0 (red). 
Tables
Table 1 Standard Parameterization in the OneDimensional Model.
Parameter  Standard value  
J  Bcd synthesis rate  84.7 s^{?1} 
D_{b}  Bcd diffusion constant  10 ?m^{2} s^{?1} 
k_{b}  Bcd degradation rate constant  7.7?10^{?4} s^{?1} 
<h_{b}>  Bcd binding rate constant to the hb enhancer  0.19 s^{?1} 
f_{b}  Bcd dissociation rate constant from the hb enhancer  7.7?10^{?3} s^{?1} 
D_{h}  Hb diffusion constant  0.1 ?m^{2} s^{?1} 
k_{h}  Hb degradation rate constant  3.85?10^{?4} s^{?1} 
<h_{h}>  Hb binding rate constant to the hb enhancer  5.8?10^{?2} s^{?1} 
f_{h}  Hb dissociation rate constant from the hb enhancer  3.9?10^{?3} s^{?1} 
g(0)  Hb synthesis rate at S?=?0  1.25?10^{?6} s^{?1} 
g(1)  Hb synthesis rate at S?=?1  1.25?10^{?6} s^{?1} 
g(2)  Hb synthesis rate at S?=?2  2.5?10^{?2} s^{?1} 
g(3)  Hb synthesis rate at S?=?3  1.25 s^{?1} 
N_{burst}  Burst size of Hb synthesis  10 molecules 
v_{ex}  Export rate constant of Hb or Bcd to interaction volume  0.4 s^{?1} 
<v_{im}>  Import rate constant of Hb or Bcd from interaction volume  5.04?10^{?4} s^{?1} 
<?x>  The average mesh size  5 ?m 
<V_{r}>  The average ratio of interaction volume to nuclear volume  10^{?2.9} 
Table 2 Standard Parameterization in the ThreeDimensional Model.
Parameter  Standard value  
J  Bcd synthesis rate  J(m) s^{?1} 
D_{b,import}  Bcd diffusion constant (import into nucleus)  1.2 ?m^{2} s^{?1} 
D_{b.export}  Bcd diffusion constant (export from nucleus)  0.3 ?m^{2} s^{?1} 
D_{b,corecotex}  Bcd diffusion constant (from core to cortex)  12.5 ?m^{2} s^{?1} 
D_{b,cortexcore}  Bcd diffusion constant (from cortex to core)  12.5?Q(m) ?m^{2} s^{?1} 
D_{b,cortexcortex}  Bcd diffusion constant (from cortex to cortex)  12.5 ?m^{2} s^{?1} 
D_{b,corecore}  Bcd diffusion constant (from core to core)  18.0 ?m^{2} s^{?1} 
k_{b, nuclear}  Bcd degradation rate constant in nuclei  7.70?10^{?4} s^{?1} 
k_{b, cortical}, k_{b, core}  Bcd degradation rate constant in cytoplasm  9.63?10^{?5} s^{?1} 
<h_{b}>  The average Bcd binding rate constant to the hb enhancer  7.70?10^{?3} s^{?1} 
f_{b}  Bcd dissociation rate constant from the hb enhancer  7.70?10^{?2} s^{?1} 
D_{h,import}  Hb diffusion coefficient (import into nucleus)  0.1 ?m^{2} s^{?1} 
D_{h,export}  Hb diffusion coefficient (export from nucleus)  0.1 ?m^{2} s^{?1} 
D_{h,corecortex}  Hb diffusion coefficient (from core to cortex)  0.1 ?m^{2} s^{?1} 
D_{h,cortexcore}  Hb diffusion coefficient (from cortex to core)  0.1 ?m^{2} s^{?1} 
D_{h,cortexcortex}  Hb diffusion coefficient (from cortex to cortex)  0.1 ?m^{2} s^{?1} 
D_{h,corecore}  Hb diffusion coefficient (from core to core)  0.1 ?m^{2} s^{?1} 
k_{h, nuclear}  Hb degradation rate constant in nuclei  3.85?10^{?4} s^{?1} 
k_{h, cortical}, k_{h, core}  Hb degradation rate constant in cytoplasm  4.81?10^{?5} s^{?1} 
<h_{h}>  The average Hb binding rate constant to the hb enhancer  1.85?10^{?2} s^{?1} 
f_{h}  Hb dissociation rate constant from the hb enhancer  3.85?10^{?2} s^{?1} 
g(0)  Hb synthesis rate at S?=?0  7.50?10^{?6} s^{?1} 
g(1)  Hb synthesis rate at S?=?1  7.50?10^{?6} s^{?1} 
g(2)  Hb synthesis rate at S?=?2  1.50?10^{?1} s^{?1} 
g(3)  Hb synthesis rate at S?=?3  7.50?10^{?1} s^{?1} 
N_{burst}  Burst size of Hb synthesis  10 molecules 
v_{ex}  Export rate constant of Hb or Bcd from interaction volume  2.04?10^{?2} s^{?1} 
<v_{im}>  The average import rate constant of Hb or Bcd to interaction volume  5.98?10^{?5} s^{?1} 
<V_{r}>  The average ratio of interaction volume to nuclear volume  10^{?2.5} 
<?>  The average side length of hexagon of each site  4.04 ?m 
<?x>  The average side length of the site along the AP axis  7 ?m 
The factor J(m) at the nuclear cycle m is J(m)?=?0.1?exp(ln(100)/14?10) for m?=?1?10 and 0.1?exp(m ln(100)/14) for m?=?11?14. The factor Q(m) at the nuclear cycle m is Q(m)?=?1 (for m?=?1?8), 0.667 (for m?=?9), 0.444 (for m?=?10), 0.198 (for m?=?11), 0.132 (for m?=?12), and 0.088 for (m?=?13?14).
Article Categories:

Previous Document: Computation of conformational coupling in allosteric proteins.
Next Document: FLORA: a novel method to predict protein function from structure in diverse superfamilies.