reaction.add_metabolites works now with SolverBasedModel