SlideShare a Scribd company logo
2
Most read
3
Most read
4
Most read
Tumor Growth Model         Method of Characteristics   Solutions for various growth rates   Future Work




                       Dynamic Tumor Growth Modelling

                                Prof. Thomas Witelski,
                     Matthew Tanzy, Ben Owens, Oleksiy Varfolomiyev


                                            NJIT, June 2011
Tumor Growth Model     Method of Characteristics   Solutions for various growth rates   Future Work




Cancer background
      Current medical understanding of some types of cancer




         1   Disease starts from a primary tumor which grows in one
             location
         2   The primary tumor may shed cancer cells which get carried to
             other parts of the body by the circulatory system
Tumor Growth Model       Method of Characteristics        Solutions for various growth rates   Future Work




      Colony size distribution of metastatic tumors with cell number x at
      time t is governed by von Foerster Equation
                                    ∂ρ   ∂
                                       +    (g (x) ρ) = 0                                      (1)
                                    ∂t   ∂x
      g (x) is tumor growth rate determined by
                                               dx
                                                  = g (x)                                      (2)
                                               dt
      Initially no metastatic tumor exists

                                               ρ(x, 0) = 0                                     (3)

      Intitial tumor birth rate
                                               ∞
                     g (1)ρ(1, t) =                β(x)ρ(x, t)dx + β(xp (t)),                  (4)
                                           1

      β(x) tumor birth rate
Tumor Growth Model   Method of Characteristics     Solutions for various growth rates   Future Work




Linear Growth Function: g (x) = k − Ex


      Assuming that in the beginning the main contribution to the tumor
      growth is given by primary tumor (i.e. neglecting the integral term
      in the BC) we have

                         ∂ρ(x, t) ∂ (g (x)ρ(x, t))
                                 +                 =0                                   (5)
                           ∂t            ∂x
                                          ρ(x, 0) = 0                                   (6)
                     g (1)ρ(1, t) = β(xp ),          β(x) = mx α                        (7)
                                  dx
                                     = k − Ex = g (x)                                   (8)
                                  dt
Tumor Growth Model   Method of Characteristics    Solutions for various growth rates    Future Work




Solution by the Method of Characteristics


                                    dρ
                                        = Eρ                                            (9)
                                    dt
      along the characteristics given by
                                        dx
                                           = k − Ex                                    (10)
                                        dt
      Therefore
                                ρ(x, t) = ρ0 (x(0)) e Et                               (11)
      along the characteristics
                                                 K
                       x(t) = x(0)e −Et +          1 − e −Et                           (12)
                                                 E
Tumor Growth Model    Method of Characteristics     Solutions for various growth rates   Future Work




Solution by the Method of Characteristics


      First we use the BC to solve for ρ0 (x(0)), i.e.
                                                                                                α
                             K                                                 K
      (K − E ) ρ0 e Et +       1 − e Et           e Et = m e −Et +               1 − e −Et
                             E                                                 E

      Finally, the colony size distribution of metastatic tumors with cell
      number x at time t is
                                                                          α
                                         mE −α         e −Et (E −k)2
                        ρ (x, t) =       k−Ex     k+      Ex−k
Tumor Growth Model      Method of Characteristics   Solutions for various growth rates   Future Work




Growth function: g (x) = Ex




                     Figure: ρ(x) at t = 10 with all parameters at 1
Tumor Growth Model      Method of Characteristics   Solutions for various growth rates   Future Work




Growth function: g (x) = k − Ex




               Figure: ρ(x) at t = 10 with E = 0.2, k = 0.3, and m = 0.1
Tumor Growth Model    Method of Characteristics   Solutions for various growth rates   Future Work




Single primary tumor size growth (no drug is active)
      Growth function: g (x) = (k − E (t)e −rt )x
      Here no drug is active (i.e. g (x) = kx)




      Figure: ρ(x) at t = 3650 with k = 0.006, α = 0.663, and m = 5.3 × 10−8
Tumor Growth Model     Method of Characteristics   Solutions for various growth rates   Future Work




Single primary tumor size growth (drug effect)
      At time t1 drug is activated
                                                   −r (t−t1 )
      In a region of constant E : x (t) = Ae kt+Ee            /r




      Figure: ρ(x) at t = 3650 with k = 0.006, α = 0.663, E = 0.0083,
      r = 1 × 10−5 and m = 5.3 × 10−8
Tumor Growth Model        Method of Characteristics       Solutions for various growth rates           Future Work




Linear Separable Growth Rate Solution



      Now we have Linear Separable Growth Rate
                                         dx
                                            = [k − E (t)] x
                                         dt
      Colony size distribution

                                                      m                m            t
                               m              −α−1− k−E (t)       α+ k−E (t)   (        k−E (s)ds )t
              ρ(x, t) =      k−E (t)      x                   e                    0
Tumor Growth Model   Method of Characteristics          Solutions for various growth rates   Future Work




Logistic Tumor Growth Rate Solution



      Now we have Tumor Growth Rate given by Logistic Equation
                                     dx
                                        = (k − Ex) x
                                     dt
      Resulting colony size distribution

                                                        e −kt
                               ρ(x, t) =                            2
                                                 (E x0 (e kt −1)+k )
Tumor Growth Model      Method of Characteristics   Solutions for various growth rates   Future Work




Future Work



             Fitting data to models, parameter estimation for multiple
             tumor and multiple patient data
             Incorporating limitations on measurable clinical data
             Comparing PDE, discrete population and polymerization
             models
             Solution of inverse problems for the birth rate
             Probability and statistics of different forms of clinical data

More Related Content

What's hot (20)

PDF
New Mathematical Tools for the Financial Sector
SSA KPI
 
PDF
A type system for the vectorial aspects of the linear-algebraic lambda-calculus
Alejandro Díaz-Caro
 
PDF
Lesson 29: Integration by Substition (worksheet solutions)
Matthew Leingang
 
PDF
Lesson 29: Integration by Substition
Matthew Leingang
 
PDF
YSC 2013
Adrien Ickowicz
 
PDF
Runtime Analysis of Population-based Evolutionary Algorithms
PK Lehre
 
PDF
Andreas Eberle
BigMC
 
PDF
Lecture on nk [compatibility mode]
NBER
 
PPTX
Systems Of Differential Equations
JDagenais
 
PDF
Speeding up the Gillespie algorithm
Colin Gillespie
 
PDF
Considerate Approaches to ABC Model Selection
Michael Stumpf
 
KEY
Trig identities
Simon Borgert
 
PDF
Prml
syou6162
 
KEY
Unit Circle - Trigonometry
Simon Borgert
 
PDF
05 history of cv a machine learning (theory) perspective on computer vision
zukun
 
PDF
Optimalpolicyhandout
NBER
 
PDF
Senior Seminar: Systems of Differential Equations
JDagenais
 
PDF
The tau-leap method for simulating stochastic kinetic models
Colin Gillespie
 
PDF
Nested sampling
Christian Robert
 
PDF
Monte Carlo Statistical Methods
Christian Robert
 
New Mathematical Tools for the Financial Sector
SSA KPI
 
A type system for the vectorial aspects of the linear-algebraic lambda-calculus
Alejandro Díaz-Caro
 
Lesson 29: Integration by Substition (worksheet solutions)
Matthew Leingang
 
Lesson 29: Integration by Substition
Matthew Leingang
 
YSC 2013
Adrien Ickowicz
 
Runtime Analysis of Population-based Evolutionary Algorithms
PK Lehre
 
Andreas Eberle
BigMC
 
Lecture on nk [compatibility mode]
NBER
 
Systems Of Differential Equations
JDagenais
 
Speeding up the Gillespie algorithm
Colin Gillespie
 
Considerate Approaches to ABC Model Selection
Michael Stumpf
 
Trig identities
Simon Borgert
 
Prml
syou6162
 
Unit Circle - Trigonometry
Simon Borgert
 
