MORE: mixed optimization for reverse engineering--an application to modeling biological networks response via sparse systems of nonlinear differential equations

A parallel metaheuristic for large mixed-integer dynamic optimization problems, with applications in computational biology

Background

We consider a general class of global optimization problems dealing with nonlinear dynamic models. Although this class is relevant to many areas of science and engineering, here we are interested in applying this framework to the reverse engineering problem in computational systems biology, which yields very large mixed-integer dynamic optimization (MIDO) problems. In particular, we consider the framework of logic-based ordinary differential equations (ODEs).

Methods

We present saCeSS2, a parallel method for the solution of this class of problems. This method is based on an parallel cooperative scatter search metaheuristic, with new mechanisms of self-adaptation and specific extensions to handle large mixed-integer problems. We have paid special attention to the avoidance of convergence stagnation using adaptive cooperation strategies tailored to this class of problems.

Results

We illustrate its performance with a set of three very challenging case studies from the domain of dynamic modelling of cell signaling. The simpler case study considers a synthetic signaling pathway and has 84 continuous and 34 binary decision variables. A second case study considers the dynamic modeling of signaling in liver cancer using high-throughput data, and has 135 continuous and 109 binaries decision variables. The third case study is an extremely difficult problem related with breast cancer, involving 690 continuous and 138 binary decision variables. We report computational results obtained in different infrastructures, including a local cluster, a large supercomputer and a public cloud platform. Interestingly, the results show how the cooperation of individual parallel searches modifies the systemic properties of the sequential algorithm, achieving superlinear speedups compared to an individual search (e.g. speedups of 15 with 10 cores), and significantly improving (above a 60%) the performance with respect to a non-cooperative parallel scheme. The scalability of the method is also good (tests were performed using up to 300 cores).

Conclusions

These results demonstrate that saCeSS2 can be used to successfully reverse engineer large dynamic models of complex biological pathways. Further, these results open up new possibilities for other MIDO-based large-scale applications in the life sciences such as metabolic engineering, synthetic biology, drug scheduling.

Reverse engineering of logic-based differential equation models using a mixed-integer dynamic optimization approach

MOTIVATION

Systems biology models can be used to test new hypotheses formulated on the basis of previous knowledge or new experimental data, contradictory with a previously existing model. New hypotheses often come in the shape of a set of possible regulatory mechanisms. This search is usually not limited to finding a single regulation link, but rather a combination of links subject to great uncertainty or no information about the kinetic parameters.

RESULTS

In this work, we combine a logic-based formalism, to describe all the possible regulatory structures for a given dynamic model of a pathway, with mixed-integer dynamic optimization (MIDO). This framework aims to simultaneously identify the regulatory structure (represented by binary parameters) and the real-valued parameters that are consistent with the available experimental data, resulting in a logic-based differential equation model. The alternative to this would be to perform real-valued parameter estimation for each possible model structure, which is not tractable for models of the size presented in this work. The performance of the method presented here is illustrated with several case studies: a synthetic pathway problem of signaling regulation, a two-component signal transduction pathway in bacterial homeostasis, and a signaling network in liver cancer cells.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

julio@iim.csic.es or saezrodriguez@ebi.ac.uk.

Detection of severe respiratory disease epidemic outbreaks by CUSUM-based overcrowd-severe-respiratory-disease-index model

A severe respiratory disease epidemic outbreak correlates with a high demand of specific supplies and specialized personnel to hold it back in a wide region or set of regions; these supplies would be beds, storage areas, hemodynamic monitors, and mechanical ventilators, as well as physicians, respiratory technicians, and specialized nurses. We describe an online cumulative sum based model named Overcrowd-Severe-Respiratory-Disease-Index based on the Modified Overcrowd Index that simultaneously monitors and informs the demand of those supplies and personnel in a healthcare network generating early warnings of severe respiratory disease epidemic outbreaks through the interpretation of such variables. A post hoc historical archive is generated, helping physicians in charge to improve the transit and future allocation of supplies in the entire hospital network during the outbreak. The model was thoroughly verified in a virtual scenario, generating multiple epidemic outbreaks in a 6-year span for a 13-hospital network. When it was superimposed over the H1N1 influenza outbreak census (2008-2010) taken by the National Institute of Medical Sciences and Nutrition Salvador Zubiran in Mexico City, it showed that it is an effective algorithm to notify early warnings of severe respiratory disease epidemic outbreaks with a minimal rate of false alerts.