Abstract:Heterogeneity in gene expression across isogenic cell populations can give rise to phenotypic diversity, even when cells are in homogenous environments. This diversity arises from the discrete, stochastic nature of biochemical reactions, which naturally arise due to the very small numbers of genes, RNA, or protein molecules in single cells. Modern measurements of single biomolecules have created a vast wealth of information about the fluctuations of these molecules, but a quantitative understanding of these complex, stochastic systems requires precise computational tools. In this article, we present modern tools necessary to model variability in biological system and to compare model results to experimental data. We review the Chemical Master Equation and approaches to solve for probability distributions of discrete numbers of biomolecules. We discuss how to fit probability distributions to single-cell data using likelihood based approaches. Finally, we provide examples of fitting discrete stochastic models to single-molecule fluorescent in-situ hybridization data that quantifies RNA levels across populations of cells in bacteria and yeast.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the phenotypic diversity in a homologous environment, even in cell populations with the same genes, which is caused by the discreteness and randomness of biochemical reactions. Specifically, the paper focuses on how to use modern computational tools to model the variability in such biological systems and compare the model results with experimental data. The paper discusses the Chemical Master Equation (CME) and its solution methods, as well as how to use the likelihood method to fit the probability distribution of single - cell data. In addition, the paper also provides specific examples of fitting discrete stochastic models to bacterial and yeast cell populations to quantify RNA levels through single - molecule fluorescence in situ hybridization (smFISH) technology. ### Core problems of the paper 1. **Phenotypic diversity**: Explain why phenotypic diversity occurs in cell populations with the same genes even in a homologous environment. 2. **Sources of randomness**: Explore how this diversity originates from the discreteness and randomness of biochemical reactions in which a very small number of genes, RNAs or protein molecules in the cell are involved. 3. **Quantitative understanding**: Propose the use of accurate computational tools to quantitatively understand these complex stochastic systems. 4. **Model verification**: Introduce how to compare model results with experimental data, especially how to use the likelihood method to fit the probability distribution to single - cell data. ### Specific problems - **Chemical Master Equation (CME)**: It details the derivation of CME and its application in describing the stochastic processes of biochemical reactions. - **Finite State Projection (FSP) method**: It discusses how the FSP method can be used to solve CME and provides an algorithm for gradually expanding the projection space to improve computational accuracy. - **Single - cell data analysis**: It shows how to use single - molecule fluorescence in situ hybridization (smFISH) data to verify and calibrate the CME model. ### Solutions - **Theoretical basis**: Provide a theoretical framework through CME for describing the randomness of biochemical reactions. - **Computational tools**: Propose computational tools such as the FSP method for solving CME and analyzing single - cell data. - **Experimental verification**: Verify the effectiveness of the model through actual single - cell data (such as smFISH data) to ensure that the model can accurately reflect the stochastic behavior of biological systems. In summary, this paper aims to deeply understand the randomness in biochemical reactions and its impact on cell phenotypic diversity through a combination of theoretical and experimental methods, and provide effective computational tools to support this research.

Stochasticity or Noise in Biochemical Reactions