Mathematical Modeling and Biology Bo Deng Introduction Examples of - - PowerPoint PPT Presentation

mathematical modeling and biology
SMART_READER_LITE
LIVE PREVIEW

Mathematical Modeling and Biology Bo Deng Introduction Examples of - - PowerPoint PPT Presentation

Mathematical Modeling and Biology Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Bo Deng Consistency Model Test Department of Mathematics Mathematical University of Nebraska Lincoln Biology Conclusion


slide-1
SLIDE 1

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical Modeling and Biology

Bo Deng

Department of Mathematics University of Nebraska – Lincoln

March 10, 2016 www.math.unl.edu/∼bdeng1

1 / 24

slide-2
SLIDE 2

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

What is modeling?

Mathematical modeling is to translate nature into mathematics

2 / 24

slide-3
SLIDE 3

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

What is modeling?

Mathematical modeling is to translate nature into mathematics to be logically consistent

2 / 24

slide-4
SLIDE 4

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

What is modeling?

Mathematical modeling is to translate nature into mathematics to be logically consistent to fit the past and to predict future

2 / 24

slide-5
SLIDE 5

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

What is modeling?

Mathematical modeling is to translate nature into mathematics to be logically consistent to fit the past and to predict future to fail against the test of time, i.e. to give way to better models

2 / 24

slide-6
SLIDE 6

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Human history has two periods – before and after calculus (1686/1687)

Issac Newton (1642-1727) is the founding father of mathematical modeling

3 / 24

slide-7
SLIDE 7

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Human history has two periods – before and after calculus (1686/1687)

Issac Newton (1642-1727) is the founding father of mathematical modeling James Clerk Maxwell (1831-1879), Albert Einstein (1879-1955), Erwin Schr¨

  • dinger (1887-1961), Claude

Shannon (1916-2001) are some of the luminary disciples

3 / 24

slide-8
SLIDE 8

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Human history has two periods – before and after calculus (1686/1687)

Issac Newton (1642-1727) is the founding father of mathematical modeling James Clerk Maxwell (1831-1879), Albert Einstein (1879-1955), Erwin Schr¨

  • dinger (1887-1961), Claude

Shannon (1916-2001) are some of the luminary disciples Calculus is the principle language of nature

3 / 24

slide-9
SLIDE 9

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Human history has two periods – before and after calculus (1686/1687)

Issac Newton (1642-1727) is the founding father of mathematical modeling James Clerk Maxwell (1831-1879), Albert Einstein (1879-1955), Erwin Schr¨

  • dinger (1887-1961), Claude

Shannon (1916-2001) are some of the luminary disciples Calculus is the principle language of nature This century is the century of mathematical biology, which is to translate Charles Darwin’s (1809-1882) theory into mathematics

3 / 24

slide-10
SLIDE 10

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model as approximation – Newton’s planetary motion

Sun Planet

  • r1
  • r2
  • r

           m1¨

  • r1

= −Gm1m2

  • r1 −

r2

  • r1 −

r23 m2¨

  • r2

= −Gm1m2

  • r2 −

r1

  • r2 −

r13

  • r

=

  • r1 −

r2

4 / 24

slide-11
SLIDE 11

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model as approximation – Newton’s planetary motion

Sun Planet

  • r1
  • r2
  • r

           m1¨

  • r1

= −Gm1m2

  • r1 −

r2

  • r1 −

r23 m2¨

  • r2

= −Gm1m2

  • r2 −

r1

  • r2 −

r13

  • r

=

  • r1 −

r2 A few calculus maneuvers lead to r(θ) = ρ 1 + ǫ cos θ with the eccentricity 0 ≤ ǫ < 1 for elliptic orbits

4 / 24

slide-12
SLIDE 12

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity – Einstein’s model of space and time

One Assumption:

The speed of light is constant for every stationary observer ¯ x ¯ y v x y K ¯ K

5 / 24

slide-13
SLIDE 13

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity – Einstein’s model of space and time

One Assumption:

The speed of light is constant for every stationary observer ¯ x ¯ y v x y K ¯ K A few calculus maneuvers lead to E = mc2, and more