05 history of cv a machine learning (theory) perspective on computer vision
zukun
 
Optimalpolicyhandout
NBER
 
Senior Seminar: Systems of Differential Equations
JDagenais
 
The tau-leap method for simulating stochastic kinetic models
Colin Gillespie
 
Nested sampling
Christian Robert
 
Monte Carlo Statistical Methods
Christian Robert
 

Similar to Dynamic Tumor Growth Modelling (10)

PDF
journal of mathematics
rikaseorika
 
PPTX
PMED Undergraduate Workshop - Modeling and Estimating Biological Heterogeneit...
The Statistical and Applied Mathematical Sciences Institute
 
PPT
American Statistical Association October 23 Minneapolis Presentation Part 2
Double Check ĆŐNSULTING
 
PPTX
ABC1 - L. Norton - Mathematical model metastases
European School of Oncology
 
PDF
Looking Inside Mechanistic Models of Carcinogenesis
Sascha Zöllner
 
PPTX
Μοντέλα διάχυσης καρκινικών όγκων εγκεφάλου
Manolis Vavalis
 
PPTX
Mathematical Models of Tumor Invasion
Sean Davis
 
PDF
2D CAT Based Modeling of Tumour Growth and Drug Transport
Editor IJMTER
 
PPTX
Principles of chemotherapy
Sheetal R Kashid
 
PDF
ThinkBayes: chapter 13  simulation
Jungkyu Lee
 
journal of mathematics
rikaseorika
 
PMED Undergraduate Workshop - Modeling and Estimating Biological Heterogeneit...
The Statistical and Applied Mathematical Sciences Institute
 
American Statistical Association October 23 Minneapolis Presentation Part 2
Double Check ĆŐNSULTING
 
ABC1 - L. Norton - Mathematical model metastases
European School of Oncology
 
Looking Inside Mechanistic Models of Carcinogenesis
Sascha Zöllner
 
Μοντέλα διάχυσης καρκινικών όγκων εγκεφάλου
Manolis Vavalis
 
Mathematical Models of Tumor Invasion
Sean Davis
 
2D CAT Based Modeling of Tumour Growth and Drug Transport
Editor IJMTER
 
Principles of chemotherapy
Sheetal R Kashid
 
ThinkBayes: chapter 13  simulation
Jungkyu Lee
 
Ad

More from Alex (Oleksiy) Varfolomiyev (7)

PDF
An Efficient Boundary Integral Method for Stiff Fluid Interface Problems
Alex (Oleksiy) Varfolomiyev
 
PDF
A non-stiff numerical method for 3D interfacial flow of inviscid fluids.
Alex (Oleksiy) Varfolomiyev
 
PDF
A non-stiff boundary integral method for internal waves
Alex (Oleksiy) Varfolomiyev
 
PDF
Optimal Finite Difference Grids for Elliptic and Parabolic PDEs with Applicat...
Alex (Oleksiy) Varfolomiyev
 
KEY
Neural Networks and Applications
Alex (Oleksiy) Varfolomiyev
 
KEY
Optimal Finite Difference Grids
Alex (Oleksiy) Varfolomiyev
 
PPT
Phase Field Method
Alex (Oleksiy) Varfolomiyev
 
An Efficient Boundary Integral Method for Stiff Fluid Interface Problems
Alex (Oleksiy) Varfolomiyev
 
A non-stiff numerical method for 3D interfacial flow of inviscid fluids.
Alex (Oleksiy) Varfolomiyev
 
A non-stiff boundary integral method for internal waves
Alex (Oleksiy) Varfolomiyev
 
Optimal Finite Difference Grids for Elliptic and Parabolic PDEs with Applicat...
Alex (Oleksiy) Varfolomiyev
 
Neural Networks and Applications
Alex (Oleksiy) Varfolomiyev
 
Optimal Finite Difference Grids
Alex (Oleksiy) Varfolomiyev
 
Phase Field Method
Alex (Oleksiy) Varfolomiyev
 
