User Tools

Site Tools


documentation:course:module1

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
documentation:course:module1 [10:38 11.01.2013] Fabian Rostdocumentation:course:module1 [16:20 13.11.2017] (current) – [Delta-Notch] Fabian Rost
Line 1: Line 1:
  
-<ditaa> 
-    +--------+   +-------+    +-------+ 
-    |        | --+ ditaa +--> |       | 
-    |  Text  |   +-------+    |diagram| 
-    |Document|   |!magic!|    |       | 
-    |     {d}|    c978 |    |       | 
-    +---+----+   +-------+    +-------+ 
-        |                         ^ 
-        |       Lots of work      | 
-        +------------+------------+ 
-                     | 
-                     | 
-                     O 
-</ditaa> 
  
  
-=== Module 1: Placeholder ===+====== Module 1: Delta-Notch (ODE systems, ODEs on a grid) ====== 
  
 Author: Fabian Rost Author: Fabian Rost
  
-Aim:+===== Aim ===== 
   * learn about ODE models (dynamics in morpheus, steady states analytically)   * learn about ODE models (dynamics in morpheus, steady states analytically)
   * develop first models   * develop first models
-Description:+ 
 +===== Description ===== 
 +==== Basic ODEs ==== 
   * get to know what students know about ODEs and adjust the module to the pre-knowledge   * get to know what students know about ODEs and adjust the module to the pre-knowledge
