Dynamique des populations et résistance aux traitements : modèles mathématiques

Size: px
Start display at page:

Download "Dynamique des populations et résistance aux traitements : modèles mathématiques"

Transcription

1 Dynamique des populations et résistance aux traitements : modèles mathématiques Alexander Lorz 1 R. Chisholm, J. Clairambault, A. Escargueil, M.E. Hochberg, T. Lorenzi, P. Markowich, B. Perthame, E. Trélat 1 Laboratoire Jacques-Louis Lions, Université Pierre et Marie Curie Team MAMBA, INRIA Rocquencourt Séminaire du Laboratoire Jacques-Louis Lions, April 11th, 2014 A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

2 Outline Outline 1 Motivation A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

3 Outline Outline 1 Motivation 2 A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

4 Outline Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

5 Outline Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells 4 Heterogeneity in the model A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

6 Motivation Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells 4 Heterogeneity in the model A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

7 Motivation Motivation: Understand Darwinian evolution A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

8 Motivation Motivation: Understand resistance to therapy Many cancers escape therapy Cells adapt and become resistance to drug(s) Tumour as an ecological system, Darwinian evolution A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

9 Motivation Resistance to therapies A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

10 Motivation Darwinian evolution of a structured population density Population models are structured by a parameter representing a phenotypical trait. The length of the giraffe s neck spac or... A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

11 Motivation Darwinian evolution of a structured population density... The shape of beaks of Darwin s finches A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

12 ... Motivation Darwinian evolution of a structured population density The activity of ABC-transporters spac A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

13 Motivation Darwinian evolution of a structured population density We study the population dynamics under selection and mutations between the traits. Limited ressources Different traits have different capacities to use the environment Competition between the traits Selection The new born can have a trait slightly different from its parent one Mutations A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

14 Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells 4 Heterogeneity in the model A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

