1
|
Mikhalychev A, Zhevno K, Vlasenko S, Benediktovitch A, Ulyanenkova T, Ulyanenkov A. Fisher information for optimal planning of X-ray diffraction experiments. J Appl Crystallogr 2021. [DOI: 10.1107/s1600576721009869] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Fisher information is a powerful mathematical tool suitable for quantification of data `informativity' and optimization of the experimental setup and measurement conditions. Here, it is applied to X-ray diffraction and an informational approach to choosing the optimal measurement configuration is proposed. The core idea is maximization of the information which can be extracted from the measured data set by the selected analysis technique, over the sets of accessible reflections and measurement geometries. The developed approach is applied to high-resolution X-ray diffraction measurements and microstructure analysis of multilayer samples, and its efficiency and consistency are demonstrated with the results of more straightforward Monte Carlo simulations.
Collapse
|
2
|
Ferreira AE, Sousa Silva M, Cordeiro C. Metabolic Network Inference from Time Series. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11347-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022] Open
|
3
|
Comprehensive experimental design for chemical engineering processes: A two-layer iterative design approach. Chem Eng Sci 2018. [DOI: 10.1016/j.ces.2018.05.047] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
4
|
Identification of optimal strategies for state transition of complex biological networks. Biochem Soc Trans 2017; 45:1015-1024. [PMID: 28733488 DOI: 10.1042/bst20160419] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 04/23/2017] [Accepted: 05/25/2017] [Indexed: 11/17/2022]
Abstract
Complex biological networks typically contain numerous parameters, and determining feasible strategies for state transition by parameter perturbation is not a trivial task. In the present study, based on dynamical and structural analyses of the biological network, we optimized strategies for controlling variables in a two-node gene regulatory network and a T-cell large granular lymphocyte signaling network associated with blood cancer by using an efficient dynamic optimization method. Optimization revealed the critical value for each decision variable to steer the system from an undesired state into a desired attractor. In addition, the minimum time for the state transition was determined by defining and solving a time-optimal control problem. Moreover, time-dependent variable profiles for state transitions were achieved rather than constant values commonly adopted in previous studies. Furthermore, the optimization method allows multiple controls to be simultaneously adjusted to drive the system out of an undesired attractor. Optimization improved the results of the parameter perturbation method, thus providing a valuable guidance for experimental design.
Collapse
|
5
|
Oguz C, Watson LT, Baumann WT, Tyson JJ. Predicting network modules of cell cycle regulators using relative protein abundance statistics. BMC SYSTEMS BIOLOGY 2017; 11:30. [PMID: 28241833 PMCID: PMC5329933 DOI: 10.1186/s12918-017-0409-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Accepted: 02/17/2017] [Indexed: 11/10/2022]
Abstract
BACKGROUND Parameter estimation in systems biology is typically done by enforcing experimental observations through an objective function as the parameter space of a model is explored by numerical simulations. Past studies have shown that one usually finds a set of "feasible" parameter vectors that fit the available experimental data equally well, and that these alternative vectors can make different predictions under novel experimental conditions. In this study, we characterize the feasible region of a complex model of the budding yeast cell cycle under a large set of discrete experimental constraints in order to test whether the statistical features of relative protein abundance predictions are influenced by the topology of the cell cycle regulatory network. RESULTS Using differential evolution, we generate an ensemble of feasible parameter vectors that reproduce the phenotypes (viable or inviable) of wild-type yeast cells and 110 mutant strains. We use this ensemble to predict the phenotypes of 129 mutant strains for which experimental data is not available. We identify 86 novel mutants that are predicted to be viable and then rank the cell cycle proteins in terms of their contributions to cumulative variability of relative protein abundance predictions. Proteins involved in "regulation of cell size" and "regulation of G1/S transition" contribute most to predictive variability, whereas proteins involved in "positive regulation of transcription involved in exit from mitosis," "mitotic spindle assembly checkpoint" and "negative regulation of cyclin-dependent protein kinase by cyclin degradation" contribute the least. These results suggest that the statistics of these predictions may be generating patterns specific to individual network modules (START, S/G2/M, and EXIT). To test this hypothesis, we develop random forest models for predicting the network modules of cell cycle regulators using relative abundance statistics as model inputs. Predictive performance is assessed by the areas under receiver operating characteristics curves (AUC). Our models generate an AUC range of 0.83-0.87 as opposed to randomized models with AUC values around 0.50. CONCLUSIONS By using differential evolution and random forest modeling, we show that the model prediction statistics generate distinct network module-specific patterns within the cell cycle network.
Collapse
Affiliation(s)
- Cihan Oguz
- Department of Biological Sciences, Virginia Tech, Blacksburg VA, 24061, USA.
| | - Layne T Watson
- Department of Computer Science, Virginia Tech, Blacksburg VA, 24061, USA.,Department of Mathematics, Virginia Tech, Blacksburg VA, 24061, USA.,Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg VA, 24061, USA
| | - William T Baumann
- Department of Electrical and Computer Engineering, Virginia Tech, Blacksburg VA, 24061, USA
| | - John J Tyson
- Department of Biological Sciences, Virginia Tech, Blacksburg VA, 24061, USA
| |
Collapse
|
6
|
Mathematical modelling of clostridial acetone-butanol-ethanol fermentation. Appl Microbiol Biotechnol 2017; 101:2251-2271. [PMID: 28210797 PMCID: PMC5320022 DOI: 10.1007/s00253-017-8137-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Revised: 01/14/2017] [Accepted: 01/16/2017] [Indexed: 12/24/2022]
Abstract
Clostridial acetone-butanol-ethanol (ABE) fermentation features a remarkable shift in the cellular metabolic activity from acid formation, acidogenesis, to the production of industrial-relevant solvents, solventogensis. In recent decades, mathematical models have been employed to elucidate the complex interlinked regulation and conditions that determine these two distinct metabolic states and govern the transition between them. In this review, we discuss these models with a focus on the mechanisms controlling intra- and extracellular changes between acidogenesis and solventogenesis. In particular, we critically evaluate underlying model assumptions and predictions in the light of current experimental knowledge. Towards this end, we briefly introduce key ideas and assumptions applied in the discussed modelling approaches, but waive a comprehensive mathematical presentation. We distinguish between structural and dynamical models, which will be discussed in their chronological order to illustrate how new biological information facilitates the ‘evolution’ of mathematical models. Mathematical models and their analysis have significantly contributed to our knowledge of ABE fermentation and the underlying regulatory network which spans all levels of biological organization. However, the ties between the different levels of cellular regulation are not well understood. Furthermore, contradictory experimental and theoretical results challenge our current notion of ABE metabolic network structure. Thus, clostridial ABE fermentation still poses theoretical as well as experimental challenges which are best approached in close collaboration between modellers and experimentalists.
Collapse
|
7
|
Efficient Optimization of Stimuli for Model-Based Design of Experiments to Resolve Dynamical Uncertainty. PLoS Comput Biol 2015; 11:e1004488. [PMID: 26379275 PMCID: PMC4574939 DOI: 10.1371/journal.pcbi.1004488] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2014] [Accepted: 08/05/2015] [Indexed: 11/19/2022] Open
Abstract
This model-based design of experiments (MBDOE) method determines the input magnitudes of an experimental stimuli to apply and the associated measurements that should be taken to optimally constrain the uncertain dynamics of a biological system under study. The ideal global solution for this experiment design problem is generally computationally intractable because of parametric uncertainties in the mathematical model of the biological system. Others have addressed this issue by limiting the solution to a local estimate of the model parameters. Here we present an approach that is independent of the local parameter constraint. This approach is made computationally efficient and tractable by the use of: (1) sparse grid interpolation that approximates the biological system dynamics, (2) representative parameters that uniformly represent the data-consistent dynamical space, and (3) probability weights of the represented experimentally distinguishable dynamics. Our approach identifies data-consistent representative parameters using sparse grid interpolants, constructs the optimal input sequence from a greedy search, and defines the associated optimal measurements using a scenario tree. We explore the optimality of this MBDOE algorithm using a 3-dimensional Hes1 model and a 19-dimensional T-cell receptor model. The 19-dimensional T-cell model also demonstrates the MBDOE algorithm’s scalability to higher dimensions. In both cases, the dynamical uncertainty region that bounds the trajectories of the target system states were reduced by as much as 86% and 99% respectively after completing the designed experiments in silico. Our results suggest that for resolving dynamical uncertainty, the ability to design an input sequence paired with its associated measurements is particularly important when limited by the number of measurements. Many mathematical models that have been developed for biological systems are limited because the complex systems are not well understood, the parameters are not known, and available data is limited and noisy. On the other hand, experiments to support model development are limited in terms of costs and time, feasible inputs and feasible measurements. MBDOE combines the mathematical models with experiment design to strategically design optimal experiments to obtain data that will contribute to the understanding of the systems. Our approach extends current capabilities of existing MBDOE techniques to make them more useful for scientists to resolve the trajectories of the system under study. It identifies the optimal conditions for stimuli and measurements that yield the most information about the system given the practical limitations. Exploration of the input space is not a trivial extension to MBDOE methods used for determining optimal measurements due to the nonlinear nature of many biological system models. The exploration of the system dynamics elicited by different inputs requires a computationally efficient and tractable approach. Our approach plans optimal experiments to reduce dynamical uncertainty in the output of selected target states of the biological system.
Collapse
|
8
|
Telen D, Van Riet N, Logist F, Van Impe J. A differentiable reformulation for E-optimal design of experiments in nonlinear dynamic biosystems. Math Biosci 2015; 264:1-7. [PMID: 25725122 DOI: 10.1016/j.mbs.2015.02.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2014] [Revised: 02/10/2015] [Accepted: 02/11/2015] [Indexed: 11/26/2022]
Abstract
Informative experiments are highly valuable for estimating parameters in nonlinear dynamic bioprocesses. Techniques for optimal experiment design ensure the systematic design of such informative experiments. The E-criterion which can be used as objective function in optimal experiment design requires the maximization of the smallest eigenvalue of the Fisher information matrix. However, one problem with the minimal eigenvalue function is that it can be nondifferentiable. In addition, no closed form expression exists for the computation of eigenvalues of a matrix larger than a 4 by 4 one. As eigenvalues are normally computed with iterative methods, state-of-the-art optimal control solvers are not able to exploit automatic differentiation to compute the derivatives with respect to the decision variables. In the current paper a reformulation strategy from the field of convex optimization is suggested to circumvent these difficulties. This reformulation requires the inclusion of a matrix inequality constraint involving positive semidefiniteness. In this paper, this positive semidefiniteness constraint is imposed via Sylverster's criterion. As a result the maximization of the minimum eigenvalue function can be formulated in standard optimal control solvers through the addition of nonlinear constraints. The presented methodology is successfully illustrated with a case study from the field of predictive microbiology.
Collapse
Affiliation(s)
- Dries Telen
- KU Leuven, Chemical Engineering Department, BioTeC & OPTEC, W. de Croylaan 46, 3001 Leuven, Belgium
| | - Nick Van Riet
- KU Leuven, Chemical Engineering Department, BioTeC & OPTEC, W. de Croylaan 46, 3001 Leuven, Belgium
| | - Flip Logist
- KU Leuven, Chemical Engineering Department, BioTeC & OPTEC, W. de Croylaan 46, 3001 Leuven, Belgium
| | - Jan Van Impe
- KU Leuven, Chemical Engineering Department, BioTeC & OPTEC, W. de Croylaan 46, 3001 Leuven, Belgium.
| |
Collapse
|
9
|
Yu H, Yue H, Halling P. Optimal Experimental Design for an Enzymatic Biodiesel Production System. ACTA ACUST UNITED AC 2015. [DOI: 10.1016/j.ifacol.2015.09.141] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
10
|
He Y, Xu Y, Geng Z, Zhu Q. Soft sensor of chemical processes with large numbers of input parameters using auto-associative hierarchical neural network. Chin J Chem Eng 2015. [DOI: 10.1016/j.cjche.2014.10.004] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
11
|
Lang M, Summers S, Stelling J. Cutting the wires: modularization of cellular networks for experimental design. Biophys J 2014; 106:321-31. [PMID: 24411264 DOI: 10.1016/j.bpj.2013.11.2960] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2013] [Revised: 11/11/2013] [Accepted: 11/12/2013] [Indexed: 01/02/2023] Open
Abstract
Understanding naturally evolved cellular networks requires the consecutive identification and revision of the interactions between relevant molecular species. In this process, initially often simplified and incomplete networks are extended by integrating new reactions or whole subnetworks to increase consistency between model predictions and new measurement data. However, increased consistency with experimental data alone is not sufficient to show the existence of biomolecular interactions, because the interplay of different potential extensions might lead to overall similar dynamics. Here, we present a graph-based modularization approach to facilitate the design of experiments targeted at independently validating the existence of several potential network extensions. Our method is based on selecting the outputs to measure during an experiment, such that each potential network extension becomes virtually insulated from all others during data analysis. Each output defines a module that only depends on one hypothetical network extension, and all other outputs act as virtual inputs to achieve insulation. Given appropriate experimental time-series measurements of the outputs, our modules can be analyzed, simulated, and compared to the experimental data separately. Our approach exemplifies the close relationship between structural systems identification and modularization, an interplay that promises development of related approaches in the future.
Collapse
Affiliation(s)
- Moritz Lang
- Department of Biosystems Science and Engineering, ETH Zürich, and Swiss Institute of Bioinformatics, Basel, Switzerland.
| | - Sean Summers
- Automatic Control Laboratory, ETH Zürich, Zurich, Switzerland
| | - Jörg Stelling
- Department of Biosystems Science and Engineering, ETH Zürich, and Swiss Institute of Bioinformatics, Basel, Switzerland
| |
Collapse
|
12
|
Borchers S, Freund S, Rath A, Streif S, Reichl U, Findeisen R. Identification of growth phases and influencing factors in cultivations with AGE1.HN cells using set-based methods. PLoS One 2013; 8:e68124. [PMID: 23936299 PMCID: PMC3732265 DOI: 10.1371/journal.pone.0068124] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2012] [Accepted: 05/30/2013] [Indexed: 01/16/2023] Open
Abstract
Production of bio-pharmaceuticals in cell culture, such as mammalian cells, is challenging. Mathematical models can provide support to the analysis, optimization, and the operation of production processes. In particular, unstructured models are suited for these purposes, since they can be tailored to particular process conditions. To this end, growth phases and the most relevant factors influencing cell growth and product formation have to be identified. Due to noisy and erroneous experimental data, unknown kinetic parameters, and the large number of combinations of influencing factors, currently there are only limited structured approaches to tackle these issues. We outline a structured set-based approach to identify different growth phases and the factors influencing cell growth and metabolism. To this end, measurement uncertainties are taken explicitly into account to bound the time-dependent specific growth rate based on the observed increase of the cell concentration. Based on the bounds on the specific growth rate, we can identify qualitatively different growth phases and (in-)validate hypotheses on the factors influencing cell growth and metabolism. We apply the approach to a mammalian suspension cell line (AGE1.HN). We show that growth in batch culture can be divided into two main growth phases. The initial phase is characterized by exponential growth dynamics, which can be described consistently by a relatively simple unstructured and segregated model. The subsequent phase is characterized by a decrease in the specific growth rate, which, as shown, results from substrate limitation and the pH of the medium. An extended model is provided which describes the observed dynamics of cell growth and main metabolites, and the corresponding kinetic parameters as well as their confidence intervals are estimated. The study is complemented by an uncertainty and outlier analysis. Overall, we demonstrate utility of set-based methods for analyzing cell growth and metabolism under conditions of uncertainty.
Collapse
Affiliation(s)
- Steffen Borchers
- Institute for Systems Theory and Automatic Control, Otto-von-Guericke University, Magdeburg, Germany
- International Max Planck Research School, Magdeburg, Germany
| | - Susann Freund
- International Max Planck Research School, Magdeburg, Germany
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
- Institute of Process Engineering, Otto-von-Guericke University, Magdeburg, Germany
| | - Alexander Rath
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
| | - Stefan Streif
- Institute for Systems Theory and Automatic Control, Otto-von-Guericke University, Magdeburg, Germany
| | - Udo Reichl
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
- Institute of Process Engineering, Otto-von-Guericke University, Magdeburg, Germany
| | - Rolf Findeisen
- Institute for Systems Theory and Automatic Control, Otto-von-Guericke University, Magdeburg, Germany
| |
Collapse
|
13
|
Maheshwari V, Rangaiah GP, Samavedham L. Multiobjective Framework for Model-based Design of Experiments to Improve Parameter Precision and Minimize Parameter Correlation. Ind Eng Chem Res 2013. [DOI: 10.1021/ie400133m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Vaibhav Maheshwari
- Department of Chemical & Biomolecular Engineering, National University of Singapore, Singapore 117576
| | - Gade Pandu Rangaiah
- Department of Chemical & Biomolecular Engineering, National University of Singapore, Singapore 117576
| | | |
Collapse
|
14
|
Chakrabarty A, Buzzard GT, Rundell AE. Model-based design of experiments for cellular processes. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2013; 5:181-203. [PMID: 23293047 DOI: 10.1002/wsbm.1204] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Ankush Chakrabarty
- School of Electrical and Computer Engineering, Purdue University, West Lafayette, IN, USA
| | | | | |
Collapse
|
15
|
Sun J, Garibaldi JM, Hodgman C. Parameter estimation using meta-heuristics in systems biology: a comprehensive review. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2012; 9:185-202. [PMID: 21464505 DOI: 10.1109/tcbb.2011.63] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
This paper gives a comprehensive review of the application of meta-heuristics to optimization problems in systems biology, mainly focussing on the parameter estimation problem (also called the inverse problem or model calibration). It is intended for either the system biologist who wishes to learn more about the various optimization techniques available and/or the meta-heuristic optimizer who is interested in applying such techniques to problems in systems biology. First, the parameter estimation problems emerging from different areas of systems biology are described from the point of view of machine learning. Brief descriptions of various meta-heuristics developed for these problems follow, along with outlines of their advantages and disadvantages. Several important issues in applying meta-heuristics to the systems biology modelling problem are addressed, including the reliability and identifiability of model parameters, optimal design of experiments, and so on. Finally, we highlight some possible future research directions in this field.
Collapse
|
16
|
Multi-objective optimisation approach to optimal experiment design in dynamic bioprocesses using ACADO toolkit. ACTA ACUST UNITED AC 2011. [DOI: 10.1016/b978-0-444-53711-9.50093-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
17
|
|
18
|
Jiménez-Hornero JE, Santos-Dueñas IM, García-García I. Optimization of biotechnological processes. The acetic acid fermentation. Part II: Practical identifiability analysis and parameter estimation. Biochem Eng J 2009. [DOI: 10.1016/j.bej.2009.01.010] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
19
|
Abstract
Mathematical models are central in systems biology and provide new ways to understand the function of biological systems, helping in the generation of novel and testable hypotheses, and supporting a rational framework for possible ways of intervention, like in e.g. genetic engineering, drug development or treatment of diseases. Since the amount and quality of experimental 'omics' data continue to increase rapidly, there is great need for methods for proper model building which can handle this complexity. In the present chapter we review two key steps of the model building process, namely parameter estimation (model calibration) and optimal experimental design. Parameter estimation aims to find the unknown parameters of the model which give the best fit to a set of experimental data. Optimal experimental design aims to devise the dynamic experiments which provide the maximum information content for subsequent non-linear model identification, estimation and/or discrimination. We place emphasis on the need for robust global optimization methods for proper solution of these problems, and we present a motivating example considering a cell signalling model.
Collapse
|
20
|
Yue H, Brown M, He F, Jia J, Kell DB. Sensitivity analysis and robust experimental design of a signal transduction pathway system. INT J CHEM KINET 2008. [DOI: 10.1002/kin.20369] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
21
|
Balsa-Canto E, Alonso AA, Banga JR. Computational procedures for optimal experimental design in biological systems. IET Syst Biol 2008; 2:163-72. [PMID: 18681746 DOI: 10.1049/iet-syb:20070069] [Citation(s) in RCA: 97] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Mathematical models of complex biological systems, such as metabolic or cell-signalling pathways, usually consist of sets of nonlinear ordinary differential equations which depend on several non-measurable parameters that can be hopefully estimated by fitting the model to experimental data. However, the success of this fitting is largely conditioned by the quantity and quality of data. Optimal experimental design (OED) aims to design the scheme of actuations and measurements which will result in data sets with the maximum amount and/or quality of information for the subsequent model calibration. New methods and computational procedures for OED in the context of biological systems are presented. The OED problem is formulated as a general dynamic optimisation problem where the time-dependent stimuli profiles, the location of sampling times, the duration of the experiments and the initial conditions are regarded as design variables. Its solution is approached using the control vector parameterisation method. Since the resultant nonlinear optimisation problem is in most of the cases non-convex, the use of a robust global nonlinear programming solver is proposed. For the sake of comparing among different experimental schemes, a Monte-Carlo-based identifiability analysis is then suggested. The applicability and advantages of the proposed techniques are illustrated by considering an example related to a cell-signalling pathway.
Collapse
|
22
|
Franceschini G, Macchietto S. Model-based design of experiments for parameter precision: State of the art. Chem Eng Sci 2008. [DOI: 10.1016/j.ces.2007.11.034] [Citation(s) in RCA: 297] [Impact Index Per Article: 18.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
23
|
Banga JR. Optimization in computational systems biology. BMC SYSTEMS BIOLOGY 2008; 2:47. [PMID: 18507829 PMCID: PMC2435524 DOI: 10.1186/1752-0509-2-47] [Citation(s) in RCA: 136] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/21/2008] [Accepted: 05/28/2008] [Indexed: 12/05/2022]
Abstract
Optimization aims to make a system or design as effective or functional as possible. Mathematical optimization methods are widely used in engineering, economics and science. This commentary is focused on applications of mathematical optimization in computational systems biology. Examples are given where optimization methods are used for topics ranging from model building and optimal experimental design to metabolic engineering and synthetic biology. Finally, several perspectives for future research are outlined.
Collapse
Affiliation(s)
- Julio R Banga
- Instituto de Investigaciones Marinas, CSIC (Spanish Council for Scientific Research), C/Eduardo Cabello 6, 36208 Vigo, Spain.
| |
Collapse
|
24
|
BALSA-CANTO E, ALONSO A, BANGA J. COMPUTING OPTIMAL DYNAMIC EXPERIMENTS FOR MODEL CALIBRATION IN PREDICTIVE MICROBIOLOGY. J FOOD PROCESS ENG 2008. [DOI: 10.1111/j.1745-4530.2007.00147.x] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
25
|
Banga JR, Balsa-Canto E, Alonso AA. Quality and Safety Models and Optimization as Part of Computer-Integrated Manufacturing. Compr Rev Food Sci Food Saf 2008. [DOI: 10.1111/j.1541-4337.2007.00023.x] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
26
|
Scheerlinck N, Berhane NH, Moles CG, Banga JR, Nicolaï BM. Optimal dynamic heat generation profiles for simultaneous estimation of thermal food properties using a hotwire probe: Computation, implementation and validation. J FOOD ENG 2008. [DOI: 10.1016/j.jfoodeng.2007.05.019] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
27
|
Identifiability and robust parameter estimation in food process modeling: Application to a drying model. J FOOD ENG 2007. [DOI: 10.1016/j.jfoodeng.2007.03.023] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
28
|
Schittkowski K. Experimental Design Tools for Ordinary and Algebraic Differential Equations. Ind Eng Chem Res 2007. [DOI: 10.1021/ie0703742] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- K. Schittkowski
- Department of Computer Science, University of Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
29
|
Balsa-Canto E, Rodriguez-Fernandez M, Banga JR. Optimal design of dynamic experiments for improved estimation of kinetic parameters of thermal degradation. J FOOD ENG 2007. [DOI: 10.1016/j.jfoodeng.2007.02.006] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
30
|
Rodriguez-Fernandez M, Mendes P, Banga JR. A hybrid approach for efficient and robust parameter estimation in biochemical pathways. Biosystems 2005; 83:248-65. [PMID: 16236429 DOI: 10.1016/j.biosystems.2005.06.016] [Citation(s) in RCA: 158] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2004] [Revised: 04/21/2005] [Accepted: 06/13/2005] [Indexed: 10/25/2022]
Abstract
Developing suitable dynamic models of biochemical pathways is a key issue in Systems Biology. Predictive models for cells or whole organisms could ultimately lead to model-based predictive and/or preventive medicine. Parameter estimation (i.e. model calibration) in these dynamic models is therefore a critical problem. In a recent contribution [Moles, C.G., Mendes, P., Banga, J.R., 2003b. Parameter estimation in biochemical pathways: a comparison of global optimisation methods. Genome Res. 13, 2467-2474], the challenging nature of such inverse problems was highlighted considering a benchmark problem, and concluding that only a certain type of stochastic global optimisation method, Evolution Strategies (ES), was able to solve it successfully, although at a rather large computational cost. In this new contribution, we present a new integrated optimisation methodology with a number of very significant improvements: (i) computation time is reduced by one order of magnitude by means of a hybrid method which increases efficiency while guaranteeing robustness, (ii) measurement noise (errors) and partial observations are handled adequately, (iii) automatic testing of identifiability of the model (both local and practical) is included and (iv) the information content of the experiments is evaluated via the Fisher information matrix, with subsequent application to design of new optimal experiments through dynamic optimisation.
Collapse
Affiliation(s)
- Maria Rodriguez-Fernandez
- Process Engineering Group, IIM-CSIC, Spanish Council for Scientific Research, C/Eduardo Cabello 6, 36208 Vigo, Spain
| | | | | |
Collapse
|
31
|
Gadkar KG, Gunawan R, Doyle FJ. Iterative approach to model identification of biological networks. BMC Bioinformatics 2005; 6:155. [PMID: 15967022 PMCID: PMC1189077 DOI: 10.1186/1471-2105-6-155] [Citation(s) in RCA: 107] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2005] [Accepted: 06/20/2005] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Recent advances in molecular biology techniques provide an opportunity for developing detailed mathematical models of biological processes. An iterative scheme is introduced for model identification using available system knowledge and experimental measurements. RESULTS The scheme includes a state regulator algorithm that provides estimates of all system unknowns (concentrations of the system components and the reaction rates of their inter-conversion). The full system information is used for estimation of the model parameters. An optimal experiment design using the parameter identifiability and D-optimality criteria is formulated to provide "rich" experimental data for maximizing the accuracy of the parameter estimates in subsequent iterations. The importance of model identifiability tests for optimal measurement selection is also considered. The iterative scheme is tested on a model for the caspase function in apoptosis where it is demonstrated that model accuracy improves with each iteration. Optimal experiment design was determined to be critical for model identification. CONCLUSION The proposed algorithm has general application to modeling a wide range of cellular processes, which include gene regulation networks, signal transduction and metabolic networks.
Collapse
Affiliation(s)
- Kapil G Gadkar
- Department of Chemical Engineering, University of California Santa Barbara, CA, USA
| | - Rudiyanto Gunawan
- Department of Chemical Engineering, University of California Santa Barbara, CA, USA
| | - Francis J Doyle
- Department of Chemical Engineering, University of California Santa Barbara, CA, USA
| |
Collapse
|
32
|
|
33
|
Hajipour AR, Malakoutikhah M. A CONVENIENT METHOD FOR PREPARATION OF TRIAZOLINEDIONES. ORG PREP PROCED INT 2004. [DOI: 10.1080/00304940409356632] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
34
|
Kremling A, Fischer S, Gadkar K, Doyle FJ, Sauter T, Bullinger E, Allgöwer F, Gilles ED. A benchmark for methods in reverse engineering and model discrimination: problem formulation and solutions. Genome Res 2004; 14:1773-85. [PMID: 15342560 PMCID: PMC515324 DOI: 10.1101/gr.1226004] [Citation(s) in RCA: 95] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2003] [Accepted: 05/18/2004] [Indexed: 11/24/2022]
Abstract
A benchmark problem is described for the reconstruction and analysis of biochemical networks given sampled experimental data. The growth of the organisms is described in a bioreactor in which one substrate is fed into the reactor with a given feed rate and feed concentration. Measurements for some intracellular components are provided representing a small biochemical network. Problems of reverse engineering, parameter estimation, and identifiability are addressed. The contribution mainly focuses on the problem of model discrimination. If two or more model variants describe the available experimental data, a new experiment must be designed to discriminate between the hypothetical models. For the problem presented, the feed rate and feed concentration of a bioreactor system are available as control inputs. To verify calculated input profiles an interactive Web site (http://www.sysbio.de/projects/benchmark/) is provided. Several solutions based on linear and nonlinear models are discussed.
Collapse
Affiliation(s)
- Andreas Kremling
- Systems Biology Group, Max-Planck-Institut für Dynamik komplexer technischer Systeme, 39106 Magdeburg, Germany.
| | | | | | | | | | | | | | | |
Collapse
|
35
|
Banga JR, Moles CG, Alonso AA. Global Optimization of Bioprocesses using Stochastic and Hybrid Methods. NONCONVEX OPTIMIZATION AND ITS APPLICATIONS 2004. [DOI: 10.1007/978-1-4613-0251-3_3] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|