Differences between revisions 1 and 2
⇤ ← Revision 1 as of 2014-08-28 11:41:17 →
Size: 1741
Editor: mark
Comment:
← Revision 2 as of 2014-09-10 07:30:51 → ⇥
Size: 1799
Editor: mark
Comment:
Deletions are marked like this. Additions are marked like this.
Line 2: Line 2:
== MCA of a 4 step pathway ==
 1.
 Create a new directory and download the [[http://mudsharkstatic.brookes.ac.uk/AccliPhot/Workshop2/Models/feedback4.spy|feedback4.spy]] model.
Line 3: Line 6:
== MCA of a 4 step pathway ==   . This is a 4 step pathway with end product inhibition on the first step.
Line 5: Line 8:
 1. Create a new directory and download the [[ http://mudsharkstatic.brookes.ac.uk/Delhi2013/Models/feedback4.spy | feedback4.spy ]] model.
    This is a 4 step pathway with end product inhibition on the first step.
Line 15: Line 16:

Line 16: Line 19:
 
 1. What relationship can you identify about the effect that any one {{{Vmax}}} value has on all reaction fluxes. Explain why this relationship exists.
 1.
 What relationship can you identify about the effect that any one {{{Vmax}}} value has on all reaction fluxes. Explain why this relationship exists.
Line 20: Line 25:
 1. The effect of end-product inhibition on control coefficients. Here we wish to scan across a very wide range of values for the inhibition constant of the first reaction, so we will increase it proportionately, and plot on a log scale. (In the code below, {{{x *= y}}} is a convenient short hand for {{{ x = x*y}}}).  1.
 The effect of end-product inhibition on control coefficients. Here we wish to scan across a very wide range of values for the inhibition constant of the first reaction, so we will increase it proportionately, and plot on a log scale. (In the code below, {{{x *= y}}} is a convenient short hand for {{{ x = x*y}}}).
Line 26: Line 33:
 m["K1_S"] *= 1.4
 K1_S = m["K1_S"]
 CJ1 = m.ScaledSensits("VM1", ["R1"])
 CJ4 = m.ScaledSensits("VM4", ["R1"])
 ds.NewRow([K1_S]+CJ1+CJ4)		         m["K1_S"] *= 1.4
        K1_S = m["K1_S"]
        CJ1 = m.ScaledSensits("VM1", ["R1"])
        CJ4 = m.ScaledSensits("VM4", ["R1"])
        ds.NewRow([K1_S]+CJ1+CJ4)
Line 35: Line 42:

2 Metabolic Control Analysis with ScrumPy

MCA of a 4 step pathway

  1. Create a new directory and download the feedback4.spy model.

    • This is a 4 step pathway with end product inhibition on the first step.
  2. Identify the Vmax parameters of each reactions.
  3. Calculate the flux control coefficient of every Vmax over every flux in the system: 
       1 >>> for v in VmaxNames:
   2         print m.ScaledSensits(v, m.sm.cnames)

(you will need to specify VmaxNames - don't just cut and paste)

  1. What relationship can you identify about the effect that any one Vmax value has on all reaction fluxes. Explain why this relationship exists.

  2. Assume that this is a metabolic engineering problem. Try and achieve a 10 fold increase in flux by manipulating only the Vmax values. Can you do this whilst minimising the impact on the concentration of the intermediates ?

  3. The effect of end-product inhibition on control coefficients. Here we wish to scan across a very wide range of values for the inhibition constant of the first reaction, so we will increase it proportionately, and plot on a log scale. (In the code below, x *= y is a convenient short hand for  x  = x*y). 

       1 >>> from Data import Datasets
   2 >>> ds = DataSets.DataSet(ItemNames  = ["K1_S", "CJ1", "CJ4"])
   3 >>> m["K1_S"] = 1e-5
   4 >>> for n in range(40):
   5         m["K1_S"] *= 1.4
   6         K1_S = m["K1_S"]
   7         CJ1 = m.ScaledSensits("VM1", ["R1"])
   8         CJ4 = m.ScaledSensits("VM4", ["R1"])
   9         ds.NewRow([K1_S]+CJ1+CJ4)
  10 
  11 >>> ds.Plotter.SetLog()
  12 >>> ds.AddToPlot(["CJ1","CJ4"])

(Note: you only need to import DataSets once per session.)

None: AccliPhot/WorkshopTwo/Practical_2 (last edited 2014-09-10 13:24:06 by mark)