Ad

Recently uploaded (20)

PDF
Therapeutic, Functional and Neutralizing Antibodies Recognize Discontinuous, ...
stavejw
 
PDF
CERT Basic Training PTT, Brigadas comunitarias
chavezvaladezjuan
 
PDF
Cultivation and goods of microorganisms-4.pdf
adimondal300
 
PDF
A Classification of Metamorphic Rocks Canada12.pdf
ErturulKanmaz
 
PPTX
animal form and function zoology miller harley
sarmadbilal3
 
PPTX
Cyclotron_Presentation_theory, designMSc.pptx
MohamedMaideen12
 
PDF
Human-to-Robot Handovers track - RGMC - ICRA 2025
Alessio Xompero
 
PPTX
Respiratory and Circulatory Sytems.pptpptx
AngeloAngeles17
 
PPTX
GDG AI for Science - Why JAX and Flax NNX
gdgforscience
 
PDF
Global Health Initiatives: Lessons from Successful Programs (www.kiu.ac.ug)
publication11
 
PDF
Thermal stratification in lakes-J. Bovas Joel.pdf
J. Bovas Joel BFSc
 
PDF
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Eviden...
Sérgio Sacani
 
PDF
HOW TO DEAL WITH THREATS FROM THE FORCES OF NATURE FROM OUTER SPACE.pdf
Faga1939
 
PPTX
PROTECTED CULTIVATION ASSIGNMENT 2..pptx
RbDharani
 
PDF
Evidence for a sub-Jovian planet in the young TWA 7 disk
Sérgio Sacani
 
PDF
RANKING THE MICRO LEVEL CRITICAL FACTORS OF ELECTRONIC MEDICAL RECORDS ADOPTI...
hiij
 
PDF
Electromagnetism 3.pdf - AN OVERVIEW ON ELECTROMAGNETISM
kaustavsahoo94
 
PDF
SCH 4103_Fibre Technology & Dyeing_07012020.pdf
samwelngigi37
 
PDF
Agentic AI: Autonomy, Accountability, and the Algorithmic Society
vs5qkn48td
 
PPTX
Philippine_Literature_Precolonial_Period_Designed.pptx
josedalagdag5
 
Therapeutic, Functional and Neutralizing Antibodies Recognize Discontinuous, ...
stavejw
 
CERT Basic Training PTT, Brigadas comunitarias
chavezvaladezjuan
 
Cultivation and goods of microorganisms-4.pdf
adimondal300
 
A Classification of Metamorphic Rocks Canada12.pdf
ErturulKanmaz
 
animal form and function zoology miller harley
sarmadbilal3
 
Cyclotron_Presentation_theory, designMSc.pptx
MohamedMaideen12
 
Human-to-Robot Handovers track - RGMC - ICRA 2025
Alessio Xompero
 
Respiratory and Circulatory Sytems.pptpptx
AngeloAngeles17
 
GDG AI for Science - Why JAX and Flax NNX
gdgforscience
 
Global Health Initiatives: Lessons from Successful Programs (www.kiu.ac.ug)
publication11
 
Thermal stratification in lakes-J. Bovas Joel.pdf
J. Bovas Joel BFSc
 
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Eviden...
Sérgio Sacani
 
HOW TO DEAL WITH THREATS FROM THE FORCES OF NATURE FROM OUTER SPACE.pdf
Faga1939
 
PROTECTED CULTIVATION ASSIGNMENT 2..pptx
RbDharani
 
Evidence for a sub-Jovian planet in the young TWA 7 disk
Sérgio Sacani
 
RANKING THE MICRO LEVEL CRITICAL FACTORS OF ELECTRONIC MEDICAL RECORDS ADOPTI...
hiij
 
Electromagnetism 3.pdf - AN OVERVIEW ON ELECTROMAGNETISM
kaustavsahoo94
 
SCH 4103_Fibre Technology & Dyeing_07012020.pdf
samwelngigi37
 
Agentic AI: Autonomy, Accountability, and the Algorithmic Society
vs5qkn48td
 
