Cellular processes critical to sustaining physiology, including growth, division and differentiation, are carefully governed by intricate control systems. Deregulations in these systems often result in complex diseases such as cancer. Hence, it is crucial to understand the interactions between molecular players of these control systems, their emergent network dynamics, and, ultimately, the overall contribution to cellular physiology. In this dissertation, we have developed a mathematical framework to understand two such cellular systems: an early checkpoint (START) in the budding yeast cell cycle (Chapter 1), and the canonical Wnt signaling pathway involved in cell proliferation and differentiation (Chapter 2). START transition is an important decision point where the cell commits to one round DNA replication followed by cell division. Several years of experimental research have gone into uncovering molecular details of this process, but a unified understanding is yet to emerge. In chapter one, we have developed a comprehensive mathematical model of START transition that incorporates several findings including information about the phosphorylation state of key START proteins and their subcellular localization. In the second chapter, we focus on modeling the canonical Wnt signaling pathway, a cellular circuit that plays a key role in cell proliferation and differentiation. The Wnt pathway is often deregulated in colon cancers. Based on some evidence of bistability in the Wnt signaling pathway, we proposed the existence of a positive feedback loop underlying the activation and inactivation of the core protein complex of the pathway. Bistability is a common feature of biological systems that toggle between ON and OFF states because it ensures robust switching back and forth between the two states. To study and explain the behavior of this dynamical system, we developed a mathematical model. Based on experimentally determined interactions, our simple model recapitulates the observed phenomena of bimodality (bistability) and hysteresis under the effects of the physiological signal (Wnt), a Wnt-mimic (LiCl), and a stabilizer of one of the key members of core complex (IWR-1). Overall, we believe that cell biologists and molecular geneticists can benefit from our work by using our model to make novel quantitative predictions for experimental verification. / Ph. D.

A complex network of genes and proteins governs the robust progression through cell cycles in the presence of inevitable noise. Stochastic modeling is viewed as a key paradigm to study the effects of intrinsic and extrinsic noise on the dynamics of biochemical networks. A detailed quantitative description of such complex and multiscale networks via stochastic modeling poses several challenges. First, stochastic models generally require extensive computations, particularly when applied to large networks. Second, the accuracy of stochastic models is highly dependent on the quality of the parameter estimation based on experimental observations. The goal of this dissertation is to address these problems by developing new efficient methods for modeling and simulation of stochastic effects in biochemical systems. Particularly, a hybrid stochastic model is developed to represent a detailed molecular mechanism of cell cycle control in budding yeast cells. In a single multiscale model, the proposed hybrid approach combines the advantages of two regimes: 1) the computational efficiency of a deterministic approach, and 2) the accuracy of stochastic simulations. The results show that this hybrid stochastic model achieves high computational efficiency while generating simulation results that match very well with published experimental measurements. Furthermore, a new hierarchical deep classification (HDC) algorithm is developed to address the parameter estimation problem in a monomolecular system. The HDC algorithm adopts a neural network that, via multiple hierarchical search steps, finds reasonably accurate ranges for the model parameters. To train the neural network in the presence of experimental data scarcity, the proposed method leverages the domain knowledge from stochastic simulations to generate labeled training data. The results show that the proposed HDC algorithm yields accurate ranges for the model parameters and highlight the potentials of model-free learning for parameter estimation in stochastic modeling of complex biochemical networks. / Doctor of Philosophy / Cell cycle is a process in which a growing cell replicates its DNA and divides into two cells. Progression through the cell cycle is regulated by complex interactions between networks of genes, transcripts, and proteins. These interactions inside the confined volume of a cell are subject to inherent noise. To provide a quantitative description of the cell cycle, several deterministic and stochastic models have been developed. However, deterministic models cannot capture the intrinsic noise. In addition, stochastic modeling poses the following challenges. First, stochastic models generally require extensive computations, particularly when applied to large networks. Second, the accuracy of stochastic models is highly dependent on the accuracy of the estimated model parameters. The goal of this dissertation is to address these challenges by developing new efficient methods for modeling and simulation of stochastic effects in biochemical networks. The results show that the proposed hybrid model that combines stochastic and deterministic modeling approaches can achieve high computational efficiency while generating accurate simulation results. Moreover, a new machine learning-based method is developed to address the parameter estimation problem in biochemical systems. The results show that the proposed method yields accurate ranges for the model parameters and highlight the potentials of model-free learning for parameter estimation in stochastic modeling of complex biochemical networks.