Estimating parameters of S-systems by an auxiliary function guided coordinate descent method

The S-system, a set of nonlinear ordinary differential equations and derived from the generalized mass action law, is an effective model to describe various biological systems. Parameters in S-systems have significant biological meanings, yet difficult to be estimated because of the nonlinearity and complexity of the model. Given time series biological data, its parameter estimation turns out to be a nonlinear optimization problem. A novel method, auxiliary function guided coordinate descent, is proposed in this paper to solve the optimization problem by cyclically optimizing every parameter. In each iteration, only one parameter value is updated and it proves that the objective function keeps nonincreasing during the iterations. The updating rules in each iteration is simple and efficient. Based on this idea, two algorithms are developed to estimate the S-systems for two different constraint situations. The performances of algorithms are studied in several simulation examples. The results demonstrate the effectiveness of the proposed method.

[1]  Jonas S. Almeida,et al.  Parameter optimization in S-system models , 2008, BMC Systems Biology.

[2]  Christine Nardini,et al.  An S-System Parameter Estimation Method (SPEM) for Biological Networks , 2012, J. Comput. Biol..

[3]  H. Sebastian Seung,et al.  Algorithms for Non-negative Matrix Factorization , 2000, NIPS.

[4]  Shuhui Li,et al.  Artificial Neural Networks for Control of a Grid-Connected Rectifier/Inverter Under Disturbance, Dynamic and Power Converter Switching Conditions , 2014, IEEE Transactions on Neural Networks and Learning Systems.

[5]  Fang-Xiang Wu,et al.  Inference of Biological S-System Using the Separable Estimation Method and the Genetic Algorithm , 2012, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[6]  Eberhard O. Voit,et al.  Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists , 2000 .

[7]  Fang-Xiang Wu,et al.  Separable Parameter Estimation Method for Nonlinear Biological Systems , 2009, 2009 3rd International Conference on Bioinformatics and Biomedical Engineering.

[8]  Chih-Hung Hsieh,et al.  An Intelligent Two-Stage Evolutionary Algorithm for Dynamic Pathway Identification From Gene Expression Profiles , 2007, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[9]  Prospero C. Naval,et al.  Parameter estimation using Simulated Annealing for S-system models of biochemical networks , 2007, Bioinform..

[10]  Lin Wu,et al.  A Fast Algorithm for Nonnegative Matrix Factorization and Its Convergence , 2014, IEEE Transactions on Neural Networks and Learning Systems.

[11]  Axel Kowald,et al.  Systems Biology in Practice: Concepts, Implementation and Application , 2005 .

[12]  Eberhard O Voit,et al.  Theoretical Biology and Medical Modelling , 2022 .

[13]  Fang-Xiang Wu,et al.  Alternating weighted least squares parameter estimation for biological S-systems , 2012, 2012 IEEE 6th International Conference on Systems Biology (ISB).

[14]  J. Collins,et al.  Inferring Genetic Networks and Identifying Compound Mode of Action via Expression Profiling , 2003, Science.

[15]  Katta G. Murty,et al.  Nonlinear Programming Theory and Algorithms , 2007, Technometrics.

[16]  James V. Beck,et al.  Parameter Estimation in Engineering and Science , 1977 .

[17]  Edward R. Dougherty,et al.  Inference of Gene Regulatory Networks using S-System: A Unified Approach , 2007, 2007 IEEE Symposium on Computational Intelligence and Bioinformatics and Computational Biology.

[18]  Masaru Tomita,et al.  Dynamic modeling of genetic networks using genetic algorithm and S-system , 2003, Bioinform..

[19]  K. Schittkowski,et al.  NONLINEAR PROGRAMMING , 2022 .

[20]  Jonas S. Almeida,et al.  Decoupling dynamical systems for pathway identification from metabolic profiles , 2004, Bioinform..