🧬 biology

Within-host HPV–immune dynamics with epithelial turnover and downstream cervical disease burden: invasion, Hopf bifurcation, and intervention modelling

This paper presents a deterministic within-host model of HPV infection and cervical disease progression that incorporates epithelial turnover and immune dynamics to prove system stability, identify a Hopf bifurcation leading to periodic oscillations, and derive optimal control strategies for therapeutic intervention, while emphasizing that the current benchmark requires further biological calibration before clinical application.

Original authors: Khalid Aldawsari, Fahad Al Basir, Konstantin B. Blyuss

Published 2026-10-20

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Khalid Aldawsari, Fahad Al Basir, Konstantin B. Blyuss

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Inside the human body, the surface of the cervix is not a static wall but a living, breathing landscape that is constantly being renewed. New cells are born at the base, push upward, and eventually shed from the surface in a continuous cycle of turnover. This process is vital for health, but it also creates a complex environment where viruses can hide and evolve. One such virus, the human papillomavirus, or HPV, is incredibly common and is the primary cause of cervical cancer. While the immune system usually clears these infections quickly, in some cases the virus persists, turning healthy cells into precancerous ones that can eventually develop into invasive cancer. The timeline for this disease is tricky: the virus replicates and the immune system reacts over days or weeks, but the slow march from a persistent infection to a full-blown tumor takes years. Understanding how these fast and slow processes interact is essential for figuring out why some infections clear up while others linger and cause harm.

A team of researchers has built a detailed mathematical model to explore this hidden world inside the body. Instead of trying to simulate every single cell, they created a simplified map that tracks the major groups: healthy cells waiting to be infected, cells currently carrying the virus, cells that have turned precancerous, and the immune cells fighting back. A key feature of their model is that it treats the tissue as a crowded space with a fixed limit. As cells multiply, they run out of room, which naturally slows down their growth. This constraint links the virus's ability to spread directly to the health of the tissue itself. The researchers also made a deliberate choice to treat cancer not as a separate, self-sustaining monster, but as a downstream result of the earlier stages. In their view, cancer is simply the accumulation of damage that flows from the precancerous stage, rather than a new system with its own independent rules. This approach allowed them to focus on the core battle between the virus, the tissue, and the immune system without getting lost in the complexities of how a tumor might grow on its own.

When the researchers ran their simulations, they discovered that the system is capable of surprising behavior. Under certain conditions, the levels of infected cells and the immune response do not settle into a steady state. Instead, they begin to rise and fall in a rhythmic pattern, oscillating over long periods. The team found a specific point where the system shifts from a calm state to this wavelike motion. In their specific example, this shift happened when the rate at which the virus is cleared from the body dropped to a certain level. At this tipping point, the system entered a state of sustained oscillation, with the infected cell population swinging between low and high levels over a cycle of roughly 759 days. It is important to note that these numbers come from a mathematical demonstration designed to show that such behavior is possible, not from a direct measurement of a specific patient. The researchers emphasize that while the model proves these oscillations can exist, real-world biology would require precise data to confirm if they happen in actual human infections.

The study also explored how medical interventions could change these outcomes. The researchers modeled three distinct types of treatment: boosting the immune system to fight the virus, removing precancerous lesions through screening, and treating established cancer. They used a mathematical strategy to determine the best way to apply these treatments over time to minimize the disease burden while keeping costs low. Their results suggest that the timing and intensity of treatment matter greatly. For instance, the model indicates that aggressive treatment for cancer is most effective when applied early in the intervention period, while immune stimulation and screening might be used at lower intensities that taper off as time goes on. When they compared a scenario with these optimal treatments against one with no treatment at all, the difference was stark. The treated scenario showed a rapid decline in cancer levels and a much lower overall cost, whereas the untreated scenario saw the disease burden remain high.

However, the researchers are careful to state that this is a theoretical framework, not a clinical guide. The numbers used in their simulations were chosen to illustrate the mathematical possibilities, not to represent a calibrated reality of human biology. They explicitly note that their model does not yet include the specific details needed to predict outcomes for real patients, such as exact rates of cell division or precise immune responses. The value of their work lies in proving that the interaction between tissue renewal, viral load, and immunity can create complex, oscillating patterns and that different interventions can be mathematically optimized. By separating the fast-moving viral dynamics from the slow accumulation of cancer, they have provided a clearer lens through which to view the long, winding path of HPV disease. This clarity offers a foundation for future studies that can eventually be tested against real-world data, potentially leading to better strategies for preventing and treating cervical cancer.

Technical Summary: Within-host HPV–immune dynamics with epithelial turnover and downstream cervical disease burden

Problem Statement
Human papillomavirus (HPV) infection is the necessary causal factor for the majority of cervical cancers. While most infections are transient, a subset persists, leading to precancerous lesions and, eventually, invasive cancer. These processes operate on markedly different time scales: viral replication and immune responses occur over days to months, while lesion establishment and progression to invasive disease occur over years. Existing mathematical models often struggle to balance biological detail with analytical tractability. Specifically, representing cancer as an autonomous tumor-growth subsystem can introduce unsupported assumptions regarding malignant proliferation and carrying capacities, potentially obscuring the dynamics of infection establishment and persistence. Furthermore, there is a need to distinguish between the rapid processes of viral infection and the slow development of cervical disease within a unified framework that incorporates epithelial turnover, immune response, and intervention strategies.

