Nonlinear biochemical signal processing via noise propagation

Systems biology–the transformative approach to integrate sciences across disciplines

Life science is the study of living organisms, including bacteria, plants, and animals. Given the importance of biology, chemistry, and bioinformatics, we anticipate that this chapter may contribute to a better understanding of the interdisciplinary connections in life science. Research in applied biological sciences has changed the paradigm of basic and applied research. Biology is the study of life and living organisms, whereas science is a dynamic subject that as a result of constant research, new fields are constantly emerging. Some fields come and go, whereas others develop into new, well-recognized entities. Chemistry is the study of composition of matter and its properties, how the substances merge or separate and also how substances interact with energy. Advances in biology and chemistry provide another means to understand the biological system using many interdisciplinary approaches. Bioinformatics is a multidisciplinary or rather transdisciplinary field that encourages the use of computer tools and methodologies for qualitative and quantitative analysis. There are many instances where two fields, biology and chemistry have intersection. In this chapter, we explain how current knowledge in biology, chemistry, and bioinformatics, as well as its various interdisciplinary domains are merged into life sciences and its applications in biological research.

The E. coli Whole-Cell Modeling Project

The Escherichia coli whole-cell modeling project seeks to create the most detailed computational model of an E. coli cell in order to better understand and predict the behavior of this model organism. Details about the approach, framework, and current version of the model are discussed. Currently, the model includes the functions of 43% of characterized genes, with ongoing efforts to include additional data and mechanisms. As additional information is incorporated in the model, its utility and predictive power will continue to increase, which means that discovery efforts can be accelerated by community involvement in the generation and inclusion of data. This project will be an invaluable resource to the E. coli community that could be used to verify expected physiological behavior, to predict new outcomes and testable hypotheses for more efficient experimental design iterations, and to evaluate heterogeneous data sets in the context of each other through deep curation.

The switch of DNA states filtering the extrinsic noise in the system of frequency modulation

There is a special node, which the large noise of the upstream element may not always lead to a broad distribution of downstream elements. This node is DNA, with upstream element TF and downstream elements mRNA and proteins. By applying the stochastic simulation algorithm (SSA) on gene circuits inspired by the fim operon in Escherichia coli, we found that cells exchanged the distribution of the upstream transcription factor (TF) for the transitional frequency of DNA. Then cells do an inverse transform, which exchanges the transitional frequency of DNA for the distribution of downstream products. Due to this special feature, DNA in the system of frequency modulation is able to reset the noise. By probability generating function, we know the ranges of parameter values that grant such an interesting phenomenon.

Noisy signal propagation and amplification in phenotypic transition cascade of colonic cells

Like genes and proteins, cells can use biochemical networks to sense and process information. The differentiation of the cell state in colonic crypts forms a typical unidirectional phenotypic transitional cascade, in which stem cells differentiate into the transit-amplifying cells (TACs), and TACs continue to differentiate into fully differentiated cells. In order to quantitatively describe the relationship between the noise of each compartment and the amplification of signals, the gain factor is introduced, and the gain-fluctuation relation is obtained by using the linear noise approximation of the master equation. Through the simulation of these theoretical formulas, the characters of noise propagation and amplification are studied. It is found that the transmitted noise is an important part of the total noise in each downstream cell. Therefore, a small number of downstream cells can only cause its small inherent noise, but the total noise may be very large due to the transmitted noise. The influence of the transmitted noise may be the indirect cause of colon cancer. In addition, the total noise of the downstream cells always has a minimum value. As long as a reasonable value of the gain factor is selected, the number of cells in colonic crypts will be controlled within the normal range. This may be a good method to intervene the uncontrollable growth of tumor cells and effectively control the deterioration of colon cancer.

Diffusion-dynamics laws in stochastic reaction networks

Many biological activities are induced by cellular chemical reactions of diffusing reactants. The dynamics of such systems can be captured by stochastic reaction networks. A recent numerical study has shown that diffusion can significantly enhance the fluctuations in gene regulatory networks. However, the universal relation between diffusion and stochastic system dynamics remains veiled. Within the approximation of reaction-diffusion master equation (RDME), we find general relation that the steady-state distribution in complex balanced networks is diffusion-independent. Here, complex balance is the nonequilibrium generalization of detailed balance. We also find that for a diffusion-included network with a Poisson-like steady-state distribution, the diffusion can be ignored at steady state. We then derive a necessary and sufficient condition for networks holding such steady-state distributions. Moreover, we show that for linear reaction networks the RDME reduces to the chemical master equation, which implies that the stochastic dynamics of networks is unaffected by diffusion at any arbitrary time. Our findings shed light on the fundamental question of when diffusion can be neglected, or (if nonnegligible) its effects on the stochastic dynamics of the reaction network.

Synthetic Biology: Engineering Living Systems from Biophysical Principles

Synthetic biology was founded as a biophysical discipline that sought explanations for the origins of life from chemical and physical first principles. Modern synthetic biology has been reinvented as an engineering discipline to design new organisms as well as to better understand fundamental biological mechanisms. However, success is still largely limited to the laboratory and transformative applications of synthetic biology are still in their infancy. Here, we review six principles of living systems and how they compare and contrast with engineered systems. We cite specific examples from the synthetic biology literature that illustrate these principles and speculate on their implications for further study. To fully realize the promise of synthetic biology, we must be aware of life's unique properties.

Bistability and Bifurcation in Minimal Self-Replication and Nonenzymatic Catalytic Networks

Bistability and bifurcation, found in a wide range of biochemical networks, are central to the proper function of living systems. We investigate herein recent model systems that show bistable behavior based on nonenzymatic self-replication reactions. Such models were used before to investigate catalytic growth, chemical logic operations, and additional processes of self-organization leading to complexification. By solving for their steady-state solutions by using various analytical and numerical methods, we analyze how and when these systems yield bistability and bifurcation and discover specific cases and conditions producing bistability. We demonstrate that the onset of bistability requires at least second-order catalysis and results from a mismatch between the various forward and reverse processes. Our findings may have far-reaching implications in understanding early evolutionary processes of complexification, emergence, and potentially the origin of life.

Guidelines for Reproducibly Building and Simulating Systems Biology Models

OBJECTIVE

Reproducibility is the cornerstone of the scientific method. However, currently, many systems biology models cannot easily be reproduced. This paper presents methods that address this problem.

METHODS

We analyzed the recent Mycoplasma genitalium whole-cell (WC) model to determine the requirements for reproducible modeling.

RESULTS

We determined that reproducible modeling requires both repeatable model building and repeatable simulation.

CONCLUSION

New standards and simulation software tools are needed to enhance and verify the reproducibility of modeling. New standards are needed to explicitly document every data source and assumption, and new deterministic parallel simulation tools are needed to quickly simulate large, complex models.

SIGNIFICANCE

We anticipate that these new standards and software will enable researchers to reproducibly build and simulate more complex models, including WC models.

Controlling E. coli Gene Expression Noise

Intracellular protein copy numbers show significant cell-to-cell variability within an isogenic population due to the random nature of biological reactions. Here we show how the variability in copy number can be controlled by perturbing gene expression. Depending on the genetic network and host, different perturbations can be applied to control variability. To understand more fully how noise propagates and behaves in biochemical networks we developed stochastic control analysis (SCA) which is a sensitivity-based analysis framework for the study of noise control. Here we apply SCA to synthetic gene expression systems encoded on plasmids that are transformed into Escherichia coli. We show that (1) dual control of transcription and translation efficiencies provides the most efficient way of noise-versus-mean control. (2) The expressed proteins follow the gamma distribution function as found in chromosomal proteins. (3) One of the major sources of noise, leading to the cell-to-cell variability in protein copy numbers, is related to bursty translation. (4) By taking into account stochastic fluctuations in autofluorescence, the correct scaling relationship between the noise and mean levels of the protein copy numbers was recovered for the case of weak fluorescence signals.

Stochastic modular analysis for gene circuits: interplay among retroactivity, nonlinearity, and stochasticity

This chapter introduces a computational analysis method for analyzing gene circuit dynamics in terms of modules while taking into account stochasticity, system nonlinearity, and retroactivity. (1) ANALOG ELECTRICAL CIRCUIT REPRESENTATION FOR GENE CIRCUITS: A connection between two gene circuit components is often mediated by a transcription factor (TF) and the connection signal is described by the TF concentration. The TF is sequestered to its specific binding site (promoter region) and regulates downstream transcription. This sequestration has been known to affect the dynamics of the TF by increasing its response time. The downstream effect-retroactivity-has been shown to be explicitly described in an electrical circuit representation, as an input capacitance increase. We provide a brief review on this topic. (2) MODULAR DESCRIPTION OF NOISE PROPAGATION: Gene circuit signals are noisy due to the random nature of biological reactions. The noisy fluctuations in TF concentrations affect downstream regulation. Thus, noise can propagate throughout the connected system components. This can cause different circuit components to behave in a statistically dependent manner, hampering a modular analysis. Here, we show that the modular analysis is still possible at the linear noise approximation level. (3) NOISE EFFECT ON MODULE INPUT-OUTPUT RESPONSE: We investigate how to deal with a module input-output response and its noise dependency. Noise-induced phenotypes are described as an interplay between system nonlinearity and signal noise. Lastly, we provide the comprehensive approach incorporating the above three analysis methods, which we call "stochastic modular analysis." This method can provide an analysis framework for gene circuit dynamics when the nontrivial effects of retroactivity, stochasticity, and nonlinearity need to be taken into account.

Stochasticity of intranuclear biochemical reaction processes controls the final decision of cell fate associated with DNA damage

A massive integrative mathematical model of DNA double-strand break (DSB) generation, DSB repair system, p53 signaling network, and apoptosis induction pathway was constructed to explore the dominant factors of unknown criteria of cell fate decision. In the proposed model, intranuclear reactions were modeled as stochastic processes and cytoplasmic reactions as deterministic processes, and both reaction sets were simulated simultaneously. The simulated results at the single-cell level showed that the model generated several sustained oscillations (pulses) of p53, Mdm2, ATM, and Wip1, and cell-to-cell variability in the number of p53 pulses depended on IR intensity. In cell populations, the model generated damped p53 oscillations, and IR intensity affected the amplitude of the first p53 oscillation. Cells were then subjected to the same IR dose exhibiting apoptosis induction variability. These simulated results are in quantitative agreement with major biological findings observed in human breast cancer epithelial MCF7, NIH3T3, and fibrosarcoma cells, demonstrating that the proposed model was concededly biologically appropriate. Statistical analysis of the simulated results shows that the generation of multiple p53 pulses is a prerequisite for apoptosis induction. Furthermore, cells exhibited considerable individual variability in p53 dynamics, which correlated with intrinsic apoptosis induction. The simulated results based on the proposed model demonstrated that the stochasticity of intranuclear biochemical reaction processes controls the final decision of cell fate associated with DNA damage. Applying stochastic simulation to an exploration of intranuclear biochemical reaction processes is indispensable in enhancing the understanding of the dynamic characteristics of biological multi-layered systems of higher organisms.