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 [13:11 06.11.2013] – [Basic ODEs] Fabian Rostdocumentation:course:module1 [16:20 13.11.2017] (current) – [Delta-Notch] Fabian Rost
Line 78: Line 78:
     * calculate steady state (do not calculate the stability, to complicated for biologists)     * 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 ==== ==== Delta-Notch ====
  
Line 86: 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 97: 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. 
 + 
 + 
 + 
 + 
  
  
Line 103: Line 132:
  
  
-<code xml |h ExponentialGrowth.xml |h>+<file xml ExponentialGrowth.xml>
  
 <?xml version='1.0' encoding='UTF-8'?> <?xml version='1.0' encoding='UTF-8'?>
-<MorpheusModel version="1">+<MorpheusModel version="3">
     <Description>     <Description>
-        <Title>Exponential Growth</Title>+        <Details></Details> 
 +        <Title></Title>
     </Description>     </Description>
     <Space>     <Space>
         <Lattice class="linear">         <Lattice class="linear">
-            <Size value="0"/>+            <Neighborhood> 
 +                <Order>1</Order> 
 +            </Neighborhood> 
 +            <Size symbol="size" value="100, 00"/>
         </Lattice>         </Lattice>
 +        <SpaceSymbol symbol="space"/>
     </Space>     </Space>
     <Time>     <Time>
         <StartTime value="0"/>         <StartTime value="0"/>
-        <StopTime value="100"/>+        <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> 
 +                <Condition boundary="x" type="periodic"/> 
 +                <Condition boundary="-x" type="periodic"/> 
 +                <Condition boundary="x" type="periodic"/> 
 +                <Condition boundary="-y" type="periodic"/> 
 +            </BoundaryConditions> 
 +        </Lattice> 
 +        <SpaceSymbol symbol="space"/> 
 +    </Space> 
 +    <Time> 
 +        <StartTime value="0"/> 
 +        <StopTime value="50"/> 
 +        <TimeSymbol symbol="time"/> 
 +        <RandomSeed value="42"/>
     </Time>     </Time>
     <CellTypes>     <CellTypes>
-        <CellType class="biological" name="cell"> +        <CellType class="biological" name="cells"> 
-            <System time-step="1.0"> +            <Property symbol="X" value="0.0"/> 
-                <DiffEqn symbol-ref="A"> +            <Constant symbol="theta" value="0.1"/> 
-                    <Expression>k*A</Expression>+            <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>                 </DiffEqn>
             </System>             </System>
-            <Property symbol="A" value="1.0"/> 
-            <Constant symbol="k" value="0.1"/> 
         </CellType>         </CellType>
     </CellTypes>     </CellTypes>
     <CellPopulations>     <CellPopulations>
-        <Population size="1" type="cell"/>+        <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>     </CellPopulations>
     <Analysis>     <Analysis>
-        <Logger interval="1"> +        <Gnuplotter time-step="50"> 
-            <Format string="A"/> +            <Plot
-            <Input+                <Cells value="X"/>
-                <Cell mapping="all" celltype="cell"/> +
-            </Input> +
-            <Plot interval="10" terminal="wxt" persist="true"> +
-                <X-axis column="1"/> +
-                <Y-axis columns="3"/>+
             </Plot>             </Plot>
-        </Logger>+            <Terminal name="png"/
 +        </Gnuplotter>
     </Analysis>     </Analysis>
 +    <Global>
 +        <Constant symbol="X" value="0.0"/>
 +    </Global>
 </MorpheusModel> </MorpheusModel>
  
- +</file>
-</code>+
  
documentation/course/module1.1383739901.txt.gz · Last modified: 13:11 06.11.2013 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki