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 [14:52 13.11.2017] – [Delta-Notch] 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 91: Line 115:
  
   * 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.1$, $n=4$,+  * 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 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>
  
 </file> </file>
  
documentation/course/module1.1510581160.txt.gz · Last modified: 14:52 13.11.2017 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki