Parameter estimation in biochemical pathways: a comparison of global optimization methods.
Here we address the problem of parameter estimation (inverse problem) of nonlinear dynamic biochemical pathways. This problem is stated as a nonlinear programming (NLP) problem subject to nonlinear differential-algebraic constraints. These problems are known to be frequently ill-conditioned and multimodal. Thus, traditional (gradient-based) local optimization methods fail to arrive at satisfactory solutions. To surmount this limitation, the use of several state-of-the-art deterministic and stochastic global optimization methods is explored. A case study considering the estimation of 36 parameters of a nonlinear biochemical dynamic model is taken as a benchmark. Only a certain type of stochastic algorithm, evolution strategies (ES), is able to solve this problem successfully. Although these stochastic methods cannot guarantee global optimality with certainty, their robustness, plus the fact that in inverse problems they have a known lower bound for the cost function, make them the best available candidates.
Moles
Carmen G
CG
Process Engineering Group, Instituto de Investigaciones Marinas (CSIC), 36208 Vigo, Spain.
Mendes
Pedro
P
Banga
Julio R
JR
D000465:Algorithms; D019295:Computational Biology / Q000379:methods* / Q000706:statistics & numerical data; D003198:Computer Simulation; D019143:Evolution, Molecular; D005818:Genetic Engineering / Q000379:methods; D016333:HIV Protease / Q000737:chemistry / Q000235:genetics / Q000378:metabolism; D017320:HIV Protease Inhibitors / Q000737:chemistry; D008956:Models, Chemical; D008967:Molecular Biology / Q000379:methods* / Q000706:statistics & numerical data; D017711:Nonlinear Dynamics; D011237:Predictive Value of Tests; D012984:Software / Q000706:statistics & numerical data; D013269:Stochastic Processes
D017320:HIV Protease Inhibitors; D016333:HIV Protease