IGEM:IMPERIAL/2006/project/Oscillator/Theoretical Analyses/Results/2D Model3a

2D Model 3: Bounded Predator and Prey Growth

Introduction

The analysis of the model with bounded growth terms for both preys and predators has hinted that in order to get oscillations it was more important to balance the respective growth and decay terms of a population than to balance the growth terms within the whole system (the idea being to achieve some form of control of the ratios of the populations). We now test the validity of this interpretation by bounding both growth and decay of the preys while leaving the terms relating to the predators identical to those in the Lotka-Volterra model.

ODE (Ordinary Differential Equation)

Physical interpretation of the equations

• All the arguments that bounded the growth of the preys still apply (production is limited by the number of promoters).
• The decay term of the preys is due to an enzymatic reaction, typically modelled with a Michaelis-Menten kinetics.
• For all the predator terms of the system (second differential equation) we make the same assumption as with Lotka-Volterra that is:

- The predator’s growth is the direct result of interacting the preys. In the present case it is proportional to the population of both preys and predators.

- The population of the predators is controlled by its natural death (proportional to the population)

Basic Results on the Steady Points

The dynamic analysis of the system yields results that differ from the previous system in a significant way:

• The system has two stationary points
• The first stationary point (0, 0) is a saddle point
• The second stationary point (coordinates) is not always stable
• The nature of the second steady point depends on the respective values of Ao and Bo:

- if Ao>Bo, the point is stable (it attracts all the trajectories in its vicinity)

- If Ao=Bo, the point is a centre (the trajectories in its vicinity spin around it)

- if Ao<Bo, the point is unstable (it attracts all the trajectories in its vicinity)

Behaviour at Infinity

The trajectories remain bounded as time goes to infinity regardless of the choice of parameters.

• We can therefore apply Poincare-Bendixson's theorem.
• Due to the saddle point nature of the origin (0,0)

- any trajectory starting on the Y-axis will be attracted to the origin

- any trajectory starting on the X-axis will be repelled to infinity

• For initial conditions (Xo,Yo) both non-null the trajectories will converge

- to the stable steady point if Ao>Bo

- else (Ao<=Bo) to a Limit Cycle: We have Oscillations

Typical Simulations

First Case: Stable System

• The first typical behaviour of the system (convergence to stable steady point) can be illustrated by the following simulation:

[[Image:]]

• The corresponding time variations of the prey and predator populations are dampened oscillations:

[[Image:]]

Second Case: Oscillations

• The second typical behaviour of the system (Oscillations around the steady point) can be illustrated by the following simulation:

[[Image:]]

• The corresponding time variations of the prey and predator populations clearly show the oscillations:

[[Image:]]

Conclusion

The bounding of the predator growth had a dramatic effect on the system and broke the balance between the growth terms of the preys and predators and their decay/degradations terms. Our attempt to re-balance the system by bounding the growth of the predators has failed to yield any oscillation. This suggests that a better way to re-balance the system is to balance the growth and decay terms of the preys and/or predators instead.

Appendices

Related documents

• Stability analysis

File:2d model 3y.PNG

File:2d model 3z.PNG