diff --git a/xbanalysis/model/xba_reaction.py b/xbanalysis/model/xba_reaction.py
index 15edd1eefa250e4603b243a2c9963a6bd10949bb..1303bf69da14f3d9b7c6f2f8d3e8a39bf11a4510 100644
--- a/xbanalysis/model/xba_reaction.py
+++ b/xbanalysis/model/xba_reaction.py
@@ -65,7 +65,7 @@ class XbaReaction:
         #     to species being consumed add the enzyme product with corresponding stoichiometry
         #     to species being produced add the enzyme product with corresponding stoichiometry, unless
         #        the species is itself the enzyme product, in which case we add the ps reaction producing it.
-        if self.sboterm == 'SBO:0000589':
+        if (self.sboterm == 'SBO:0000589') or (self.sboterm == 'SBO:0000205') :
             enz_sid = ''
             for sid in self.products:
                 if (species[sid].sboterm == 'SBO:0000252') or (species[sid].sboterm == 'SBO:0000250'):
diff --git a/xbanalysis/problems/gba_problem.py b/xbanalysis/problems/gba_problem.py
index 0488eae2e111152a011d2ac0acef55eba2183fc2..1f242fbf915040da7cadad3cc6c565fc2ca77f19 100644
--- a/xbanalysis/problems/gba_problem.py
+++ b/xbanalysis/problems/gba_problem.py
@@ -41,7 +41,8 @@ class GbaProblem:
                         if ((r.sboterm == 'SBO:0000167') or
                             (r.sboterm == 'SBO:0000179') or
                             (r.sboterm == 'SBO:0000182'))]
-        self.ps_rids = [r.id for r in xba_model.reactions.values() if r.sboterm == 'SBO:0000589']
+        self.ps_rids = [r.id for r in xba_model.reactions.values()
+                        if (r.sboterm == 'SBO:0000589') or (r.sboterm == 'SBO:0000205')]
         spont_rids = [rid for rid in self.mr_rids if xba_model.reactions[rid].enzyme == '']
 
         # calculate stoichiometric for mass balance equality constraints