🧬 biology

Migration Model of Lung Cancer Cell Dynamics and Immune System Interplay

This paper presents a partial differential equation-based mathematical model integrating lung cancer cell migration, growth, and tumor-immune interactions to analyze disease dynamics and treatment responses, while also exploring fluid dynamics in obstructed airways and identifying a "Go-or-Grow" mechanism as a key descriptor for specific cell lines.

Original authors: Willie Kudakwashe Chidaushe, Thomas Musora, Steady Mushayabasa

Published 2026-09-11

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Willie Kudakwashe Chidaushe, Thomas Musora, Steady Mushayabasa

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

Lung cancer remains one of the most formidable challenges in modern medicine, a disease where cells in the airways grow out of control, damaging the very organ that keeps us alive. While the causes are well known, ranging from the toxins in cigarette smoke to invisible radioactive gases seeping from the ground, the way these tumors spread and interact with the body's defenses is a complex dance of biology and physics. When cancer cells break away from a primary tumor, they do not simply wander aimlessly; they move through tissues and fluids, often hitching a ride on the body's own circulatory systems to reach distant organs. At the same time, the immune system sends out specialized cells to hunt down these invaders, creating a constant, invisible battle within the lungs. Understanding this struggle requires more than just looking at cells under a microscope; it demands a way to map how these cells move, how heat spreads through tissue, and how drugs travel through the bloodstream.

A team of researchers from Zimbabwe has tackled this complexity by building a detailed mathematical map of lung cancer dynamics. Instead of relying solely on physical experiments, they constructed a set of equations that act like a virtual laboratory, simulating how tumor cells, immune cells, and treatment drugs interact over time and space. Their work focuses on the physical forces at play: how cancer cells migrate, how the immune system responds to infection, and how the flow of blood and air in the lungs affects the spread of the disease. By treating the tumor and the surrounding fluids as a continuous system, they were able to model the movement of particles, the transfer of heat, and the pressure changes that occur when a tumor grows and blocks airways. This approach allows scientists to visualize processes that are too small or too fast to see directly, offering a new lens through which to view the disease.

The researchers began by formulating a model that tracks three main groups: the cancer cells themselves, the immune cells activated to fight them, and the therapeutic particles used in treatment. They incorporated the physical reality of the lungs, accounting for how blood flows, how temperature changes, and how gases like radon from cigarette smoke might influence the environment. In their simulations, they found that the movement of these cells is not random but follows specific physical laws, much like how water flows through a pipe or heat spreads through a metal rod. They discovered that the density of the tumor and the speed of blood flow significantly alter how quickly cancer cells can spread or how effectively immune cells can reach them. For instance, when the tumor becomes very dense, it can restrict the movement of water and other fluids around it, creating a bottleneck that changes how the disease progresses.

A key finding in their work involves the concept of "Go-or-Grow," a behavior where cells must choose between moving to a new location or staying put to multiply. The researchers' model suggests that certain types of lung cancer cell lines fit this pattern, where the cells prioritize migration over growth or vice versa, depending on the conditions. This distinction is crucial because it helps explain why some tumors spread quickly while others remain localized. The team also simulated how heat moves through the lung tissue, noting that tumors often generate more heat than healthy tissue due to their high metabolic activity. They calculated that the heat capacity of the lung tissue can change dramatically in the presence of a tumor, reaching values over 1,000 joules per kelvin, which could serve as a detectable signal for diagnosis.

The study also explored the role of fluid dynamics in drug delivery. By modeling the velocity of blood and the movement of therapeutic particles, the researchers showed that the flow rate of fluids in the lungs directly impacts how much heat builds up around a tumor. Faster fluid flow tends to cool the area, reducing thermal buildup, while slower flow allows heat to accumulate. This relationship is vital for treatments that rely on temperature, such as thermal therapies, where controlling the heat distribution is essential for killing cancer cells without damaging healthy tissue. The simulations revealed that the speed at which immune cells and drugs move is heavily influenced by the physical properties of the lung environment, including the viscosity of the blood and the pressure differences across the airway walls.

To ensure their model was accurate, the researchers compared their mathematical predictions against real-world data, fitting their equations to observed cell counts and tumor sizes. They used a statistical method to determine which version of their model best explained the data, finding that a specific subset of cell lines was indeed best described by the "Go-or-Grow" behavior. This validation step gave them confidence that their virtual simulations reflected real biological processes. They also calculated specific physical values, such as the diffusion coefficient for immune cells, which they found to be extremely high, indicating that these cells move rapidly and spontaneously through the lung tissue. In contrast, the diffusion coefficient for tumor cells in some scenarios appeared negative in their calculations, a mathematical sign that the cells are so densely packed that they restrict the movement of water around them, effectively trapping themselves in a crowded mass.

The researchers emphasized that while their model provides a powerful tool for understanding lung cancer, it is a simulation based on mathematical principles and does not replace clinical judgment. The complexity of the human body means that many variables, such as the exact rate of drug degradation or the specific behavior of immune cells in different patients, are still being refined. However, by bridging the gap between fluid dynamics, heat transfer, and cell biology, this work offers a clearer picture of how lung cancer evolves. It suggests that treating the disease may require not just attacking the cells, but also managing the physical environment they inhabit, from the flow of blood to the temperature of the tissue. As the researchers noted, these models can help clinicians design better treatment plans, predict how a tumor might respond to therapy, and ultimately improve the chances of a successful outcome for patients.

Technical Summary: Migration Model of Lung Cancer Cell Dynamics and Immune System Interplay

