Optimizing defence, counter-defence and counter-counter defence in parasitic and trophic interactions -- A modelling study

Jan Ewald; Stefan Schuster; Sybille D\"uhring; Thomas Dandekar

arxiv: 1907.04820 · v1 · pith:PVR4CBQGnew · submitted 2019-07-10 · 🧬 q-bio.SC

Optimizing defence, counter-defence and counter-counter defence in parasitic and trophic interactions -- A modelling study

Stefan Schuster , Jan Ewald , Thomas Dandekar , Sybille D\"uhring This is my paper

Pith reviewed 2026-05-24 23:17 UTC · model grok-4.3

classification 🧬 q-bio.SC

keywords defence mechanismscounter-counter defenceenzyme inhibitionHaber's rulehost-pathogen interactionsbeta-lactamasesmathematical modellingevolutionary biology

0 comments

The pith

A model of host-pathogen defence shows that counter-counter defence is advantageous only when the inhibitor binds the inactivating enzyme strongly.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models a three-level interaction in which a host produces a toxin, a parasite counters with an enzyme that degrades it, and the host may counter that with an inhibitor of the enzyme. Using an objective function based on the time-integral of toxin concentration according to Haber's rule, the authors derive an analytical threshold for the inhibition constant. Below this threshold, allocating resources to the inhibitor rather than solely to the toxin is the better evolutionary strategy. This result accounts for why counter-counter defences are not universal across biological systems.

Core claim

The central discovery is that the optimal resource allocation between primary defence (toxin production) and counter-counter defence (inhibitor production) includes a positive investment in the inhibitor only when the dissociation constant of the inhibitor-enzyme complex falls below a specific threshold value calculated from the model parameters.

What carries the argument

The threshold value for the inhibition constant (dissociation constant) that determines whether investing in the counter-counter defence improves the objective function over intensifying the primary defence alone.

Load-bearing premise

Fitness depends only on minimizing the time integral of toxin concentration, with no additional costs or constraints on producing the inhibitor or allocating resources.

What would settle it

Finding a natural system where the inhibitor has a dissociation constant above the model's threshold but counter-counter defence is still maintained through evolution would falsify the claim.

read the original abstract

In host-pathogen interactions, often the host (attacked organism) defends itself by some toxic compound and the parasite, in turn, responds by producing an enzyme that inactivates that compound. In some cases, the host can respond by producing an inhibitor of that enzyme, which can be considered as a counter-counter defence. An example is provided by cephalosporins, beta-lactamases and clavulanic acid (an inhibitor of beta-lactamases). Here, we tackle the question under which conditions it pays, during evolution, to establish a counter-counter defence rather than to intensify or widen the defence mechanisms. We establish a mathematical model describing this phenomenon, based on enzyme kinetics for competitive inhibition. We use an objective function based on Haber's rule, which says that the toxic effect is proportional to the time integral of toxin concentration. The optimal allocation of defence and counter-counter defence can be calculated in an analytical way despite the nonlinearity in the underlying differential equation. The calculation provides a threshold value for the dissociation constant of the inhibitor. Only if the inhibition constant is below that threshold, that is, in the case of strong binding of the inhibitor, it pays to have a counter-counter defence. This theoretical prediction accounts for the observation that not for all defence mechanisms, a counter-counter defence exists. Our results should be of interest for computing optimal mixtures of beta-lactam antibiotics and beta-lactamase inhibitors such as sulbactam, as well as for plant-herbivore and other molecular-ecological interactions and to fight antibiotic resistance in general.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives a threshold on inhibitor dissociation constant below which hosts should evolve counter-counter defense instead of stronger primary toxin, using Haber's rule on a competitive-inhibition model.

read the letter

The core result is a threshold value for the inhibitor's dissociation constant: below it, allocating resources to counter-counter defense beats just increasing the primary defense. They model this with standard enzyme kinetics for competitive inhibition and optimize an objective that integrates toxin concentration over time per Haber's rule. The claim is that this explains why layered defenses appear in some systems like beta-lactamase/clavulanic acid but not others. That threshold is new as stated and comes directly from the equations without fitted parameters. The analytic optimum despite the nonlinear ODEs is the part that stands out if the steps hold up. The model is simple and the prediction is testable against real dissociation constants in known interactions. The main limitation is the objective function. It treats fitness as solely the toxin integral with no metabolic or production cost for the inhibitor enzyme and no total protein budget. Adding even a modest linear cost for enzyme synthesis would likely move or remove the threshold. The abstract says the optimum is analytic, but the derivation, boundary conditions, and handling of the time integral are not shown, so the central math needs checking. The model is also a one-shot allocation rather than explicit multi-generation evolution. This is for readers working on molecular arms races, plant-herbivore systems, or antibiotic combinations. It is worth sending to referees because the threshold is a concrete, falsifiable claim even with the modeling simplifications. Ask for the full derivation and a sensitivity check on added costs.

