Difference between revisions of "The Magrathea Manual: Building Coordinated Agent Models"

From irefindex
Line 550: Line 550:
  
 
{|border="1" cellpadding="2"  
 
{|border="1" cellpadding="2"  
| align="center" style="background:grey" colspan="2"|'''if line format'''
+
| align="center" style="background:grey" colspan="3"|'''if line format'''
 
|-
 
|-
 
| position||Meaning and possible values||Example
 
| position||Meaning and possible values||Example
Line 569: Line 569:
  
 
mean that the rule applies for each copy of molecule type 2, with feature 1 in state 1.  And the following 2.4 (positions 4 - 5) would indicate that the molecule must have feature 2 in state 4.
 
mean that the rule applies for each copy of molecule type 2, with feature 1 in state 1.  And the following 2.4 (positions 4 - 5) would indicate that the molecule must have feature 2 in state 4.
 +
 +
{|border="1" cellpadding="2"
 +
| align="center" style="background:yellow" colspan="2"|'''[rule] competitive-with [/rule]'''
 +
|-
 +
| key||competitive-with
 +
|-
 +
| format||competitive-with [integer]
 +
|-
 +
| description||Used to point out that the present rule and the rule having the specified id can not occur simultaneously.
 +
|-
 +
| required||no
 +
|-
 +
| default value||there are no default values
 +
|-
 +
| example||competitive-with 1
 +
|}
 +
Notes: The competitive-with line specifies the identifier for a rule that is competitive to the current rule. That means that the rule including this line and the rule having the specified id can not occur simultaneously, even if the specified conditions in the if line were met, so the algorithm will choose one of them or neither of them, but never both. This is useful when having molecules that can follow different reactions.
 +
 +
{|border="1" cellpadding="2"
 +
| align="center" style="background:yellow" colspan="2"|'''[rule] then [/rule]'''
 +
|-
 +
| key||then
 +
|-
 +
| format||then [string|integer] [string] [string|integer]
 +
|-
 +
| description||Used to define events that occur to a given molecule
 +
|-
 +
| required||no
 +
|-
 +
| default value||there are no default events
 +
|-
 +
| example||then 1 is-destroyed
 +
|}
 +
Notes:The then line is used to define events that may occur to a given molecule if the condition defined by the if line is met.  There are several different possibilities for this line that are differentiated by the second string value.  They are:
 +
 +
<pre>
 +
then 6.1 is-created-at-location-of 6
 +
</pre>
 +
 +
In this case, 1 instance of molecule type 6 is created at the same location of the molecule.
 +
 +
<pre>
 +
then 1 is-destroyed
 +
</pre>
 +
 +
In this case, the current molecule (1) is destroyed.
 +
 +
<pre>
 +
then 1 changes-state-to 1.2
 +
</pre>
 +
 +
In this case, the current molecule (tag 1) will change state according to the state description 1.2 (feature 1 changes state to state 2). Multiple feature–state changes may be specified on one line; i.e., 1.1.2.4.5.6.33 means that feature 1 assumes state 2, feature 4 assumes state 5 and feature 6 assumes state 33. 
 +
 +
<pre>
 +
then 1 changes-colour-to .1:.32:.11
 +
</pre>
 +
 +
In this case, the molecule will change to the colour specified by .1:.32:.11 (the colons separate three float values in the RGB colour description; i.e. .1 .32 .11).
 +
 +
Currently, translocations are not supported as events for molecule rules (as they are for coordinator rules).

Revision as of 12:39, 9 March 2009

Making models

The model.txt file

Magrathea opens a file called “model.txt” at the beginning of a simulation. Information is read from the file to populate the simulation with molecules and coordinators. The properties of different types of molecules and coordinators are specified in the model.txt file along with the number and positions of all of the starting molecules and coordinators. The sections below describe the format of the model.txt file.


A simple example

An example model.txt file is shown below in Figure 1.

The model.txt file is a text file and must be called “model.txt”. Lines beginning with # are comments and are ignored. Blank lines are ignored. Any text in the file after the line

EOF

will be ignored.

There are two types of sections in a model file. A section used to describe a type of molecule and create instances of that molecule type at the beginning of the simulation begins with the line:

[moleculeType]

and ends with the line

[/moleculeType]

A section that is used to describe a coordinator type and create instances of that coordinator at the beginning of the simulation begins with the line:

[coordinatorType]

and ends with the line

[/coordinatorType]

An example model.txt file for Magrathea

###########
#ab complex
###########

[moleculeType]
id 1
label a
color .81 .60 .21
shape Sphere 1
neighborhood-size 5
coordinator 1 3
instantiate 10 random spherical volume 10 30 10 10
[/moleculeType]

[moleculeType]
id 2        
label b
color 0.67 0.08 0.37
shape Sphere 1
neighborhood-size 5
coordinator 1 3
instantiate 4 random spherical volume 10 15 10 10
[/moleculeType]

[coordinatorType]
id 3
label a-b-complex
shape Cube 1 1 1
molecule 1 1
molecule 2 1
[/coordinatorType]

EOF


Key-value format

Each line in these two sections types encodes a separate piece of information as a key-value pair. So, for instance, members of the first molecule type are labeled “a” as indicated by the line

label a

In this line “label” is the key and “a” is the value. Keys and their values are separated by white-spaces. There may be several values specified for a given key. The key is always the first block of text on a line. A brief description of the example model.txt file

The above example model file, specifies a simulation that has two molecule types and one coordinator type.

Each molecule of type “1” has a label of “a”, has some specified color, is drawn as a sphere, can detect neighboring molecules within a 5 unit radius and can be coordinated by some coordinator type 3 on channel 1. These ideas are encoded by the lines

id 1
label a
color .81 .60 .21
shape Sphere 1
neighborhood-size 5
coordinator 1 3

The simulation begins with 10 instances of this molecule type distributed randomly throughout some specified spherical volume. This is encoded by the line:

instantiate 10 random spherical volume 10 30 10 10

The second molecule type has a label of “b”, is also drawn as a sphere with a different color. The simulation begins with only four instances of this molecule type that are also distributed in a spherical volume at a separate location in the simulation. Both molecules a and b can be coordinated by coordinator type 3 on channel 1 (as specified in the line):

coordinator 1 3

for both molecule a and b.


There is only one type of coordinator in this simulation. It has an id of 3 and is labeled “a-b-complex”. The simulation does not begin with any instances of this coordinator type. Instead, an instance of this coordinator type will be created once an instance of molecule “a” and an instance of molecule “b” diffuse within one another’s “neighborhood distance” during the simulation. The coordinator will then appear between the two molecules as a cube. The coordinator has two molecule type members (each with a stoichiometry of 1) as indicated by the lines:

molecule 1 1
molecule 2 1

Once the two molecules become a part of this coordinator, they will move towards the coordinator (at some position between where the two molecules started) and remain there.

A screen shot of the resulting simulation is shown below in the Figure below. The remainder of this section of the manual shows exactly how molecule and coordinator types can be specified.


Screen shot of a simulation resulting from the example model.txt file Insert screen shot here

Describing molecule types

An example molecule type description is given below in Figure 3. Each of the key-value entry lines are described in a separate section below. For each key-value pair, a format is given. Square brackets indicate types of values that must are to be inserted into the line by the user (i.e., either a float (like 0.1), an integer (like 200) or a string value (like myComplex).


An example molecule type description

[moleculeType]
id 1
label molecule_a
color .81 .60 .21
shape Sphere 1
draw-as-point 0
opacity 1
show-neighbor-lines 0
neighborhood-size 5
dynamic-neighborhood-sizing 0
feature 234 Tyr234
feature-state 234 0 unphosphorylated
feature-state 234 1 phosphorylated
coordinator 1 3
resident-compartment-type 4
instantiate 10 random spherical volume 10 30 10 10
[/moleculeType]

Molecule parameters

id
key id
format id [integer]
description A unique identifier for this molecule type
required yes
default value not applicable
example id 1

Notes: id is usually the first key-value pair specified for a moleculeType. This id is used to refer to this molecule type in other key-value pairs of the model file.


label
key label
format label [string]
description A label for the molecule type.
required no
default value molecule
example label molecule_a

Notes: The label line allows the user to give a name to all molecules of this molecule type. The name will not appear during the simulation unless show-label is set to TRUE (see below). Including spaces in a molecule type name will truncate the name at the space.


show-label
key show-label
format show-label [1 or 0]
description Boolean to determine if label should be shown.
required no
default value 0
example show-label 1

Notes: See label (above). All other values besides “1” will be interpreted as FALSE.


color
key color
format color [0…1] [0…1] [0…1]
description A description of the color used to draw this molecule type.
required no
default value some randomly chosen color
example color .81 .60 .21

Notes: A color may be specified by the user for drawing molecules of this molecule type.


shape
key shape
format shape Sphere [float] [float]
or
shape Cube [float] [float] [float]
description A description of the shape used to draw this molecule type.
required no
default value Sphere 1.0
example shape Sphere 0.5
example shape Sphere 0.5 2
example shape Cube 1.0 2.0 3.5

Notes:

Two shapes may be used to draw a molecule, a sphere or a cube.

If a Sphere is chosen, at least one additional float value must be given to specify the radius of the sphere. A second, optional float value may be given to specify the “effective radius” of the molecule. “Effective radius” is used to maintain distances between “molecules” that act as compartment markers. Compartment markers are used in place of representing all of the molecules that are really present in a membrane. They can be thought of as representing a single molecule at the center of a circular collection of molecules where the radius of the circle is the “effective radius”. Compartment markers will maintain a distance from one another equivalent to (er1 + er2) where er1 and er2 are the effective radii of compartment markers 1 and 2. Other molecules (non compartment markers) will maintain a distance from compartment markers of (r1 + r2) where r1 and r2 are the radii of the compartment marker and the molecule respectively. A second float value must be specified after the shape key in order for the molecule to behave as a compartment marker.


draw-as-point
key draw-as-point
format draw-as-point [1 or 0]
description A Boolean flag indicating that the molecule is to be drawn as a point.
required no
default value draw-as-point 0
example draw-as-point 1

Notes: If draw–as-point is set to 1 (TRUE) then the molecule will appear in the simulation as a point. The radius and effective radius will still be as specified in the shape line. This flag is used to speed up the simulation by reducing the amount of required rendering. It may also be used to hide molecule types that are not of interest to the user. Values other than 1 are interpreted as 0 (FALSE).

opacity
key opacity
format opacity [0…1]
description A description of how opaque the molecule should appear (where 0 is transparent and 1 is totally opaque).
required no
default value opacity 1
example opacity 0.1

Notes: Opacity can be set to some value below 1 in order to make the molecule somewhat transparent. This allows the user to see through molecules of a certain type that might otherwise block views to other parts of the simulation.


show-neighbor-lines
key show-neighbor-lines
format show-neighbor-lines [1,0]
description A Boolean flag indicating that lines are to be drawn between the molecule and its neighbor (as determined by the molecule’s neighborhood-size).
required no
default value show-neighbor-lines 0
example show-neighbor-lines 1

Notes: If show-neighbor-lines is set to 1 (TRUE) then lines will appear between the molecule and its neighbors. This can be useful when debugging a simulation. Values other than 1 are interpreted as 0 (FALSE).


neighborhood-size
key neighborhood-size
format neighborhood-size [float]
description A description of the molecule’s neighborhood size.
required no
default value neighborhood-size 1.0
example neighborhood-size 2.5

Notes: The neighborhood-size is a float value that specifies the radius of a sphere that has the same centre as the molecule. This value specifies how far the molecule “looks” when finding other neighbor molecules. When a molecule is looking for neighboring coordinators, it will only consider coordinators within this radius AND any coordinators of neighboring molecules.


dynamic-neighborhood-sizing
key dynamic-neighborhood-sizing
format dynamic-neighborhood-sizing [1,0]
description A Boolean flag indicating if this molecule is to use dynamic-neighborhood-sizing.
required no
default value dynamic-neighborhood-sizing 0
example dynamic-neighborhood-sizing 1

Notes: If dynamic-neighborhood-sizing is set to 1 (TRUE) then the molecule’s neighborhood size will increase by 1 unit per simulation iteration until all of its channels are occupied by some coordinator. Values other than 1 are interpreted as 0 (FALSE).


feature
key feature
format feature [integer] [string]
description A feature is described using a feature identifier and a label. Each feature may have one of multiple possible states. These possible states are described using the feature-state line.
required no
default value There are no default values
example feature 234 Tyr234

Notes: A molecule may have multiple features. Each feature is described using a separate feature line. Each feature is described using a feature id (integer value following the feature key word) followed by a label (string value with no spaces). The feature is not associated with any state information. For example, the line

feature 234 Tyr234

simply indicates that there is some feature (with id 234 and label Tyr234) that can exist in one of many multiple states (see feature-state below). For protein molecules, amino acids numbers could be used as feature identifiers for their corresponding amino acids; however, this is not required. Any integer may be used for a feature identifier so long as that integer is used only once for any given molecule.

N.B. features and feature-values are not used by Magrathea right now, but may be used in future to present information to the user. They do not have to be specified as a prerequisite to using state information in a coordinator rule (see below).


feature-state
key feature-state
format feature-state [integer] [integer] [string]
description Each feature described by a feature line may have one of multiple possible states. These possible states are described using the feature-state line.
required no
default value There are no default values
example feature-state 234 0 unphosphorylated
feature-state 234 1 phosphorylated

Notes: A molecule may have multiple features. Each feature is described using a separate feature line (above). Each feature can exist in one of many multiple states. Each of these possible states is described by a separate feature-state line. The first integer indicates the feature identifier that this state applies to. The second integer is a state-identifier indicating one possible state for this feature. An integer identifier should be unique to a feature on a molecule. The string value provides a label that describes this particular state.

Note: features and feature-values are not used by Magrathea, but may be used in future to present information to the user. They do not have to be specified as a prerequisite to using state information in a coordinator rule (see below).


coordinator
key coordinator
format coordinator [integer] [integer]
description A description of the coordinators that a molecule can “listen” to.
required no
default value there are no default coordinator assignments
example coordinator 1 3

Notes: Molecules of a given molecule type can be coordinated by a coordinator of a given coordinator type. Molecules may have multiple coordinators. A molecule joins a coordinator using one of its “channels”. A channel may be occupied by only one coordinator at a time. However, the same channel may be occupied by one of several different coordinator types. Channels could be thought of as binding sites and coordinators may be thought of as interactions (though the analogy is not applicable to all molecule coordinator pairings, hence the use of the “coordinator and channel” terminology). A binding site may be involved in one of several different interactions (i.e. with different binding partners). Therefore, multiple coordinator lines may be listed for a single molecule type. For instance, the lines:

coordinator 1 3
coordinator 1 4
coordinator 2 5

would be used to specify a molecule type that could pair (mutually exclusively) with either coordinator 3 or 4 using channel 1 and (at the same time) could pair with coordinator 5 using channel 2. The abstract ideas of channels and coordinator are more thoroughly explained in previous sections of this manual.


resident-compartment-type
key resident-compartment-type
format resident-compartment-type [integer]
description Specifies that a molecule is a resident of some coordinator that represents a compartment
required no
default value there are no default assignments
example resident-compartment-type 4

Notes: Spherical point coordinators may be used to represent compartments. At the beginning of a simulation, any molecules found inside a coordinator compartment will be kicked out unless this flag has been set indicating that the molecule may reside in that compartment type (the coordinator id is used to specify the compartment type). If a molecule is a resident of a compartment, it will move randomly within the space of that compartment until it encounters the inside face of the compartment. Collisions with the inside face of a compartment will result in the molecule being repelled (contained) back towards the center of the compartment.

It is very important that the isSizeDynamic flag be set to 1 for coordinators that are used to represent compartments. If not, the size and shape of the coordinator rendered in the simulation is meaningless (it is drawn only to mark the point-location of the coordinator). When isSizeDynamic is TRUE, the coordinator is drawn as a sphere and its size reflects the radius of the imaginary sphere around the coordinator centre that maintains the offset distance of member molecules from the coordinator. This offset distance is also used to determine whether molecules are inside or outside the compartment.


[rule] …[/rule]
key See below. Multiple keys are associated with this subsection including id, if, then, and parameters
format
[rule]
id [integer]
if [string]
.
.
competitive-with [integer]
.
.
then [string|integer] [string] [string|integer]
.
.
parameters [string] [float]
.
.
[/rule]
description Used to define events that occur to a molecule.
required no
default value there are no default rules
example
[rule]
id 1
if 1.1.1
competitive-with 2
then 1 is-destroyed
parameters p 0.001
[/rule]

Notes: The [rule] … [/rule] subsection of molecules is used to define events that occur to a single molecule. Multiple rule subsections may be associated with one molecule. Each rule is defined using multiple lines with keywords id, if, competitive-with, then and parameters. These keywords and their line formats are described below.

In general, each rule has:

1) an identifier that is unique to the molecule

2) a set of prerequisite conditions (specified by if lines) that must be met in order for

3) some set of events to be carried out (specified by then lines) according to

