Integrated Modeling and Optimization (Fundamental)

Location
Wageningen Campus, The Netherlands
Forum Building (building nr. 102)
December 11, room B0314,
December 12, room B0326,
December 13, 14, and 15, room B0217
Course coordinator

Prof. Natal van Riel, Eindhoven University of Technology
Dr. Dragan Bosnacki, Eindhoven University of Technology

Lecturers

Software: Matlab (with following toolboxes: Optimization toolbox, Statistics toolbox)

Course description

Living organisms are characterized by an amazing degree of hierarchical complexity. Although our ability to collect measurements at different spatial levels and time-scales has grown dramatically, it has become clear that only measurements cannot provide the answer to unravelling biological complexity. This is because the dynamical behavior of complex systems cannot be reduced to the linear sum of the functions of their parts. Hence, computational modelling is an absolute requisite to gain understanding of the mechanisms underlying patterns observed in experimental data, in particular when studying dynamic phenomena. Mathematical models allow in a relatively cheap way to generate and test hypotheses about these mechanisms. However, given the huge complexity and peculiar features of biological systems, it is necessary to carefully understand the specific modelling requirements they pose, in order to define what a good model should look like. In this way one could say that modelling is a craftsmanship, that can only be learned via intense exercising and ‘learning by doing’. In this course we offer the participants the possibility to learn and exercise the modeling process.

In validating models one always meet with the need to fit models to data. So, the parameters that are present in any realistic model have to be chosen based on comparison of model predictions with data. In this matching process optimization techniques are indispensable. That’s why a considerable part of this course is spent on getting you acquainted with the optimization techniques that are nowadays available and widely used. Numerical optimization also is the basis for so-called flux balance analysis (FBA), commonly used to study large metabolic networks. This type of models and their analysis and simulation is also introduced in the course.

The course is a mixture of theory sessions and computer practicals. The course is completed with an assignment to be finalized afterwards, for those who want to acquire 3 ECTS.

 

 

Daily schedule

9:30-12:30 lecture
12:30-14:00 lunch + meet-the-teacher
14:00-16:00 pen & paper exercises and/or computer practical
16:15-17:00 research talk + wrap up

Except Friday: 9.30-12.30 lecture & practical; 14.00 – 15.00 wrap up

Day 1 (11 December)

Introduction Systems Biology, math tools, optimization
Prof. Jaap Molenaar

Topics:
– brief intro participants at the start of the course
– introduction of modelling with (ordinary) differential equations
– introduction of optimization concepts and methods
– practical (1): phase plane analysis
– practical (2): Matlab intro (Hans Stigter, Dragan Bošnački)
Research talk: Dragan Bosnacki discusses the topic linked to the assignment

Day 2 (12 December)

Constraint-based modeling and analysis
Dr. Rob Smith

Topics:
– reaction kinetics
– Flux Balance Analysis
– linear programming optimisation
– Dynamic FBA
– practical using the COBRA package in MATLAB

Research talk: Fabian Frohlich (Crick Institute, UK) ‘The Hitchhiker’s Guide to Complex Models of Intracellular Signalling’

Day 3 (13 December)

Networks, Network motifs, Network inference
Dr. Aalt-Jan van Dijk

Topics:
– network structure, motifs, dynamics
– network inference
– bias/variance: overfitting
Practical (Matlab): network inference with thickening – thinning

Research talk: Monica Garcia Gomez (Theoretical Biology and Bioinformatics, Utrecht University) `A multi-scale model of asymmetric cell fate decisions in the root stem cell niche’

Day 4 (14 December)

Numerical optimization, Parameter estimation, Sensitivity analysis
Dr. Hans Stigter

Topics:
– Maximum Likelihood Estimation (MLE)
– linear and nonlinear least squares
– optimization algorithms (local and global), steepest descent, Gauss-Newton, Levenberg-Marquardt
– sensitivity (local)
– Fisher Information Matrix
– confidence intervals
– experimental design
– practical (Matlab)

Research talk: Dr. Henk van Lingen (WU) ‘Bayesian calibration of dynamic models’

17:30  Social event with YoungCB (details to be announced)

Day 5 (15 December)

Developing credible simulation models of biological systems
Prof. Dr. Natal van Riel

Topics:
– concepts on uncertainty and uncertainty quantification
– profile likelihood analysis
– bias-variance trade-off (regularization)
– multi-start gradient search
– experimental design

Discussion about the assignment (Dragan Bošnački)

 

Learning objectives

The students will be provided with a theoretical basis, a variety of methods, and a computational hands-on experience to set-up systems biology models and handle numerical optimization.

In the course the students will learn:

  • To understand the common ground and the differences for applications of dynamic modeling in metabolic, regulatory, signaling, and multi-scale biological processes
  • How to set-up a dynamic model to represent biological networks using different interaction mechanisms
  • To implement, simulate and analyze dynamic network models
  • To understand the wide variety of problems in modelling that can be solved with optimization
  • To apply different types of numerical optimization methods
  • The combination of dynamic modeling and optimization to integrate experimental data in modelling, estimate model parameters and design experiments.
  • To understand how numerical optimization (linear programming) works in flux balance analysis to simulate metabolic network models.

Techniques include:

  • Nonlinear differential equations, numerical simulation, parameter sensitivity analysis.
  • Parameter estimation, identifiability, uncertainty quantification, experimental design, regularization.
  • Global and local search methods: steepest descent, Levenberg-Marquardt, genetic algorithms, linear programming.

Target audience

The course is aimed at PhD students with a background in bioinformatics, systems biology, computer science or a related field, and life sciences. Participants from the private sector are also welcome. A working knowledge of mathematics, especially differential equations, is recommendable, but we will distribute preparation material to be studied by students missing the required background. Furthermore, at the start we offer a math refresher to help those participants who are not (yet) involved in modelling on a daily basis.

Software

Examples and computer practical make use of Matlab. A computer with a working version of Matlab is needed and some programming experience and knowledge of Matlab are required to take the course. A short introductory training in Matlab will be made available (online) for those without Matlab skills.

Registration

Fees for this 5-day course are:

  • PhD student: 480 euro (excl. VAT)
  • Academic researcher (PI/Postdoc): 720 euro (excl. VAT)
  • Industry: 1080 euros (excl.VAT)

The course fee includes:

  • Course material
  • Catering: coffee, tea and lunch will be provided.

You can register for the course by filling out this registration form. Find general enrollment information here.