Parameter estimation using Simulated Annealing for S-system models of biochemical networks

Abstract

Motivation: High-throughput technologies now allow the acquisition of biological data, such as comprehensive biochemical time-courses at unprecedented rates. These temporal profiles carry topological and kinetic information regarding the biochemical network from which they were drawn. Retrieving this information will require systematic application of both experimental and computational methods. Results: S-systems are non-linear mathematical approximative models based on the power-law formalism. They provide a general framework for the simulation of integrated biological systems exhibiting complex dynamics, such as genetic circuits, signal transduction and metabolic networks. We describe how the heuristic optimization technique simulated annealing (SA) can be effectively used for estimating the parameters of S-systems from time-course biochemical data. We demonstrate our methods using three artificial networks designed to simulate different network topologies and behavior. We then end with an application to a real biochemical network by creating a working model for the cadBA system in Escherichia coli. Availability: The source code written in C++ is available at http://www. engg.upd.edu.ph/~naval/bioinformcode.html. All the necessary programs including the required compiler are described in a document archived with the source code. Contact:

Biochemical Systems Theory: A Review

“…Also using evolutionary methods, Mendoza's group developed a method for estimating parameters in GMA systems, based on particle-swarm optimization [279][280][281]; see also [282,283]. e group also tested ant colony optimization [284,285] and simulated annealing [286,287] (see also [288,289]) and provided a benchmark system for comparing different approaches [290]. McKinney and Tian proposed an arti�cial immune system for the same estimation purposes and compared various underlying models [291].…”

Reverse engineering of drug induced DNA damage response signalling pathway reveals dual outcomes of ATM kinase inhibition

“…Over the last two decades the systems biology community has witnessed rapid advances in system identification methods development including the following: linear, loglinear and nonlinear differential equations; artificial neural networks [53]; genetic algorithm [54]; evolutionary optimization with data collocation [55]; interval analysis [56]; alternating regression [57]; parameter estimation for noisy metabolic profiles using newton-flow analysis [58]; simulated annealing [59]; ant colony optimization algorithm for parameter estimation and network inference [60]; substitution of slopes for differentials; dynamic-flux estimation [61]; eigen-vector optimization [62]; transposive and repressive regression method [63][64][65]; and so on. For instance, the following presented ODE-based methods that use dynamic data and steady-state measurements to capture and identify complex systems: Sorribas and Cascante et al [66]; Irvine [67]; Savageau et al [68]; Tominaga and Okamoto [69]; etc.…”

Using Restricted Additive Tree Model for Identifying the Large-Scale Gene Regulatory Networks

“…Kikuchi et al [9] proposed a method to predict not only the network structure but also its dynamics using a unified extension of a Genetic Algorithm and an S-system formalism. Gonzalez et al [10] described how the heuristic optimization technique simulated annealing could be effectively used for estimating the parameters of S-system from time-course biochemical data, and used three artificial networks designed to simulate different network topologies and behavior to demonstrate the method. Palafox et al [11] implemented a variation of particle swarm optimization, called dissipative PSO (DPSO) to optimize the parameters of the popular non-linear differential equation model named S-system in order to infer small-scale gene regulatory networks.…”