This website provides video tutorials to help interdisciplinary scientists understand biological systems using mathematical reasoning, physical creativity, and visual art. Use this curriculum to prepare for introductory courses in undergraduate math, science, and systems biology. These videos and slides are open course ware made available under a Creative Commons license so you can use them to develop customized lessons. If you find these tutorials helpful, these external resources and textbooks might interest you as well.

Pre-algebra, algebra, geometry, and pre-calculus

We want to emphasize the following ideas about numbers and functions. (1) The number line is a sequenced collection of distinct addresses like the arrangment of house numbers along a street. (2) Infinity is not a number. (3) A variable, such as "x," is a stand in for a number from a set of interest, at once arbitrary in the sense that x could be any number from that set, but at the same time specific and particular in that x is considered at any one time only one number. (4) f(x) is the single result of applying function f to a particular x. f(x) can be plotted for a collection of values of x, as well as expressed in writing as an "association rule" using arithmetic and other operations on variable x. (5) Functions can be composed. (6) Plots of inverse functions are mirror images reflected about the y = x diagonal. If you are familiar with these ideas you can probably comfortably skip to the first set of blue videos.

Topic Slides Video Description
"What is a number?" PPTX slides MP4 video Street numbers, money in bank accounts, points on number lines, quantum particles
as contrasted with distinguishable manipulatives
For our purposes, infinity is not a number
Algebra PPTX slides MP4 video Variables: At once arbitrary, yet specific and particular (a.s.a.p.)
Functions, composition, and inverse
f(x) is not a function. f(x) is not the function f. f(x) is the particular value associated, in the big picture of the function f, with x, a number that is at once arbitrary, yet particular and specific
Inverse functions do not always exist
First glimpse at the complex plane and i := √ -1 

Additional activity for this module

The following four videos are required for describing simple dynamic systems. As we will see much later, solving for the dynamics of systems with more than one moving part often reduces to solving for the roots of a polynomial expression. The simplest toy models solved this way involve polynomials of order 2. This is why we must understand how to solve a quadratic equation. We need trigonometry so we can use sinusoidal functions to describe oscillations. The videos on summation and combinatorics provide the underpinnings of the binomial theorem, which we later use to talk about the exponential function ex, which appears when dynamics involve repetitive proportional scalings.

Topic Slides Video Description
Quadratic equation PPTX slides MP4 video Linear combination of terms in a polynomial
Zeroes or "roots" of a function
Completing the square
Euclidean geometry and trigonometry PPTX slides MP4 video Flat space, curved space, non-embedded curved space
Pythagorean theorem (ca. 300 BCE)
The unit-radius circle, the unit-hypotenuse triangle, jya-ardha (sine), koti-jya (cosine) (ca. 510)
Geometric construction for approximating π
Angle-sum identities

Additional activity for this module

Summation PPTX slides MP4 video Geometric series
Harmonic series; sums do not always exist
Gauss summation trick
Enumerative combinatorics PPTX slides MP4 video Permutations and factorials
Combinations
Binomial theorem
Small parameter expansion

Calculus (ca. 1700s)

In this unit we define constructions in calculus (derivatives and integrals) so that we can eventually identify them, using explicit and concrete figures and words, with the processes of (1) describing the kinetic processes and (2) inferring long-term time-course outcomes for collections of chemically reacting biological molecules. This line of discussion presents differential equation models of biology as potentially useful approximations. Biochemistry is not always idiomatic for calculus; in fact, cellular biochemistry is an example of where applying calculus can be very dicey.

