ISRES+: an improved evolutionary strategy for function minimization to estimate the free parameters of systems biology models

Abstract

Motivation Mathematical models in systems biology help generate hypotheses, guide experimental design, and infer the dynamics of gene regulatory networks. These models are characterized by phenomenological or mechanistic parameters, which are typically hard to measure. Therefore, efficient parameter estimation is central to model development. Global optimization techniques, such as evolutionary algorithms (EA), are applied to estimate model parameters by inverse modeling, i.e., calibrating models by minimizing a function that evaluates a measure of the error between model predictions and experimental data. EAs estimate model parameters “fittest individuals” by generating a large population of individuals using strategies like recombination and mutation over multiple “generations”. Typically, only a few individuals from each generation are used to create new individuals in the next generation. Improved Evolutionary Strategy by Stochastic Ranking (ISRES), proposed by Runnarson and Yao, is one such EA that is widely used in systems biology to estimate parameters. ISRES uses information at most from a pair of individuals in any generation to create a new population to minimize the error. In this paper we propose an efficient evolutionary strategy, ISRES+, which builds on ISRES by combining information from all individuals across the population and across all generations to develop a better understanding of the fitness landscape. Results ISRES+ uses the additional information generated by the algorithm during evolution to approximate the local neighborhood around the best-fit individual using linear least squares fits in one and two dimensions, enabling efficient parameter estimation. ISRES+ outperforms ISRES and results in fitter individuals with a tighter distribution over multiple runs, such that a typical run of ISRES+ estimates parameters with a higher goodness-of-fit compared to ISRES. Availability Algorithm and implementation: Github - https://github.com/gtreeves/isres-plus-bandodkar-2022. Supplementary information Supplementary data are available at Bioinformatics online.

Spatiotemporal dynamics of NF-κB/Dorsal inhibitor IκBα/Cactus inDrosophilablastoderm embryos

“…Therefore, we developed a material balance model of Cact-LT/GFP dynamics and constrained it with our quantitative measurements of the nuclear import/export rates, mol fraction of Cact-LT/GFP and the NCR of intensities in 1x1x, 2xGFP, 2xLT, and control embryos, as well as an estimate of the total GFP concentration in the embryo (see Supplementary Information). We used an evolutionary optimization strategy, ISRES+ (Bandodkar et al, 2023; see Supplementary Information), which results in distributions of inferred parameters (Fig. 5).…”

Optimal Performance Objectives in the Highly Conserved Bone Morphogenetic Protein Signaling Pathway

“…We implemented Improved Stochastic Ranking Evolutionary Strategy plus (ISRES+) [51] to estimate the 10 variable parameters in the model. The model was fit to nuclear EGFP-Smad1 concentrations in HaCat cells in the presence of TGF-β signaling and later when an inhibitor of the Smad pathway was externally added.…”

Identifying effective evolutionary strategies for uncovering reaction kinetic parameters under the effect of measurement noises

“…Indeed, the efficacy of ESs has caught the attention of systems biologists. For example, they have been used as to the explanative modelling of circadian clock [39], development and pattern formation [40–43], and iron metabolism [44]. More broadly, they have been employed to help discern network mechanisms [45], reverse engineer pathways [46, 47], and infer signaling dynamics [48].…”