Mathematical models of biological reactions at the system-level lead to a large set of ordinary differential equations with many unknown parameters that need to be inferred using relatively few experimental measurements. Having a reliable and robust algorithm for parameter inference and prediction of the hidden dynamics has been one of the core subjects in systems biology, and is the focus of this study. We have developed a novel systems-biology-informed deep learning algorithm that incorporates the system of ordinary differential equations into the neural networks. Enforcing these equations effectively adds constraints to the optimization procedure that manifests itself as an imposed structure on the observational data. Using few scattered and noisy measurements, we are able to infer the dynamics of unobserved species, systematic forcing and the unknown model parameters. We have successfully tested the algorithm for three different benchmark problems.
Author summaryThe dynamics of systems biological processes are usually modeled using ordinary differential equations (ODEs), which introduce various unknown parameters that need to be estimated efficiently from noisy measurements of concentration for a few species only. In this work, we present a novel "systems-informed neural network" to infer the dynamics of experimentally unobserved species as well as the unknown parameters in the system of equations. By incorporating the system of ODEs into the neural networks we effectively add constraints to the optimization algorithm, which makes the method robust to measurement noise and few scattered observations. 2 holistic approach to deciphering the complexity of biological systems. To understand 3 the biological systems, we must understand the structures of the systems (both their 4 components and structural relationships), and their dynamics [1]. The dynamics of 5 systems biological processes are usually modeled using ordinary differential equations 6 (ODEs) that describe the time evolution of chemical and molecular species 7 concentrations. Once the pathway structure of chemical reactions is known, the 8 corresponding equations can be derived using widely accepted kinetic laws, such as the 9 law of mass action or the Michaelis-Menten kinetics [2]. 10 November 28, 2019 1/15Most system-level biological models introduce various unknown parameters, which 11 need to be estimated efficiently. Thus, the central challenge in computational modeling 12 of these systems could be the prediction of model parameters such as rate constants or 13 initial concentrations, and model trajectories such as time evolution of experimentally 14 unobserved concentrations. Due to the importance of parameter estimation, a lot of 15 attention has been given to this problem in the systems biology community. A lot of 16 research has been conducted on the applications of several optimization techniques, such 17 as linear and nonlinear least-squares fitting [3], genetic algorithms [4], and evolutionary 18 computation [5]. Considerable interest has a...