21 tumor sizes Flashcards
Here, we shall focus on two types of
non-dangerous tumours that are found in human bodies:
Tumours consisting entirely of living, but abnormal, cells, see
- Tumours with a “necrotic core”. In the core, cells are dead because the concentration of
nutrients is too low to sustain them
tumors dont produce oxygen does not transport oxygen, but rely on that of the healthy tissues which do transport it and then oxygen permeates the tissues by diffusion
(distance means less oxygen cells die)
instead of modelling the tumour cells, we model the..
model variables
assumptions
instead of modelling the tumour cells, we model the
nutrient concentration c(r, t) inside the tumor as a function of the distance from the core and time
r_2 larger radius includes living cells
r_1 radius of necrotic core
nutrient concentration in healthy tissue: c_2 (outside the tumor) and is CONSTANT
nutrients diffuse inside the tumour (avascular tumour: no production/transport inside the tumour)
tumour living cells (red) assumed consume at constant rate k > 0
minimal nutrient concentration to sustain cells: c₁ < c₂
if concentration goes below this the cells die
necrotic cells (grey) do not consume nutrients
Tumor models:
Hence, the diffusion(-like) equation obeyed by the concentration c(r) is…
∂c(r, t)/∂t =
{D∇²c(r, t) in zones with no living cells
{D∇²c(r, t) − k in zones with living cells .
(no consumption)
Tumours with no necrotic core: Modelling
Under what conditions can the tumour be formed by all living cells?
BCs?
∂C(r,t)/∂t =
diffusion- consumption
here as no necrotic core r∈[0,r₂]
(diffusion modelled with laplacian where D is the diffusion coefficient inside the tumor; consumption assumed constant)
∂C(r,t)/∂t =
** D∇²c(r,t) **- K
PDE 2nd order we need BCs:
∂c(0,t)/∂r = 0 (no nutrients exchange at the core, no flux)
c(r₂,t)=c₂ concentration in healthy tissue outside is the same
Tumor model:
∂C(r,t)/∂t =
** D∇²c(r,t) **- K
PDE 2nd order we need BCs:
∂c(0,t)/∂r = 0
c(r₂,t)=c₂
meaning
BCs
∂c(0,t)/∂r = 0 (no nutrients exchange at the core, rate of change is 0 and nutrients are constant/steady
c(r₂,t)=c₂
concentration in healthy tissue at boundary of tumor is the same as the nutrient concentration in the healthy tissue
oxygen concentration diffuses inside the tumor
Tumor model:
∂C(r,t)/∂t =
** D∇²c(r,t) **- K
PDE 2nd order we need BCs:
∂c(0,t)/∂r = 0
c(r₂,t)=c₂
meaning
Model terms
oxygen concentration diffuses inside the tumor
similar to diffusion eq
∂C(r,t)/∂t =
diffusion- consumption
is a reaction diffusion eq
diffusion from the concentration higher outside
consumption process as the cells are living
Tumours with no necrotic core: Modelling
spherical?
Tumors assumed to be spherical so spacial dependence of c(r,t) varies only with distance r from the centre with laplacian
∇²c(r,t) = (1/r²)(∂/∂r)(r² ∂c(r,t)/∂r))
Tumours with no necrotic core: stationary solution assumptions
We assume that the tumour is stationary
(i.e., stable in time, not expanding):
Under what conditions can we have a stationary tumour formed by living cells only?
stationary (time-independent) concentration c(r) = lim_t→∞ c(r, t)
Even if tumor growing in the body, the time scale at which this growth happens is much larger than time scale of diffusion of oxygen. so its slow enough to ignore at times.
condition
assumption that the tumor is stationary gives
∂c(r,t)/∂t=0
so it becomes: that c(r,t)=c(r) as it does not dep on time
concentration only depends on r
0= D∇²c(r) -k
spherical coords
∇²c(r) = (1/r²)(d/dr)(r² dc(r,t)/dr))
D(1/r²)(d/dr)(r² dc(r)/dr)) - k = 0
time not involved!
with bcs
c’(0)=0 c(r_2)=c_2
gives a solution…
Tumours with no necrotic core: stationary solution
stationary (time-independent) D(1/r²)(d/dr)(r² dc(r)/dr)) - k = 0
time not involved!
with bcs
c’(0)=0 c(r_2)=c_2
gives a solution…
valid for?
c(r) = c₂ - k/6D
quadratic function
minimal cocnentration at centre r=0
c(0)=cₘᵢₙ = c₂ - r₂²[k/6D]
graph against r
parabola increases
bounded above by c₂ at r₂
intercepts y-axis at r=0cₘᵢₙ
only valid for r in[0,r₂]
as expected from diffusion: you expect to decrease to the core inversely proportional to r².
bounded below by c₁? Yes bc we are considering healthy tissue in the tumor, all concentration inside must be above c₁.
Ie we require cₘᵢₙ≥c₁
c₂ - r₂²[k/6D]≥c₁
[6D/k] (c₂ -c₁) ≥ r₂² =R²𝒸ᵣᵢₜ
sqrt of this gives us the critical values
R²𝒸ᵣᵢₜ
tumor model:
We want healthy tissue, using our solution under stationary assumptions solution is
c(r) = c₂ - k/6D
bounded below by c₁? Yes bc we are considering healthy tissue in the tumor, all concentration inside must be above c₁.
Ie we require cₘᵢₙ≥c₁
c₂ - r₂²[k/6D]≥c₁
[6D/k] (c₂ -c₁) ≥ r₂² =R²𝒸ᵣᵢₜ
sqrt of this gives us the critical values for which the healthy tumor is still healthy (viable)
Maximal viable size
of a tumor relates to R𝒸ᵣᵢₜ
Which gives us the maximal size for a tumour composed of living cells only, given the nutrient concentration parameters:
c1, c2, D and the uptake k
What happens if the tumor r_2 is larger than R𝒸ᵣᵢₜ
we will have a tumor with a necrotic core.
Then diffusion of nutrients is not sufficient to guarantee enough nutrients at the core
Tumours with necrotic core
model: when is it true?
model
tumor r_2 is larger than R𝒸ᵣᵢₜ
then we have
∂c(r, t)/∂t =
{D∇²c(r, t) − k if r ∈ [r1, r2] (living cells)
{D∇²c(r, t) if r ∈ [0, r1) (necrotic core)
notice when modelling for the outer healthy tissue part we have the original eq, necrotic core is just diffusion no consumption
Can we say anything about the size r1 of the necrotic core for a
‘large’ tumour (r_2 > R𝒸ᵣᵢₜ)?
.
Tumours with necrotic core
model:
BCs
stationary solution
stationary so assumed….
c(r₂) = c₂ (concentration in healthy tissue outiside)
c(r₁) = c₁ (minimum concentration at the necrotic boundary so need this at the start of the red region assumed to be = as starts dying as soon as youre less than r₁)
c′(r₁) = 0 (no nutrients exchange with the necrotic core)
stationary sol: no time dep so set time derivative to 0
c(r, t) ≡ c(r)
0 = ∂c/∂t = D∇²c(r, t) − k, r ∈ [r₁, r₂]
(as working in the red healthy region only)