5 / 24

slide-14
SLIDE 14

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity — Einstein’s model of space and time

One Assumption:

The speed of light is constant for every stationary observer ¯ x ¯ y x y K ct vt L [ √ c2 − v2]t = c¯ t ¯ K

6 / 24

slide-15
SLIDE 15

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity — Einstein’s model of space and time

One Assumption:

The speed of light is constant for every stationary observer ¯ x ¯ y x y K ct vt L [ √ c2 − v2]t = c¯ t ¯ K Prediction: Time dilation for K-frame observer t = L c

  • 1 − (v/c)2 > L

c = ¯ t

6 / 24

slide-16
SLIDE 16

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

General Relativity — Model of space and time in acceleration

x y ¯ x ¯ y c∆t v0∆t c∆t v1∆t v1 = a∆t + v0

7 / 24

slide-17
SLIDE 17

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

General Relativity — Model of space and time in acceleration

x y ¯ x ¯ y c∆t v0∆t c∆t v1∆t v1 = a∆t + v0 Prediction: Light beam bends under acceleration or near massive bodies

7 / 24

slide-18
SLIDE 18

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical model need not be mathematical

Gregor Johann Mendel (1822-1884) found the first mathematical model in biology, leading to the discovery of gene Parent Genotype Offspring Genotype mrr × frr mrD × frD mDD × fDD mrr × frD or mrD × frr mrr × fDD or mDD × frr mrD × fDD or mDD × frD z′

rr

1 1/4 1/2 z′

rD

1/2 1/2 1 1/2 z′

DD

1/4 1 1/2

8 / 24

slide-19
SLIDE 19

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Structure of DNA by modeling

Rosalind Franklin and Maurice Wilkins had the data, but James D. Watson and Francis Crick had the frame of mind to model the data (1953)

9 / 24

slide-20
SLIDE 20

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Another More – Predation in Ecology

The mathematical model was discovered by Crawford Stanley (Buzz) Holling (1930- ) in 1959 Td — average time a predator takes to discover a prey

10 / 24

slide-21
SLIDE 21

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Another More – Predation in Ecology

The mathematical model was discovered by Crawford Stanley (Buzz) Holling (1930- ) in 1959 Td — average time a predator takes to discover a prey Tk — average time a predator takes to kill a prey

10 / 24

slide-22
SLIDE 22

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Another More – Predation in Ecology

The mathematical model was discovered by Crawford Stanley (Buzz) Holling (1930- ) in 1959 Td — average time a predator takes to discover a prey Tk — average time a predator takes to kill a prey Td,k = Td + Tk — average time a predator takes to discovery and kill a prey

10 / 24

slide-23
SLIDE 23

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Another More – Predation in Ecology

The mathematical model was discovered by Crawford Stanley (Buzz) Holling (1930- ) in 1959 Td — average time a predator takes to discover a prey Tk — average time a predator takes to kill a prey Td,k = Td + Tk — average time a predator takes to discovery and kill a prey Rd = 1 Td — rate of discovery, i.e. number of preys a predator would find in a unit time Rk = 1 Tk — rate of killing, i.e. number of preys a predator would kill in a unit time Rd,k = 1 Td,k = 1 Td + Tk — rate of discovery and killing

10 / 24

slide-24
SLIDE 24

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model of Predation in Ecology

And Holling’s predation function form: Rd,k = 1 Td + Tk = 1/Td 1 + Tk(1/Td) = Rd 1 + TkRd

11 / 24

slide-25
SLIDE 25

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model of Predation in Ecology

And Holling’s predation function form: Rd,k = 1 Td + Tk = 1/Td 1 + Tk(1/Td) = Rd 1 + TkRd Prediction: Assume the discovery rate is proportional to the prey population X, Rd = aX. Then the Holling Type II predation rate must saturate as X → ∞ lim

X→∞ Rd,k = lim X→∞

aX 1 + TkaX = 1 Tk X Rd,k

1 Tk

11 / 24

slide-26
SLIDE 26

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Consistency

