|
171 | 171 | "output_type": "stream", |
172 | 172 | "text": [ |
173 | 173 | "build status: reconstruction complete\n", |
174 | | - "Inc. reactions: 33/101\n", |
| 174 | + "Inc. reactions: 20/60\n", |
175 | 175 | " - unclear: 0/0\n", |
176 | | - " - exclude: 32/100\n", |
| 176 | + " - exclude: 19/59\n", |
177 | 177 | " - low and medium: 0/0\n", |
178 | 178 | " - high: 1/1\n", |
179 | 179 | "\n" |
|
192 | 192 | "cell_type": "markdown", |
193 | 193 | "metadata": {}, |
194 | 194 | "source": [ |
195 | | - "The metric you see are for irreversible reactions. We can obtain the reversible reconstruction using:" |
| 195 | + "The metric you see are for reversible reactions. We can obtain the irreversible reconstruction metrics by using:" |
196 | 196 | ] |
197 | 197 | }, |
198 | 198 | { |
|
203 | 203 | }, |
204 | 204 | "outputs": [ |
205 | 205 | { |
206 | | - "data": { |
207 | | - "text/plain": [ |
208 | | - "33" |
209 | | - ] |
210 | | - }, |
211 | | - "execution_count": 6, |
212 | | - "metadata": {}, |
213 | | - "output_type": "execute_result" |
| 206 | + "name": "stdout", |
| 207 | + "output_type": "stream", |
| 208 | + "text": [ |
| 209 | + "build status: reconstruction complete\n", |
| 210 | + "Inc. reactions: 32/101\n", |
| 211 | + " - unclear: 0/0\n", |
| 212 | + " - exclude: 31/100\n", |
| 213 | + " - low and medium: 0/0\n", |
| 214 | + " - high: 1/1\n", |
| 215 | + "\n" |
| 216 | + ] |
214 | 217 | } |
215 | 218 | ], |
216 | 219 | "source": [ |
217 | | - "rec = opt.cobra_model(name=\"minimum\")\n", |
218 | | - "len(rec.reactions)" |
| 220 | + "print(opt.info(reversible=False))" |
219 | 221 | ] |
220 | 222 | }, |
221 | 223 | { |
|
245 | 247 | "output_type": "stream", |
246 | 248 | "text": [ |
247 | 249 | "build status: reconstruction complete\n", |
248 | | - "Inc. reactions: 39/102\n", |
| 250 | + "Inc. reactions: 24/61\n", |
249 | 251 | " - unclear: 0/0\n", |
250 | | - " - exclude: 37/100\n", |
| 252 | + " - exclude: 22/59\n", |
251 | 253 | " - low and medium: 0/0\n", |
252 | 254 | " - high: 2/2\n", |
253 | 255 | "\n" |
|
278 | 280 | "name": "stdout", |
279 | 281 | "output_type": "stream", |
280 | 282 | "text": [ |
281 | | - "38\n", |
| 283 | + "37\n", |
282 | 284 | "2pg <=> pep\n", |
283 | 285 | "gtp + oaa <=> gdp + pep\n" |
284 | 286 | ] |
285 | 287 | } |
286 | 288 | ], |
287 | 289 | "source": [ |
288 | | - "rec2 = opt.cobra_model(\"plus_pep\")\n", |
289 | | - "print(len(rec2.reactions))\n", |
290 | | - "use = rec2.metabolites.pep.reactions\n", |
| 290 | + "rec = opt.cobra_model(\"plus_pep\")\n", |
| 291 | + "print(len(rec.reactions))\n", |
| 292 | + "use = rec.metabolites.pep.reactions\n", |
291 | 293 | "for r in use: print(r.reaction)" |
292 | 294 | ] |
293 | 295 | }, |
294 | 296 | { |
295 | 297 | "cell_type": "markdown", |
296 | 298 | "metadata": {}, |
297 | 299 | "source": [ |
298 | | - "As we can see the algorithm actually included two different ways of producing pep, from 3pg and from oaa and several other reactions. This is due to automatic inclusion of redundant pathways (which is good since it gives your model some robustness). If we do not want that feature we can modify the parameter n in the CORDA initializer." |
| 300 | + "By default CORDA uses redundancy. This means, in case there are several minimal pathways to reach your objective, CORDA will include several of those (which is good since it gives your model some robustness). If we do not want that feature we can modify the parameter n in the CORDA initializer which denotes the maximum number of redundant pathways to include." |
299 | 301 | ] |
300 | 302 | }, |
301 | 303 | { |
|
309 | 311 | "name": "stdout", |
310 | 312 | "output_type": "stream", |
311 | 313 | "text": [ |
312 | | - "33\n", |
| 314 | + "31\n", |
313 | 315 | "gtp + oaa <=> gdp + pep\n" |
314 | 316 | ] |
315 | 317 | } |
|
318 | 320 | "opt = CORDA(mod, conf, met_prod=\"pep\", n=1)\n", |
319 | 321 | "opt.build()\n", |
320 | 322 | "\n", |
321 | | - "rec3 = opt.cobra_model(\"plus_pep_nored\")\n", |
322 | | - "print(len(rec3.reactions))\n", |
323 | | - "use = rec3.metabolites.pep.reactions\n", |
| 323 | + "rec_min = opt.cobra_model(\"plus_pep_nored\")\n", |
| 324 | + "print(len(rec_min.reactions))\n", |
| 325 | + "use = rec_min.metabolites.pep.reactions\n", |
324 | 326 | "for r in use: print(r.reaction)" |
325 | 327 | ] |
326 | 328 | }, |
|
0 commit comments