Problem Statement
Lung cancer remains a leading cause of global cancer-related mortality, driven by genetic damage often linked to smoking, radon exposure, and air pollution. While the biological mechanisms of tumor migration, invasion, and immune evasion are well-documented, there is a need for rigorous mathematical frameworks that integrate these complex dynamics. Specifically, the interplay between tumor cell migration, immune system activation, fluid dynamics within the lungs, and the effects of therapeutic interventions requires a unified modeling approach. Current literature often addresses these factors in isolation; this study aims to formulate a comprehensive system of Partial Differential Equations (PDEs) to model the spatiotemporal dynamics of lung cancer cells, immune responses, and the physical transport phenomena (heat and fluid flow) associated with tumor progression and treatment.

Methodology
The authors developed a mathematical framework based on a system of coupled Partial Differential Equations (PDEs) incorporating migration, growth, and interaction terms. The model utilizes a spherical coordinate system to represent the lung environment, assuming radial symmetry.

Governing Equations: The core of the study consists of five coupled PDEs (Eq. 2.1–2.5) describing:
- Tumor Cell Concentration ( $C$ ): Modeled with a logistic growth term, a diffusion-advection term, radiative heat transfer effects (Stefan-Boltzmann law), a term representing toxic gas volume (radon from cigarettes), and a clearance term by immune cells.
- Immune Cell Density ( $n$ ): Modeled with proliferation, clearance by tumor cells, radiative heat transfer, and sensitivity to chemokines.
- Therapeutic Particle Concentration ( $H$ ): Includes terms for direct immune boosting by particles, washout, activation/inhibition via PD-1 pathways, and metabolic degradation.
- Blood Velocity ( $v$ ) and Temperature ( $T$ ): Derived from the continuity equation and bio-heat equations, incorporating blood perfusion, pressure gradients, and thermal conduction.
Physical Principles: The model integrates the continuity equation (mass conservation), Starling's principle for fluid movement across semi-permeable capillary walls, and the Stefan-Boltzmann law for radiative heat transfer. It also considers mixed convection flow, accounting for both forced convection (breathing) and natural convection (heat-induced density differences).
Analytical and Numerical Approaches:
- Continuity Solutions: The authors derived general solutions for the continuity equation in spherical coordinates, transforming them into standard convection-diffusion equations.
- Parameter Estimation: Diffusion coefficients ( $\delta_C, \delta_n, \delta_T$ ) and other parameters were estimated using experimental data ranges, logistic growth assumptions, and binary diffusion formulas (Fuller-Schechter-Giddings method for radon and oxygen).
- Numerical Simulation: Implicit finite difference schemes (FDS) were employed to solve the non-linear PDEs. The study utilized Python and MATLAB for graphical analysis and numerical stability.
- Model Selection: The Akaike Information Criterion (AIC) was used to compare model performance, balancing goodness-of-fit (Sum of Squares Error) against model complexity to avoid overfitting.

Key Contributions and Results

Unified Mathematical Framework: The study presents a novel system of PDEs that simultaneously models tumor-immune interactions, fluid dynamics, and thermal regulation in the lungs. This allows for the analysis of how physical factors (temperature, blood flow) influence biological processes (migration, drug delivery).
Diffusion Coefficient Analysis: The authors calculated diffusion coefficients for tumor cells, immune cells, and temperature. Notably, the model identified that a negative diffusion coefficient for tumor cells ( $\delta_C < 0$ ) mathematically corresponds to extremely high tumor cellularity and restricted water movement, rather than a physical impossibility. This "negative diffusion" was shown to lead to catastrophic instability in the model, contrasting with the stable Gaussian spreading observed when $\delta_C > 0$ .
Thermal and Fluid Dynamics: Numerical results demonstrated that increasing fluid flow rates ( $Q_w$ ) inversely scales local thermal buildup ( $T$ ), suggesting that enhanced fluid exchange can limit thermal accumulation in tumor microenvironments. Similarly, shear rates and skin friction coefficients were analyzed, showing their dependence on mixed convection parameters ( $\lambda$ ).
"Go-or-Grow" Classification: The study applied a "Go-or-Grow" hypothesis to classify cell line behaviors. By normalizing migration and proliferation coefficients, the authors established criteria to determine whether a specific cell population is dominated by migration or proliferation.
Wave Speed and Proliferation: Simulations of wave propagation revealed that smaller diffusion coefficients accelerate tumor population decline in specific steady-state scenarios, while immune cell proliferation can exhibit explosive growth (blow-up singularities) depending on threshold parameters.
Data Fitting: The models were fitted to datasets, revealing that a subset of cell lines is best described by the "Go-or-Grow" type, a special case of the formulated model.

Significance and Claims
The paper claims that its primary significance lies in providing a quantitative framework for understanding the complex, multi-scale dynamics of lung cancer. By integrating fluid dynamics and thermal physics with biological tumor-immune interactions, the model offers a tool for:

Predictive Analysis: Simulating tumor growth, therapy response, and the impact of drug delivery under varying physiological conditions.
Therapeutic Optimization: Assisting in the design of thermal therapies and drug delivery strategies by predicting heat capacity and fluid flow interactions within the tumor.
Mechanistic Insight: Offering a deeper understanding of how physical barriers (like airway obstruction) and thermal properties influence metastasis and immune evasion.

The authors emphasize that while the model bridges theoretical pharmacology and patient care, it is intended to support, not replace, clinical judgment. They acknowledge that the complexity of the model presents challenges in parameter estimation due to the scarcity of precise biological data and the dynamic nature of the tumor microenvironment. Consequently, the study advocates for continued development of high-quality experimental data and advanced parameter estimation methods to refine these models further. The work positions mathematical modeling as a critical component in the move toward personalized, data-driven lung cancer management.

More like this