Optimal Design of Single-Cell Experiments within Temporally Fluctuating Environments

Modern biological experiments are becoming increasingly complex, and designing these experiments to yield the greatest possible quantitative insight is an open challenge. Increasingly, computational models of complex stochastic biological systems are being used to understand and predict biological behaviors or to infer biological parameters. Such quantitative analyses can also help to improve experiment designs for particular goals, such as to learn more about specific model mechanisms or to reduce prediction errors in certain situations. A classic approach to experiment design is to use the Fisher information matrix (FIM), which quantifies the expected information a particular experiment will reveal about model parameters. The Finite State Projection based FIM (FSP-FIM) was recently developed to compute the FIM for discrete stochastic gene regulatory systems, whose complex response distributions do not satisfy standard assumptions of Gaussian variations. In this work, we develop the FSP-FIM analysis for a stochastic model of stress response genes in S. cerevisae under time-varying MAPK induction. We verify this FSP-FIM analysis and use it to optimize the number of cells that should be quantified at particular times to learn as much as possible about the model parameters. We then extend the FSP-FIM approach to explore how different measurement times or genetic modifications help to minimize uncertainty in the sensing of extracellular environments, and we experimentally validate the FSP-FIM to rank single-cell experiments for their abilities to minimize estimation uncertainty of NaCl concentrations during yeast osmotic shock. This work demonstrates the potential of quantitative models to not only make sense of modern biological data sets, but to close the loop between quantitative modeling and experimental data collection.

[1]  S. Kay Fundamentals of statistical signal processing: estimation theory , 1993 .

[2]  Michael P H Stumpf,et al.  An information-theoretic framework for deciphering pleiotropic and noisy biochemical signaling , 2018, Nature Communications.

[3]  Jakob Ruess,et al.  Shaping bacterial population behavior through computer-interfaced control of individual cells , 2017, Nature Communications.

[4]  M. Peter,et al.  Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings , 2013, Nature Methods.

[5]  J. Timmer,et al.  Systems biology: experimental design , 2009, The FEBS journal.

[6]  W. Marsden I and J , 2012 .

[7]  Nicolas E. Buchler,et al.  BayFish: Bayesian inference of transcription dynamics from population snapshots of single-molecule RNA FISH in single cells , 2017, Genome Biology.

[8]  Jared E. Toettcher,et al.  Stochastic Gene Expression in a Lentiviral Positive-Feedback Loop: HIV-1 Tat Fluctuations Drive Phenotypic Diversity , 2005, Cell.

[9]  Jeffrey J. Tabor,et al.  Optogenetic control of Bacillus subtilis gene expression , 2019, Nature Communications.

[10]  João Pedro Hespanha,et al.  Approximate Moment Dynamics for Chemically Reacting Systems , 2011, IEEE Transactions on Automatic Control.

[11]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[12]  Scott A. Rifkin,et al.  Imaging individual mRNA molecules using multiple singly labeled probes , 2008, Nature Methods.

[13]  Christophe Zimmer,et al.  smiFISH and FISH-quant – a flexible single RNA detection approach with super-resolution capability , 2016, Nucleic acids research.

[14]  Guoliang Li,et al.  Generating kinetic environments to study dynamic cellular processes in single cells , 2019, Scientific Reports.

[15]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[16]  Brian Munsky,et al.  Integrating single-molecule experiments and discrete stochastic models to understand heterogeneous gene transcription dynamics. , 2015, Methods.

[17]  Brian Munsky,et al.  The finite state projection based Fisher information matrix approach to estimate information and optimize single-cell experiments , 2019, PLoS Comput. Biol..

[18]  Francesc Posas,et al.  Response to Hyperosmotic Stress , 2012, Genetics.

[19]  Heinz Koeppl,et al.  Parallel feedback loops control the basal activity of the HOG MAPK signaling cascade. , 2015, Integrative biology : quantitative biosciences from nano to macro.

[20]  John Lygeros,et al.  Designing experiments to understand the variability in biochemical reaction networks , 2013, Journal of The Royal Society Interface.

[21]  Robert Gentleman,et al.  Software for Computing and Annotating Genomic Ranges , 2013, PLoS Comput. Biol..