Referee Report

2 major / 1 minor

Summary. The paper develops a mathematical model of host-parasite defense interactions based on competitive enzyme inhibition kinetics. Using an objective function given by the time integral of toxin concentration (Haber's rule), it derives an analytic threshold value of the inhibitor dissociation constant below which it is evolutionarily advantageous for the host to allocate resources to a counter-counter defense (enzyme inhibitor) rather than to intensify primary defense.

Significance. If the analytic threshold result holds after verification, the work supplies a parameter-free, falsifiable condition explaining why counter-counter defenses are not universal across defense systems. This has direct relevance to antibiotic-inhibitor combinations and to molecular ecology. The attempt to obtain a closed-form optimum from nonlinear ODEs is a methodological strength, though the absence of shown steps limits immediate utility.

major comments (2)

Abstract: the central claim that an analytic optimum and threshold are obtained despite nonlinearity is load-bearing, yet the derivation steps, boundary conditions, and explicit handling of the time-integral objective are not provided. This prevents verification of the reported threshold on the inhibition constant.
Model formulation (implicit in abstract and skeptic note): the objective function is defined solely by the time-integral of toxin concentration with no metabolic or production costs for the inhibitor enzyme and no constraint on total protein synthesis. Because this single-scalar fitness measure directly determines the threshold, the paper must show how the condition changes when linear or nonlinear costs are added.

minor comments (1)

The abstract cites cephalosporins, beta-lactamases and clavulanic acid as motivation; a short paragraph mapping model parameters (e.g., dissociation constants) to measured biochemical values would improve the applied section.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will revise the manuscript to improve clarity and robustness.

read point-by-point responses

Referee: [—] Abstract: the central claim that an analytic optimum and threshold are obtained despite nonlinearity is load-bearing, yet the derivation steps, boundary conditions, and explicit handling of the time-integral objective are not provided. This prevents verification of the reported threshold on the inhibition constant.

Authors: We agree that the explicit derivation steps were omitted from the manuscript, which hinders verification. In the revised version we will insert the full analytical derivation, beginning from the system of nonlinear ODEs describing the competitive inhibition kinetics, stating the initial and boundary conditions, and showing how the time-integral objective (Haber's rule) is minimized to obtain the closed-form threshold on the inhibitor dissociation constant. revision: yes
Referee: [—] Model formulation (implicit in abstract and skeptic note): the objective function is defined solely by the time-integral of toxin concentration with no metabolic or production costs for the inhibitor enzyme and no constraint on total protein synthesis. Because this single-scalar fitness measure directly determines the threshold, the paper must show how the condition changes when linear or nonlinear costs are added.

Authors: The present model deliberately omits explicit metabolic costs to isolate the kinetic threshold. We accept that demonstrating robustness to costs is necessary. In revision we will add a supplementary analysis for linear production costs, showing the resulting shift in the threshold value, and will discuss that nonlinear costs generally require numerical methods while preserving the qualitative prediction that sufficiently tight binding favors counter-counter defense. revision: partial

Circularity Check

0 steps flagged

No significant circularity; analytical threshold derived from kinetic ODEs and Haber's-rule objective

full rationale

The paper constructs a system of nonlinear ODEs from competitive enzyme inhibition kinetics and defines an objective function as the time-integral of toxin concentration (Haber's rule). It then derives an analytical threshold on the inhibitor dissociation constant by optimizing resource allocation under that single scalar measure. No fitted parameters are invoked, no self-citations are load-bearing for the threshold result, and the derivation does not reduce to a renaming or redefinition of its inputs. The central claim therefore remains independent of the target prediction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard mass-action enzyme kinetics and the Haber's-rule objective; no new entities are postulated and the threshold is derived rather than fitted.

axioms (2)

domain assumption Toxin effect is proportional to the time integral of toxin concentration (Haber's rule)
Invoked to define the objective function that is optimized.
standard math Competitive inhibition follows standard reversible mass-action kinetics
Used to write the rate equations for the three-species interaction.

pith-pipeline@v0.9.0 · 5827 in / 1391 out tokens · 17097 ms · 2026-05-24T23:17:56.861827+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Only if the inhibition constant is below that threshold... it pays to have a counter-counter defence. ... objective function based on Haber’s rule... time integral of toxin concentration.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

dTox/dt = −Vmax Tox / (KM(1+Inh/KI)+Tox)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.