Not every piece of mathematics can be a physical law or

  • model. Logical consistency is the first and necessary

constraint

Time Invariance Principle (TIP)

A model must has the same functional form for every time independent observation

12 / 24

slide-27
SLIDE 27

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Consistency

Not every piece of mathematics can be a physical law or

  • model. Logical consistency is the first and necessary

constraint

Time Invariance Principle (TIP)

A model must has the same functional form for every time independent observation Newtonian mechanics is TIP-consistent: s t x0 x(s, x0) x(t + s, x0) = x(t, x(s, x0))

12 / 24

slide-28
SLIDE 28

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity is self-consistent

Let P be a point, having K = (x, y, z, t) coordinate in the K-frame and ¯ K = (¯ x, ¯ y, ¯ z, ¯ t) coordinate in the ¯ K-frame. Then they are exchangeable via a linear transformation depending the speed v: ¯ K = KL(v)

13 / 24

slide-29
SLIDE 29

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity is self-consistent

Let P be a point, having K = (x, y, z, t) coordinate in the K-frame and ¯ K = (¯ x, ¯ y, ¯ z, ¯ t) coordinate in the ¯ K-frame. Then they are exchangeable via a linear transformation depending the speed v: ¯ K = KL(v) Let ˜ K = (˜ x, ˜ y, ˜ z, ˜ t) be the coordinate of the same point in a ˜ K-frame moving at speed u with respect to the ¯ K-frame. Then we have ˜ K = ¯ KL(u) = KL(v)L(u) = KL(w) with w = u + v 1 + uv

c2

13 / 24

slide-30
SLIDE 30

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Special Relativity is self-consistent

Let P be a point, having K = (x, y, z, t) coordinate in the K-frame and ¯ K = (¯ x, ¯ y, ¯ z, ¯ t) coordinate in the ¯ K-frame. Then they are exchangeable via a linear transformation depending the speed v: ¯ K = KL(v) Let ˜ K = (˜ x, ˜ y, ˜ z, ˜ t) be the coordinate of the same point in a ˜ K-frame moving at speed u with respect to the ¯ K-frame. Then we have ˜ K = ¯ KL(u) = KL(v)L(u) = KL(w) with w = u + v 1 + uv

c2

The operation u ⊕ v = u + v 1 + uv

c2

for elements u, v ∈ (−c, c) defines a commutative group

13 / 24

slide-31
SLIDE 31

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Holling’s predation model is consistent

Tc — average time to consume a prey

14 / 24

slide-32
SLIDE 32

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Holling’s predation model is consistent

Tc — average time to consume a prey Td,k,c = Td + Tk + Tc — average time to discover, kill, and consume a prey

14 / 24

slide-33
SLIDE 33

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Holling’s predation model is consistent

Tc — average time to consume a prey Td,k,c = Td + Tk + Tc — average time to discover, kill, and consume a prey Then the rate of predation is self-consistent: Rd,k,c = 1 Td,k,c = 1 Td + Tk + Tc = Rd,k 1 + TcRd,k = Rd 1 + (Tk + Tc)Rd

14 / 24

slide-34
SLIDE 34

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Pay the TIP, or else

All differential equation models are TIP-consistent

15 / 24

slide-35
SLIDE 35

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Pay the TIP, or else

All differential equation models are TIP-consistent Most mapping models in ecology are TIP-inconsistent

15 / 24

slide-36
SLIDE 36

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Pay the TIP, or else

All differential equation models are TIP-consistent Most mapping models in ecology are TIP-inconsistent Example: Logistic map xn+1 = Qλ(xn) = λxn(1 − xn) cannot be a model for which n represents time

15 / 24

slide-37
SLIDE 37

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Pay the TIP, or else

All differential equation models are TIP-consistent Most mapping models in ecology are TIP-inconsistent Example: Logistic map xn+1 = Qλ(xn) = λxn(1 − xn) cannot be a model for which n represents time The time n + 2 observation yields a different functional form: xn+2 = Qλ(xn+1) = Qλ(Qλ(xn)) = Qµ(xn) for any value µ. Strike one on the logistic map

