/rb_vs_ode_model_of_darpp-32

Repository storing molecular rule-based models in Kappa language and Python-based figures published as the BioRxiv preprint: https://www.biorxiv.org/content/10.1101/2022.03.26.485918v1

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network

Authors: Emilia M. Wysocka, Matthew Page, James Snowden, T. Ian Simpson

Repository storing molecular rule-based models in Kappa language and Python-based figures published as the BioRxiv preprint.

Abstract

Dynamic modelling has significantly improved our understanding of the complex molecular mechanisms. The most detailed and popular approach to model dynamics of molecular systems are ordinary diffferential equations (ODEs). However, their application in signalling networks, characterised by multi-state molecular complexes, can be forbidding. These issues have been addressed by contemporary methods, such as rule-based (RB) modelling. Although advantages of RB modelling over ODE modelling have been presented and discussed in numerous reviews, no direct comparison of time courses of one molecular system encoded in the two frameworks have been presented before. To perform such comparison, a set of reactions that underlie an ODE model by Fernandez et al. 2006 [1] was manually encoded into one of RB frameworks, the Kappa language. Closer comparison of models was performed on the level of model specification and results acquired through model simulations. We established that rule patterns reduce reaction number in a certain type of model components. Furthermore, the Kappa model recapitulated general dynamics represented by its ODE counterpart with minor discrepancies. They occurred whenever a molecules have multiple sites binding the same interactor. Notation of such rules requires a complete listing of all possible binding configurations. Furthermore, activation of such molecules is slower than in the ODE model. Therefore, rate constants that define this type of rules ought to be revised. We conclude that the RB representation offers a more expressive and flexible syntax that eases access to fine-grain details of the model, facilitating model reuse. Alongside these analyses, this manuscript reports a refactored model of DARPP-32 interaction network that can serve as a canvas for development of more complex interaction network to study this particular molecular system.

List of figures reproduced in main.ipynb:

  • Figure 7
  • Figure 8
  • Figure 10
  • Figure 11
  • Figure 12
  • Figure S1
  • Figure S2
  • Figure S3
  • Figure S4

References

[1] Fernandez E, Schiappa R, Girault JA, Le Novère N. DARPP-32 is a robust integrator of dopamine and glutamate signals. PLoS Comput Biol. 2006 Dec 22;2(12):e176. doi: 10.1371/journal.pcbi.0020176. Epub 2006 Nov 6. PMID: 17194217; PMCID: PMC1761654.