Practical 1

Part 1 - Exploring Python

Before staring with the exercises, have a look at the python tutorial, here. Try to repeat some of the code in your ScrumPy window.

(a) Start a ScrumPy session, and make the following assignments:

   1 >>> num_type_1 = 8
   2 >>> num_type_2 = 9.0
   3 >>> str_type = 'a'

What do you think the value and type of the following expressions are?

num_type_1/2
num_type_1/2.0
num_type_2/3
5 + 2*10
str_type*4

Control your answers using Python. Any deviations from expected? Why?

(b) Create a list with six items, e.g.:

   1 >>> test_list= ['a', 'b', 'c', 'd', 'e', 'f']

Make another list (reversed_list) that is a copy of the first, e.g.:

   1 >>> reversed_list=test_list[:]

(why will this assignment result in a copy of test_list? Hint: this is a slicing operation.)

Now, reverse the order of the items in reversed_list by repeating the following assignment with suitable choice of n and m:

   1 >>> reversed_list[n]=test_list[m]

(c) Define two empty dictionaries, one called ara the other rad. Both will be used for collection information about plants, ara about Arabidopsis thaliana, and rad about radish.

For both dictionaries define the keys name with a empty dictionary as value, usedFor and relatedTo (both with empty lists as values). For the nested dictionaries (name) define the keys Latin (with another empty dictionary as value) and English with the English name of the plant as value (for A. thaliana this is thale cress). For the Latin dictionary define the keys genus (with the value being the genus name) and sp. (with the value being the species name). For A. thaliana, genus is Arabidopsis, sp. is thaliana. The Latin name for radish is Raphanus sativus.

For both ara and rad define the key reletedTo, with an empty list as value. Since the two plants are in the same family, for each dictionary append the other dictionary to the relatedTo list.

How would you find, using one single command, the English name of an organism that A. thaliana is related to?

Explore what happens if you try to access the nested relatedTo organism, starting from any of the two organisms (i.e. start from one organism, access the organism it is related to, then continue down one level to the organism the second organism is related to).

(d) The volume of a sphere can be calculated as 4/3*π*r³ given a radius r. Import the math module to obtain an approximation of π and calculate the volume of a sphere with radius of 4. The method pow(...) can be used to calculate x^y. Note that 201.06 is not the correct answer.

Part 2 - Creating Models

Have a look at this tutorial before starting.

(a) Create and load toy_model.spy as described in the tutorial.

(b) See what happens if you remove the tilde sign (~) from any of the reactions in the model file. Implement the modification either by selecting ScrumPy > Compile in the model file, or do:

   1 >>> m.Reload()

in the ScrumPy window (after saving the modification in the model file).

(c) Modify the reaction stoichiometry in any of the reactions in toy_model.spy so that all reactants are identical, e.g.:

   1 A_tx:
   2     A -> A
   3     ~

Reload or compile the model as above. What happens?

Part 3 - The Stoichiometry Matrix and Null-space

(a) Familiarize yourself with the stoichiometry matrix of the toy model, e.g. repeat the examples in the section of the tutorial relating to the stoichiometry matrix.

(b) Generate a kernel of the stoichiometry matrix (using the .NullSpace() method).

(c) Multiply the stoichiometry matrix by the kernel matrix, using the .Mul() method of stoichiometry matrix.

(d) Repeat the multiplication above using the external stoichiometry matrix instead.

(e) Create a kernel of the external stoichiometry matrix - what is different?

(f) Multiply the internal and external stoichiometry matrices by your new kernel.

(g) Try adding or removing reactions from the model, and repeat the exercises above, explain the new results.