Topic Slides Video Description
Limits PPTX slides MP4 video ε-δ definition of limit, notion of "arbitrarily close"
Example of calculating a limit
Limits do not always exist
Differentiation PPTX slides MP4 video Slopes and derivatives
Derivatives do not always exist
Common derivatives: power law, chain rule, product rule, quotient rule, sine, and cosine
Infinitessimals describe processes; in this course they are not numbers
Partial differentiation PPTX slides MP4 video When a function depends on multiple independent variables, the ∂ denotes slopes calculated by jiggling only one independent variable at a time
Taylor expansions PPTX slides MP4 video Second derivative and curvature
Local approximations and Taylor series
(Successful) power-series representations do not always exist
Power-series expansion of sine and cosine, iterative calculation of π
Taylor expansions II l'Hôpital's rule ("0/0" version)
Newton-Raphson method for approximating zeros of a function
Integration PPTX slides MP4 video Anti-derivatives, Riemann sums, and integrals
Example kosher calculation of a simple integral
Deductive inference of integral by definition as anti-derivative
"Backwards chain rule"--u substitution
Integration II Integration by parts
Integration III Integrals do not always exist
Triangle/trigonometric substitution
Partial fractions
Separation of variables PPTX slides MP4 video Two wrongs make a right
Tear two differentials apart as though they retained meaning in isolation
Slap on the smooth S integral sign as though it were a unit of meaning itself, even without a differential
You get the same integral expression you would obtain long-hand using u-substitution or "change of variables" in integrals

Applications of calculus

Topic Slides Video Description
Euler's number I PPTX slides MP4 video Euler's number 1a: Compound interest
Compounding interest with arbitrarily short compounding periods
Power series representation of ex
MP4 video Euler's number 1b: e to the zero
e0 = 1
MP4 video Euler's number 1c: Exponent multiplication identity
(ex)p = epx
MP4 video Euler's number 1d: Exponent addition identity
exey = ex+y
MP4 video Euler's number 1e: Andrew Jackson
Mnemonic for memorizing e = 2.718281828459045...
MP4 video Euler's number 1f: Natural logarithm
The natural logarithm is the inverse of the exponential
ln(ex) = eln(x) = 1
MP4 video Euler's number 1g: Integral of 1/x
∫(1/x)dx = ln(x) + C
Projects day Make an Euler's number chart
Make a slide rule
Stochasticity PPTX slides MP4 video (This section is a prerequisite for "protein dynamics 101"). Many of the homework problems from undergraduate calculus and differential equations involve notions of stochasticity.
* For-real stochasticity: Fundamental indeterminism
* Fake stochasticity: Periodic, deterministic hidden variables
* Fake stochasticity: Aperiodic, deterministic (chaos)
Markov models

Additional activity for this module

Protein dynamics 101 PPTX slides MP4 video This is a canonical worked problem from introductory systems biology. We will explain one way to fantasize about the classic protein dynamics equation dx/dt = β - αx and analytically demonstrate that protein "rise time" depends on degradation rate only.
Additional activity: See textbook presentation by Alon, An Introduction to Systems Biology: Design Principles of Biological Circuits, Boca Raton: Chapman & Hall/CRC, 2007 (p. 18-22).
Mass action PPTX slides MP4 video Mass action a: Law of mass action
Collision picture
MP4 video Mass action b: Cooperativity
Cooperativity of the simple kind and Hill functions
MP4 video Mass action c: Bistability
Combining molecular production rates with nonlinear dose-dependence with unimolecular degradation can generate systems with multiple stable steady states
Additional activity: See textbook presentation by Alon, An Introduction to Systems Biology: Design Principles of Biological Circuits, Boca Raton: Chapman & Hall/CRC, 2007 (sections 2.3-2.3.4, p. 7-16).
Evolutionary game theory I PPTX slides MP4 video Collisional population dynamics and tabular game theory
An outcome of the prisoner's dilemma is simultaneous survival of the relatively most fit with decrease in overall fitness

Additional activity: Access McKenzie, A.J., "Evolutionary Game Theory", The Stanford Encyclopedia of Philosophy (Fall 2009 Edition), Zalta, E.N. (ed.) (online) and compare the replicator dynamics described there with the collisional population dynamics in this tutorial. Watch Deborah Gordon talk about colony expansion, task allocation, and organization without central control in ant colonies (TED-talk video online).

Probability, statistics, and stochastic processes

Topic Slides Video Description
Statistics 101 PPTX slides MP4 video Distributions, averages, variances, useful identities
Statistical independence
Routinely-exploited expressions: Covariances vanish and variances of sums are sums of variances
Probability 101 PPTX slides MP4 video Probability 101a: Binomial distribution
Bernoulli (ca. 1700s) coin-toss process
Independent events
PPTX slides MP4 video Probability 101b: Poisson limit
Independent + "rare" events
PPTX slides MP4 video Probability 101c: Stirling's approximation
Comparison with integral of natural logarithm
PPTX slides MP4 video Probability 101d: Central limit theorem
Many independent events
Binomial distribution in limit of many coin tosses
Gaussian distribution
PPTX slides MP4 video Probability 101e: Stories that relate central limit theorem to physics and biology
Physics: Taylor expansion in the context of tightly-controlled, narrow instrument noise
Biology: Logarithm of product of fluctuating factors: Log-normal distributions
Uncertainty propagation PPTX slides MP4 video Uncertainty propagation a: Quadrature
Quadrature formula is a result of Taylor expanding functions of multiple fluctuating variables, assuming that fluctuations are independent, and then applying the identity "variances of sums are sums of variances"
PPTX slides MP4 video Uncertainty propagation b: Sample estimates
Standard deviation vs. sample standard deviation
Mean vs. sample mean
Standard deviation of the mean vs. standard error of the mean
Rule of thumb for thinking about whether error bars overlap
MP4 video Uncertainty propagation c: Illusory sample size
"I quantitated staining intensity for 1 million cells from 5 patients, everything I measure is statistically significant!" It is quite possible that you need to use n = 5, instead of 5 million, for the √ n  factor in the standard error.
PPTX slides MP4 video Uncertainty propagation d: Sample variance curve fitting
Reduced chi-square χ2 fitting

Do not assume that parameter fit uncertainties from black-box software packages are appropriate to interpret in a "covariance = zero" context (Gutenkunst, Sorger)

Additional resource: Web page on data fitting from the Harvey Mudd College physics kiosk (online)

Describing heterogeneity Why on earth have we multiple measurements of distribution heterogeneity?
Some measures increase, some flatline, and some decrease with "effective number" size
  • Standard deviation
  • Coefficient of variation
  • Fano factor
  • Entropy, Shannon information, diversity

Stochastic dynamics

Topic Slides Video Description
Basic stochastic simulation PPTX slides MP4 video Basic stochastic simulation a: Master equation
MP4 video Basic stochastic simulation b: Stochastic simulation algorithm
Derivation of exponential distribution of waiting times
Luria-Delbrück fluctuation analysis P0 method
Numerical stochastic simulation PPTX slides Pseudo-random number generators
Example stochastic birth-death process script
Genetic drift Moran birth-death model
Langevin integration method PPTX slides Stochastic differential equations
Gaussian white noise
Path integrals The probability of ending up at a final condition starting out from an initial condition can be written as a sum of the probabilities of distinct ways to travel between the two conditions.

Linear algebra

Topic Slides Video Description
LA I PPTX slides MP4 video LA 1a: Teaser
Motivating example: Modeling dynamics of web start-up company customer base
MP4 video LA 1b: Vectors
Vectors, vector spaces, and coordinate systems
MP4 video LA 1c: Operators
Linear operators, matrix representation, matrix multiplication
MP4 video LA 1d: Solution of teaser
Using eigenvalue-eigenvector analysis to solve for the dynamics of the demographics of the web-startup customer base
Additional activity for this module
Quasispecies PPTX slides MP4 video Simple quasispecies eigendemographics and eigenrates
Additional activity: Read the green box on p. 0454 from Bull, Meyers, and Lachmann, "Quasispecies made simple," PLoS Comp Biol, 1(6):e61 (2005) (online).
Euler's number II PPTX slides MP4 video Euler's formula: Expanding the exponential function in terms of sine and cosine
Complex exponentials in the complex plane
Euler's identity e = -1
LA II PPTX slides MP4 video Rotation matrix
Complex eigenvectors and eigenvalues
LA III Orthogonality, dot product, and length
Determinants
Invertibility

