Guido Sanguinetti, Magnus Rattray and Neil D. Lawrence

1. A Probabilistic Dynamical Model for Quantitative Inference of the Regulatory Mechanism of Transcription Guido Sanguinetti, Magnus Rattray and Neil D. Lawrence

2. Talk plan • Overview of the problem • Extending regression • Introducing dynamics • Modelling separately concentrations • What next?

3. The problem • The Central Dogma Genes mRNA Proteins Life Transcription Easy to measure Translation Hard to measure Protein interactions COMPLEX!

4. Specific problem • Transcription factors produce proteins that promote or repress transcription of other genes; they play a fundamental role in gene networking • Deduce the activity of the transcription factors’ proteins (in an experimental condition) from the mRNA expression data.

5. Why not use the TFs expressions? TFs are often low expressed, noisy TFs are post-transcriptionally regulated TFs interact non-trivially with each other

6. Current approaches • Integrate with ChIP-on-chip data • ChIP-on-chip gives a binary matrix X of transcription factors binding genes (connectivity matrix) • Regress microarray expression data on X bmt is the transcription factor activity (TFA) of TF m at time t, monotonically linked to protein concentrations (Liao et al, Boulesteix and Strimmer, Gao et al,...)

7. Problems • All genes bound by the TF contribute equally to the estimate of the TFA, regardless of the regulation type. • TFAs are gene-independent, but the influence of a transcription factor varies from gene to gene (and according to condition) • The model is linear (inevitable)

8. Extending Regression Modify the regression model to allow different TFAs for different genes and experiments Reduce the number of parameters by placing a prior distribution over the gene-specific TFAs. The choice of the prior distribution depends on the situation we model. E.g., for independent samples we may assume TFAs at different time points to be independent

9. Introducing dynamics • To model time series data, we choose a Kalman filter prior on the rows of B where This is equivalent to assuming TFAs vary smoothly

10. Likelihood function • Given the model and the prior, we can obtain a likelihood The likelihood can be estimated efficiently using the sparsity of the covariance and recursion relations.

11. Estimating the TFAs TFAs can be estimated a posteriori using Bayes’s Theorem and moment matching Error bars associated with each TFA are given by the squared root of the diagonal entries in the posterior covariance. Mean TFAs can be obtained by averaging gene-specific TFAs over the target genes.

12. Testing the model • We compared our averaged TFAs with the ones obtained by regression for the Spellman dataset (Mol.Biol.Cell,1998), ChIP data from Lee et al. (Science 2002). The diagrams show the TFA for ACE2p.

13. ...but we also get... TFA for CTS1 TFA for SCW11 TFA for YER124C TFA for YKL151C

14. ...and we can do more! • Error bars allow to determine which regulations are significant • Correlations among TFs can be obtained from Σ

15. Decoupling action and concentration • It is not clear in the model whether a high gene-specific TFA is the result of a high affinity or of a high protein concentration • We modify the model to distinguish the effects of protein concentration and affinity • Specifically, we model

16. Estimating the parameters • The model is no longer exact. • Approximate inference is performed using a variational EM algorithm • This exploits Jensen’s inequality to get a bound on the log likelihood Under a factorization assumption on the approximating distribution q, the E-step becomes exactly solvable via fixed point equations.

17. Results The left hand picture shows the expression level of ACE2 in the yeast cell cycle, the middle shows the inferred protein concentration and right shows the significance of the activities.

18. Problems • ChIP data is notoriously noisy; for example the same transcription factor (MSN4) in the same conditions (rich medium) is found to bind 32 genes in Lee et al. and 57 genes in Harbison et al. (the intersection is 20 genes). • Posterior estimation helps with false positives, not with false negatives. • The model is additive (in log space) and doesn’t model combinatorial effects.

19. What next? • Collaborate with biologists to validate our predictions on novel data • Microarray and ChIP data from same lab should be more consistent • Use the model results as a starting point for systems biology modeling • Introduce combinatorial effects