The ‘finite number effect’ in chemical reactions, known as chemical noise, is the main cause of creating heterogeneity in a population of genetically identical cells. With the help of tools from physical chemistry, statistical physics and nonlinear dynamics, we investigate theoretically and computationally chemical noise propagation in networks of chemical reactions consisting of various types of feedback and feed-forward regulatory motifs. We are developing theoretical/computational methods to analyze the steady state properties of the dynamical systems that regulate cell fate decision making processes in living cells. Our goal is to uncover the topological requirements of chemical reaction networks in producing robust response in the face of chemical noise inside a living cell. A particular interest is to understand the underlying potential energy landscapes of the cellular differentiation pathways in multicellular organisms.