# Hydrobates

A mathematician thinks aloud

## Archive for January, 2012

### Modelling Dictyostelium aggregation, yet again

January 29, 2012

In the last post I discussed the second chapter of the book ‘La Vie Oscillatoire’ which is concerned with glycolytic oscillations. The third chapter is on calcium oscillations, a theme which I have written about more than once in recent posts. Here I will say something about the subject of the fourth chapter, signalling by the cellular slime mould Dictyostelium. I have written some things about Dictyostelium in previous posts. It has fascinated many scientists by its ability to go from a state where the cells are independent to a state which looks like a multicellular organism as a reaction to a scarcity of food. The cells gather by means of chemotaxis. It is usual for mathematical talks about chemotaxis to start with nice pictures portraying the life cycle of Dictyostelium discoideum, the most famous organism of this type. It gave rise to the formulation of the Keller-Segel model which is rather popular among mathematicians. As I have mentioned in previous posts it is not so clear to what extent the Keller-Segel model, as attractive as it is, is really relevant to the life of D. discoideum. In his book Goldbeter seems to share this sceptical point of view while leaving open the possibility that the Keller-Segel model might be a reasonable model for aggregation in a less famous relative of D. discoideum, Dictyostelium minutum.

Returning to D. discoideum, it is a known fact that the cells of this organism signal to their neighbours by producing cAMP (cyclic adenosine monophosphate). This process can be modelled by a system of ODE (or a similar system with diffusion) called the Martiel-Goldbeter model. Experimentally it is seen that cultures of starving D. discoideum develop circular or spiral waves centered on certain pacemaker cells. In the book a decription is given of how this process can be understood on the basis of the Martiel-Goldbeter model. It is useful to draw a diagram of the dynamical properties of this model as a function of the activities of two key enzymes. There is a region where solutions of the MG model converge to a stationary solution, a region where they show excitable behaviour and a region where there is a limit cycle. As a cell develops following the beginning of a period of hunger it moves around in this parameter space. It starts from constant production (1), becomes excitable (2), produces pulsations (3) and then becomes excitable again (4). This is a statement which relates to the ODE system but when diffusion is added the idea is that the pacemaker cells are in stage (3) of their development while cells in stage (4) can then lead to waves which propagate away from the pacemakers. Because of the variability of the cells (of their intrinsic properties or their life histories) the different stages of development can be present at the same time.

I now want to say some more about the MG model itself. The variables in the system are the fraction $\rho_T$ of receptor in the active state and the intracellular ($\beta$) and extracellular ($\gamma$) concentrations of cAMP. The evolution equation for $\rho_T$ is of the form $\frac{d\rho_T}{dt}=-f_1(\gamma)\rho_T+f_2(\gamma)(1-\rho_T)$ where $f_1$ and $f_2$ are ratios of linear functions. The evolution equation for $\gamma$ is $\frac{d\gamma}{dt}=\frac{k_t\beta}{h}-k_e\gamma$ where the quantities other then the unknowns are constants. The most complicated evolution equation is that for $\beta$ which is $\frac{d\beta}{dt}=q\sigma\Phi(\rho_T,\gamma,\alpha)-(k_i+k_t)\beta$. Here $\alpha$ is the concentration of ATP, taken to be constant in this model, and all other quantities except the unknowns are constants. The function $\Phi$ is complicated and will not be written down here. It is obtained from a system with more equations by a quasi-steady state assumption and in its dependence on $\gamma$ it is a ratio of two quadratic functions. A further quasi-steady state assumption leads to a simplified system for $\rho_T$ and $\gamma$ alone which is more tractable for analytical considerations.

### Albert Goldbeter and glycolytic oscillations

January 21, 2012

This Christmas, at my own suggestion, I was given the book ‘La Vie Oscillatoire’ by Albert Goldbeter as a present. This book is concerned with oscillatory phenomena in biological systems and how they can be explained and modelled mathematically. After the introduction the second chapter is concerned with glycolytic oscillations. I had a vague acquaintance with this subject but the book has given me a much better picture. The chapter treats both the theoretical and experimental aspects of this subject.

If yeast cells are fed with glucose they convert it into alcohol. Those of us who appreciate alcoholic beverages can be grateful to them for that. In the presence of a supply of glucose with a small constant rate alcohol is produced at a constant rate. When the supply rate is increased something more interesting happens. The output starts to undergo periodic oscillations although the input is constant. It is not that the yeast cells are using some kind of complicated machine to produce these. If the cells are broken down to make yeast extract the effect persists. In fact for yeast extract the oscillations go away again for very high concentrations of glucose, an effect not seen for intact cells. This difference is not important for the basic mechanism of production of oscillations. The breakdown of sugar in living organisms takes place via a process called glycolysis consisting of a sequence of chemical reactions. By replacing the input of glucose by an input of each of the intermediate products it was possible to track down the place where the oscillations are generated. The enzyme responsible is phosphofructokinase (PFK), which converts fructose-6-phosphate into fructose-1,6-bisphosphate while converting ATP to ADP to obtain energy. Now ADP itself increases the activity of PFK, thus giving a positive feedback loop. This is what leads to the oscillations. The process can be modelled by a two-dimensional dynamical system called the Higgins-Selkov oscillator. Let $S$ and $P$ denote the concentrations of substrate and product respectively. The substrate concentration satisfies an equation of the form $\dot S=k_0-k_1SP^2$. The substrate is supplied at a constant rate and used up at a rate which increases with the concentration of the product. (Here we are thinking of ADP as the product and ignoring other possible effects.) The product concentration correspondingly satisfies $\dot P=k_1 SP^2-k_2 P$.

The Higgins-Selkov oscillator gives rise to a limit cycle by means of a Hopf bifurcation. The ODE system is similar to the Brusselator. There are two clear differences. The substance which is being supplied from ouside occurs linearly in the nonlinear term in the Higgins-Selkov system and quadratically in the Brusselator. In the Higgins-Selkov system the nonlinear term occurs with a negative sign in the evolution equation for the substance being supplied from outside while in the Brusselator it occurs with a positive sign. In the book of Goldbeter the Higgins-Selkov oscillator seems to play the role of a basic example to illustrate the nature of biological oscillations.
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 20, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0}, "config": {"markdown_headings": true, "markdown_code": false, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "github"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "english_confidence": 0.9382222294807434, "perplexity_flag": "head"}