Differences between revisions 1 and 7 (spanning 6 versions)
Revision 1 as of 2014-08-28 11:41:17
Size: 1741
Editor: mark
Comment:
Revision 7 as of 2014-09-10 13:21:42
Size: 2025
Editor: mark
Comment:
Deletions are marked like this. Additions are marked like this.
Line 2: Line 2:
{{attachment:Slides.pdf|Slides.pdf|title=""}}

If the slides cannot been seen above, please download them [[http://mudsharkstatic.brookes.ac.uk/AccliPhot/Workshop2/Slides/ScrumPyMCA/Slides.pdf|here]]
Line 4: Line 7:
 1.
 Create a new directory and download the [[http://mudsharkstatic.brookes.ac.uk/AccliPhot/Workshop2/Models/feedback4.spy|feedback4.spy]] model.
Line 5: Line 10:
 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.
  . This is a 4 step pathway with end product inhibition on the first step.
Line 15: Line 20:

Line 16: Line 23:
 
 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 29:
 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 22: Line 33:
>>> from Data import Datasets >>> from Data import DataSets
Line 26: Line 37:
 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 32: Line 43:
>>> ds.Plotter.SetLog() >>> ds.Plotter.SetLog()# or ds.SetLogAxis()
Line 35: Line 46:

2 Metabolic Control Analysis with ScrumPy

Slides.pdf

If the slides cannot been seen above, please download them here

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()# or ds.SetLogAxis()
      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)