Computational Systems Biology has naturally emerged as a branch of computational models applied to biological systems due to large-scale biological data production in the past decades. Although several models have been developed in order to simulate cellular subsystems, only a few tackle the problem of integrating several cellular processes. In previous works, an extended biochemical network modeling framework was proposed aiming at integration of cellular subsystems into a single biochemical network [1]. Additionally, constraint-based (CB) methods have been successfully employed to study the dynamics of metabolism based on its respective network [2]. In this work, we prose to extend the current CB methods to a wider cellular context using the extended biochemical networks as the constraint source. For example, Fig. 1 depicts the E. coli’s glucose metabolic network from which the dynamics was simulated using the proposed CB method. The obtained results and validation against experimental data will be presented. This work is a first step towards the simulation of integrated biochemical systems based mainly on its respective networks.