-  * give them very simple sketches of biomolecular models, which they should translate into ODEs+  * give them very simple sketches of biomolecular models, which they should translate into ODEs, e.g.
  
  
-  example ODEs:+ 
 + 
 +<graphviz> 
 +digraph finite_state_machine { 
 + rankdir=LR; 
 + size="9,5" 
 + 
 + node [shape = circle]; 
 +        a[label="", fixedsize="false", width=0, height=0, shape=none]; 
 + a -> A [ label = "k1" ]; 
 +
 +</graphviz> 
 +<graphviz> 
 +digraph finite_state_machine { 
 + rankdir=LR; 
 + size="9,5" 
 + 
 + node [shape = circle]; 
 + A -> A [ label = "k2" ]; 
 +
 +</graphviz> 
 +<graphviz> 
 +digraph finite_state_machine { 
 + rankdir=LR; 
 + size="9,5" 
 + 
 + node [shape = circle]; 
 +        a[label="", fixedsize="false", width=0, height=0, shape=none]; 
 + A -> a [ label = "k3" ]; 
 +
 +</graphviz> 
 +<graphviz> 
 +digraph finite_state_machine { 
 + rankdir=LR; 
 + size="9,5" 
 + 
 + node [shape = circle]; 
 +        a[label="", fixedsize="false", width=0, height=0, shape=none]; 
 +        b[label="", fixedsize="false", width=0, height=0, shape=none]; 
 +        b -> A [ label = "k4" ]; 
 + A -> a [ label = "k5" ]; 
 +         
 +
 +</graphviz> 
 + 
 + 
 + 
 +  could be translated to the following ODEs:
  
 \begin{align} \begin{align}
-\dot A &\\ +\dot A &k_1 \\ 
-\dot A & A \\ +\dot A &k_2  A \\ 
-\dot A &k * +\dot A &= - k_3 A \\ 
 +\dot A &= k_4 - k_5 
 \end{align} \end{align}
  
   * discuss those ODEs by    * discuss those ODEs by 
-    * steady state analysis+    * calculate steady state (do not calculate the stability, to complicated for biologists)
     * simulate in morpheus     * simulate in morpheus
 +  * typical models for biochemistry/systems biology use hill kinetics:
 +
 +<graphviz>
 +digraph {
 + size="9,5"
 +
 + node [shape = circle];
 +  a[label="", fixedsize="false", width=0, height=0, shape=none];
 +  ab[label="", fixedsize="false", width=0, height=0, shape=none];
 +  c[label="", fixedsize="false", width=0, height=0, shape=none];
 +
 +  a -> ab[arrowhead=None];
 +  ab -> A;
 +  A -> c
 +  B -> ab[arrowhead=tee];
 +
 +  {rank=same; a; ab; A; c};
 +}
 +
 +</graphviz>
 +
 +\begin{align}
 +\dot A &= \frac{\theta^n}{\theta^n + B^n} - k A
 +\end{align}
 +==== Delta-Notch ====
 +
   * then discuss the delta-notch sketch with two species   * then discuss the delta-notch sketch with two species
     * start with the Collier model     * start with the Collier model
Line 45: Line 110:
     * they could come up with something like:     * they could come up with something like:
 \begin{align} \begin{align}
-  \dot X_1 &\frac{\theta^n}{\theta^n + X_2^n} - k X_1 \\ +  \dot A_1 &= \frac{\theta^n}{\theta^n + A_2^n} - A_1 \\ 
-  \dot X_2 &\frac{\theta^n}{\theta^n + X_1^n} - k X_2+  \dot A_2 &= \frac{\theta^n}{\theta^n + A_1^n} - A_2
 \end{align} \end{align}
  
   * this system is bistable for certain parameter ranges, if the students are advanced they might find this out themselves   * this system is bistable for certain parameter ranges, if the students are advanced they might find this out themselves
-  * bistable e.g. for $\theta=0.5$, $n=4$, $c=k=1$+  * bistable e.g. for $\theta=0.5$, $n=4$,
   * if they have this system running in morpheus go spatial and let them simulate the system on a square and hexagonal grid   * if they have this system running in morpheus go spatial and let them simulate the system on a square and hexagonal grid
   * then you could also move to shaped cpm cells or even moving cells   * then you could also move to shaped cpm cells or even moving cells
Line 56: Line 121:
   * students won't do so much on their own in this session, it is a lot teaching on ODEs (don't be theoretical here, not enough time!) and introducing morpheus   * students won't do so much on their own in this session, it is a lot teaching on ODEs (don't be theoretical here, not enough time!) and introducing morpheus
 Paper: Paper:
-  * Collier, JR., McInerney, D., Schnell, S., MainiP. K.GavaghanD. J.HoustonP.& SternCD(2000)A cell cycle model for somitogenesis: mathematical formulation and numerical simulationJTheor. Biol., 207(3), 305-316.+  * Joanne R. Collier, Nicholas A.MMonkPhilip KMainiJulian HLewisPattern Formation by Lateral Inhibition with Feedback: a Mathematical Model of Delta-Notch Intercellular SignallingJournal of Theoretical BiologyVolume 183Issue 421 December 1996Pages 429-446ISSN 0022-5193http://dx.doi.org/10.1006/jtbi.1996.0233.
  
  
  
-Morpheus models: 
  
-[[examples:miscellaneous#French Flag: Morphogen gradient|Link to Example page]] 
  
-<code xml |h MorpheusModel.xml |h> 
  
-<MorpheusModel version="1">+ 
 +===== Morpheus models ===== 
 + 
 + 
 +<file xml ExponentialGrowth.xml> 
 + 
 +<?xml version='1.0' encoding='UTF-8'?> 
 +<MorpheusModel version="3">
     <Description>     <Description>
-        <Title>Example-Model</Title>+        <Details></Details> 
 +        <Title></Title>
     </Description>     </Description>
     <Space>     <Space>
         <Lattice class="linear">         <Lattice class="linear">
-            <Size value="100 0 0"/>+            <Neighborhood> 
 +                <Order>1</Order> 
 +            </Neighborhood> 
 +            <Size symbol="size" value="1000, 0"/> 
 +        </Lattice> 
 +        <SpaceSymbol symbol="space"/> 
 +    </Space> 
 +    <Time> 
 +        <StartTime value="0"/> 
 +        <StopTime value="5"/> 
 +        <TimeSymbol symbol="time"/> 
 +    </Time> 
 +    <Global> 
 +        <Constant symbol="k1" value="1.0"/> 
 +        <Variable symbol="A" value="1.0"/> 
 +        <System solver="heun" time-step="0.01"> 
 +            <DiffEqn symbol-ref="A"> 
 +                <Expression>k1 * A</Expression> 
 +            </DiffEqn> 
 +        </System> 
 +    </Global> 
 +    <Analysis> 
 +        <Logger time-step="0.01"> 
 +            <Input> 
 +                <Symbol symbol-ref="A"/> 
 +            </Input> 
 +            <Output> 
 +                <TextOutput/> 
 +            </Output> 
 +            <Plots> 
 +                <Plot time-step="5"> 
 +                    <Style style="points"/> 
 +                    <Terminal terminal="png"/> 
 +                    <X-axis> 
 +                        <Symbol symbol-ref="time"/> 
 +                    </X-axis> 
 +                    <Y-axis> 
 +                        <Symbol symbol-ref="A"/> 
 +                    </Y-axis> 
 +                </Plot> 
 +            </Plots> 
 +        </Logger> 
 +    </Analysis> 
 +</MorpheusModel> 
 + 
 +</file> 
 + 
 +<file xml | LateralInhibition_2cell.xml> 
 + 
 +<?xml version='1.0' encoding='UTF-8'?> 
 +<MorpheusModel version="3"> 
 +    <Description> 
 +        <Details></Details> 
 +        <Title></Title> 
 +    </Description> 
 +    <Space> 
 +        <Lattice class="linear"> 
 +            <Neighborhood> 
 +                <Order>1</Order> 
 +            </Neighborhood> 
 +            <Size symbol="size" value="100, 0, 0"/> 
 +        </Lattice> 
 +        <SpaceSymbol symbol="space"/> 
 +    </Space> 
 +    <Time> 
 +        <StartTime value="0"/> 
 +        <StopTime value="8"/> 
 +        <TimeSymbol symbol="time"/> 
 +    </Time> 
 +    <Global> 
 +        <Variable symbol="A" value="0.51"/> 
 +        <Variable symbol="A_neighbour" value="0.5"/> 
 +        <System solver="heun" time-step="0.01"> 
 +            <DiffEqn symbol-ref="A"> 
 +                <Expression>theta^n / (theta^n + A_neighbour^n) 
 +- k * A</Expression> 
 +            </DiffEqn> 
 +            <DiffEqn symbol-ref="A_neighbour"> 
 +                <Expression>theta^n / (theta^n + A^n) 
 +- k * A_neighbour</Expression> 
 +            </DiffEqn> 
 +        </System> 
 +        <Constant symbol="theta" value="0.1"/> 
 +        <Constant symbol="n" value="6"/> 
 +        <Constant symbol="k" value="1.0"/> 
 +    </Global> 
 +    <Analysis> 
 +        <Logger time-step="0.1"> 
 +            <Input> 
 +                <Symbol symbol-ref="A"/> 
 +            </Input> 
 +            <Output> 
 +                <TextOutput/> 
 +            </Output> 
 +            <Plots> 
 +                <Plot time-step="8"> 
 +                    <Style style="points"/> 
 +                    <Terminal terminal="png"/> 
 +                    <X-axis> 
 +                        <Symbol symbol-ref="time"/> 
 +                    </X-axis> 
 +                    <Y-axis> 
 +                        <Symbol symbol-ref="A"/> 
 +                        <Symbol symbol-ref="A_neighbour"/> 
 +                    </Y-axis> 
 +                </Plot> 
 +            </Plots> 
 +        </Logger> 
 +    </Analysis> 
 +</MorpheusModel> 
 + 
 + 
 +</file> 
 + 
 +<file xml | LateralInhibition_space.xml> 
 + 
 +<?xml version='1.0' encoding='UTF-8'?> 
 +<MorpheusModel version="3"> 
 +    <Description> 
 +        <Details></Details> 
 +        <Title></Title> 
 +    </Description> 
 +    <Space> 
 +        <Lattice class="hexagonal"> 
 +            <Neighborhood> 
 +                <Order>1</Order> 
 +            </Neighborhood> 
 +            <Size symbol="size" value="10,10,0"/>
             <BoundaryConditions>             <BoundaryConditions>
                 <Condition boundary="x" type="periodic"/>                 <Condition boundary="x" type="periodic"/>
 +                <Condition boundary="-x" type="periodic"/>
 +                <Condition boundary="x" type="periodic"/>
 +                <Condition boundary="-y" type="periodic"/>
             </BoundaryConditions>             </BoundaryConditions>
         </Lattice>         </Lattice>
 +        <SpaceSymbol symbol="space"/>
     </Space>     </Space>
     <Time>     <Time>
         <StartTime value="0"/>         <StartTime value="0"/>
-        <StopTime value="100"/>+        <StopTime value="50"/> 
 +        <TimeSymbol symbol="time"/> 
 +        <RandomSeed value="42"/>
     </Time>     </Time>
 +    <CellTypes>
 +        <CellType class="biological" name="cells">
 +            <Property symbol="X" value="0.0"/>
 +            <Constant symbol="theta" value="0.1"/>
 +            <Constant symbol="n" value="6"/>
 +            <Property symbol="X_neighbour" value="0.0"/>
 +            <NeighborhoodReporter>
 +                <Input scaling="cell" value="X"/>
 +                <Output symbol-ref="X_neighbour" mapping="average"/>
 +            </NeighborhoodReporter>
 +            <System solver="heun" time-step="0.01">
 +                <DiffEqn symbol-ref="X">
 +                    <Expression>theta^n/(theta^n+X_neighbour^2)
 +- X</Expression>
 +                </DiffEqn>
 +            </System>
 +        </CellType>
 +    </CellTypes>
 +    <CellPopulations>
 +        <Population size="1" type="cells">
 +            <InitCellLattice/>
 +            <InitProperty symbol-ref="X">
 +                <!--    <Disabled>
 +        <Expression>if(cell.id == 1, 1, 0)</Expression>
 +    </Disabled>
 +-->
 +                <Expression>rand_uni(0, 0.01)</Expression>
 +            </InitProperty>
 +        </Population>
 +    </CellPopulations>
 +    <Analysis>
 +        <Gnuplotter time-step="50">
 +            <Plot>
 +                <Cells value="X"/>
 +            </Plot>
 +            <Terminal name="png"/>
 +        </Gnuplotter>
 +    </Analysis>
 +    <Global>
 +        <Constant symbol="X" value="0.0"/>
 +    </Global>
 </MorpheusModel> </MorpheusModel>
  
-</code>+</file>
  
documentation/course/module1.1357897130.txt.gz · Last modified: 10:38 11.01.2013 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki