Cancer is caused by genetic abnormalities, such as mutations of oncogenes or tumor suppressor genes, which alter downstream signal transduction pathways and protein-protein interactions. Comparisons of the interactions of proteins in cancerous and normal cells can shed light on the mechanisms of carcinogenesis.

Lung cancer is the leading cause of cancer deaths worldwide. Many studies have investigated the carcinogenic process and identified the biomarkers for signature classification. However, based on the research dedicated to this field, there is no highly sensitive network-based method for carcinogenesis characterization and diagnosis from the systems perspective.

With the accumulation of increasing omics data, a key goal of systems biology is to construct networks at different cellular levels to investigate cellular machinery of the cell. However, there is currently no satisfactory method to construct an integrated cellular network that combines the gene regulatory network and the signaling regulatory pathway. In this study, we integrated different kinds of omics data and developed a systematic method to construct the integrated cellular network based on coupling dynamic models and statistical assessments. The proposed method was applied to S. cerevisiae stress responses, elucidating the stress response mechanism of the yeast. From the resulting integrated cellular network under hyperosmotic stress, the highly connected hubs which are functionally relevant to the stress response were identified. Beyond hyperosmotic stress, the integrated network under heat shock and oxidative stress were also constructed and the crosstalks of these networks were analyzed, specifying the significance of some transcription factors to serve as the decision-making devices at the center of the bow-tie structure and the crucial role for rapid adaptation scheme to respond to stress. In addition, the predictive power of the proposed method was also demonstrated. We successfully construct the integrated cellular network which is validated by literature evidences. The integration of transcription regulations and protein-protein interactions gives more insight into the actual biological network and is more predictive than those without integration. The method is shown to be powerful and flexible and can be used under different conditions and for different species. The coupling dynamic models of the whole integrated cellular network are very useful for theoretical analyses and for further experiments in the fields of network biology and synthetic biology.

Motivation: Transcription factor binding sites are known to co-occur in the same gene owing to cooperativity of the transcription factors (TFs) that bind to them. Genome-wide location data can help us understand how an individual TF regulates its target gene. Nevertheless, how TFs cooperate to regulate their target genes still needs further study. In this study, genome-wide location data and expression profiles are integrated to reveal how TFs cooperate to regulate their target genes from the stochastic system perspective.

Motivation: The robustness of a biochemical network is defined as the tolerance of variations in kinetic parameters with respect to the maintenance of steady state. Robustness also plays an important role in the fail-safe mechanism in the evolutionary process of biochemical networks. The purposes of this paper are to use the synergism and saturation system (S-system) representation to describe a biochemical network and to develop a robustness measure of a biochemical network subject to variations in kinetic parameters. Since most biochemical networks in nature operate close to the steady state, we consider only the robustness measurement of a biochemical network at the steady state. Results: We show that the upper bound of the tolerated parameter variations is related to the system matrix of a biochemical network at the steady state. Using this upper bound, we can calculate the tolerance (robustness) of a biochemical network without testing many parametric perturbations. We find that a biochemical network with a large tolerance can also better attenuate the effects of variations in rate parameters and environments. Compensatory parameter variations and network redundancy are found to be important mechanisms for the robustness of biochemical networks. Finally, four biochemical networks, such as a cascaded biochemical network, the glycolytic–glycogenolytic pathway in a perfused rat liver, the tricarboxylic acid cycle in Dictyostelium discoideum and the cAMP oscillation network in bacterial chemotaxis, are used to illustrate the usefulness of the proposed robustness measure. Supplementary information:http://www.ee.nthu.edu.tw/~bschen/robustness_bio-networks/ Contact:bschen@ee.nthu.edu.tw

Background: Noise has many important roles in cellular genetic regulatory functions at the nanomolar scale. At present, no good theory exists for identifying all possible mechanisms of genetic regulatory networks to attenuate the molecular noise to achieve regulatory ability or to amplify the molecular noise to randomize outcomes to the advantage of diversity. Therefore, the noise filtering of genetic regulatory network is an important topic for gene networks under intrinsic fluctuation and extrinsic noise. Results: Based on stochastic dynamic regulation equation, the intrinsic fluctuation in reaction rates is modeled as a state-dependent stochastic process, which will influence the stability of gene regulatory network, especially, with low concentrations of reacting species. Then the mechanisms of genetic regulatory network to attenuate or amplify extrinsic fluctuation are revealed from the nonlinear stochastic filtering point of view. Furthermore, a simple measure of attenuation level or amplification level of extrinsic noise for genetic regulatory networks is also introduced by nonlinear robust filtering method. Based on the global linearization scheme, a convenient method is introduced to measure noise attenuation or amplification for each gene of the nonlinear stochastic regulatory network by solving a set of filtering problems, which correspond to a set of linearized stochastic regulatory networks. Finally, by the proposed methods, several simulation examples of genetic regulatory networks are given to measure their robust stability under intrinsic fluctuations, and to estimate the genes' attenuation and amplification levels under extrinsic noises. Conclusion: In this study, a stochastic nonlinear dynamic model is developed for genetic regulatory networks under intrinsic fluctuation and extrinsic noise. By the method we proposed, we could determine the robust stability under intrinsic fluctuations and identify the genes that are significantly affected by extrinsic noises, which we call the weak structure of the network. This method will be potential for robust gene circuit design in future, on which a drug design could be based.