Introducing a new reaction to a model

[Maria Carolina Sisco]

Hello,

I'm adding a new reaction to an existing model. I already added the metabolites and the reaction with Cobrapy., no problem there. My question is how to properly add these new metabolites (consumed and produced) to the biomass reaction?

Also, my model has a reaction named GROWTH and when I checked, this is the output:

model.reactions.get_by_id('GROWTH')

Reaction identifier
GROWTH
Name
Artificial reaction
Memory address
0x7e768fd8ae90
Stoichiometry

61.0 ATP + 0.38 CELLWALL + 0.00244 DNA + 61.0 H2O + LIPIDS + POOL + 0.263 Protein + 0.01814 RNA --> 61.0 ADP + BIOMASS + 61.0 PI

61.0 ATP + 0.38 Cellwall (1 g dry cellwall) + 0.00244 Deoxyribonucleic acid + 61.0 Water + Artificial metabolite. The metabolites found in the biomass. Measured in grams + Artificial metabolite. A...

Lower bound
0.0
Upper bound
1000.0

I appreciate any help. Thanks!