15 / 24

slide-38
SLIDE 38

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Finding the Best Fit

¯ x1, . . . , ¯ xn — Observed states at time t1, . . . , tn for a natural process which are modeled by competing models y(t; y0, p) and z(t; z0, q), respectively, with parameter p, q, and initial state y0, z0

16 / 24

slide-39
SLIDE 39

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Finding the Best Fit

¯ x1, . . . , ¯ xn — Observed states at time t1, . . . , tn for a natural process which are modeled by competing models y(t; y0, p) and z(t; z0, q), respectively, with parameter p, q, and initial state y0, z0 Model selection criterion: All else being equal whichever has a smaller error is the benchmark model by default: Ey = min

(y0,p) n

  • i=1

[y(ti; y0, p) − ¯ xi]2 Ez = min

(z0,q) n

  • i=1

[z(ti; z0, q) − ¯ xi]2

16 / 24

slide-40
SLIDE 40

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Finding the Best Fit

¯ x1, . . . , ¯ xn — Observed states at time t1, . . . , tn for a natural process which are modeled by competing models y(t; y0, p) and z(t; z0, q), respectively, with parameter p, q, and initial state y0, z0 Model selection criterion: All else being equal whichever has a smaller error is the benchmark model by default: Ey = min

(y0,p) n

  • i=1

[y(ti; y0, p) − ¯ xi]2 Ez = min

(z0,q) n

  • i=1

[z(ti; z0, q) − ¯ xi]2 A parameter value is only meaningful to its model, and it can only be derived by best-fitting the observed data to the model

16 / 24

slide-41
SLIDE 41

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Fit the past, predict the future

Edmond Halley (1656-1742) used Newtonian mechanics to predict the 1758 return of Halley’s Comet, giving the comet its name

17 / 24

slide-42
SLIDE 42

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Fit the past, predict the future

Edmond Halley (1656-1742) used Newtonian mechanics to predict the 1758 return of Halley’s Comet, giving the comet its name Arthur Eddington (1882-1944) used the total solar eclipse

  • f May 29, 1919 to confirm general relativity’s prediction

for the bending of starlight by the Sun, making Einstein an instant world celebrity

17 / 24

slide-43
SLIDE 43

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Fit the past, predict the future

Edmond Halley (1656-1742) used Newtonian mechanics to predict the 1758 return of Halley’s Comet, giving the comet its name Arthur Eddington (1882-1944) used the total solar eclipse

  • f May 29, 1919 to confirm general relativity’s prediction

for the bending of starlight by the Sun, making Einstein an instant world celebrity Gregor Mendel’s Laws of Inheritance (1866) was rediscovered in 1900, ushering in the science of modern genetics

17 / 24

slide-44
SLIDE 44

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Model Test – Fit the past, predict the future

Edmond Halley (1656-1742) used Newtonian mechanics to predict the 1758 return of Halley’s Comet, giving the comet its name Arthur Eddington (1882-1944) used the total solar eclipse

  • f May 29, 1919 to confirm general relativity’s prediction

for the bending of starlight by the Sun, making Einstein an instant world celebrity Gregor Mendel’s Laws of Inheritance (1866) was rediscovered in 1900, ushering in the science of modern genetics Holling’s model of predation is ubiquitous in theoretical ecology

17 / 24

slide-45
SLIDE 45

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical Biology — To Translate Evolution to Mathematics

Example: One Life Rule

Every organism lives only once and must die in any finite time in the presence of infinite population density

18 / 24

slide-46
SLIDE 46

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical Biology — To Translate Evolution to Mathematics

Example: One Life Rule

Every organism lives only once and must die in any finite time in the presence of infinite population density In math translation: Let xt be the population at time t. Then the per-capita change must satisfy xt − x0 x0 = xt x0 − 1 ≥ −1

18 / 24

slide-47
SLIDE 47

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical Biology — To Translate Evolution to Mathematics

Example: One Life Rule

