reproduction of fig 4D from the publication the integration was performed using SBMLodeSolver (20080507)