What is systems biology?

1. What is systems biology? Being a mathematician in a biologist’s and a bioinformatician’s world Zofia Jones, PhD Mathematical Sciences

2. Outline • Intro • Mass action / ODE example. • Metabolic modelling example.

3. A mathematicalbiologist is… by definition interdisciplinary. My main motivation: a chance to use a broader range of mathematical skills than I could in physics and fewer rules. (but that’s just me, lots of people love biology for its own sake) However, progress is somewhat dependent on effective collaboration and communication.

4. First: What is the computational focus of this group? • Something like…

5. First: What is the computational focus of this group? • Identification of species at OTU level. • Fitting diversity distributions to hypotheses – rarefication, neutral theory with immigration from one or more metacommunities • Accounting for diversity using environmental variables. • Lots of bioinformatics – de-noising, assembly, visualisation, phyologenic trees, identification • Metagenomics – takes the bioinformatics challenge up a notch in complexity and looks at function as well as phylogeny

6. But more about diversity than function • But to understand diversity we need to ask how microorganisms are competitive • And for that we need to understand how they function • How does the available energy in the environment translate to fitness? • What are the principles underlying diversity? • Then match to data.

7. The big goals • Pharmaceuticals £££ • Alternative fuels – bacteria to produce ethanol from waste biomass • Bioremediation, phytoremediation, mycoremediation • Use microorganisms to grow building materials and cellulose based “plastics” – check out BioMason, Ecovative… • Want to engineer specific metabolic pathways and their efficiency

8. Scientific Method- hypothesis and evidence Deduction / Induction

9. Need predictive models • on a cell scale • or/and on community / ecological scale This covers with a LOT of science and expertise

10. What do predictive models do? • Need to integrate information on metabolic pathways, regulation, kinetics... • See if we can reproduce what we observe in experiments on a computer • Predict growth/no growth, specific pathways, co-existence, inhibition factors • Help plan experiments. • Help save money. • Help save time. • (Inspire funding)

11. The most common tasks ofsystems biologistswhen modelling a standardgene regulatory network- an example systems biology model

12. Figure 2. Network models, derived from the heuristic MIM shown in Figure 1, for simulation. Detailed Diagram -lots of interaction. Lots of research and Biologist’s input on this -many years of work. Kim S, Aladjem MI, McFadden GB, Kohn KW (2010) Predicted Functions of MdmX in Fine-Tuning the Response of p53 to DNA Damage. PLoSComputBiol 6(2): e1000665. doi:10.1371/journal.pcbi.1000665 http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1000665

13. Research Similar Models with elements similar to yours.

14. Make assumptions on rate constants - group them by type, evidence of relative magnitude

15. Write down some equations -here we have simple mass action (rate proportional to concentration) -though more complicated with less information

16. Make assumptions on Initial conditions

17. Fit to experimental data -some well developed data here

18. Sensitivity Analysis – how sensitive are outputs to parameters?

19. Figure 10. Prediction of late response to DNA damage. Make some predictions on dynamics -qualitative statements are best eg. oscillation or decay? Give this back to the Biologists. Kim S, Aladjem MI, McFadden GB, Kohn KW (2010) Predicted Functions of MdmX in Fine-Tuning the Response of p53 to DNA Damage. PLoSComputBiol 6(2): e1000665. doi:10.1371/journal.pcbi.1000665 http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1000665

20. Figure 7. Bifurcation diagram of the effects of MdmX on p53 oscillatory behavior. • So some stability analysis • where in parameter • space do these behaviours • occur? Kim S, Aladjem MI, McFadden GB, Kohn KW (2010) Predicted Functions of MdmX in Fine-Tuning the Response of p53 to DNA Damage. PLoSComputBiol 6(2): e1000665. doi:10.1371/journal.pcbi.1000665 http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1000665

21. What might you think of when you think about mathematical modelling? Experimental biologists are pragmatic people…

22. What might you think of when you think about mathematical modelling? • Something complicated or time consuming. • Memories of tediousmaths lessons at school. • Michaelis-Menten enzyme kinetics. • Don’t know where to start. • Something to be used sparingly for specific questions or problems. • There are interdisciplinary courses available.

23. Actually need a lot of qualitative knowledge. I spend a lot of time reading the small print. This is the real research bottleneck for numerical people in biology. I miss sums … Bad luck Charlie Brown.

24. What I think of … • Asking lots of questions • What is the science I need to learn? • What is the experimental evidence? • What do we want to explain? • What is the fundamental cause/effect • behind these results and discussion papers? • How can I account for these relations • quantitatively?

25. What I think of … • I need to know a lot about the biology: details on what drives syntrophy, evolutionary trade-offs, information on primary and secondary metabolites, phylogeny of enzymes, thermodynamic limitations, thermodynamic gradients, co-regulation, regulation of pH, limiting factors, summary of open questions, good experimental data. • Need specific, precise information and clearly expressed ideas and theories. • I need to know a lot about the relevant mathematicswhich fits this biological information. • The maths is not driven by complexity or difficulty but relevance. • ODEs, graph theory, flux balance analysis with further constraints….

26. Again… • Need specific, precise information and clearly expressed ideas and theories. (and biology is big on vagueness … ) • The maths is not driven by complexity or difficulty but relevanceand context. • Should be helpful and accessible to both mathematicians and biologist’s specialised in that area. • If it is confusing they haven’t done their job.

27. Let’s look at some infographics“The Virtuous Cycle”

28. to constrain models to be more precise • infer regulatory networks • better kinetic data, information on interactions • to put the bigger picture theories • on evolutionary biology or cell function to the test Modern research : Different people involved in different steps. • microarray • transcriptomics • diversity data • environmental variables • metabolic reconstruction • mass action kinetics • directed graphs

29. Eg. One lady Jeanette Johnson • Who works with mathematicians in Oxford on diffusion problems. • Said that she often find things in her experiments which the modellers then can explain in theory. • And they can predict things they find in experiment. • So that’s the ideal situation.

30. Requires lots of skills and lots of people… Requires communication, teamwork and patience.

31. After all, each cell can be viewed as a tiny computer with a core program modified byexperimental evidence.

32. Scientific Deduction / Induction

33. Use your imagination!

34. Need predictive models • on a cell scale • or/and on community / ecological scale This covers with a LOT of science and expertise

35. Another example.

36. Metabolic Modelling ie. Only need genome to get stoichiometry of network and some estimates of parameters which can be provided by Kbase. Then can improve annotation as required.

37. Solution of metabolic model == • Net flux at each node = 0 • Flux is concentration x rate • Extracellular source terms for substrates etc. • Sink terms for biomass. • Assign a vector to be optimisedeg. biomass flux • Standard linear discrete optimisation problem. • Many alternative solutions are usually possible. • Finding the biologically meaningful ones…

38. How to build a metabolic model • Genome: get it, annotate it. • All you need is databases. • Get it: sequencing then assembly, NCBI, arCOG • Initial annotation: curated genomes NCBI, RAST • Additional annotation: comparative genomics with MAUVE, literature review

39. How to build a metabolic model • Draft a model: join the dots and create a sbml file. • SEED or Kbase – same software different packages. • SEED: slow and over-subscribed. • Kbase: command line and faster. • Draft built on RAST annotation then can add additional missing reactions manually.

40. How to build a metabolic model • Edit model: add or delete reactions • Fit model to growth or no growth data. • This data is usually specified by choice of media or gene deletions. • False growth positive requires more gaps. • False growth negative requires more reactions. • Can be done efficiently on KBase. • Put model in paper.

41. Some Kbase Commands kbws-login zofia1 -pl******1 kbws-urlhttp://kbase.us/services/workspace_service kbws-workspace glasgow kbws-listobj kbws-urlhttp://kbase.us/services/fba_model_services kbfba-getmediaacetate_minimalglasgow -e >> acetate_minimal.txt annotate_genomecontigs.fastaconcilii kbfba-buildfbamodelconciliimethanosaeta kbfba-gapfillmethanosaeta –macetate_minimal = gapfilled model

42. Curate from getting gapfilling results.

43. = good enough model Metabolic models mainly used to show we have correct understanding of network via growth/ no growth data.

44. Metabolic model building can be used to check understanding as we make guesses about which pathways are present or active in a given environment. • Set constraints to reflect presence/absence of substrates. • Run flux balance analysis to get steady state solution. • See which pathways are active in solution. • Use solution to help interpret transcriptomics, microarray data, metagenomics or PCR. • Adjust understanding of organism or model as necessary.

45. How to use a metabolic model • Check understanding of topology: the most complicated bit • Done mainly by referring to growth or no growth data. • Databases are light on info on archaea and non-pathogenic micro-organisms. • Difficult is microorganisms can only grow on a limited range of media – fewer experimental options. • Challenge to chose the correct edits.

46. How to use a metabolic model • Curate detail: learn about your pet organisms • Some detail to add… proton/electron ratio = getting the net ATP produced correct

47. How to use a metabolic model • Curate detail: learn about your pet organisms • Some detail to add… accurate rate of ATP production -> accurate growth rate Minimal ATP requirement. ATP requirement increases linearly with growth.

48. How to use a metabolic model • Curate detail: learn about your pet organisms • Some detail to add… biomass composition

49. How to use a metabolic model • Curate detail: learn about your pet organisms • Some detail to add… Thermodynamics: is reaction reversible or not?

50. How to use a metabolic model • Add kinetics: model works out yield • Just add Michaelis-Menten kinetics for carbon and nitrogen sources.