Differential equations

Topic Slides Video Description
DEs I PPTX slides MP4 video Direction fields, quiver plots, and integral curves
Numerical integration of systems of differential equations
DEs II Integrating factor method
Solution by power-series expansion
DEs III PPTX slides MP4 video DEs IIIa: Transcription-translation
Canonical mRNA-protein system from systems biology 101
Additional activity: See textbook presentation by Alon, An Introduction to Systems Biology: Design Principles of Biological Circuits, Boca Raton: Chapman & Hall/CRC, 2007 (problem 2.2, p. 23).
MP4 video DEs IIIb: Eigenvector-eigenvalue analysis
Determine the directions of "unbending" trajectories for a more precise hand sketch of the phase portrait
MP4 video DEs IIIc: The cribsheet of linear stability analysis
Use eigenvalue-eigenvector analysis to find analytic solutions for linear systems and describe the qualitative features of trajectories approaching, side-swiping, or departing from steady state.
DEs IV PPTX slides MP4 video DEs IVa: Adaptation
Adaptation is not absence of change; instead it is the presence of eventually compensatory changes
Additional activity: Read Ma, Trusina, El-Samad, Lim, and Tang, "Defining network topologies that can achieve biochemical adaptation," Cell, 138: 760-773 (2008) (online).
MP4 video DEs IVb: Cribsheet of almost linear stability analysis
Linear analysis of nonlinear systems
Local linearization: Jacobian
Additional activity: See "Sketching non-linear systems," from Differential Equations: Unit IV First-Order Systems, MIT OpenCourseWare (online) and Harris, K., "Perturbations in linear systems (2008 November 12)," Math 216: Differential Equations, University of Michigan (online).
DEs V PPTX slides MP4 video Heuristic picture of oscillations in 2-d
Intuitive introduction to 2-d oscillations (Romeo and Juliet)
Twisting nullclines
Time-delays
Stochastic resonance
Additional activity: You may skim Ferrell, Jr., Tsai, and Yang, "Modeling the cell cycle: Why do certain circuits oscillate?" Cell, 144: 874-885 (2011)(online). Comment on how the positive-feedback term in Eqtn. 25 (pg. 882) contributes to the difference between the phase portraits in Fig. 4B (pg. 878) and Fig. 8B (pg. 883). The article describes the positive-feedback in terms of a time delay. Please describe the contribution of the positive-feedback term to stable oscillations instead in terms of "twisting nullclines" from the video tutorial.
Additional homework TBA

Spatially-resolved models, cellular automata, and preview of agent-based models

Topic Slides Video Description
Develop expectations Seeing what computers can do
In this activity, you will play against the computer in Blizzard's StarCraft for 2 hrs and in Sid Meier's Civilization for 2 hrs. WARNING: This activity might require rehabilitation and video game addiction treatment (PubMed).
PPTX slides MP4 video Cellular automata a
Deterministic cellular automata
In this video, we see that limiting dispersal of seeds of annual plants can increase the proportion of the copper-colored subpopulation, whereas thorough mixing instead allows the denim plant subpopulation to dominate quickly.

Additional activities: Refer to a similar model in Nowak and May, "Evolutionary games and spatial chaos," Nature 359:826-829 (1992) (online). Watch Athena Aktipis talk about the walk-away model, which can contribute to the evolution of cooperation in highly-mobile populations (University of California, Los Angeles, Center for Behavior, Evolution, and Culture 2009, 1-hr video online)

Cellular automata b
Stochastic cellular automata
Toy agent-based model
You will program a simple ABM
For more extensive discussion, see Athena Aktipis's page on agent-based modeling (online).
Fast-Fourier transform
Efficient computation of local linear interactions FFT convolution trick

Newtonian physics

