17.4 Biological Oscillators: Negative Feedback as the Driving Force

17.4 Biological Oscillators: Negative Feedback as the Driving Force

The Rayleigh oscillator demonstrated that stable oscillations require two essential ingredients: a mechanism that generates periodic motion and a mechanism that stabilizes its amplitude.

Biological systems obviously do not contain springs or air flows. Instead, they rely on regulatory networks composed of genes, proteins, metabolites, and signaling molecules.

This raises an important question:

Which biological mechanisms play the role of the nonlinear friction in the Rayleigh oscillator?

The answer is surprisingly simple.

Stable biological oscillations are almost always generated by negative feedback loops.

However, not every negative feedback loop oscillates. Most biological feedback systems produce stable homeostasis. Understanding why some remain stationary while others oscillate is one of the central questions of systems biology.

In this section, we investigate one of the best-known examples of biological oscillations: the hormonal regulation of the hypothalamic–pituitary–gonadal (HPG) axis.

The hypothalamic–pituitary–gonadal axis

The HPG axis controls reproductive physiology in vertebrates through the interaction of three endocrine organs:

• the hypothalamus,
• the pituitary gland, and
• the gonads.

The hypothalamus secretes gonadotropin-releasing hormone (GnRH), which stimulates the pituitary gland.

The pituitary releases luteinizing hormone (LH) and follicle-stimulating hormone (FSH), which stimulate the gonads.

Finally, the gonads produce sex hormones such as testosterone, estradiol, or progesterone.

These hormones inhibit hormone release from the hypothalamus, thereby closing the regulatory loop.

The overall architecture therefore consists of two positive regulatory interactions followed by one negative interaction, forming a classical negative feedback loop.

Building the model

As in previous chapters, we begin by constructing the simplest model that captures the essential biological interactions.

To avoid unnecessary biological complexity, we represent each endocrine organ by a single variable:

• (H): hormone concentration released by the hypothalamus,
• (P): hormone concentration released by the pituitary,
• (G): hormone concentration released by the gonads.

The stimulatory interactions are assumed to be proportional to the concentration of the upstream hormone.

Each hormone is also degraded continuously.

The pituitary and gonadal dynamics can therefore be written as

$$P' = H-k_1P$$

and

$$G' = P-k_2G.$$

These equations should already look familiar.

They have exactly the same structure as many of the dynamical models developed throughout this book: production increases one variable, while degradation removes it.

The interesting part of the model is the negative feedback acting on the hypothalamus.

Two different kinds of negative feedback

At first sight, it might seem natural to model the inhibition simply by subtracting the gonadal hormone concentration.

However, this would not accurately represent the biology.

It is therefore useful to distinguish between two fundamentally different kinds of negative feedback.

The first type is concentration feedback.

In this case, one molecular species directly removes or degrades another.

Examples include enzymatic degradation, molecular sequestration, or protein degradation.

Here, the feedback acts directly on the concentration of a molecule.

The second type is rate-change feedback.

In endocrine regulation, the gonadal hormones do not directly remove GnRH molecules from the bloodstream.

Instead, they reduce the rate at which the hypothalamus produces new hormone.

The feedback therefore changes the production rate rather than the concentration itself.

This distinction is subtle but extremely important.

One type of feedback modifies a state variable directly, whereas the other modifies the differential equation governing its dynamics.

Recognizing this difference is an important modelling skill because many biological regulatory systems operate through changes in production rates rather than direct removal of molecules.

Choosing an appropriate mathematical function

Once again, we encounter a recurring modelling problem.

We understand the biological mechanism qualitatively, but we still need to translate it into mathematics.

How should the inhibition depend on the hormone concentration?

Rather than searching for an equation immediately, we first ask what biological properties the function should possess.

The function should

• decrease continuously as the inhibitor concentration increases,
• never become negative,
• approach saturation for very large concentrations,
• and allow us to control how sharply the response changes.

Notice that we are not choosing a function because it is mathematically convenient.

Instead, we are selecting a function whose qualitative behavior matches our biological knowledge.

A decreasing sigmoidal function naturally satisfies all of these requirements.

This illustrates an important general principle that will appear repeatedly throughout this book:

Mathematical functions are chosen because they encode biological assumptions.

Different biological assumptions require different mathematical functions.

The mathematics is therefore not arbitrary—it is a formal language for expressing biological knowledge.

Sensitivity as a model parameter

One particularly important property of the sigmoidal feedback function is its steepness.

A shallow curve produces only a weak response to changes in hormone concentration.

The controller reacts gradually.

A steep curve behaves very differently.

Even a small change in hormone concentration produces a large change in regulatory activity.

In other words, the controller becomes highly sensitive.

Mathematically, this sensitivity is controlled by a single parameter, often denoted by the exponent (n).

Rather than representing a specific biochemical quantity, (n) summarizes how strongly the regulatory system amplifies deviations from its operating point.

As we will see in the next section, this seemingly simple parameter determines whether the entire system converges to a stable equilibrium or develops sustained oscillations.

Key concepts

• Stable biological oscillators are typically based on negative feedback loops.
• The HPG axis is a classical example of an oscillatory endocrine feedback system.
• Mathematical models begin with qualitative biological assumptions.
• Negative feedback can act on concentrations or on production rates.
• Biological assumptions determine the mathematical functions used in a model.
• Feedback sensitivity is controlled by the steepness of the regulatory function.

Summary

The HPG axis illustrates how biological oscillators can be modeled using negative feedback. Constructing such a model begins by identifying the underlying biological interactions and translating them into mathematical relationships. An important modelling principle is that mathematical functions are selected because they encode specific biological assumptions. In endocrine regulation, negative feedback acts primarily by changing production rates rather than directly removing molecules. The sensitivity of this feedback will prove to be a key determinant of oscillatory behavior.

Self-check questions

1. Why is the HPG axis a negative feedback loop?
2. Why is each endocrine organ represented by a single variable in the model?
3. What is the difference between concentration feedback and rate-change feedback?
4. Why is a decreasing sigmoidal function an appropriate model for hormonal inhibition?
5. Why should mathematical functions be chosen based on biological properties rather than mathematical convenience?
6. What biological property is represented by the steepness of the feedback function?