User Tools

Site Tools


documentation:course:module6

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:module6 [11:53 11.01.2013] – [Aim] Fabian Rostdocumentation:course:module6 [10:06 17.11.2017] (current) – [Description] Fabian Rost
Line 12: Line 12:
 ===== Description ===== ===== Description =====
  
-  * Step by step description of module +  * good description on what has been done can be found in the report by the students (see section Papers) 
-  * Written such that other tutors can teach the same stuff+  * it's been two days of teaching: 
 +    * first explain CPM and let them play with that a little (e.g. cell sorting) 
 +    * explain PDE system we are going to use 
 +      * I had to explain about phase plane analysis quite a bit 
 +  * then on day 3 start to couple PDE & CPM 
 +    * start with coupling in one direction only 
 +      * e.g. first give a fixed gradient and let CPM cells perform chemotaxis 
 +      * then formulate a PDE where a substance is only produced where CPM cells are 
 +    * then couple in both directions (implement the Savill model) 
 +      * perform clustering and slug movement experiment
    
 ===== Paper ===== ===== Paper =====
Line 20: Line 29:
   * Stan Marée (2002), From pattern formation to morphogenesis: multicellular coordination in //Dictyostelium discoideum// [[http://www-binf.bio.uu.nl/stan/Thesis/|link]]   * Stan Marée (2002), From pattern formation to morphogenesis: multicellular coordination in //Dictyostelium discoideum// [[http://www-binf.bio.uu.nl/stan/Thesis/|link]]
   * Alberto Quintero, Mirko Myllykoski, Anna Igolkina, Alexandra Rohde O’Sullivan Freltoft, Nitya Dixit, Fabian Rost (2012), Morphogenesis and Dynamics of Multicellular Systems, [[http://www.mafy.lut.fi/EcmiNL/issues.php?action=downart&ID=280|link]]   * Alberto Quintero, Mirko Myllykoski, Anna Igolkina, Alexandra Rohde O’Sullivan Freltoft, Nitya Dixit, Fabian Rost (2012), Morphogenesis and Dynamics of Multicellular Systems, [[http://www.mafy.lut.fi/EcmiNL/issues.php?action=downart&ID=280|link]]
 +
  
  
 ===== Morpheus models ===== ===== Morpheus models =====
  
 +<code xml |h ClusterFormation.xml |h>
 +<?xml version='1.0' encoding='UTF-8'?>
 +<MorpheusModel version="1">
 +    <Description>
 +        <Title>Cluster Formation</Title>
 +    </Description>
 +    <Space>
 +        <Lattice class="square">
 +            <Size symbol="size" value="200 200 0"/>
 +            <BoundaryConditions>
 +                <Condition boundary="x" type="constant"/>
 +                <Condition boundary="y" type="constant"/>
 +                <Condition boundary="-x" type="constant"/>
 +                <Condition boundary="-y" type="constant"/>
 +            </BoundaryConditions>
 +            <NodeLength value="0.37"/>
 +        </Lattice>
 +        <SpaceSymbol symbol="space"/>
 +    </Space>
 +    <Time>
 +        <StartTime value="0"/>
 +        <StopTime symbol="end" value="5000"/>
 +        <TimeSymbol symbol="time"/>
 +    </Time>
 +    <CellTypes>
 +        <CellType class="medium" name="medium">
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +        <CellType class="biological" name="ameba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <Chemotaxis>
 +                <Layer symbol-ref="c"/>
 +                <Strength symbol-ref="mu"/>
 +            </Chemotaxis>
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <PDEReporter>
 +                <Input symbol-ref="c"/>
 +                <Maximum symbol-ref="max_c"/>
 +            </PDEReporter>
 +            <Property symbol="phase" value="1"/>
 +            <Property symbol="phaseTime" value="0"/>
 +            <Equation symbol-ref="phaseTime">
 +                <Expression>if(phase == 2, phaseTime+MCStime, 0)</Expression>
 +            </Equation>
 +            <Equation symbol-ref="phase" name="Rule_2_3">
 +                <Expression>if(phase == 1 and max_c > 0.1, 2,
 +if(phase == 2 and phaseTime>phase2duration, 3,
 +if(phase == 3 and max_c &lt; 0.05, 1, 
 +phase
 +))) </Expression>
 +            </Equation>
 +            <Property symbol="phase2duration" value="2"/>
 +            <Function symbol="mu">
 +                <Expression>if(phase == 2, 300, 0)</Expression>
 +            </Function>
 +        </CellType>
 +        <CellType class="biological" name="autoAmeba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <!--
 +  <Disabled>
 +    <Chemotaxis>
 +        <Layer symbol-ref="c"/>
 +        <Strength symbol-ref="mu"/>
 +    </Chemotaxis>
 +  </Disabled>  -->
 +            <Property symbol="b" value="0.5"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +    </CellTypes>
 +    <CPM>
 +        <Interaction>
 +            <Contact type1="ameba" type2="ameba" value="4"/>
 +            <Contact type1="ameba" type2="medium" value="2"/>
 +            <Contact type1="ameba" type2="autoAmeba" value="4"/>
 +            <Contact type1="autoAmeba" type2="medium" value="3"/>
 +            <Contact type1="autoAmeba" type2="autoAmeba" value="-15"/>
 +        </Interaction>
 +        <MCSDuration symbol="MCStime" value="0.1"/>
 +        <MetropolisKinetics temperature="1.0" yield="0.1" stepper="edgelist">
 +            <Neighborhood>
 +                <Order>2</Order>
 +            </Neighborhood>
 +        </MetropolisKinetics>
 +    </CPM>
 +    <PDE>
 +        <Layer symbol="c">
 +            <Diffusion rate="1"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <Layer symbol="r">
 +            <Diffusion rate="0"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <System solver="euler" time-step="0.1">
 +            <DiffEqn symbol-ref="c">
 +                <Expression>alpha*(-f -r)</Expression>
 +            </DiffEqn>
 +            <DiffEqn symbol-ref="r">
 +                <Expression>alpha*epsilon*(3.5*c-b-r)</Expression>
 +            </DiffEqn>
 +            <Function symbol="epsilon">
 +                <Expression>if(c &lt; 0.0065, 0.5,
 +if(c &lt; 0.841, 0.0589,
 +0.5))</Expression>
 +            </Function>
 +            <Function symbol="f">
 +                <Expression>if(c &lt; 0.0065, 20*c, 
 +if(c &lt; 0.841, -3*c+0.15,
 +15*(c-1)))</Expression>
 +            </Function>
 +            <Function symbol="alpha">
 +                <Expression>if(cell.id != 0, 1, 0)</Expression>
 +            </Function>
 +        </System>
 +    </PDE>
 +    <CellPopulations>
 +        <Population size="1" type="autoAmeba">
 +            <InitRectangle cells="1" type="regular">
 +                <Dimensions size="10 10 0" center="100 100 0"/>
 +            </InitRectangle>
 +        </Population>
 +        <Population size="1" type="ameba">
 +            <InitRectangle cells="500" type="random">
 +                <Dimensions size="200 200 0" origin="0 0 0"/>
 +            </InitRectangle>
 +        </Population>
 +    </CellPopulations>
 +    <Analysis>
 +        <Gnuplotter interval="10">
 +            <Terminal opacity="0.7" name="png"/>
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +            <Cells symbol-ref="phase" min="0" max="3"/>
 +            <Cells symbol-ref="cell.volume" min="15" max="25"/>
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +        </Gnuplotter>
 +    </Analysis>
 +</MorpheusModel>
 +</code>
  
-[[examples:miscellaneous#French Flag: Morphogen gradient|Link to Example page]] 
  
-<code xml |h MorpheusModel.xml |h>+  * {{documentation:course:clusters_new_parameters.mp4|Movie: Cluster Formation & cAmp Spirals}} 
 +  * for the movie the system size was 500x500 with 4000 cells randomly placed, everything else was like in the example above
  
 +
 +<code xml |h CrawlingSlug.xml |h>
 +
 +<?xml version='1.0' encoding='UTF-8'?>
 <MorpheusModel version="1"> <MorpheusModel version="1">
     <Description>     <Description>
-        <Title>Example-Model</Title>+        <Title>test</Title>
     </Description>     </Description>
     <Space>     <Space>
-        <Lattice class="linear"> +        <Lattice class="square"> 
-            <Size value="100 0"/>+            <Size symbol="size" value="100 200 0"/>
             <BoundaryConditions>             <BoundaryConditions>
-                <Condition boundary="x" type="periodic"/>+                <Condition boundary="x" type="constant"/> 
 +                <Condition boundary="y" type="constant"/> 
 +                <Condition boundary="-x" type="constant"/> 
 +                <Condition boundary="-y" type="constant"/>
             </BoundaryConditions>             </BoundaryConditions>
 +            <NodeLength value="0.37"/>
         </Lattice>         </Lattice>
 +        <SpaceSymbol symbol="space"/>
     </Space>     </Space>
     <Time>     <Time>
-        <StartTime value="0"/> +        <StartTime value="1.0"/> 
-        <StopTime value="100"/>+        <StopTime symbol="end" value="5000"/> 
 +        <TimeSymbol symbol="time"/>
     </Time>     </Time>
 +    <CellTypes>
 +        <CellType class="medium" name="medium">
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +        <CellType class="biological" name="ameba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <Chemotaxis>
 +                <Layer symbol-ref="c"/>
 +                <Strength symbol-ref="mu"/>
 +            </Chemotaxis>
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <PDEReporter>
 +                <Input symbol-ref="c"/>
 +                <Maximum symbol-ref="max_c"/>
 +            </PDEReporter>
 +            <Property symbol="phase" value="1"/>
 +            <Property symbol="phaseTime" value="0"/>
 +            <Equation symbol-ref="phaseTime">
 +                <Expression>if(phase == 2, phaseTime+MCStime, 0)</Expression>
 +            </Equation>
 +            <Equation symbol-ref="phase" name="Rule_2_3">
 +                <Expression>if(phase == 1 and max_c > 0.1, 2,
 +if(phase == 2 and phaseTime>phase2duration, 3,
 +if(phase == 3 and max_c &lt; 0.05, 1, 
 +phase
 +))) </Expression>
 +            </Equation>
 +            <Property symbol="phase2duration" value="0.5"/>
 +            <Function symbol="mu">
 +                <Expression>if(phase == 2, 50, 0)</Expression>
 +            </Function>
 +        </CellType>
 +        <CellType class="biological" name="autoAmeba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <!--
 +  <Disabled>
 +    <Chemotaxis>
 +        <Layer symbol-ref="c"/>
 +        <Strength symbol-ref="mu"/>
 +    </Chemotaxis>
 +  </Disabled>  -->
 +            <Property symbol="b" value="0.5"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +    </CellTypes>
 +    <CPM>
 +        <Interaction>
 +            <Contact type1="ameba" type2="ameba" value="3"/>
 +            <Contact type1="ameba" type2="medium" value="2"/>
 +            <Contact type1="ameba" type2="autoAmeba" value="4"/>
 +            <Contact type1="autoAmeba" type2="medium" value="3"/>
 +            <Contact type1="autoAmeba" type2="autoAmeba" value="-2"/>
 +        </Interaction>
 +        <MCSDuration symbol="MCStime" value="0.01"/>
 +        <MetropolisKinetics temperature="1.0" yield="2" stepper="edgelist">
 +            <Neighborhood>
 +                <Order>2</Order>
 +            </Neighborhood>
 +        </MetropolisKinetics>
 +    </CPM>
 +    <PDE>
 +        <Layer symbol="c">
 +            <Diffusion rate="1"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <Layer symbol="r">
 +            <Diffusion rate="0"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <!--
 +  <Disabled>
 +    <Function symbol="b">
 +        <Expression>0</Expression>
 +    </Function>
 +  </Disabled>  -->
 +        <System solver="euler" time-step="0.1">
 +            <DiffEqn symbol-ref="c">
 +                <Expression>alpha*(-f -r)</Expression>
 +            </DiffEqn>
 +            <DiffEqn symbol-ref="r">
 +                <Expression>alpha*epsilon*(3.5*c-b-r)</Expression>
 +            </DiffEqn>
 +            <Function symbol="f">
 +                <Expression>if(c &lt; 0.0065, 20*c, 
 +if(c &lt; 0.841, -3*c+0.15,
 +15*(c-1)))</Expression>
 +            </Function>
 +            <Function symbol="alpha">
 +                <Expression>if(cell.id != 0, 1, 0)</Expression>
 +            </Function>
 +            <Function symbol="epsilon">
 +                <Expression>if(c &lt; 0.0065, 0.5,
 +if(c &lt; 0.841, 0.0589,
 +0.5))</Expression>
 +            </Function>
 +        </System>
 +        <!--
 +  <Disabled>
 +    <Equation symbol-ref="r">
 +        <Expression>if(cell.id != 0, r, 0)</Expression>
 +    </Equation>
 +  </Disabled>  -->
 +    </PDE>
 +    <CellPopulations>
 +        <Population size="2" type="autoAmeba">
 +            <InitRectangle cells="2" type="regular">
 +                <Dimensions size="5 5 0" center="50 50 0"/>
 +            </InitRectangle>
 +        </Population>
 +        <Population size="1" type="ameba">
 +            <InitRectangle cells="50" type="regular">
 +                <Dimensions size="35 50 0" center="50 25 0"/>
 +            </InitRectangle>
 +        </Population>
 +        <!--
 +  <Disabled>
 +    <Population size="2" type="autoAmeba">
 +        <InitRectangle cells="2" type="regular">
 +            <Dimensions size="5 5 0" center="55 50 0"/>
 +        </InitRectangle>
 +    </Population>
 +  </Disabled>  -->
 +    </CellPopulations>
 +    <Analysis>
 +        <Gnuplotter interval="10">
 +            <Terminal opacity="0.7" name="png"/>
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +            <Cells symbol-ref="phase" min="1" max="3">
 +                <ColorMap>
 +                    <Color value="1" color="yellow"/>
 +                    <Color value="2" color="green"/>
 +                    <Color value="3" color="red"/>
 +                </ColorMap>
 +            </Cells>
 +            <!--
 +  <Disabled>
 +    <Cells symbol-ref="cell.volume" min="15" max="25"/>
 +  </Disabled>  -->
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +            <!--
 +  <Disabled>
 +    <PDE symbol-ref="r" min="-0.5" max="3.0">
 +        <ColorMap>
 +            <Color value="-0.5" color="dark-violet"/>
 +            <Color value="0" color="white"/>
 +            <Color value="3.0" color="red"/>
 +        </ColorMap>
 +    </PDE>
 +  </Disabled>  -->
 +        </Gnuplotter>
 +        <!--
 +  <Disabled>
 +    <Logger interval="0.1">
 +        <Filename name="log"/>
 +        <Format string="c r"/>
 +        <Input>
 +            <PDE mapping="all" slice="1"/>
 +        </Input>
 +        <Plot interval="0.1" terminal="png">
 +            <X-axis column="1"/>
 +            <Y-axis columns="2 3"/>
 +        </Plot>
 +    </Logger>
 +  </Disabled>  -->
 +    </Analysis>
 </MorpheusModel> </MorpheusModel>
 +
  
 </code> </code>
 +
 +  * {{documentation:course:slug_new_parameters.mp4|Movie: Crawling Slug}}
  
documentation/course/module6.1357901598.txt.gz · Last modified: 11:53 11.01.2013 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki