All cellular phenomena are implemented by chemical reactions of diverse molecules with small numbers being encapsulated within a microscopic cellular environment. How should we describe and model stochastic behaviors of the molecules? How can the small number of molecules induce qualitative transition? What kind of thermodynamic law does an equilibrium / non-equilibrium intracellular reaction system follow? We address these fundamental questions and develop methods for their verification using quantitative data.