The Structure, Function, and Evolution of Biological Systems

1. The Structure, Function, and Evolution of Biological Systems Instructor: Van Savage Spring 2010 Quarter 4/27/2010

2. Outline Finite size populations and genetic drift Coalescence Understanding directional and nondirectional forces 3. General Diffusion Equation 4. Biology a. Bacterium’s use of 3 (physical constraint) b. Population genetics--combining selection and drift in evolution (Conceptual analogy) Derive equations in this context

3. Genetic Drift--Nondirectional force Drift--by random, transient, non-genetic events, individuals that would be highly reproductive in a “repeated” experiment, are lost to death. Imagine a population of identical individuals that are being chosen for mating. Random chance that a given individual will be chosen.

4. Genetic Drift Model this by randomly sampling from entire population (Wright-Fisher model). Population size, N, is constant, and individuals are randomly selected for mating Each individual has 1/N chance of reproducing. We get a binomial tree that depends on frequency, p, and total population size, N. Generation 0

5. Genetic Drift frequency, p 1 time So, rate of spread of the width of distribution is~p(1-p)/2N

6. Coalescence: Look backwards in time

7. Quick results • P(fixation)=p • P(fixation new individual mutant)=1/2N • Total fixation rate of mutations=2Nμ(1/2N)=μ • Probability of coalescing t+1 generations back is (1/2N)e-t/2N

8. Using last result • <T>=2N • StDev(T)=2N

9. General diffusion equationandcombining selection and drift

10. Directional Forces: Go with the flow

11. t1 t2 t3 t4 t0 Directional Forces Position Crowd running towards a celebrity or away from a fire. Pushing or rolling any ball or object A river flowing towards the sea or ocean

12. Nondirectional forces: No flow

13. v Nondirectional forces Lost in a crowded intersection Drop of dye in water Smoke Choosing each step after flipping a coin

14. equilibrium boundary(wall) Net Flow--Directional Forces f f(t,i) is abundance of or probability of being in bin i at time, t. v(i) is speed of flow out of bin i. x i-1 i i+1 Net Flow=Flow In-Flow Out =Flow from left(i-1->i) - Flow to right(i->i+1) ->f(t+1,i)-f(t,i)=v(i-1)f(i-1,t)-v(i)f(i,t) Continuum limit: ->df/dt=-d(vf)/dx (i.e., distance=velocity*time)

15. Net Flow--Nondirectional Forces f Equilibrium (D=constant) x i-1 i i+1 Net Flow=Flow In-Flow Out =Flow from left (i-1->i)+Flow from right (i+1->i) -Flow to right (i->i-1)-Flow to left (i->i+1) ->f(t+1,i)-f(t,i)=D(i-1)f(i-1,t)+D(i+1)f(i,t)-2*D(i)f(i,t) D(i) is the diffusion rate Continuum limit: ->df/dt=d2(Df)/dx2 (Second derivative) Local process and affects width of distribution, not mean

16. Global signature of diffusion Random walk x(t+1)=x(t)±1 ->x2(t+1)=x2(t)+2x(t)+1(1/2 of time) =x2(t)-2x(t)+1(1/2 of time) <x2(t+1)>=(1/2)* <(x2(t)+2x(t)+1)> +(1/2)*<(x2(t)-2x(t)+1)> =<x2(t)>+1=<x2(t-1)>+2 Iterating this gives: <x2(t+1)>=Number of time steps~t

17. Diffusion properties and simulation http://web.mit.edu/course/3/3.091/www/diffusion/ Internal dynamics of diffusion are not immediately obvious from global behaviors, whereas it is for directional forces Nondirectional force (diffusion) affects width of distribution, and directional forces affect the mean. Diffusion process is NOT time reversible. Initial conditions are forgotten.

18. 2. Combined Effects Person trying to walk north (directional) through a busy intersection (nondirectional) Net Flow=Directional Flow+Nondirectional Flow Diffusion Equation (Also known as Kolmogorov forward equation)

19. Often, v and D are constant, so:

20. 3. Physics--Brownian Motion Molecule in glass of water is analogous to our person walking through a crowd. Since molecule is so small (mass is so little), gravity’s effect (directional force) is negligible. Hence, (Heat Equation)

21. 4. Biology—Magnetotactic Bacteria

22. Magnetotactic Bacteria Better conditions at bottom (oxygen pressure) Weigh so little that gravity is negligible, and they do not know which way is down. They have internalized enough magnetite particles so that earth’s magnetic field can just overcome nondirectional forces of Brownian Motion. Since magnetic fields go into earth, they can now sense down. So, they “solve” problem from previous slide! Bacteria in north and south are polarized differently.

