SBML2Julia: interfacing SBML with efficient nonlinear Julia modelling and solution tools for parameter optimization

Motivation: Estimating model parameters from experimental observations is one of the key challenges in systems biology and can be computationally very expensive. While the Julia programming language was recently developed as a high-level and high-performance language for scientific computing, systems biologists have only started to realise its potential. For instance, we have recently used Julia to cut down the optimization time of a microbial community model by a factor of 140. To facilitate access of the systems biology community to the efficient nonlinear solvers used for this optimisation, we developed SBML2Julia. SBML2Julia translates optimisation problems specified in SBML and TSV files (PEtab format) into Julia for Mathematical Programming (JuMP), executes the optimization and returns the results in tabular format. Availability and implementation: SBML2Julia is freely available under the MIT license. It comes with a command line interface and Python API. Internally, SBML2Julia calls the Julia LTS release v1.0.5 for optimisation. All necessary dependencies can be pulled from Docker Hub (this https URL). Source code and documentation are available at this https URL.

[1]  Sungho Shin,et al.  Scalable nonlinear programming framework for parameter estimation in dynamic biological system models , 2018, bioRxiv.

[2]  Iain Dunning,et al.  JuMP: A Modeling Language for Mathematical Optimization , 2015, SIAM Rev..

[3]  Alan Edelman,et al.  Julia: A Fresh Approach to Numerical Computing , 2014, SIAM Rev..

[4]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

[5]  Julio R. Banga,et al.  PEtab -- interoperable specification of parameter estimation problems in systems biology , 2020 .

[6]  Julio R. Banga,et al.  Benchmarking optimization methods for parameter estimation in large kinetic models , 2018, bioRxiv.

[7]  WächterAndreas,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006 .

[8]  Julio R. Banga,et al.  Robust and efficient parameter estimation in dynamic models of biological systems , 2015, BMC Systems Biology.

[9]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[10]  Jens Timmer,et al.  Benchmark problems for dynamic modeling of intracellular processes , 2018, bioRxiv.

[11]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[12]  Fabian J. Theis,et al.  Data2Dynamics: a modeling environment tailored to parameter estimation in dynamical systems , 2015, Bioinform..

[13]  B. Novák,et al.  Model scenarios for switch‐like mitotic transitions , 2015, FEBS letters.

[14]  Eva Balsa-Canto,et al.  AMIGO2, a toolbox for dynamic modeling, optimization and control in systems biology , 2016, Bioinform..