[22]  John Lygeros,et al.  Iterative experiment design guides the characterization of a light-inducible gene expression circuit , 2015, Proceedings of the National Academy of Sciences.

[23]  L. Tsimring,et al.  Accurate information transmission through dynamic biochemical signaling networks , 2014, Science.

[24]  J. Timmer,et al.  Experimental Design for Parameter Estimation of Gene Regulatory Networks , 2012, PloS one.

[25]  Eric J. Deeds,et al.  Fundamental trade-offs between information flow in single cells and cellular populations , 2017, Proceedings of the National Academy of Sciences.

[26]  M. Khammash,et al.  Systematic Identification of Signal-Activated Stochastic Gene Regulation , 2013, Science.

[27]  Fabian Rudolf,et al.  Transient Activation of the HOG MAPK Pathway Regulates Bimodal Gene Expression , 2011, Science.

[28]  Aviv Regev,et al.  Deconstructing transcriptional heterogeneity in pluripotent stem cells , 2014, Nature.

[29]  Gregory Batt,et al.  Balancing a genetic toggle switch by real-time feedback control and periodic forcing , 2017, Nature Communications.

[30]  J. Lygeros,et al.  Moment-based inference predicts bimodality in transient gene expression , 2012, Proceedings of the National Academy of Sciences.

[31]  Brian Munsky,et al.  Distribution shapes govern the discovery of predictive models for gene regulation , 2017, Proceedings of the National Academy of Sciences.

[32]  Michael P. H. Stumpf,et al.  Maximizing the Information Content of Experiments in Systems Biology , 2013, PLoS Comput. Biol..

[33]  M. Khammash,et al.  The finite state projection algorithm for the solution of the chemical master equation. , 2006, The Journal of chemical physics.

[34]  A. van Oudenaarden,et al.  Using Gene Expression Noise to Understand Gene Regulation , 2012, Science.

[35]  Brian Munsky,et al.  Finite state projection based bounds to compare chemical master equation models using single-cell data. , 2016, The Journal of chemical physics.

[36]  Brian Munsky,et al.  The finite state projection based Fisher information matrix approach to estimate information and optimize single-cell experiments , 2018, bioRxiv.

[37]  E. Gilles,et al.  Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors , 2002, Nature Biotechnology.

[38]  M. Khammash,et al.  An Optogenetic Platform for Real-Time, Single-Cell Interrogation of Stochastic Transcriptional Regulation , 2018, Molecular cell.

[39]  David Baker,et al.  Modular and tunable biological feedback control using a de novo protein switch , 2019, Nature.

[40]  Aleksandra M Walczak,et al.  Information transmission in genetic regulatory networks: a review , 2011, Journal of physics. Condensed matter : an Institute of Physics journal.

[41]  E. Klipp,et al.  Integrative model of the response of yeast to osmotic shock , 2005, Nature Biotechnology.

[42]  Ilya Nemenman,et al.  Universal Properties of Concentration Sensing in Large Ligand-Receptor Networks. , 2019, Physical review letters.

[43]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[44]  Alexander Thiemicke,et al.  Generating kinetic environments to study dynamic cellular processes in single cells , 2019 .

[45]  Michael P H Stumpf,et al.  Sensitivity, robustness, and identifiability in stochastic chemical kinetics models , 2011, Proceedings of the National Academy of Sciences.

[46]  I. Nemenman,et al.  Information Transduction Capacity of Noisy Biochemical Signaling Networks , 2011, Science.

[47]  Bruce Tidor,et al.  Sloppy models, parameter uncertainty, and the role of experimental design. , 2010, Molecular bioSystems.

[48]  Christoph Zimmer,et al.  Experimental Design for Stochastic Models of Nonlinear Signaling Pathways Using an Interval-Wise Linear Noise Approximation and State Estimation , 2016, PloS one.

[49]  F S Fay,et al.  Visualization of single RNA transcripts in situ. , 1998, Science.

[50]  Akademii︠a︡ medit︠s︡inskikh nauk Sssr Journal of physics , 1939 .

[51]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.

[52]  A. Oudenaarden,et al.  A Systems-Level Analysis of Perfect Adaptation in Yeast Osmoregulation , 2009, Cell.