23. Apply magnetic field to diffusing particles This gives a directional force, and since magnetic force is much stronger than gravity, this is not negligible. Must return to full equation. v depends on strength of magnetic field

24. Types of multidisciplinary influences Physical Process Conceptual and Mathematical Analogy Magnetotactic bacteria Evolution and Population Genetics -> Combine natural selection and genetic drift

25. 4. B. Combining Selection and Drift We can also understand process of evolution by means of diffusion equation. Requires different sort of extension to biology. It’s not just understanding how biological organisms use and are constrained by physics, but it’s using analogies to mathematical physics to understand biological problems.

26. Selection--Directional Force Let a population (wild type) suddenly have a few individuals with a mutation that forms a new allele. If fitness (as measured by growth rate--number of offspring per individual per generation that survive to next generation) of wild type is normalized to 1, and mutants have fitness 1+s Population Size Mutant Wild type time

27. If population size is fixed (finite resources), only a matter of time, until mutant takes over (fixation). frequency of mutants, p 1 time, t Position space, x, is replaced by frequency space, p, for frequency of mutants. Velocity of selection force is~ p(1-p)s

28. Genetic Drift frequency, p 1 time So, rate of spread of the width of distribution is~p(1-p)/N

29. Strong Analogy 1. Selection Gravity Pushed by marathon 2. Drift Brownian motion Crowd at intersection 3. Small organisms Small populations (Brownian>>Gravity) (Drift>>Selection) 4. Large organisms Large populations (Gravity>>Brownian) (Selection>>Drift)

30. Equation for Population Genetics p0 is initial frequency of mutants in the population. Questions we can answer using this equation. 1. Is mutant population likely to go extinct or take over population (fixation)? 2. How long does it take before extinction or fixation occurs? 3. For a given N and s, how large does p0 need to be before mutants are likely to take over.

31. More proper derivation Probability of moving from p-ε to p probability density of having frequency p at time t+dt Taylor expand in p around epsilon to get Kolmogorov forward equations

32. More proper derivation Looking backward in time, as or coalescence gives Kolmogorov backward equation Sign of directional term flips because now going backwards in times and is time reversible. Non-directional term does not flip sign because non-reversible.

33. Probability of Fixation Solve equation at and impose boundary condition for p0=0 and p0=1. Probability of Fixation of mutants

34. Investigate some limits Large population, strong selection: e-2Nsp<<1 -> u(p0)~1 (guaranteed to fix) 2. Large population, really weak selection such that 2Nsp0<<1: e-2Nsp~1-2Nsp0-> u(p0)~2Nsp0 When one mutant,p0=1/N,and u(1/N)~2s (fixation probability increases linearly with s) 3. Under very weak selection (s->0): e-2Nsp~1-2Nsp0-> u(p0)~ p0, When one mutant, p0=1/N,and u(1/N)~1/N (same as for pure drift) All limits check out.

35. 5. Economics Black-Scholes model Price of stock is like position space (physics) or frequency space (population genetics). Directional force--general increase in worth of the market, represented by interest rate. Nondirectional force--random forces in market. Individual stocks or groups of stock will wander randomly in price. (major insight of this model! Also, because it shows value of volatility and how to make money from it.)

36. Assumptions of Black-Scholes • Price follows Brownian motion • It is possible to short sell stock (options) • No arbitrage is possible (no asymmetry of which to take advantage) • Trading is continuous • No transactions costs or taxes • Stock’s price is continuous and can be arbitrarily small • Risk-free interest rate is constant

37. Results from Black-Scholes Provides method for calculating fair cost of an option. Provides method for hedging “bets” and getting risk-free investment that allows one to make money according to the overall growth of the market.

38. Black-Scholes PDE S-stock price V-option cost r-interest rate -variance of random process

39. Impact of this work One equation, similar concepts, applications to multiple fields with its own set of insights • Fourier developed the heat equation • Brownian motion, Einstein’s greatest achievement? • Applied in cosmology, particle physics, etc. • Huge advance in population genetics, used to study molecular motors, and lots of intracellular processes • 1997 Nobel prize in economics for Black-Scholes

40. Conclusions • Diffusion equations describe directional and nondirectional forces. (Could also have forces on higher-order moments by extending this.) 2. Because of generality of 1, we can apply them to many different types of problems in many different fields. (Multidisciplinary) 3. Diffusion equations have already proved very useful in physics, biology, and economics. 4. Examples of two types of multidisciplinary science: a. Results from one field directly place constraints on or are utilized by agents in the other field. b. By the correct choice of analogy between fields, mathematical treatments and results can be used to draw new conclusions and insights within another field.