Effect of Linearization in a WNT Signaling Model

A nonlinear model consisting of a system of coupled ordinary differential equations (ODE), describing a biological process linked with cancer development, is linearized using Taylor series and tested against different magnitudes of input perturbations, in order to investigate the extent to which the linearization is accurate. The canonical wingless/integrated (WNT) signaling pathway is considered. The linearization procedure is described, and special considerations for linearization validity are analyzed. The analytical properties of nonlinear and linearized systems are studied, including aspects such as existence of steady state and initial value sensitivity. Linearization is a useful tool for speeding up drug response computations or for providing analytical answers to problems such as required drug concentrations. A Monte Carlo-based error testing workflow is employed to study the errors introduced by the linearization for different input conditions and parameter vectors. The deviations between the nonlinear and the linearized system were found to increase in a polynomial fashion w.r.t. the magnitude of tested perturbations. The linearized system closely followed the original one for perturbations of magnitude within 10% of the base input vector which yielded the state-space fixed point used for the linearization.