Systems biology rests on the idea that biological complexity can be better unraveled through the interplay of modeling and experimentation. However, the success of this approach depends critically on the informativeness of the chosen experiments, which is usually unknown a priori. Here, we propose a systematic scheme based on iterations of optimal experiment design, flow cytometry experiments, and Bayesian parameter inference to guide the discovery process in the case of stochastic biochemical reaction networks. To illustrate the benefit of our methodology, we apply it to the characterization of an engineered light-inducible gene expression circuit in yeast and compare the performance of the resulting model with models identified from nonoptimal experiments. In particular, we compare the parameter posterior distributions and the precision to which the outcome of future experiments can be predicted. Moreover, we illustrate how the identified stochastic model can be used to determine light induction patterns that make either the average amount of protein or the variability in a population of cells follow a desired profile. Our results show that optimal experiment design allows one to derive models that are accurate enough to precisely predict and regulate the protein expression in heterogeneous cell populations over extended periods of time. T he use of quantitative mathematical models to investigate biochemical reaction networks is nowadays common practice. Typically, models are built based on the available biological knowledge and used to generate hypotheses, which are then refined or invalidated through experimentation. For this process to be successful, it is of paramount importance to design and perform experiments that yield the information required to identify the model under consideration. Optimal experiment design techniques have been extensively studied for ordinary differential equation models (1-5), which are typically used to describe the average behavior of cell populations (6-8). With the development of high-throughput measurement techniques, such as flow cytometry, it has, however, become evident that restricting the attention only to the average population behavior neglects the potentially valuable information contained in the full population distribution (9-11). This additional information can be captured by stochastic models. Recently, methods for parameter inference (12-15) and optimal experiment design (16, 17) for stochastic models have been developed and applied to a number of biological systems (12, 18). However, a systematic characterization procedure that exploits the information gained from each performed experiment has not yet been fully developed or experimentally validated.Here, we provide the first study, to our knowledge, in which a noisy biochemical reaction network is characterized and ultimately also controlled through iterations of optimally designed flow cytometry experiments and stochastic modeling. Specifically, we consider a gene expression circuit in yeast that has been eng...