Every organism lives only once and must die in any finite time in the presence of infinite population density In math translation: Let xt be the population at time t. Then the per-capita change must satisfy xt − x0 x0 = xt x0 − 1 ≥ −1 Lead to One Life Rule ⇐ ⇒ lim

x0→∞

xt − x0 x0 = −1 and to the logistic equation ˙ x(t) = rx(t)[1 − x(t)/K] with x(t) = xt, r the max per-capita growth rate, and K the carrying capacity

18 / 24

slide-48
SLIDE 48

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Footnote: model or no model, generalization or relativism is often the problem

Strike two on the logistic map x1 = λx0(1 − x0): lim

x0→∞

x1 − x0 x0 = lim

x0→∞[λ(1 − x0) − 1] = −∞ = −1

19 / 24

slide-49
SLIDE 49

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Footnote: model or no model, generalization or relativism is often the problem

Strike two on the logistic map x1 = λx0(1 − x0): lim

x0→∞

x1 − x0 x0 = lim

x0→∞[λ(1 − x0) − 1] = −∞ = −1

While the logistic equation, x′(t) = rx(t)(1 − x(t)/K), dogged another consistency bullet lim

x0→∞

x(t; x0) − x0 x0 = lim

x0→∞

  • K

x0 + (K − x0)e−rt − 1

  • = −1

19 / 24

slide-50
SLIDE 50

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Footnote: model or no model, generalization or relativism is often the problem

Strike two on the logistic map x1 = λx0(1 − x0): lim

x0→∞

x1 − x0 x0 = lim

x0→∞[λ(1 − x0) − 1] = −∞ = −1

While the logistic equation, x′(t) = rx(t)(1 − x(t)/K), dogged another consistency bullet lim

x0→∞

x(t; x0) − x0 x0 = lim

x0→∞

  • K

x0 + (K − x0)e−rt − 1

  • = −1

There should be no different versions of the same reality, but refined approximations

19 / 24

slide-51
SLIDE 51

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

20 / 24

slide-52
SLIDE 52

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

The AT pair has one weak O-H bond but the GC pair has two O-H bonds. Hence, the GC pair takes longer to complete binding than the AT pair does

20 / 24

slide-53
SLIDE 53

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

21 / 24

slide-54
SLIDE 54

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

Start with a conceptual model: DNA replication is a communication channel

21 / 24

slide-55
SLIDE 55

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

Start with a conceptual model: DNA replication is a communication channel Every communication is characterized by the transmission data rate in bits per second, i.e. the information entropy per second

21 / 24

slide-56
SLIDE 56

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

For 2n paired bases, the replication rate is R2n = log2(2n)

τ12+τ34+···+τ(2n−1)(2n) n

in bits per time

22 / 24

slide-57
SLIDE 57

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

For 2n paired bases, the replication rate is R2n = log2(2n)

τ12+τ34+···+τ(2n−1)(2n) n

in bits per time If 5

3 ≤ τGC τAT ≤ 2.7, then:

max

n

R2n = R4

22 / 24

slide-58
SLIDE 58

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

One More Example: Why DNA is coded in 4 bases?

For 2n paired bases, the replication rate is R2n = log2(2n)

τ12+τ34+···+τ(2n−1)(2n) n

in bits per time If 5

3 ≤ τGC τAT ≤ 2.7, then:

max

n

R2n = R4 Punch Line: Life is a reality show on your DNA channel

22 / 24

slide-59
SLIDE 59

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Closing Comments

Mathematics is driven by open problems, but science is driven by existing solutions Mathematical modeling is to find the equation to which nature fits as a solution Mathematics is to create more hays but modeling is to find the needle in haystack Mathematical biology is not to solve mathematical problems of models but to find mathematical models for biological problems Training to be a mathematical modeler does need to solve mathematical problems of reasonable models.

23 / 24

slide-60
SLIDE 60

Mathematical Modeling and Biology Bo Deng Introduction Examples of Models Consistency Model Test Mathematical Biology Conclusion

Mathematical modeling is to construct the picture so that the consequence of the picture is the picture of the consequence. – Anonymous or by Heinrich Hertz (1857-1894)

24 / 24