Abstract
Many biological processes are described with coupled non-linear systems of ordinary differential equations that contain a plethora of parameters. The goal is to understand these systems and to predict the effect of different influences. This asks for a dynamical systems approach where numerical continuation methods and bifurcation analysis are used to detect the solutions and their stability as a function of the parameters. We developed PyNCT – Python Numerical Continuation Toolbox – an open source Python package that implements numerical continuation methods and can perform bifurcation analysis based on sparse linear algebra. The software gives the user the choice of different solvers (direct and iterative) and allows the use of preconditioners to reduce the number of iterations and guarantee the convergence when working with complex non-linear models.
In this paper we demonstrate the usefulness of the toolbox with a class of models pertaining to auxin transport between cells in plant organs.We show how easy it is to compute the steady state solutions for different parameter values, to calculate how they depend on each other and to map parts of the solution landscape.
An interactive model development and discovery cycle is key in bio-systems research. It allows one to investigate and compare different model parameter settings and even different models and gauge the model’s usefulness. Our toolbox allows for such quick experimentation and has a low entry barrier for non-technical users.
Although PyNCT was developed particularly for the study of transport models in biology, its implementation is generic and extensible, and can be used in many other dynamical system applications.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Allgower, E., Georg, K.: Numerical Path Following. Springer, Berlin (1994)
Ascher, D.: Numpy, Numerical Python (2001). http://www.numpy.org/. Accessed June 2015
Beemster, G., Baskin, T.: Analysis of cell division and elongation underlying the developmental acceleration of root growth in arabidopsis thaliana. Plant Physiol. 116, 515–526 (1998)
Brunoud, G., Wells, D., Oliva, M., Larrieu, A., Mirabet, V., Burrow, A., Beeckman, T., Kepinski, S., Traas, J., Bennett, M., et al.: A novel sensor to map auxin response and distribution at high spatio-temporal resolution. Nature 482, 103–106 (2012)
Dhooge, A., Govaerts, W., Kuznetsov, Y.A.: Matcont: a matlab package for numerical bifurcation analysis of odes. ACM Trans. Math. Softw. (TOMS) 29(2), 141–164 (2003)
Doedel, E.J.: Auto: A program for the automatic bifurcation analysis of autonomous systems. Congr. Numer. 30, 265–284 (1981)
Draelants, D., Avitabile, D., Vanroose, W.: Localized auxin peaks in concentration-based transport models of the shoot apical meristem. J. R. Soc. Interface 12(106) (2015). http://rsif.royalsocietypublishing.org/content/12/106/20141407.full
Jones, E., Oliphant, T., Peterson, P., et al.: SciPy: open source scientific tools for Python (2001). http://www.scipy.org/. Accessed June 2015
Kelley, C.T.: Iterative methods for linear and nonlinear equations. SIAM: Frontiers in Applied Mathematics 16 (1995)
Kelley, C.T.: Solving Nonlinear Equations with Newton’s Method, vol. 1. Society for Industrial and Applied Mathematics, Philadelphia (2003)
Kłosiewicz, P.: Pypts, python plant tissue simulations (2015). https://pypi.python.org/pypi/PyPTS/. Accessed June 2015
Krauskopf, B., Osinga, H., Galán-Vioque, J.: Numerical Continuation methods for Dynamical Systems. Springer, Netherlands (2007)
Merks, R., Guravage, M., Inzé, D., Beemster, G.: Virtualleaf: an open-source framework for cell-based modeling of plant tissue growth and development. Plant Physiol. 155(2), 656–666 (2011)
Salinger, A.G., Bou-Rabee, N.M., Pawlowski, R.P., Wilkes, E.D., Burroughs, E.A., Lehoucq, R.B., Romero, L.A.: Loca 1.0 library of continuation algorithms: theory and implementation manual. Sandia National Laboratories, Albuquerque, NM, Technical Report No. SAND2002-0396 (2002)
Schmidt, H., Jirstrand, M.: Systems biology toolbox for matlab: A computational platform for reasearch in systems biology. Bioinformatics Advance Access (2005)
Seydel, R.: Practical Bifurcation and Stability Analysis: From Equilibrium to Chaos, vol. 5. Springer, New York (1994)
Siso-Nadal, F., Ollivier, J.F., Swain, P.S.: Facile: a command-line network compiler for systems biology. BMC Syst. Biol. 1(1), 36 (2007)
Smith, R., Guyomarc’h, S., Mandel, T., Reinhardt, D., Kuhlemeier, C., Prusinkiewicz, P.: A plausible model of phyllotaxis. PNAS 103(5), 1301–1306 (2006)
SymPy Development Team: SymPy: Python library for symbolic mathematics (2014). http://www.sympy.org
Acknowledgements
DD acknowledges financial support from the Department of Mathematics and Computer Science of the University of Antwerp. This work is part of the Geconcerteerde Onderzoeksactie (G.O.A.) research grant “A System Biology Approach of Leaf Morphogenesis” granted by the research council of the University of Antwerp. We acknowledge Giovanni Samaey for sharing basic version of a continuation code.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2015 Springer International Publishing Switzerland
About this paper
Cite this paper
Draelants, D., Kłosiewicz, P., Broeckhove, J., Vanroose, W. (2015). Solving General Auxin Transport Models with a Numerical Continuation Toolbox in Python: PyNCT. In: Abate, A., Šafránek, D. (eds) Hybrid Systems Biology. HSB 2015. Lecture Notes in Computer Science(), vol 9271. Springer, Cham. https://doi.org/10.1007/978-3-319-26916-0_12
Download citation
DOI: https://doi.org/10.1007/978-3-319-26916-0_12
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-26915-3
Online ISBN: 978-3-319-26916-0
eBook Packages: Computer ScienceComputer Science (R0)