470 likes | 573 Views
Analyzing the systemic function of genes and proteins. Rui Alves. Organization of the talk. From networks to physiological behavior Network representations Mathematical formalisms Studying a mathematical model. In silico networks are limited as predictors of physiological behavior.
E N D
Analyzing the systemic function of genes and proteins Rui Alves
Organization of the talk • From networks to physiological behavior • Network representations • Mathematical formalisms • Studying a mathematical model
In silico networks are limited as predictors of physiological behavior Probably a very sick mutant? What happens?
How to predict behavior from network? • Build mathematical models!!!!
Organization of the talk • From networks to physiological behavior • Network representations • Mathematical formalisms • Studying a mathematical model
A B A B A B A B B A A B Function Function Function Function Function Network representation is fundamental for clarity of analysis What does this mean? Possibilities:
Defining network conventions C - + A B Full arrow represents a flux between A and B Dashed arrow with a plus sign represents positive modulation of a flux Dashed arrow with a minus sign represents negative modulation of a flux Dashed arrow represents modulation of a flux
Organization of the talk • From networks to physiological behavior • Network representations • Mathematical formalism • Studying a mathematical model
Flux Linear A or C Saturating Sigmoid What is the form of the function? C + A B
What if the form of the function is unknown? C + A B Taylor Theorem: f(A,C) can be written as a polynomial function of A and C using the function’s mathematical derivatives with respect to the variables (A,C)
Are all terms needed? C + A B f(A,C) can be approximated by considering only a few of its mathematical derivatives with respect to the variables (A,C)
Linear approximation C + A B Taylor Theorem: f(A,C) is approximated with a linear function by its first order derivatives with respect to the variables (A,C)
What if system is non-linear? • Use a first order approximation in a non-linear space.
Logarithmic space is non-linear C + A B Use Taylor theorem in Log space g<0 inhibits flux g=0 no influence on flux g>0 activates flux
Why log space? • Intuitive parameters • Simple, yet non-linear • Linearizes exponential space • Many biological processes are close to exponential → Linearizes mathematics
Why is formalism important? • Reproduction of observed behavior • Tayloring of numerical methods to specific forms of mathematical equations
Organization of the talk • From networks to physiological behavior • Network representations • Mathematical formalism • Studying a mathematical model
A model of a biosynthetic pathway _ Constant X0 X1 X2 X3 + X4 Protein using X3
What can you learn? • Steady state response • Long term or homeostatic systemic behavior of the network • Transient response • Transient of adaptive systemic behavior of the network
What else can you learn? Sensitivity of the system to perturbations in parameters or conditions in the medium Stability of the homeostatic behavior of the system Understand design principles in the network as a consequence of evolution
How is homeostasis of the flux affected by changes in X0? Increases in X0 always increase the homeostatic values of the flux through the pathway Log[V] Large g10 Medium g10 Low g10 Log[X0]
How is flux affected by changes in demand for X3? Log[V] Large g13 Medium g13 Low g13 Log[X4]
How is homeostasis affected by changes in demand for X3? Log[X3] Large g13 Medium g13 Low g13 Log[X4]
What to look for in the analysis? • Steady state response • Long term or homeostatic systemic behavior of the network • Transient response • Transient of adaptive systemic behavior of the network
Transient response analysis Solve numerically
Specific adaptive response Get parameter values Get concentration values Substitution [X3] Solve numerically Change in X4 Time
General adaptive response Normalize Unstable system, uncapable of homeostasis if feedback is strong!! Solve numerically with comprehensive scan of parameter values [X3] High g13 Increasing g13 Low g13 Threshold g13 Increase in X4 Time
Sensitivity analysis • Sensitivity of the system to changes in environment • Increase in demand for product causes increase in flux through pathway • Increase in strength of feedback increases response of flux to demand • Increase in strength of feedback decreases homeostasis margin of the system
Stability analysis • Stability of the homeostatic behavior • Increase in strength of feedback decreases homeostasis margin of the system
How to do it • Download programs/algorithms and do it • PLAS, GEPASI, COPASI SBML suites, MatLab, Mathematica, etc. • Use an on-line server to build the model and do the simulation • V-Cell, Basis
Design principles • Why is a given pathway design prefered over another? • Overall feedback in biosynthetic pathways • Why are there alternative designs of the same pathway? • Dual modes of gene control
Why regulation by overall feedback? _ Overall feedback X0 X1 X2 X3 + _ _ _ X4 X0 X1 X2 X3 Cascade feedback + X4
Overall feedback improves functionality of the system [C] [C] Overall Overall Cascade Cascade Overall Cascade Stimulus Spurious stimulation Proper stimulus Time
Demand theory of gene control High demand for gene expression→ Positive Regulation Low demand for gene expression → Negative mode of regulation Wall et al, 2004, Nature Genetics Reviews
How to do it • Download programs/algorithms and do it • BST Lab, Mathematica, Maple
Summary • From networks to physiological behavior • Network representations • Mathematical formalism • Studying a mathematical model
Papers to present • Vasquez et al, Nature • Alves et al. Proteins
Computational tools in Molecular Biology • Predictions & Analysis • Identification of components • Organization of components • Conectivity of components • Behavior of systems • Evolution & Design • Prioritizing wet lab experiments • Most likely elements to test • Most likely processes to test
The Taylor theorem f(C) ith order 1st order f(C) ith + jth order 0 order 2nd order C
Are all terms needed? C + A B f(A,C) can be approximated by considering only a few of its mathematical derivatives with respect to the variables (A,C)
Linear approximation C + A B Taylor Theorem: f(A,C) is approximated with a linear function by its first order derivatives with respect to the variables (A,C)
What if flux is non linear? C + A B Use Taylor theorem with large number of terms or Use Taylor theorem in Non-Linear space!
Y Y X X How does the transformation between spaces work?
How does the Taylor approximation work in another space? Variables: A, B, C, … Variables: A, B, C, … Transform to new space f(A,B,…) f(A,B,…) ~f(A,B,…) ~f(A,B,…) Taylor theorem Return to original space