FBA(flux balance analysis) with Boolean rules for representing regulatory events has correctly predicted cellular behaviors, such as optimal flux distribution, maximal growth rate, metabolic by-product, and substrate concentration changes, with various environmental conditions. However, until now, since FBA has not taken into account a hierarchical regulatory network, it has limited the representation of the whole transcriptional regulation mechanism and interactions between specific regulatory proteins and genes. In this paper, in order to solve these problems, we describe the construction of hierarchical regulatory network with defined symbols and the introduction of a weight for representing interactions between symbols. Finally, the whole cellular behaviors with time were simulated through the linkage of a hierarchical regulatory network module and dynamic simulation module including FBA. The central metabolic network of E. coli was chosen as the basic model to identify our suggested modeling method.