Team:Imperial College London/Project Auxin Modelling

From 2011.igem.org




Module 2: Auxin Xpress

Auxin, or Indole 3-acetic acid (IAA), is a plant growth hormone which is produced by several soil bacteria. We have taken the genes encoding the IAA-producing pathway from Pseudomonas savastanoi and expressed them in Escherichia coli. Following chemotaxis towards the roots and uptake by the Phyto Route module, IAA expression will promote root growth with the aim of improving soil stability.




Modelling

Collapse all | Expand all

1. Introduction

Adequate production of the plant hormone auxin or indole-3-acetic acid (IAA) from genetically modified (GM) Escherichia coli could efficiently promote plant root growth. However, IAA is toxic to plants if its concentration is too high. Therefore it is important to be able to predict the IAA expression level for a given promoter, and then adjust the promoter strength to ensure the IAA produced by our E.coli will optimally increase root growth.

IAA increases root growth in terms of root length and number of branches. In order to study how different concentrations of IAA affect root growth patterns, modelling tools were combined with wet lab results to predict and visualise Arabidopsis root growth.

Collapse

2. Modelling of IAA synthesis

2.1 Objective

Determine the IAA expression level of a single E. coli cell with IAA promoter strength 4.536 RNA/min/μg substrate DNA, and then predict the number of bacteria to be placed in the seed coat to induce optimal root growth.


2.2 Description

The genetically engineered IAA pathway involves two genes, IaaM and IaaH, both of which are constitutively expressed. The IaaM gene encodes tryptophan-2-monooxygenase (T-2-monase) that catalyses the conversion of tryptophan (Trp) to indole-3-acetamide (IAM), which is then hydrolysed to release indole-3-acetic acid (IAA) by the hydrolase iaaH[3]. At the same time, the synthesised IAM and IAA will competitively inhibit the enzyme activity of tryptophan-2-monooxygenase, thereby inducing a negative feedback loop on the expression of IAA. The enzymatic reactions involved in the pathway are illustrated in Figure 1.


Figure 1: IAA synthesis pathway. (Diagram by Imperial College London iGEM team).


In addition, based on research carried out by mathematical biology research group from the University of Colorado[4], tryptophan is also negatively controlled inside bacteria. Therefore the tryptophan synthesis pathway should be integrated into the above model. Furthermore, in order to reduce the numbers of parameters, as most of the parameters in Figure 1 are not available to us, the IAA pathway is then simplified into two Michaelis Menten equations, which are combined with the tryptophan pathway and constitutive gene expression for T-2-monase and iaaH. The whole tryptophan IAA pathway model is described in Equation 1[5], and the parameters are defined in the parameters section below.

In this model, we made the following assumptions:

(1) We neglect the short time delay due to synthesis of Trp-T-2-monase (substrate-enzyme (ES) complex), IAM-iaaH (ES complex), IAM-T-2-monase (inhibitor-enzyme (EI) complex) and IAA- T-2-monase (inhibitor-enzyme (EI) complex) and assume that these species reach their equilibrium almost instantaneously.

(2) The degradation rate of IAA in the dark is extremely low compared to bacterial growth rate. Therefore we used bacterial growth rate as the IAA degradation rate in this model.

(3) During the modelling of the IAA pathway, we found that the rate-determining species for IAA synthesis is IAM, not the enzyme iaaH, since the production of IAM is inhibited by itself and IAA.


2.3 Results and discussion

Fig.2 below represents the output of our model.

1. How the concentration of each protein species varies with time.

2. The simulation of the enzymatic reaction for each of the species, with initial concentrations of OF = 1.54×10-4 µM, MF = 3.78×10-4 µM , E = 0.378 µM, Trp = 4.1 µM and all others = 0 [4].

Also, Fig.2 shows that the IAA expression level is 72.35 µM, which means that each bacterium produces 7.24×10-14 µmole per bacteria at steady state with bacterial volume equals to 10-15 dm3. From wet lab experiments we know that the optimal concentration of IAA to promote root growth is 0.1 nM, and the volume of an Arabidopsis seed coat approximately equals 3.6×10-9 m3 [7]. Therefore the number of bacteria required to be present in seed coat to maximally increase root growth of roots is 4.97 x 106.

The details of the calculation are listed below:

1. Number of moles produced by a single bacterium is 72.35 µm * 10-18 m3 (7.235 * 10-17µmol).

2. The number of moles of IAA need to be in one seed coat is 3.6 * 10-9 m3 * 0.1 nmol.

3. The number of bacteria need to be place in one seed coat is (3.6 * 10-10 µmol) / (7.235 * 10-17) = 4.97 * 106.

Figure 2 (a): The evolution of IAM vs. time. Figure 2 (b): The evolution of IAA vs. time.


In conclusion, we obtained the concentration of IAA produced by the IAA DNA construct (72.25 μM). From this value, we calculated the number of bacteria that would need to enter the root under ideal growth conditions and ignoring death and division of bacteria, which was found to be 4.97x106 bacteria for Arabidopsis. This value varies due to the variation in seed size of different of plants. Since bacterial cells can be lost from the plant and not all the bacteria from the seed coat will enter the plant, the number of bacteria required in the seed coat will be higher than this. Furthermore, the modelling of root growth in the next part helped us visualise the root growth.

Collapse

3. IAA uptake and root growth

3.1 Objective

1. Create a graphic program to demonstrate the growing process of the Arabidopsis root system, based on the principles of Lindenmayer system and plant physiology.

2. Use Matlab data fitting tools to develop relations between IAA concentration and growth rate and number of branches.


3.2 Description

We used computational tools to visulise the root growth phenomena (primary rootlength, branching, root density, etc) in different enviromental conditions, and specifically considered root order and root length. Root order describes the branching “generation” of a root system; a root without branching is defined as a zero-order root. Root growth depends on environmental factors, such as gravitation and soil heterogeneities.


3.2.1 Tropisms

A root system starts with a single root tip of a zero-order root. Then the root grows away from the plant stem in a conical way[10].

Root growth depends on environmental factors, such as gravitation and soil heterogeneities. Therefore, two more variables are defined to describe the plant adaptation[10]:

σ:-

- how strong the roots direction changes per 1 cm growth

- a larger value indicates a more deflected root and a more twisted root system

N:-

- the number of trials for the roots to find the optimal angles α and β for the rotation

of the downward movement

- N can be any real number. If N = 1.5, if means that N can be either 1 or 2.


Figure 3: Differences between root systems with different values of N and σ. (Modelling based on Daniel Leitner et al from BOKU, and codes were modified by Imperial College London iGEM team 2011).


3.2.2 Lindenmayer system and root growth modelling

An L-system is a parallel rewriting system, namely a variant of a formal grammar, most famously used to model the growth processes of plant development, but also able to model the morphology of a variety of organisms, due to the two main properties: recursiveness and self-similarity[11]. Plant models and natural-looking organic forms are easy to define, as by increasing the recursion level the form slowly 'grows' and becomes more complex.

Using L-systems for generating graphical images requires that the symbols in the model refer to elements of a drawing on the computer screen. It interprets each constant in an L-system model as a turtle command.

3.2.3 Root growth

The modelling of IAA uptake can give predictive answers to questions such as these:

"What is the primary root growth rate?"

"What does the root system look like after a certain period of time?"

"How does Arabidopsis respond to different IAA concentrations?"

We modified a MATLAB program developed by Daniel Leitner's research group[6] to demonstrate the 3D root system based on the principles of Lindenmayer system (turtle commands) and the root growth modelling toolbox developed by Daniel Leitner et al from BOKU (Universität für Bodenkultur Wien, University of Natural Resources and Life Sciences, Vienna)[12].

3.2.4 Data fitting

We took growth rate parameters from the literature and analysed the raw data from our wet lab to give more accurate and suitable parameters for our project. Arabidopsis plants were planted and root length and number of branches recorded every two days from day 0 to day 9. Root length, daily root growth rate and number of branches were plotted against time and IAA concentration.

3.3 Results and discussion


3.3.1 Visualisation of a root system

The values from the literature give the relationship between external IAA concentration and elongation of the roots as 5*10-5 mol/L → 200 µm elongation in 30 min. The modelling parameter of growth speed is therefore 9.6*10-3 m/day[9].

We observed real root growth patterns and modified our simulation to give a more reliable and accurate prediction of root growth. Arabidopsis has a primary root with zeroth order and it is thicker than the branches. Arabidopsis normally grows to the depth of 20-30 cm inside the soil and branches only once. The 3D picture shown below predicts the root growth with different elongation rate (with IAA = 0.96 cm/day; without IAA = 0.46 cm/day[10]). They can be compared with the photo of a real root system.

Figure 4: 3D visalisation of the arabidopsis root system. The graph Figure 4(a)(b) shows the demonstration of the Arabidopsis root system in different IAA concentration conditions with a growth time of 20 days as our simulation results. Figure 4(c)(d) are the actual photos of real Arabidopsis plants taken from the literature and our wet lab. (Modelling based on Daniel Leitner et al from BOKU, and codes were modified by Imperial College London iGEM team 2011).





Figure 5: The graphic analysis of the 3D root system. By taking cross and longitudinal sections of the root system at different positions, Figure 5(a) is the root length distribution (by fraction) against the root depth. Figure 5(c) is the surface area distribution map. This graph will give guidance for placing GM E. coli inside the soil. E. coli should be placed at the depth with the maximal root surface area (6 cm) inside the soil for the maximal uptake efficiency. (Modelling based on Daniel Leitner et al from BOKU, and codes were modified by Imperial College London iGEM team 2011).