Philippine_Literature_Precolonial_Period_Designed.pptx
josedalagdag5
 

Dynamic Tumor Growth Modelling

  • 1. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Dynamic Tumor Growth Modelling Prof. Thomas Witelski, Matthew Tanzy, Ben Owens, Oleksiy Varfolomiyev NJIT, June 2011
  • 2. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Cancer background Current medical understanding of some types of cancer 1 Disease starts from a primary tumor which grows in one location 2 The primary tumor may shed cancer cells which get carried to other parts of the body by the circulatory system
  • 3. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Colony size distribution of metastatic tumors with cell number x at time t is governed by von Foerster Equation ∂ρ ∂ + (g (x) ρ) = 0 (1) ∂t ∂x g (x) is tumor growth rate determined by dx = g (x) (2) dt Initially no metastatic tumor exists ρ(x, 0) = 0 (3) Intitial tumor birth rate ∞ g (1)ρ(1, t) = β(x)ρ(x, t)dx + β(xp (t)), (4) 1 β(x) tumor birth rate
  • 4. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Linear Growth Function: g (x) = k − Ex Assuming that in the beginning the main contribution to the tumor growth is given by primary tumor (i.e. neglecting the integral term in the BC) we have ∂ρ(x, t) ∂ (g (x)ρ(x, t)) + =0 (5) ∂t ∂x ρ(x, 0) = 0 (6) g (1)ρ(1, t) = β(xp ), β(x) = mx α (7) dx = k − Ex = g (x) (8) dt
  • 5. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Solution by the Method of Characteristics dρ = Eρ (9) dt along the characteristics given by dx = k − Ex (10) dt Therefore ρ(x, t) = ρ0 (x(0)) e Et (11) along the characteristics K x(t) = x(0)e −Et + 1 − e −Et (12) E
  • 6. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Solution by the Method of Characteristics First we use the BC to solve for ρ0 (x(0)), i.e. α K K (K − E ) ρ0 e Et + 1 − e Et e Et = m e −Et + 1 − e −Et E E Finally, the colony size distribution of metastatic tumors with cell number x at time t is α mE −α e −Et (E −k)2 ρ (x, t) = k−Ex k+ Ex−k
  • 7. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Growth function: g (x) = Ex Figure: ρ(x) at t = 10 with all parameters at 1
  • 8. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Growth function: g (x) = k − Ex Figure: ρ(x) at t = 10 with E = 0.2, k = 0.3, and m = 0.1
  • 9. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Single primary tumor size growth (no drug is active) Growth function: g (x) = (k − E (t)e −rt )x Here no drug is active (i.e. g (x) = kx) Figure: ρ(x) at t = 3650 with k = 0.006, α = 0.663, and m = 5.3 × 10−8
  • 10. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Single primary tumor size growth (drug effect) At time t1 drug is activated −r (t−t1 ) In a region of constant E : x (t) = Ae kt+Ee /r Figure: ρ(x) at t = 3650 with k = 0.006, α = 0.663, E = 0.0083, r = 1 × 10−5 and m = 5.3 × 10−8
  • 11. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Linear Separable Growth Rate Solution Now we have Linear Separable Growth Rate dx = [k − E (t)] x dt Colony size distribution m m t m −α−1− k−E (t) α+ k−E (t) ( k−E (s)ds )t ρ(x, t) = k−E (t) x e 0
  • 12. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Logistic Tumor Growth Rate Solution Now we have Tumor Growth Rate given by Logistic Equation dx = (k − Ex) x dt Resulting colony size distribution e −kt ρ(x, t) = 2 (E x0 (e kt −1)+k )
  • 13. Tumor Growth Model Method of Characteristics Solutions for various growth rates Future Work Future Work Fitting data to models, parameter estimation for multiple tumor and multiple patient data Incorporating limitations on measurable clinical data Comparing PDE, discrete population and polymerization models Solution of inverse problems for the birth rate Probability and statistics of different forms of clinical data