pyFOOMB: Python Framework for Object Oriented Modelling of Bioprocesses

Quantitative characterization of biotechnological production processes requires the determination of different key performance indicators (KPIs) such as titer, rate and yield. Classically, these KPIs can be derived by combining black-box bioprocess modelling with non-linear regression for model parameter estimation. The presented pyFOOMB package enables a guided and flexible implementation of bioprocess models in the form of ordinary differential equation systems (ODEs). By building on Python as powerful and multi-purpose programming language, ODEs can be formulated in an object-oriented manner, which facilitates their modular design, reusability and extensibility. Once the model is implemented, seamless integration and analysis of the experimental data is supported by various Python packages that are already available. In particular, for the iterative workflow of experimental data generation and subsequent model parameter estimation we employed the concept of replicate model instances, which are linked by common sets of parameters with global or local properties. For the description of multi-stage processes, discontinuities in the right-hand sides of the differential equations are supported via event handling using the freely available assimulo package. Optimization problems can be solved by making use of a parallelized version of the generalized island approach provided by the pygmo package. Furthermore, pyFOOMB in combination with Jupyter notebooks also supports education in bioprocess engineering and the applied learning of Python as scientific programming language. Finally, the applicability and strengths of pyFOOMB will be demonstrated by a comprehensive collection of notebook examples.

[1]  Dragan D. Nikolic Parallelisation of equation-based simulation programs on heterogeneous computing systems , 2018, PeerJ Comput. Sci..

[2]  Eva Balsa-Canto,et al.  Optimality and identification of dynamic models in systems biology: an inverse optimal control framework , 2018, Bioinform..

[3]  Dragan D. Nikolic DAE Tools: equation-based object-oriented modelling, simulation and optimisation software , 2016, PeerJ Comput. Sci..

[4]  Joel Nothman,et al.  SciPy 1.0-Fundamental Algorithms for Scientific Computing in Python , 2019, ArXiv.

[5]  José Ignacio Hidalgo,et al.  Parallel Architectures and Bioinspired Algorithms , 2010, Parallel Comput..

[6]  Dario Izzo,et al.  A parallel global multiobjective framework for optimization: pagmo , 2020, J. Open Source Softw..

[7]  Fabian J. Theis,et al.  Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks , 2016, bioRxiv.

[8]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[9]  K. Jarrod Millman,et al.  Array programming with NumPy , 2020, Nat..

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

[11]  Johan Åkesson,et al.  Assimulo: A unified framework for ODE solvers , 2015, Math. Comput. Simul..

[12]  Fabian J. Theis,et al.  Parameter estimation for dynamical systems with discrete events and logical operations , 2016, Bioinform..

[13]  Jan Marienhagen,et al.  Miniaturized and automated adaptive laboratory evolution: Evolving Corynebacterium glutamicum towards an improved d-xylose utilization. , 2017, Bioresource technology.

[14]  Carole J. Proctor,et al.  PyCoTools: a Python toolbox for COPASI , 2018, Bioinform..

[15]  Wolfgang Wiechert,et al.  Corynebacterium glutamicum Chassis C1*: Building and Testing a Novel Platform Host for Synthetic Biology and Industrial Biotechnology. , 2017, ACS synthetic biology.

[16]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..