Topic Slides Video Description
Vectors Review of vector addition and vector calculus
Momentum Velocity v := dx/dt
Momentum p := mv
Spatial translational invariance and momentum conservation
Momentum exchange for an isolated pair of bodies
Force Superposition of pairwise momentum exchange
FNET = ma
Newtonian gravitation r2-law
Approximate uniform gravitational acceleration for small relative changes in radii
Uniformly accelerated motion
Non-fundamental forces Normal force
Tension
Kinetic and static friction

Simple-harmonic-oscillation (pendulum)

Lab: Measuring the gravitational field with a pendulum

Elastic materials (Hookean)
Elastic modulus
Propagation of vibrations on a taut string
Standing waves
Lab: Measure the speed of sound using a flute
Lab: Build a graduated cylinder water harmonica (equal-tempered tuning)

Derivation of the shape of a suspension bridge
Lab: Calculate the linear mass-density of the Golden Gate bridge

Momentum conservation
Work and energy conservation

Newton's law of gravitation
Rocket science
Electricity and magnetism

Modern physics

Topic Slides Video Description

Light has wave-like properties

Lab: Location of phantom laser spots seen while wearing diffraction-grating glasses

Postulates of quantum mechanics
Heisenberg uncertainty relationship (not a principle, but a result): You postulate that eigenvalues of Hermitian operators represent experimental measurement outcomes; you define some operators in terms of other operators in ways such that they necessarily cannot share eigenvectors, and then you are surprised that dispersion in some observables can only be squished at the expense of precision in other observables
Ehrenfest theorem and correspondence principle, Newtonian physics from quantum mechanics
Particle in a box

Hydrogen atom Single-electron solution
Multi-electron atom Taylor-expansion of interaction potential
See also: http://arxiv.org/abs/physics/0407126
Multi-electron molecules

Symmetries, conserved quantities, and extremized quantities
Path-integral formulation of quantum mechanics and wave optics, extremizing functions of path
Entanglement, Einstein-Podolsky-Rosen paradox, and Bell's inequalities

General relativity

Sine-Gordon equation and solitons

Statistical physics

Topic Slides Video Description
Modern statistical mechanics Heuristic justification of equal a priori probability postulate
Second "law" of thermodynamics
Boltzmann distribution PPTX slides Ways of sharing energy between a small system and a large reservoir
Entropy
Free energy
Ising model Phase transition
Hysteresis
Bose-Einstein statistics Bose-Einstein systems obey Boltzmann distributions! A quantum index need not be interpreted as referring to a collection of distinct particles.
Lagrange multiplier

Blackbody radiation
Diffusion, Brownian motion, and fluctuation-dissipation (Smoluchowski treatment)
Self-organized criticality
Chaos is not a lack of determinism

Biological physics and biophysics

Topic Slides Video Description

Optical coherence microscopy
Fluorescence correlation spectroscopy and auto-correlation function
Lateral displacement mechanical flow cell sorter (Loutherback-Sturm)

Introduction to population dynamics topics PPTX slides See table below
Fancy name Descriptive name Example insight
Fixation time of mutants Allele proportions fluctuate stochastically in finite populations of true-breeding replicators Rapid loss of heterogeneity is possible in small populations
Quasispecies Interplay between dispersive forces that generate genotypic heterogeneity and populating forces that favor replication of fit genotypes determine whether populations expand, maintain heterogeneity, or disperse throughout genotypic space and lose macroscopic population numbers Fitness, in the sense of rapid net production of progeny, is not sufficient to ensure survival. It also matters whether the progeny actually look like their parents.
Evolutionary game theory Collections of (mostly) true-breeding cell subpopulations interact by modulating each other's net replication rates in a pairwise, linear manner In a population where population composition influences the absolute fitness of both relatively unfit and relatively fit subpopulations, survival of the relatively fit can lead to decline in fitness for all individuals. The rich and the poor both get poorer.
Cellular automata Probabilistic or deterministic prescription for changing the discrete state of a lattice point into the future is a function of the current configuration of its local spatial neighborhood Cell subpopulations that would otherwise die out in a well-mixed population can survive by huddling together and minimizing their interactions with hostile neighbors. Spatiality can help sustain diversity.