NEW SINCE EUROPE JAMBOREE

3.3.2 Data fitting

We designed an experiment to support the data fitting and tested it before the European Regional Jamboree. The results are integrated into Section 2. The effect of IAA on Arabidopsis roots. Click here to see the details of the first data-fitting experiment.

After the European Regional Jamboree, we improved the protocol to include the following considerations and improvements:-

a) The degradation of IAA:

The wet lab results shows that the IAA concentration in phytogel media is negligible after 9 days, the degradation of IAA is an important factor when root growth is modelled. The IAA concentration can be calculated everyday using the equation: remaining IAA (%) = 95.75 - 2.9x + 0.05x2, where x is the number of days after the auxin has been added.[14]

b) The number of samples:

We prepared 20 replica for each auxin concentration to reduce the error.

c) The measurements technique:

Since the uptake rate and the IAA degradation rate changes rapidly, we measured Arabidopsis roots in phytogel everyday instead of once every three days. Therefore we have correlated the daily growth rate with IAA concentration.

Figure 6: Daily growth rate of Arabidopsis against the logarithms of the additional IAA concentration. The relation can be described by the Gaussian equation.


Then, the mathematic model we got from the data fitting experiments was integrated with our 3D root simulation system. This produced a MATLAB toolbox with IAA concentration as the input and the 3D root model as the output. We defined a valid IAA concentration range from 10-3 mol/L (1 mM) to 10-14 mol/L (0.01 pM) for the simulation based on the following reasons:

a) When the IAA concentration is higher than 10-3 mol/L (1 mM), the Arabidopsis will be killed.

b) When the IAA concentration is lower than 10-14 mol/L (0.01 pM), the Arabidopsis will not be affected since the IAA produced by Arabidopsis itself is much more. Also, it is reasonable to define a limitation of IAA dilution to keep the percentage error low.

c) The Gaussian equation can only describe the relation between the growth rate and the IAA concentration within this range accurately.

Our wetlab results shows two different optimal IAA concentrations for the Arabidopsis root growth, 0.1nM in liquid media and 1pM in phytogel. The data fitting experiment is done inside phytogel.

Figure 7: The results of the simulation with different input.



The video below shows how the MATLAB toolbox predicting the root system after 25 days given a certain IAA concentration.The results simulated are supported by our wetlab experiments.(Please select the 720HD version when you play the video to see the details clearly.)

Collapse

4. Parameters
5. Matlab code
6. References

[1] Hutcheson SW, Kosuge T (1985) Regulation of 3-indole acetic acid production in Pseudomonas syringae pv savastanoi: Purification and properties of tryptophan 2-monooxygenase. The Journal of Biological Chemistry 260(10): 6281-6287.

[2] Saepen S, Vanderleyden J, Remans R (2007) Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiology Reviews 31(4) 425-448.

[3] Cheng Y, Dai C, Zhao Y (2006), Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis. Genes & Dev 1790-199. Doi: 10.1101/gad.1415106

[4] http://mathbio.colorado.edu/mediawiki/index.php/MBW:Modelling_the_Tryptophan_Operon

[5] Santillan M, Mackey MC (2001) Dynamic regulation of the tryptophan operon: a modeling study and comparison with experimental data. PNAS 98(4) 1364-1369.

[6] http://www.ccdb.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi

[7] http://www.seedgenes.org/Tutorial.html

[8] Leitner D, Klepsch S, Bodner G, Schnepf A (2009) A dynamic root system growth model based on L-Systems: tropisms and coupling to nutrient uptake from soil. Plant and Soil 332(1-2): 177-192.

[9] Rayle DL, Evans ML, Hertel L (1970) Action of auxin on cell elongation. Proc Nat Acd Sci USA 65(1): 184-191.

[10] Beemster GTS, Baskin TI (1998) Analysis of cell division and elongation underlying the developmental acceleration of root growth in Arabidopsis thaliana. Plant Physiology 116(4): 1515-1526.

[11] Hu H, Cao X, Lin B (2003) Three dimensional Lindenmayer system. The Journal of Engineering Graphics.

[12] Leitner D, Schnepf A, Klepsch S, Roose T (2009) Comparison of nutrient uptake between 3-dimensional simulation and an averaged root system model. Plant Biosystems 144(2): 443-447.

[13] Brenda: The Comprehensive Enzyme Information System (http://www.brenda-enzymes.org/).

[14] James D, Stephen K, Robert M (1986) The Effect of Salt Concentration on Auxin Stability in Culture Media. Plant Physiology 81: 934-936.

Collapse

M2: Design M2: Assembly