4) a set of parameters specified on one or more parameter lines.


See this example for an example model file that uses molecule rules.


[rule] id [/rule]
key id
format id [integer]
description Used to identify a rule for a molecule
required yes
default value there are no default rules
example id 1

Notes: The id line defines a single integer identifier for a rule. It must be the first line of the rule subsection and it must be unique to the rules defined for the molecule to which belongs to.


[rule] if [/rule]
key if
format if [string]
description Used to define a prerequisite condition for some event. The precondition is expressed as the molecule id with or without some feature-state specification
required no
default value there are no default values
example if 2.1.1

Notes: The if line is used to define a prerequisite condition for some event. The precondition is expressed as the molecule id with or without some feature-state specification. The line contains a single value that is a dot delimited list of integers. For example,

if 2.1.1.2.4

This dot-delimited condition is of the form:


if line format
position Meaning and possible values Example
1 Required molecule type (an integer corresponding to a molecule type id). See “id” key under [moleculeType] 2
2 + 2n Feature identifier specified by some integer. 1
3 + 2n Required state of feature (previous position) specified by some integer. 1

Note: n is an integer from 0, 1, 2… specifying a feature-state condition (see explanation below).

So the first three numbers in the example above

if 2.1.1

mean that the rule applies for each copy of molecule type 2, with feature 1 in state 1. And the following 2.4 (positions 4 - 5) would indicate that the molecule must have feature 2 in state 4.

