Biological systems consist of many components and interactions between them. In Systems Biology the principal problem is modeling complex biological systems and reconstructing interactions between their building blocks. Symbolic machine learning approaches have the power to model structured domains and relations among objects. However biological domains require uncertainty handling due to their hidden complex nature. Statistical machine learning approaches have the potential to model uncertainty in a robust manner. In this paper we apply a hybrid symbolic-statistical framework to modeling metabolic pathways and show through experiments that complex phenomenon such as biochemical reactions in cell’s metabolic networks can be modeled and simulated in the proposed framework.