Methodology
The authors formulate a deterministic within-host model consisting of six state variables: susceptible epithelial cells ( $E_S$ ), infected epithelial cells ( $E_I$ ), effective viral load ( $V$ ), precancerous epithelial fractions ( $P$ ), downstream cervical cancer burden ( $C$ ), and local cellular immune activity ( $Z$ ).

Model Structure: The model incorporates a shared tissue-capacity constraint ( $H = 1 - E_S - E_I - P$ ) to regulate epithelial proliferation, coupling cellular turnover directly to infection dynamics. Crucially, cancer ( $C$ ) is modeled not as an autonomous tumor-growth subsystem but as a downstream progression/removal compartment. It receives input from precancerous lesions but does not feed back into the epithelial–viral–immune subsystem.
Intervention Controls: Three control mechanisms are integrated: therapeutic immune stimulation ( $u_1$ ), screen-linked treatment of precancerous lesions ( $u_2$ ), and oncologic therapy for cancer ( $u_3$ ).
Mathematical Analysis:
- Qualitative Analysis: The authors prove the positivity and uniform boundedness of solutions. They derive the infection-free, infected-boundary, and positive coexistence equilibria.
- Invasion Threshold: A basic reproduction number ( $R_{HPV}$ ) is derived, incorporating both viral-load-mediated infection and infected-cell proliferation.
- Stability and Bifurcation: Local stability is analyzed using the Routh–Hurwitz criteria on the characteristic polynomial of the five-dimensional core subsystem (excluding the cancer variable). Numerical continuation is employed to investigate Hopf bifurcations using the viral clearance rate ( $\mu_V$ ) as the bifurcation parameter.
- Optimal Control: An optimal control problem is formulated to minimize precancerous lesions, cancer burden, and intervention costs. Necessary conditions for optimality are derived using Pontryagin's Maximum Principle, resulting in a two-point boundary value problem solved via a forward-backward sweep algorithm.
- Sensitivity Analysis: Normalized sensitivity indices are calculated for the invasion threshold to identify influential parameters.

Key Results

Dynamics and Stability: The model demonstrates that the stability of the system is governed by the five-dimensional epithelial–viral–immune–lesion subsystem. The cancer direction is automatically stable provided the cancer removal rate is positive.
Hopf Bifurcation: For a specific illustrative benchmark (not calibrated to clinical data), numerical continuation identifies a supercritical Hopf bifurcation at $\mu_V \approx 0.08722$ . At this point, the system transitions from a stable equilibrium to a stable periodic attractor. Long-term integrations (up to $10^5$ days) confirm sustained oscillations with a period of approximately 759 days, though convergence is exceptionally slow due to weak instability near the bifurcation point.
Invasion Threshold: The derived threshold $R_{HPV}$ explicitly accounts for two mechanisms: the proliferation of infected cells and the infection of susceptible cells via viral load. Sensitivity analysis reveals that the transmission rate ( $\beta$ ) and viral production rate ( $\pi_V$ ) have the largest positive influence on invasion, while viral clearance ( $\mu_V$ ) and immune killing ( $a$ ) have the strongest negative influence.
Optimal Control Simulations: Numerical simulations of the optimal control problem indicate that combined interventions (immune stimulation, screening-linked treatment, and oncologic therapy) significantly reduce the burden of infection, precancerous lesions, and cancer compared to uncontrolled scenarios. The optimal strategy involves varying intensities: oncologic treatment is applied at maximum intensity initially and tapers off, while immune stimulation and screening treatments start at low intensity and decrease over time.

Significance and Claims
The paper claims to provide a unified mathematical framework that separates the mechanisms of HPV establishment and persistence from the downstream accumulation of cancer burden. By avoiding autonomous tumor-growth assumptions, the model allows for a clearer investigation of invasion thresholds and oscillatory instabilities driven by epithelial turnover and immune interactions.

The authors emphasize that the results, particularly the Hopf bifurcation and the specific parameter values used, serve as a mathematical demonstration rather than a quantitatively calibrated representation of human HPV infection. The benchmark parameters are "search-selected" to demonstrate a mathematically admissible regime with substantial tissue occupancy and are not derived from clinical data or the PAPCLEAR study (though the study's time scales motivated the model structure).

The study concludes that while the model offers insights into the potential for sustained temporal variability (oscillations) in HPV dynamics and the theoretical benefits of multi-stage interventions, quantitative biological calibration—particularly of epithelial proliferation and immune/viral parameters—is necessary before any clinical interpretation can be made. The sensitivity analysis serves to identify parameters that are mathematically influential and candidates for future calibration, rather than providing direct evidence for specific therapeutic interventions in vivo.

More like this