[rule] competitive-with [/rule]
key competitive-with
format competitive-with [integer]
description Used to point out that the present rule and the rule having the specified id can not occur simultaneously.
required no
default value there are no default values
example competitive-with 1

Notes: The competitive-with line specifies the identifier for a rule that is competitive to the current rule. That means that the rule including this line and the rule having the specified id can not occur simultaneously, even if the specified conditions in the if line were met, so the algorithm will choose one of them or neither of them, but never both. This is useful when having molecules that can follow different reactions.

[rule] then [/rule]
key then
format integer] [string] [string|integer]
description Used to define events that occur to a given molecule
required no
default value there are no default events
example then 1 is-destroyed

Notes:The then line is used to define events that may occur to a given molecule if the condition defined by the if line is met. There are several different possibilities for this line that are differentiated by the second string value. They are:

then 6.1 is-created-at-location-of 6

In this case, 1 instance of molecule type 6 is created at the same location of the molecule.

then 1 is-destroyed

In this case, the current molecule (1) is destroyed.

then 1 changes-state-to 1.2

In this case, the current molecule (tag 1) will change state according to the state description 1.2 (feature 1 changes state to state 2). Multiple feature–state changes may be specified on one line; i.e., 1.1.2.4.5.6.33 means that feature 1 assumes state 2, feature 4 assumes state 5 and feature 6 assumes state 33.

then 1 changes-colour-to .1:.32:.11

In this case, the molecule will change to the colour specified by .1:.32:.11 (the colons separate three float values in the RGB colour description; i.e. .1 .32 .11).

Currently, translocations are not supported as events for molecule rules (as they are for coordinator rules).