15 Darwinian evolution of a structured population density selection proliferation competition {}}{{}}{ t n(t,x) = n(t,x)[ p(x) d(x)i(t)] }{{ ) } x R d : phenotypic trait, n(t,x): density of trait x, R ( x,i(t) I(t) := n(t,x)dx R d I(t): total number of individuals, R(x,I): growth and death rates of trait x, A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

16 Example R = 4(x 0.3) 2 I n(t = 0,x) and R(x,I) t A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

17 Example R = 4(x 0.3) 2 I n(t = 1,x) and R(x,I) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

18 Example R = 4(x 0.3) 2 I n(t = 2,x) and R(x,I) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

19 Example R = 4(x 0.3) 2 I n(t = 3,x) and R(x,I) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

20 Example R = 4(x 0.3) 2 I n(t = 0,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

21 Example R = 4(x 0.3) 2 I n(t = 3,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

22 Example R = 4(x 0.3) 2 I n(t = 10,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

23 Example R = 4(x 0.3) 2 I n(t = 40,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

24 Example R = 4(x 0.3) 2 I n(t = 100,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

25 Surface plot of n(t, x) from the side t A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

26 Surface plot of n(t, x) from the top x A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

27 BV-bound: I(t) t I t t n(t,x) = t n(t,x)r ( x,i(t) ) + n(t,x)r I ( x,i(t) )İ(t) t t n(t,x) = n(t,x)r 2( x,i(t) ) ( )İ(t) + n(t,x)r I x,i(t) d İ(t) = n(t,x)r 2( x,i(t) ) dx + dt }{{} İ(t) ( ) n(t,x)r I x,i(t) dx =: J(t) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

28 BV-bound: I(t) t I t t n(t,x) = t n(t,x)r ( x,i(t) ) + n(t,x)r I ( x,i(t) )İ(t) t t n(t,x) = n(t,x)r 2( x,i(t) ) ( )İ(t) + n(t,x)r I x,i(t) d n(t,x)r 2( x,i(t) ) dx + dt İ(t) ( ) n(t,x)r I x,i(t) dx İ(t) = }{{} =: J(t) J(t) J(t) d dt d dt J(t) J(t) K J(t) J(t = 0) e Kt n(t,x)r I ( x,i(t) ) dx A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

29 BV-bound: I(t) t I t t n(t,x) = t n(t,x)r ( x,i(t) ) + n(t,x)r I ( x,i(t) )İ(t) t t n(t,x) = n(t,x)r 2( x,i(t) ) ( )İ(t) + n(t,x)r I x,i(t) d n(t,x)r 2( x,i(t) ) dx + dt İ(t) ( ) n(t,x)r I x,i(t) dx İ(t) = }{{} =: J(t) J(t) J(t) d dt d dt J(t) J(t) K J(t) J(t = 0) e Kt J(t) = İ(t) + 2J(t) n(t,x)r I ( x,i(t) ) dx J(t) dt I M + 2J(t = 0) e Kt dt A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

30 Idea of proof 0 < I m I(t) = n(t,x)dx I M & BV-bound I(t) t I 1) max x R(x,I ) > 0 x x 2) max x R(x,I ) < 0 3) max x R(x,I ) = 0 x n (x) = I δ(x x ) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

31 Generalisation of BV-bound t n(t,x) = n(t,x)r ( x,i(t) ), x R d, t 0, β d dt I(t) = Q( I(t),ρ(t) ), ρ(t) := n(t,x)dx. with β small A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

32 Darwinian evolution of a structured population density selection proliferation competition {}}{{}}{ t n(t,x) = n(t,x)[ p(x) d(x)i(t)] }{{ ) } x R d : phenotypic trait, n(t;x): density of trait x, R ( x,i(t) I(t) := n(t,x)dx R d I(t): total number of individuals, R(x,I): growth and death rates of trait x, A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

33 Darwinian evolution of a structured population density selection + mutations proliferation competition {}}{{}}{ t n(t,x) = n(t,x)[ p(x) d(x)i(t)] + n(t,x) }{{ ) } x R d : phenotypic trait, n(t;x): density of trait x, R ( x,i(t) I(t) := n(t,x)dx R d I(t): total number of individuals, R(x,I): growth and death rates of trait x, A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

34 Darwinian evolution of a structured population density We study the large time behavior of the population density for rare or small mutations We expect: discrete set of traits, extinction and branching A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

35 Rescaled equation for mutations with small effects ǫ t n ǫ = n ǫ R ( x,i ǫ (t) ) + ǫ 2 n ǫ, t > 0, x R d, with I ǫ (t) = n ǫ (x,t)dx. R d we observe very small mutations for a long time A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

36 Rescaled equation for mutations with small effects 4 ǫ t n ǫ = n ǫ R ( x,i ǫ (t) ) +ǫ 2 n ǫ, }{{}}{{} concentrates expands A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

37 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t = 0,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

38 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t =.125, x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

39 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t =.25,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

40 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t =.5,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

41 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t =.75,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

42 Example R = 4(x 0.3) 2 I + 0.5, ǫ = n(t = 1,x) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

43 Rescaled equation for ǫ = 0 0 = nr ( x,i(t) ) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

44 Rescaled equation for ǫ = 0 0 = nr ( x,i(t) ) 1 d and x R(x,I) monotonic R(x,I) x A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

45 Rescaled equation for ǫ = 0 0 = nr ( x,i(t) ) 1 d and x R(x,I) monotonic R(x,I) x n = ρ(t)δ(x x) 0 = R ( x,i(t) ) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

46 WKB ansatz δ(x x) 1 2πǫ e x x 2 /2ǫ = e ( x x 2 ǫln(2πǫ))/2ǫ. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

47 WKB ansatz δ(x x) 1 2πǫ e x x 2 /2ǫ = e ( x x 2 ǫln(2πǫ))/2ǫ. n ǫ (t,x) = e uǫ(t,x)/ǫ. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

48 WKB ansatz δ(x x) 1 2πǫ e x x 2 /2ǫ = e ( x x 2 ǫln(2πǫ))/2ǫ. n ǫ (t,x) = e uǫ(t,x)/ǫ. { t u ǫ = u ǫ 2 + R(x,I ǫ (t)) + ǫ u ǫ, x R d, t 0, u ǫ (t = 0) = ǫln(n 0 ǫ) := u 0 ǫ. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

49 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t =.0625,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

50 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t =.125,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

51 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t =.25,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

52 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t =.5,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

53 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t =.75,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

54 Example R = 4(x 0.3) 2 I ln(max(10 11,n(t = 1,x)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

55 Passing to the limit t u = u 2 + R(x,I(t)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

56 Passing to the limit t u = u 2 + R(x,I(t)) max R d u(t,x) = 0 supp n(t, ) {u(t, ) = 0} A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

57 Passing to the limit t u = u 2 + R(x,I(t)) max R d u(t,x) = 0 supp n(t, ) {u(t, ) = 0} open question: uniqueness A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

58 Challenges x(t) and I(t) are not continuous u(t,x) for t = t 1 < t 2 < t 3 A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

59 Results Theorem (A.L., S. Mirrahimi, B. Perthame) In R d, under concavity assumptions we obtain weak convergence in the sense of measures for a subsequence n ǫ : n ǫ (t,x) ǫ 0 ρ(t) δ ( x x(t) ), (initially one Dirac one Dirac for all times). There are many by many authors including: O. Diekmann, P.-E. Jabin, L. Desvillettes, St. Mischler, G. Barles, S. Genieys, M. Gauduchon, S. Cuadrado, J. Carillo, S. Mirrahimi, P. E. Souganidis A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

60 Can we get more on x(t)? A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

61 Canonical Equation u ǫ (t, x ǫ (t)) = 0 t u ǫ(t, x ǫ (t)) + D 2 xu ǫ (t, x ǫ (t)) x ǫ (t) = 0 D 2 xu ǫ (t, x ǫ (t)) x ǫ (t) = t u ǫ(t, x ǫ (t)) = x R( x ǫ (t),i ǫ (t)) ǫ x u ǫ A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

62 Canonical Equation Passing to the limit d dt x(t) = ( D2 1 u) R( x(t),ī(t)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

63 Canonical Equation Passing to the limit d dt x(t) = ( D2 1 u) R( x(t),ī(t)) Let s illustrate that A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

64 Canonical Equation Passing to the limit d dt x(t) = ( D2 1 u) R( x(t),ī(t)) Let s illustrate that A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

65 Interplay of Canonical Equation and R ( x(t), Ī(t)) = 0 d dt x(t) = ( D2 1 u) R( x(t),ī(t)) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

66 with cut-off: Motivation The tail problem A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

67 with cut-off: Motivation The tail problem A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

68 with cut-off: Motivation The tail problem A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

69 with cut-off: Motivation The tail problem Further out: number of individuals < 1 should be zero A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

70 with cut-off: Motivation penalise the tails t n σ (t,x) = n σ (t,x)r ( x,i(t) ) + n σ (t,x) n σ σ 21 {n σ µ σ α } seems to cut exponential tails A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

71 The tail problem with cut-off: simplified version t n σ (t,x) = n σ (t,x)r ( x,i(t) ) + n σ (t,x) n σ σ 21 {n σ µ σ α } f = n σ (t,x) + n σ (t,x) n σ σ 21 {n σ µ σ α } solutions very non-unique variational solution maximal solution minimal solution on the right scale: interesting free boundary problem A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

72 Back to cancer Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells 4 Heterogeneity in the model A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

73 Back to cancer : cancer cells Biological problem a resistant sub-populations is selected can be viewed as Darwinian evolution this can be formulated in the frame work above t n = nr ( x,i(t) ) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

74 Back to cancer : cancer cells 1d model: cancer cells with therapy t n(x,t) = [ p(x) d(x)i(t) ] n(x, t) n(x,t): number of cells with trait x resistance phenotype x [0,1], x larger more resistant more resistant less proliferative (cost of resistance) two possible targets for therapy: slowing down proliferation or increasing cell death A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

75 Back to cancer : cancer cells 1d model: cancer cells with therapy t n(x,t) = [ ] p(x) d(x)i(t) µ(x)c n(x, t) n(x,t): number of cells with trait x resistance phenotype x [0,1], x larger more resistant more resistant less proliferative (cost of resistance) two possible targets for therapy: slowing down proliferation or increasing cell death cytotoxic drug c A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

76 8 2 Back to cancer : cancer cells 1d model: cancer cells with therapy t n(x,t) = [ ] p(x) d(x)i(t) µ(x)c n(x, t) level-sets of n(x,t)/ n(x,t)dx x t A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

77 Back to cancer : healthy and cancer cells How to avoid resistance? Drug resistance is a major obstacle for anti-cancer therapy using different drugs to reduce resistance. Cytotoxic drugs: killing cells Cytostatic drugs: slowing down cell proliferation Cytotoxic drugs + cytostatic drugs A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

78 Back to cancer : healthy and cancer cells Polychemotherapy cancer cells n C compare with effects on healthy cells n H t n C (x,t) = [ ] p C (x) d C (x)i C (t) c 1 µ C (x) n C (x,t), where I C = ρ c 1 : cytotoxic therapy A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

79 Back to cancer : healthy and cancer cells Polychemotherapy cancer cells n C compare with effects on healthy cells n H t n C (x,t) = [ ] 1 p C (x) d C (x)i C (t) c 1 µ C (x) n C (x,t), 1 + α C c 2 where I C = ρ c 1 : cytotoxic therapy c 2 : cytostatic therapy A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

80 Back to cancer : healthy and cancer cells Polychemotherapy compare with effects on cancer n C and healthy cells n H t n C (x,t) = [ ] 1 p C (x) d C (x)i C (t) c 1 µ C (x) n C (x,t), 1 + α C c 2 where space c 1 : cytotoxic therapy c 2 : cytostatic therapy A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

81 Back to cancer : healthy and cancer cells Polychemotherapy compare with effects on cancer n C and healthy cells n H t n C (x,t) = [ ] 1 p C (x) d C (x)i C (t) c 1 µ C (x) n C (x,t), 1 + α C c 2 t n H (x,t) = [ ] 1 p H (x) d H (x)i H (t) c 1 µ H (x) n H (x,t), 1 + α H c 2 where space c 1 : cytotoxic therapy c 2 : cytostatic therapy A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

82 Back to cancer : healthy and cancer cells Polychemotherapy compare with effects on cancer n C and healthy cells n H t n C (x,t) = [ ] 1 p C (x) d C (x)i C (t) c 1 µ C (x) n C (x,t), 1 + α C c 2 t n H (x,t) = [ ] 1 p H (x) d H (x)i H (t) c 1 µ H (x) n H (x,t), 1 + α H c 2 where I C = a CH ρ H + a CC ρ C c 1 : cytotoxic therapy c 2 : cytostatic therapy I H = a HH ρ H + a HC ρ C A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

83 Back to cancer : healthy and cancer cells Parameter functions: proliferation/death Cost of resistance Cytotoxic drugs Adapted to cancer cells p > 0, p < 0. µ > 0, µ < 0. µ H < µ C α H < α C. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

84 Back to cancer : healthy and cancer cells Effect of therapy on R: p C (x) c 1 µ C (x) space c 1 = A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

85 Back to cancer : healthy and cancer cells Effect of therapy on R: p C (x) c 1 µ C (x) space c 1 = A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

86 Back to cancer : healthy and cancer cells Effect of therapy on R: p C (x) c 1 µ C (x) space c 1 = A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

87 Back to cancer : healthy and cancer cells Effect of therapy on R: p C (x) c 1 µ C (x) space c 1 = A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

88 Back to cancer : healthy and cancer cells Preliminaries: initial conditions 0.5 n H (t=0,x) x A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

89 Back to cancer : healthy and cancer cells Cytotoxic agents only (c 1 = 0, 1.75, 3.5 with c 2 = 0) n C (t=2000,x) n H (t=2000,x) 60 c 1 =0, c 2 =0 25 c 1 =0, c 2 =0 c 1 =1.75, c 2 =0 c 1 =1.75, c 2 =0 c 1 =3.5, c 2 =0 c 1 =3.5, c 2 = x cancer cells x healthy cells c 1 = 0 selection for strong proliferative potential. c 1 increases both the number of cancer cells and healthy cells become smaller + selection for resistance. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

90 Back to cancer : healthy and cancer cells Cytostatic agents only (c 2 = 1, 3, 7 with c 1 = 0) n C (t=2000,x) n H (t=2000,x) 30 c 1 =0, c 2 =1 30 c 1 =0, c 2 =1 c 1 =0, c 2 =3 c 1 =0, c 2 =3 c 1 =0, c 2 =7 c 1 =0, c 2 = x cancer cells x healthy cells c 2 increases the dynamics of healthy cells is kept unaltered while the proliferation of cancer cells is reduced + the dynamics of cancer cells is slowed down. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

91 Back to cancer : healthy and cancer cells Simultaneous action of cytostatic and cytotoxic agents (c 1 = c 2 = 0, 1, 1.5, 2) cancer cells A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

92 Back to cancer : healthy and cancer cells Simultaneous action of cytostatic and cytotoxic agents (c 1 = c 2 = 0, 1, 1.5, 2) healthy cells A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

93 Back to cancer : healthy and cancer cells Action of cytostatic and cytotoxic agents A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

94 Back to cancer : healthy and cancer cells Simultaneous action of cytostatic and cytotoxic agents cytostatic agents slow down the growth of resistant cells cytotoxic agents kill the cancer cells A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

95 Heterogeneity in the model Outline 1 Motivation 2 3 Back to cancer : cancer cells : healthy and cancer cells 4 Heterogeneity in the model A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

96 Heterogeneity in the model Resistance to therapies A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

97 Heterogeneity in the model Resistance to therapies In ecology, competitive exclusion principle : In a ecological niche only one single species can exist Why is there such a large phenotypical heterogeneity in tumour cells? A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

98 Heterogeneity in the model Resistance to therapies In ecology, competitive exclusion principle : In a ecological niche only one single species can exist Why is there such a large phenotypical heterogeneity in tumour cells? A large phenotypical heterogeneity in tumour cells can make them more resistant. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

99 Heterogeneity in the model for a spherical tumor Collaboration with J. Clairambault, A. Escargueil, T. Lorenzi, B. Perthame Nutrients and cytotoxic agents diffuse in the tumour through the boundary. A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

100 Heterogeneity in the model for spherical tumour t n(t,r,x) = [ { proliferation death }}{{}}{ p(x)s(t,r) d(x)ρ(t) ] n(t,r,x), n(t,r,x): cancer cells density at radius r [0,1] with resistance phenotype x [0, 1] ρ(t): total number of cancer cells s(t, r): nutrients c(t, r): cytotoxic therapy p(x),d(x),µ(x) given, σ s,σ c,γ s,γ c constants A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

101 Heterogeneity in the model for spherical tumour t n(t,r,x) = [ { proliferation death }}{{}}{ p(x)s(t,r) d(x)ρ(t) effect of therapy {}}{ µ(x)c(t,r) ] n(t,r,x), n(t,r,x): cancer cells density at radius r [0,1] with resistance phenotype x [0, 1] ρ(t): total number of cancer cells s(t, r): nutrients c(t, r): cytotoxic therapy p(x),d(x),µ(x) given, σ s,σ c,γ s,γ c constants A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

102 Heterogeneity in the model for spherical tumour t n(t,r,x) = [ { proliferation death }}{{}}{ p(x)s(t,r) d(x)ρ(t) [ σ s s(t,r) + }{{} γ }{{} s + p(x)n(t,r,x)dx diffusion degradation } {{ } consumption effect of therapy {}}{ µ(x)c(t,r) ] n(t,r,x), ] s(t,r) = 0, n(t,r,x): cancer cells density at radius r [0,1] with resistance phenotype x [0, 1] ρ(t): total number of cancer cells s(t, r): nutrients c(t, r): cytotoxic therapy p(x),d(x),µ(x) given, σ s,σ c,γ s,γ c constants A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

103 Heterogeneity in the model for spherical tumour t n(t,r,x) = [ { proliferation death }}{{}}{ p(x)s(t,r) d(x)ρ(t) [ σ s s(t,r) + }{{} γ }{{} s + p(x)n(t,r,x)dx diffusion degradation [ σ c c(t,r) + γ c + } {{ } consumption µ(x)n(t,r,x)dx effect of therapy {}}{ µ(x)c(t,r) ] n(t,r,x), ] s(t,r) = 0, ] c(t,r) = 0. n(t,r,x): cancer cells density at radius r [0,1] with resistance phenotype x [0, 1] ρ(t): total number of cancer cells s(t, r): nutrients c(t, r): cytotoxic therapy p(x),d(x),µ(x) given, σ s,σ c,γ s,γ c constants A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

104 Heterogeneity in the model for spherical tumour t n(t,r,x) = [ { proliferation death }}{{}}{ p(x)s(t,r) d(x)ρ(t) [ σ s s(t,r) + }{{} γ }{{} s + p(x)n(t,r,x)dx diffusion degradation [ σ c c(t,r) + γ c + } {{ } consumption µ(x)n(t,r,x)dx effect of therapy {}}{ µ(x)c(t,r) ] n(t,r,x), ] s(t,r) = 0, ] c(t,r) = 0. boundary data: s(t,r = 1),c(t,r = 1) n(t,r,x): cancer cells density at radius r [0,1] with resistance phenotype x [0, 1] ρ(t): total number of cancer cells s(t, r): nutrients c(t, r): cytotoxic therapy p(x),d(x),µ(x) given, σ s,σ c,γ s,γ c constants A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

105 Heterogeneity in the model nutrient profile s(t,r) r A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

106 Heterogeneity in the model Without cytotoxic therapy phenotype averaged over radius r phenotype at radius r n(t,r,x)dr/ρ(t) n(t,r,x)/ n(t,r,x)dx x x decreases of heterogeneity r A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

107 Heterogeneity in the model With cytotoxic therapy phenotype averaged over radius r phenotype at radius r n(t,r,x)dr/ρ(t) n(t,r,x)/ n(t,r,x)dx x x decreases of heterogeneity r A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

108 Heterogeneity in the model Definition of heterogeneity averaged dominant trait X(t) := 1 ρ(t) xn(t,r,x)dxdr heterogeneity or trait variance xn(t,r,x)dx h(t) := ρ(t,r) X(t) 2 dr where ρ(t,r) = n(t,r,x)dx A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

109 Heterogeneity in the model Future Work analytic results e.g. reduction of heterogeneity by chemotherapy include growing or shrinking tumour model validation with data therapy optimisation: polychemotherapy, adaptive therapy A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

110 Heterogeneity in the model THANK YOU Gatenby R.A., A change of strategy in the war on cancer, Nature, 459 (2009) A. Lorz (LJLL, UPMC & INRIA MAMBA) Dynamique des populations Apr 11th, / 62

Effects of space structure and combination therapies on phenotypic heterogeneity and drug resistance in solid tumors

Effects of space structure and combination therapies on phenotypic heterogeneity and drug resistance in solid tumors Effects of space structure and combination therapies on phenotypic heterogeneity and drug resistance in solid tumors Alexander Lorz 123 Tommaso Lorenzi 123 Jean Clairambault 312 Alexandre Escargueil 45

More information

Why Is Evolution Important in Cancer and What Mathematics Should Be Used to Treat Cancer? Focus on Drug Resistance

Why Is Evolution Important in Cancer and What Mathematics Should Be Used to Treat Cancer? Focus on Drug Resistance Why Is Evolution Important in Cancer and What Mathematics Should Be Used to Treat Cancer? Focus on Drug Resistance Luis Almeida, Rebecca Chisholm, Jean Clairambault, Tommaso Lorenzi, Alexander Lorz, Camille

More information

Mathematics Meets Oncology

Mathematics Meets Oncology .. Mathematics Meets Oncology Mathematical Oncology Philippe B. Laval Kennesaw State University November 12, 2011 Philippe B. Laval (Kennesaw State University)Mathematics Meets Oncology November 12, 2011

More information

Modeling Multi-Mutation And Drug Resistance: A Case of Immune-Suppression

Modeling Multi-Mutation And Drug Resistance: A Case of Immune-Suppression Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 14, Number 6 (2018), pp. 787-807 Research India Publications http://www.ripublication.com/gjpam.htm Modeling Multi-Mutation And Drug

More information

Cancer Treatment Using Multiple Chemotheraputic Agents Subject to Drug Resistance

Cancer Treatment Using Multiple Chemotheraputic Agents Subject to Drug Resistance Cancer Treatment Using Multiple Chemotheraputic Agents Subject to Drug Resistance J. J. Westman Department of Mathematics University of California Box 951555 Los Angeles, CA 90095-1555 B. R. Fabijonas

More information

Modeling the Effects of HIV on the Evolution of Diseases

Modeling the Effects of HIV on the Evolution of Diseases 1/20 the Evolution of Mentor: Nina Fefferman DIMACS REU July 14, 2011 2/20 Motivating Questions How does the presence of HIV in the body changes the evolution of pathogens in the human body? How do different

More information

Monte Carlo Analysis of Univariate Statistical Outlier Techniques Mark W. Lukens

Monte Carlo Analysis of Univariate Statistical Outlier Techniques Mark W. Lukens Monte Carlo Analysis of Univariate Statistical Outlier Techniques Mark W. Lukens This paper examines three techniques for univariate outlier identification: Extreme Studentized Deviate ESD), the Hampel

More information

A mathematical model for the primary tumor of mcrc

A mathematical model for the primary tumor of mcrc A mathematical model for the primary tumor of mcrc Marta Leocata Joint work with F.Flandoli, C. Ricci, M.C. Polito, V. De Mattei December 2, 2016 University of Pisa Plan of the talk General Project; A

More information

Maria-Athina Altzerinakou1, Xavier Paoletti2. 9 May, 2017

Maria-Athina Altzerinakou1, Xavier Paoletti2. 9 May, 2017 An adaptive design for the identification of the optimal dose using joint modelling of efficacy and toxicity in phase I/II clinical trials of molecularly targeted agents Maria-Athina Altzerinakou1, Xavier

More information

Mathematical modelling of spatio-temporal glioma evolution

Mathematical modelling of spatio-temporal glioma evolution Papadogiorgaki et al. Theoretical Biology and Medical Modelling 213, 1:47 RESEARCH Open Access Mathematical modelling of spatio-temporal glioma evolution Maria Papadogiorgaki 1*, Panagiotis Koliou 2, Xenofon

More information

Modeling origin and natural evolution of low-grade gliomas

Modeling origin and natural evolution of low-grade gliomas Modeling origin and natural evolution of low-grade gliomas Mathilde Badoual Paris Diderot University, IMNC lab 2nd HTE workshop: Mathematical & Computer Modeling to study tumors heterogeneity in its ecosystem,

More information

Epidemiological Model of HIV/AIDS with Demographic Consequences

Epidemiological Model of HIV/AIDS with Demographic Consequences Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 5, Number 1 (2014), pp. 65-74 International Research Publication House http://www.irphouse.com Epidemiological Model of HIV/AIDS with

More information

Structured models for dengue epidemiology

Structured models for dengue epidemiology Structured models for dengue epidemiology submitted by Hannah Woodall for the degree of Doctor of Philosophy of the University of Bath Department of Mathematical Sciences September 24 COPYRIGHT Attention

More information

Mathematics and Physics of Cancer: Questions. Robijn Bruinsma, UCLA KITP Colloquium May 6, ) Cancer statistics and the multi-stage model.

Mathematics and Physics of Cancer: Questions. Robijn Bruinsma, UCLA KITP Colloquium May 6, ) Cancer statistics and the multi-stage model. Mathematics and Physics of Cancer: Questions Robijn Bruinsma, UCLA KITP Colloquium May 6, 2009 1) Cancer statistics and the multi-stage model. 2) Cancer microevolution and clonal expansion. 3) Metastasis:

More information

Stochastic and Numerical Modeling of HIV/AIDS. Spread in a Complex System and Its Application to. the HIV Cases in Indonesia

Stochastic and Numerical Modeling of HIV/AIDS. Spread in a Complex System and Its Application to. the HIV Cases in Indonesia Applied Mathematical Sciences, Vol. 9, 2015, no. 122, 6095-6106 HIKAI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.59564 Stochastic and Numerical Modeling of HIV/AIDS Spread in a Complex System

More information

ON CHARACTERIZING THE BIVARIATE EXPONENTIAL

ON CHARACTERIZING THE BIVARIATE EXPONENTIAL Ann. Inst. Statist. Math. Vol. 40, No. 2, 267-271 (1988) ON CHARACTERIZING THE BIVARIATE EXPONENTIAL AND GEOMETRIC DISTRIBUTIONS K. R. MURALEEDHARAN NAIR AND N. UNNIKRISHNAN NAIR Department of Mathematics

More information

arxiv: v2 [math.oc] 15 Nov 2016

arxiv: v2 [math.oc] 15 Nov 2016 Optimizing chemoradiotherapy to target multi-site metastatic disease and tumor growth arxiv:1603.00349v2 [math.oc] 15 Nov 2016 Hamidreza Badri, 1 Ehsan Salari, 2 Yoichi Watanabe, 3 Kevin Leder 1 1 Department

More information

Conditional Distributions and the Bivariate Normal Distribution. James H. Steiger

Conditional Distributions and the Bivariate Normal Distribution. James H. Steiger Conditional Distributions and the Bivariate Normal Distribution James H. Steiger Overview In this module, we have several goals: Introduce several technical terms Bivariate frequency distribution Marginal

More information

Modeling of cancer virotherapy with recombinant measles viruses

Modeling of cancer virotherapy with recombinant measles viruses Modeling of cancer virotherapy with recombinant measles viruses Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Collaborators: Krešimir Josić Dept. of Mathematics, University

More information

Mathematical Model Solid Tumor at the Stage of Angiogenesis with Immune Response

Mathematical Model Solid Tumor at the Stage of Angiogenesis with Immune Response Mathematical Model Solid Tumor at the Stage of Angiogenesis with Immune esponse Deep Shikha Dixit, Deepak Kumar, Sanjeev Kr,ajesh Johri. Department of Mathematics, CET IILM, Gr. Noida.. Department of Mathematics,

More information

The Chemostat: Stability at Steady States. Chapter 5: Linear & Non-Linear Interaction Models. So, in dimensional form, α 1 > 1 corresponds to

The Chemostat: Stability at Steady States. Chapter 5: Linear & Non-Linear Interaction Models. So, in dimensional form, α 1 > 1 corresponds to Introduction & Simple Models Logistic Growth Models The Chemostat: Stability at Steady States 1 So, in dimensional form, α 1 > 1 corresponds to K max < V F. As K max is max bacterial repro rate with unlimited

More information

Optimisation of cancer drug treatments using cell population dynamics

Optimisation of cancer drug treatments using cell population dynamics Optimisation of cancer drug treatments using cell population dynamics Frédérique Billy 1, Jean Clairambault 1, and Olivier Fercoq 2 1 INRIA BANG team, BP 105, F78153 Rocquencourt, and LJLL, UPMC, 4 Place

More information

Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy

Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy Hang Xie 1, Yang Jiao 2, Qihui Fan 3, Miaomiao Hai 1, Jiaen Yang 1, Zhijian Hu 1, Yue Yang 4,

More information

Simulating the Tumor Growth with Cellular Automata Models

Simulating the Tumor Growth with Cellular Automata Models Simulating the Tumor Growth with Cellular Automata Models S. Zouhri Université Hassan II- Mohammédia, Faculté des Sciences Ben M'sik Département de Mathématiques, B.7955, Sidi Othmane, Casablanca, Maroc

More information

A Method for Utilizing Bivariate Efficacy Outcome Measures to Screen Agents for Activity in 2-Stage Phase II Clinical Trials

A Method for Utilizing Bivariate Efficacy Outcome Measures to Screen Agents for Activity in 2-Stage Phase II Clinical Trials A Method for Utilizing Bivariate Efficacy Outcome Measures to Screen Agents for Activity in 2-Stage Phase II Clinical Trials Michael W. Sill 1 and Greg Yothers 2 1 Adjunct Instructor; Research Assistant

More information

Volume Effects in Chemotaxis

Volume Effects in Chemotaxis Volume Effects in Chemotaxis Thomas Hillen University of Alberta supported by NSERC with Kevin Painter (Edinburgh), Volume Effects in Chemotaxis p.1/48 Eschirichia coli Berg - Lab (Harvard) Volume Effects

More information

Number 263 October 2015 THE HEALTH HUMP. Holger Strulik ISSN:

Number 263 October 2015 THE HEALTH HUMP. Holger Strulik ISSN: Number 263 October 2015 THE HEALTH HUMP Holger Strulik ISSN: 1439-2305 The Health Hump Holger Strulik October 2015. Abstract. In this paper, I suggest a novel explanation for a hump-shaped ageconsumption

More information

Hypothesis testing: comparig means

Hypothesis testing: comparig means Hypothesis testing: comparig means Prof. Giuseppe Verlato Unit of Epidemiology & Medical Statistics Department of Diagnostics & Public Health University of Verona COMPARING MEANS 1) Comparing sample mean

More information

Empirical assessment of univariate and bivariate meta-analyses for comparing the accuracy of diagnostic tests

Empirical assessment of univariate and bivariate meta-analyses for comparing the accuracy of diagnostic tests Empirical assessment of univariate and bivariate meta-analyses for comparing the accuracy of diagnostic tests Yemisi Takwoingi, Richard Riley and Jon Deeks Outline Rationale Methods Findings Summary Motivating

More information

Defining kinetic properties of HIV-specific CD8 + T-cell responses in acute infection

Defining kinetic properties of HIV-specific CD8 + T-cell responses in acute infection Defining kinetic properties of HIV-specific CD8 + T-cell responses in acute infection Yiding Yang 1 and Vitaly V. Ganusov 1,2,3 1 Department of Microbiology, University of Tennessee, Knoxville, TN 37996,

More information

Optimal Control of The HIV-Infection Dynamics

Optimal Control of The HIV-Infection Dynamics Australian Journal of Basic and Applied Sciences, 5(9): 1051-1059, 2011 ISSN 1991-8178 Optimal Control of The HIV-Infection Dynamics 1 H. Ghiasi zadeh, 1 H. Chahkandi nejad, 1 R. Jahani, 2 H. Roshan Nahad

More information

Mathematical Model of Colorectal Cancer with Monoclonal Antibody Treatments

Mathematical Model of Colorectal Cancer with Monoclonal Antibody Treatments British Journal of Mathematics & Computer Science X(X): XX-XX, 2013 SCIENCEDOMAIN international www.sciencedomain.org Mathematical Model of Colorectal Cancer with Monoclonal Antibody Treatments Research

More information

Mathematical biology From individual cell behavior to biological growth and form

Mathematical biology From individual cell behavior to biological growth and form Mathematical biology From individual cell behavior to biological growth and form Lecture 8: Multiscale models Roeland Merks (1,2) (1) Centrum Wiskunde & Informatica, Amsterdam (2) Mathematical Institute,

More information

Avian influenza in migratory birds and basic reproductive ratio

Avian influenza in migratory birds and basic reproductive ratio Avian influenza in migratory birds and basic reproductive ratio Xiang-Sheng Wang Mprime Centre for Disease Modelling York University, Toronto (joint work with Jianhong Wu) Avian influenza in migratory

More information

MATHEMATICAL MODELS OF THE ADAPTIVE IMMUNE RESPONSE IN A NOVEL CANCER IMMUNOTHERAPY

MATHEMATICAL MODELS OF THE ADAPTIVE IMMUNE RESPONSE IN A NOVEL CANCER IMMUNOTHERAPY MATHEMATICAL MODELS OF THE ADAPTIVE IMMUNE RESPONSE IN A NOVEL CANCER IMMUNOTHERAPY by Bryan A. Dawkins A dissertation submitted to the faculty of The University of Central Oklahoma in partial fulfillment

More information

Supplementary Materials: Identifying critical transitions and their leading biomolecular networks in complex diseases

Supplementary Materials: Identifying critical transitions and their leading biomolecular networks in complex diseases Supplementary Materials: Identifying critical transitions and their leading biomolecular networks in complex diseases Rui Liu, Meiyi Li, Zhiping Liu, Jiarui Wu, Luonan Chen, Kazuyuki Aihara Contents A

More information

Liénard Oscillator Modeling of Bipolar Disorder

Liénard Oscillator Modeling of Bipolar Disorder Liénard Oscillator Modeling of Bipolar Disorder Jessica Snyder 5 August 23 Abstract Bipolar disorder, also known as manic-depression, is a disorder in which a person has mood swings out of proportion or

More information

Bayesian Models for Combining Data Across Subjects and Studies in Predictive fmri Data Analysis

Bayesian Models for Combining Data Across Subjects and Studies in Predictive fmri Data Analysis Bayesian Models for Combining Data Across Subjects and Studies in Predictive fmri Data Analysis Thesis Proposal Indrayana Rustandi April 3, 2007 Outline Motivation and Thesis Preliminary results: Hierarchical

More information

Optimal Delivery of Chemotherapeutic Agents in Cancer

Optimal Delivery of Chemotherapeutic Agents in Cancer 16th European Symposium on Computer Aided Process Engineering and 9th International Symposium on Process Systems Engineering W. Marquardt, C. Pantelides (Editors) 2006 Published by Elsevier B.V. 1643 Optimal

More information

Extrapolating Tumor Invasion Margins for Physiologically Determined Radiotherapy Regions

Extrapolating Tumor Invasion Margins for Physiologically Determined Radiotherapy Regions Extrapolating Tumor Invasion Margins for Physiologically Determined Radiotherapy Regions Ender Konuko glu 1, Olivier Clatz 1, Pierre-Yves Bondiau 2, Hervé Delingette 1, and Nicholas Ayache 1 1 Asclepios

More information

How Genetic Algorithms Can Improve a Pacemaker Effciency

How Genetic Algorithms Can Improve a Pacemaker Effciency How Genetic Algorithms Can Improve a Pacemaker Effciency Laurent Dumas Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie 75252 Paris Cedex 05, France dumas@ann.jussieu.fr Linda El Alaoui

More information

IMMUNOLOGICAL SURVEILLANCE AND NEOPLASTIC DEVELOPMENT

IMMUNOLOGICAL SURVEILLANCE AND NEOPLASTIC DEVELOPMENT ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 9, Number 1, Winter 1979 IMMUNOLOGICAL SURVEILLANCE AND NEOPLASTIC DEVELOPMENT GEORGE W. SWAN 1. Introduction. In 1908 Ehrlich suggested that one of the functions

More information

Part I. Boolean modelling exercises

Part I. Boolean modelling exercises Part I. Boolean modelling exercises. Glucose repression of Icl in yeast In yeast Saccharomyces cerevisiae, expression of enzyme Icl (isocitrate lyase-, involved in the gluconeogenesis pathway) is important

More information

Density dependent diffusion and spread of epidemics in a metapopulation model.

Density dependent diffusion and spread of epidemics in a metapopulation model. Erice (Italy) August 30 2015 September 5 2015. Density dependent diffusion and spread of epidemics in a metapopulation model. Marta Pellicer UNIVERSITAT DE GIRONA Catalunya (Spain) Joint work with A. Avinyó,

More information

Analyzing Immunotherapy and Chemotherapy of Tumors through Mathematical Modeling Summer Student-Faculty Research Project

Analyzing Immunotherapy and Chemotherapy of Tumors through Mathematical Modeling Summer Student-Faculty Research Project Analyzing Immunotherapy and Chemotherapy of Tumors through Mathematical Modeling Summer Student-Faculty Research Project by Lisette de Pillis, Weiqing Gu, William Chang, Lindsay Crowl, Eric Malm, Katherine

More information

Motivation Empirical models Data and methodology Results Discussion. University of York. University of York

Motivation Empirical models Data and methodology Results Discussion. University of York. University of York Healthcare Cost Regressions: Going Beyond the Mean to Estimate the Full Distribution A. M. Jones 1 J. Lomas 2 N. Rice 1,2 1 Department of Economics and Related Studies University of York 2 Centre for Health

More information

10-1 MMSE Estimation S. Lall, Stanford

10-1 MMSE Estimation S. Lall, Stanford 0 - MMSE Estimation S. Lall, Stanford 20.02.02.0 0 - MMSE Estimation Estimation given a pdf Minimizing the mean square error The minimum mean square error (MMSE) estimator The MMSE and the mean-variance

More information

Extrapolating Glioma Invasion Margin in Brain Magnetic Resonance Images: Suggesting New Irradiation Margins

Extrapolating Glioma Invasion Margin in Brain Magnetic Resonance Images: Suggesting New Irradiation Margins Extrapolating Glioma Invasion Margin in Brain Magnetic Resonance Images: Suggesting New Irradiation Margins Ender Konukoglu, Olivier Clatz, Pierre-Yves Bondiau, Hervé Delingette, Nicholas Ayache To cite

More information

MECHANISMS AND PATTERNS OF EVOLUTION

MECHANISMS AND PATTERNS OF EVOLUTION MECHANISMS AND PATTERNS OF EVOLUTION Evolution What is it again? Evolution is the change in allele frequencies of a population over generations Mechanisms of Evolution what can make evolution happen? 1.

More information

How Relevant is the Asymptomatic Population in Dengue Transmission?

How Relevant is the Asymptomatic Population in Dengue Transmission? Applied Mathematical Sciences, Vol. 12, 2018, no. 32, 1699-1708 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.810150 How Relevant is the Asymptomatic Population in Dengue Transmission?

More information

A mathematical model of tumor dynamics during stereotactic body radiation therapy for non-small cell lung cancer

A mathematical model of tumor dynamics during stereotactic body radiation therapy for non-small cell lung cancer A mathematical model of tumor dynamics during stereotactic body radiation therapy for non-small cell lung cancer Russell Injerd, Emma Turian Abstract Image guided radiotherapy allows tumor volume dynamics

More information

Pharmacokinetics and pharmacodynamics Peter Dalgaard April/May Introduction to metabolic models

Pharmacokinetics and pharmacodynamics Peter Dalgaard April/May Introduction to metabolic models Pharmacokinetics and pharmacodynamics Peter Dalgaard April/May 2003 Introduction to metabolic models 1 16 Ideas 1 Ideas Mainly follow Carson, Cobelli, and Finkelstein, Ch.2+5 1 Ideas Mainly follow Carson,

More information

0.1 Immunology - HIV/AIDS. 0.2 History & biology of HIV

0.1 Immunology - HIV/AIDS. 0.2 History & biology of HIV 0.1 mmunology - HV/ADS n our previous models we assumed a homogeneous population where everyone was susceptible and infectious to the same degree. n contrast, the dynamics of STDs is affected by the general

More information

A Theory of Rational Junk-Food Consumption

A Theory of Rational Junk-Food Consumption University of Wollongong Research Online Faculty of Business - Economics Working Papers Faculty of Business 2002 A Theory of Rational Junk-Food Consumption Amnon Levy University of Wollongong, levy@uow.edu.au

More information

Problem for the Optimal Control of Cigarette Addiction

Problem for the Optimal Control of Cigarette Addiction Contemporary Engineering Sciences, Vol. 12, 2019, no. 1, 41-49 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2019.914 Problem for the Optimal Control of Cigarette Addiction Anibal Muñoz L.,

More information

Modeling the Diffusion of TGF-β1 from a Fibrin Scaffold through Alveolar Bone Martin Barrio and Sindhu Raghunandan 10/28/2011 BENG221

Modeling the Diffusion of TGF-β1 from a Fibrin Scaffold through Alveolar Bone Martin Barrio and Sindhu Raghunandan 10/28/2011 BENG221 10/28/2011 BENG221 I. Background Periodontitis is defined as infection and inflammation of the primary tissues surrounding one or more teeth [5]. While the initial stages involve the infection and inflammation

More information

Statistical Methods for Improving Biosurveillance System Performance

Statistical Methods for Improving Biosurveillance System Performance Statistical Methods for Improving Biosurveillance System Performance Associate Professor Ron Fricker Naval Postgraduate School Army Conference on Applied Statistics October 17, 27 1 The New Status Quo?

More information

SCHOOL OF MATHEMATICS AND STATISTICS

SCHOOL OF MATHEMATICS AND STATISTICS Data provided: Tables of distributions MAS603 SCHOOL OF MATHEMATICS AND STATISTICS Further Clinical Trials Spring Semester 014 015 hours Candidates may bring to the examination a calculator which conforms

More information

UNIVERSITY OF CALIFORNIA SANTA CRUZ A STOCHASTIC DYNAMIC MODEL OF THE BEHAVIORAL ECOLOGY OF SOCIAL PLAY

UNIVERSITY OF CALIFORNIA SANTA CRUZ A STOCHASTIC DYNAMIC MODEL OF THE BEHAVIORAL ECOLOGY OF SOCIAL PLAY . UNIVERSITY OF CALIFORNIA SANTA CRUZ A STOCHASTIC DYNAMIC MODEL OF THE BEHAVIORAL ECOLOGY OF SOCIAL PLAY A dissertation submitted in partial satisfaction of the requirements for the degree of BACHELOR

More information

Real-time computational attention model for dynamic scenes analysis

Real-time computational attention model for dynamic scenes analysis Computer Science Image and Interaction Laboratory Real-time computational attention model for dynamic scenes analysis Matthieu Perreira Da Silva Vincent Courboulay 19/04/2012 Photonics Europe 2012 Symposium,

More information

Limiting the Development of Anti-Cancer Drug Resistance in a Spatial Model of Micrometastases

Limiting the Development of Anti-Cancer Drug Resistance in a Spatial Model of Micrometastases Limiting the Development of Anti-Cancer Drug Resistance in a Spatial Model of Micrometastases arxiv:1601.03412v2 [q-bio.to] 2 Mar 2016 September 1, 2018 Ami B. Shah Department of Biology The College of

More information

MEA DISCUSSION PAPERS

MEA DISCUSSION PAPERS Inference Problems under a Special Form of Heteroskedasticity Helmut Farbmacher, Heinrich Kögel 03-2015 MEA DISCUSSION PAPERS mea Amalienstr. 33_D-80799 Munich_Phone+49 89 38602-355_Fax +49 89 38602-390_www.mea.mpisoc.mpg.de

More information

MMSE Interference in Gaussian Channels 1

MMSE Interference in Gaussian Channels 1 MMSE Interference in Gaussian Channels Shlomo Shamai Department of Electrical Engineering Technion - Israel Institute of Technology 202 Information Theory and Applications Workshop 5-0 February, San Diego

More information

Parameter Optimization for FEM based modeling of singlet oxygen during PDT using COMSOL

Parameter Optimization for FEM based modeling of singlet oxygen during PDT using COMSOL Presented at the COMSOL Conference 21 Boston COMSOL Conference, Boston, 1/8/21 Parameter Optimization for FEM based modeling of singlet oxygen during PDT using COMSOL Xing Liang, Ken Kang-Hsin Wang, and

More information

Challenges in Synthesis of Optimal Controls For Problems Arising in Biomedicine

Challenges in Synthesis of Optimal Controls For Problems Arising in Biomedicine Challenges in Synthesis of Optimal Controls For Problems Arising in Biomedicine ITN SADCO meeting January 21-25, 2013, Funchal, Portugal Urszula Ledzewicz Dept. of Mathematics and Statistics Southern Illinois

More information

A Brief Introduction to Bayesian Statistics

A Brief Introduction to Bayesian Statistics A Brief Introduction to Statistics David Kaplan Department of Educational Psychology Methods for Social Policy Research and, Washington, DC 2017 1 / 37 The Reverend Thomas Bayes, 1701 1761 2 / 37 Pierre-Simon

More information

Modeling Consequences of Reduced Vaccination Levels on the Spread of Measles

Modeling Consequences of Reduced Vaccination Levels on the Spread of Measles Bridgewater State University Virtual Commons - Bridgewater State University Honors Program Theses and Projects Undergraduate Honors Program 5-2016 Modeling Consequences of Reduced Vaccination Levels on

More information

Mathematical Modelling of Pulmonary and Extra-pulmonary Tuberculosis

Mathematical Modelling of Pulmonary and Extra-pulmonary Tuberculosis Mathematical Modelling of Pulmonary and xtra-pulmonary Tuberculosis ita H. Shah # and Jyoti Gupta * # Professor, Department of Mathematics, Gujarat University, Ahmedabad, Gujarat, India, Research Scholar,

More information

Local Stability of AIDS Epidemic Model Through Treatment and Vertical Transmission with Time Delay

Local Stability of AIDS Epidemic Model Through Treatment and Vertical Transmission with Time Delay Journal of Physics: Conference Series PAPER OPEN ACCESS Local Stability of AIDS Epidemic Model Through Treatment and Vertical Transmission with Time Delay To cite this article: Cascarilla Novi W and Dwi

More information

E 490 FE Exam Prep. Engineering Probability and Statistics

E 490 FE Exam Prep. Engineering Probability and Statistics E 490 FE Exam Prep Engineering Probability and Statistics Dispersion, Mean, Median, Mode 1. The population standard deviation of the data points 2,1,6 is: (A) 1.00 (B) 1.52 (C) 2.16 (D) 2.56 2. A certain

More information

Reinforcement Learning in Steady-State Cellular Genetic Algorithms

Reinforcement Learning in Steady-State Cellular Genetic Algorithms Reinforcement Learning in Steady-State Cellular Genetic Algorithms Cin-Young Lee and Erik K. Antonsson Abstract A novel cellular genetic algorithm is developed to address the issues of good mate selection.

More information

Abhimanyu Khan, Ronald Peeters. Cognitive hierarchies in adaptive play RM/12/007

Abhimanyu Khan, Ronald Peeters. Cognitive hierarchies in adaptive play RM/12/007 Abhimanyu Khan, Ronald Peeters Cognitive hierarchies in adaptive play RM/12/007 Cognitive hierarchies in adaptive play Abhimanyu Khan Ronald Peeters January 2012 Abstract Inspired by the behavior in repeated

More information

Statistical Tolerance Regions: Theory, Applications and Computation

Statistical Tolerance Regions: Theory, Applications and Computation Statistical Tolerance Regions: Theory, Applications and Computation K. KRISHNAMOORTHY University of Louisiana at Lafayette THOMAS MATHEW University of Maryland Baltimore County Contents List of Tables

More information

Toward Comparison-based Adaptive Operator Selection

Toward Comparison-based Adaptive Operator Selection Toward Comparison-based Adaptive Operator Selection Álvaro Fialho, Marc Schoenauer, Michèle Sebag Orsay, France Evolutionary Algorithms cycle and parameters Performance very sensitive to the parameter

More information

Mathematical Formulation and Numerical Simulation of Bird Flu Infection Process within a Poultry Farm

Mathematical Formulation and Numerical Simulation of Bird Flu Infection Process within a Poultry Farm Mathematical Formulation and Numerical Simulation of Bird Flu Infection Process within a Poultry Farm Arrival Rince Putri, Tertia Delia Nova and M. Watanabe, Graduate School Environmental and Life Science,

More information

Identification of influential proteins in the classical retinoic acid signaling pathway

Identification of influential proteins in the classical retinoic acid signaling pathway Ghaffari and Petzold Theoretical Biology and Medical Modelling (2018) 15:16 https://doi.org/10.1186/s12976-018-0088-7 RESEARCH Open Access Identification of influential proteins in the classical retinoic

More information

Analysis of Vaccine Effects on Post-Infection Endpoints Biostat 578A Lecture 3

Analysis of Vaccine Effects on Post-Infection Endpoints Biostat 578A Lecture 3 Analysis of Vaccine Effects on Post-Infection Endpoints Biostat 578A Lecture 3 Analysis of Vaccine Effects on Post-Infection Endpoints p.1/40 Data Collected in Phase IIb/III Vaccine Trial Longitudinal

More information

Australian Journal of Basic and Applied Sciences. Stability Analysis and Qualitative Behavior of Giving up Smoking Model with Education Campaign

Australian Journal of Basic and Applied Sciences. Stability Analysis and Qualitative Behavior of Giving up Smoking Model with Education Campaign Australian Journal of Basic and Applied Sciences, 9(17) Special 215, Pages: 46-53 ISSN:1991-8178 Australian Journal of Basic and Applied Sciences Journal home page: www.ajbasweb.com Stability Analysis

More information

Bayesian Statistics Estimation of a Single Mean and Variance MCMC Diagnostics and Missing Data

Bayesian Statistics Estimation of a Single Mean and Variance MCMC Diagnostics and Missing Data Bayesian Statistics Estimation of a Single Mean and Variance MCMC Diagnostics and Missing Data Michael Anderson, PhD Hélène Carabin, DVM, PhD Department of Biostatistics and Epidemiology The University

More information

The Time Course of Negative Priming

The Time Course of Negative Priming The Time Course of Negative Priming Hendrik Degering Bernstein Center for Computational Neuroscience Göttingen University of Göttingen, Institute for Nonlinear Dynamics 11.12.2009, Disputation Aging Effects

More information

Application of act and wait control to oscillatory network desynchronization

Application of act and wait control to oscillatory network desynchronization Application of act and wait control to oscillatory network desynchronization Center for Physical Sciences and Technology A. Gostauto 11, LT-01108 Vilnius LITHUANIA XXXIII Dynamics Days Europe, Madrid,

More information

Data Analysis of Manufacturing Test Results for DRAM Module

Data Analysis of Manufacturing Test Results for DRAM Module Data Analysis of Manufacturing Test Results for DRAM Module White Paper March 2018 Copyright 2017 Wiwynn. All rights reserved Contents Introduction... 3 Quality Analysis of DIMMs, SSDs and HDDs by Vendor

More information

INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY

INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY A PATH FOR HORIZING YOUR INNOVATIVE WORK A SLIDING MODE CONTROL ALGORITHM FOR ARTIFICIAL PANCREAS NIRLIPTA RANJAN MOHANTY

More information

Review. Imagine the following table being obtained as a random. Decision Test Diseased Not Diseased Positive TP FP Negative FN TN

Review. Imagine the following table being obtained as a random. Decision Test Diseased Not Diseased Positive TP FP Negative FN TN Outline 1. Review sensitivity and specificity 2. Define an ROC curve 3. Define AUC 4. Non-parametric tests for whether or not the test is informative 5. Introduce the binormal ROC model 6. Discuss non-parametric

More information

Physiologically structured cell population dynamic models with applications to combined drug delivery optimisation in oncology

Physiologically structured cell population dynamic models with applications to combined drug delivery optimisation in oncology Physiologically structured cell population dynamic models with applications to combined drug delivery optimisation in oncology Jean Clairambault, Olivier Fercoq To cite this version: Jean Clairambault,

More information

3.1 The Lacker/Peskin Model of Mammalian Ovulation

3.1 The Lacker/Peskin Model of Mammalian Ovulation Mammalian Ovulation 3.1 The Lacker/Peskin Model of Mammalian Ovulation Reference: Regulation of Ovulation Number in Mammals A large reserve pool of follicles is formed before birth in many female mammals.

More information

Assessing drug synergy in combination therapies. Adam Palmer, Ph.D. Laboratory of Systems Pharmacology Harvard Medical School

Assessing drug synergy in combination therapies. Adam Palmer, Ph.D. Laboratory of Systems Pharmacology Harvard Medical School Assessing drug synergy in combination therapies Adam Palmer, Ph.D. Laboratory of Systems Pharmacology Harvard Medical School Outline Introduction Pharmacokinetic and pharmacodynamic drug interactions Drug

More information

Evolution of pathogens: a within-host approach

Evolution of pathogens: a within-host approach / 52 Evolution of pathogens: a within-host approach Vitaly V. Ganusov Theoretical Biology Utrecht University, Utrecht, The Netherlands Outline Introduction evolution of virulence 2 Evolution of infectious

More information

Optimal dose selection considering both toxicity and activity data; plateau detection for molecularly targeted agents

Optimal dose selection considering both toxicity and activity data; plateau detection for molecularly targeted agents Optimal dose selection considering both toxicity and activity data; plateau detection for molecularly targeted agents Maria-Athina Altzerinakou1,2,3 and Xavier Paoletti3,1,2 1. CESP OncoStat, Inserm, Villejuif,

More information

MODELING GLUCOSE-INSULIN METABOLIC SYSTEM AND INSULIN SECRETORY ULTRADIAN OSCILLATIONS WITH EXPLICIT TIME DELAYS. Yang Kuang

MODELING GLUCOSE-INSULIN METABOLIC SYSTEM AND INSULIN SECRETORY ULTRADIAN OSCILLATIONS WITH EXPLICIT TIME DELAYS. Yang Kuang MODELING GLUCOSE-INSULIN METABOLIC SYSTEM AND INSULIN SECRETORY ULTRADIAN OSCILLATIONS WITH EXPLICIT TIME DELAYS Yang Kuang (joint work with Jiaxu Li and Clinton C. Mason) Department of Mathematics and

More information

Modelling the H1N1 influenza using mathematical and neural network approaches.

Modelling the H1N1 influenza using mathematical and neural network approaches. Biomedical Research 2017; 28 (8): 3711-3715 ISSN 0970-938X www.biomedres.info Modelling the H1N1 influenza using mathematical and neural network approaches. Daphne Lopez 1, Gunasekaran Manogaran 1*, Jagan

More information

The Role of Implicit Motives in Strategic Decision-Making: Computational Models of Motivated Learning and the Evolution of Motivated Agents

The Role of Implicit Motives in Strategic Decision-Making: Computational Models of Motivated Learning and the Evolution of Motivated Agents Games 2015, 6, 604-636; doi:10.3390/g6040604 Article OPEN ACCESS games ISSN 2073-4336 www.mdpi.com/journal/games The Role of Implicit Motives in Strategic Decision-Making: Computational Models of Motivated

More information

Article Epidemic Analysis and Mathematical Modelling of H1N1 (A) with Vaccination

Article Epidemic Analysis and Mathematical Modelling of H1N1 (A) with Vaccination Article Epidemic Analysis and Mathematical Modelling of H1N1 (A) with Vaccination Jagan Mohan Jonnalagadda and Kartheek Gaddam Department of Mathematics, Birla Institute of Technology & Science Pilani,

More information

Analysis of the basic reproduction number from the initial growth phase of the outbreak in diseases caused by vectors

Analysis of the basic reproduction number from the initial growth phase of the outbreak in diseases caused by vectors Analysis of the basic reproduction number from the initial growth phase of the outbreak in diseases caused by vectors University of São Paulo Medicine School rosangelasanches@usp.br November,2013 Introduction

More information

In Silico Tumor Growth: Application to Glioblastomas

In Silico Tumor Growth: Application to Glioblastomas In Silico Tumor Growth: Application to Glioblastomas Olivier Clatz 1, Pierre-Yves Bondiau 1, Hervé Delingette 1,Grégoire Malandain 1, Maxime Sermesant 1, Simon K. Warfield 2, and Nicholas Ayache 1 1 Epidaure

More information

Equation Development of Tumor-Immune ODE System

Equation Development of Tumor-Immune ODE System Equation Development of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 22, 2002 ThisworkwassupportedinpartbyagrantfromtheW.M.KeckFoundation 0-0 TUMOR-IMMUNE EQUATION DEVELOPMENT Overview

More information

Modeling the HIV Transmission with Reinfection

Modeling the HIV Transmission with Reinfection Applied Mathematical Sciences, Vol. 9, 2015, no. 87, 4323-4330 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.52174 Modeling the HIV Transmission with Reinfection Juan C. Castillo Paz,

More information

NEOADJUVANT TRADE-OFFS IN ER+ BREAST CANCER: A GAME THEORETIC APPROACH

NEOADJUVANT TRADE-OFFS IN ER+ BREAST CANCER: A GAME THEORETIC APPROACH Wednesday, 25th, 10:30 NEOADJUVANT TRADE-OFFS IN ER+ BREAST CANCER: A GAME THEORETIC APPROACH Jeffrey West jeffrey.west@moffitt.org H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida,

More information