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 [12:42 13.11.2017] – [Morpheus models] 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 108: 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'?>
Line 163: Line 187:
 </MorpheusModel> </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>
 +    <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>
  
-</code>+</file>
  
documentation/course/module1.1510573338.txt.gz · Last modified: 12:42 13.11.2017 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki