Parameter Estimation for Linear Compartmental Models—A Sensitivity Analysis Approach

Linear compartmental models are useful, explanatory tools, that have been widely used to represent the dynamic behavior of complex biological systems. This paper addresses the problem of the numerical identification of such models, i.e., the estimation of the parameter values that will generate predictions closest to experimental observations. Traditional local optimization techniques find it difficult to arrive at satisfactory solutions to such a parameter estimation problem, especially when the number of parameters is large and/or few data are available from experiments. We present herewith a method based on a prior sensitivity analysis, which enables division of a large optimization problem into several smaller and simpler subproblems, on which only sensitive parameters are estimated, before the whole optimization problem is tackled from starting points that are already close to the optimum values. This method has been applied successfully to a linear 13-compartment, 21-parameter model describing the postprandial metabolism of dietary nitrogen in humans. The effectiveness of the method has been demonstrated using simulated and real data obtained in the intestine, blood and urine of healthy humans after the ingestion of a [15N]-labeled protein meal.

[1]  M. Elia,et al.  Amino acid metabolism in muscle and in the whole body of man before and after ingestion of a single mixed meal. , 1989, The American journal of clinical nutrition.

[2]  C. Cobelli,et al.  A model to measure insulin effects on glucose transport and phosphorylation in muscle: a three-tracer study. , 1996, The American journal of physiology.

[3]  Klaus Schittkowski,et al.  Numerical Data Fitting in Dynamical Systems: A Practical Introduction with Applications and Software , 2002 .

[4]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[5]  J. Huneau,et al.  A high-protein meal exceeds anabolic and catabolic capacities in rats adapted to a normal protein diet. , 2000, The Journal of nutrition.

[6]  R. Benamouzig,et al.  Increasing habitual protein intake accentuates differences in postprandial dietary nitrogen utilization between protein sources in humans. , 2003, The Journal of nutrition.

[7]  D. Tomé,et al.  DEVELOPMENT AND CALIBRATION OF A MODELING TOOL FOR THE ANALYSIS OF CLINICAL DATA IN HUMAN NUTRITION , 2005 .

[8]  J. Jacquez Compartmental analysis in biology and medicine , 1985 .

[9]  W. Griffith,et al.  An optimization strategy for a biokinetic model of inhaled radionuclides. , 1991, Fundamental and applied toxicology : official journal of the Society of Toxicology.

[10]  R. Fletcher Practical Methods of Optimization , 1988 .

[11]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[12]  N. Krebs,et al.  Human zinc metabolism: advances in the modeling of stable isotope data. , 1998, Advances in experimental medicine and biology.

[13]  Maria Rodriguez-Fernandez,et al.  A hybrid approach for efficient and robust parameter estimation in biochemical pathways. , 2006, Bio Systems.

[14]  James R. Wilson,et al.  A revised simplex search procedure for stochastic simulation response-surface optimization , 1998, WSC '98.

[15]  N. Deutz,et al.  Infusion of soy and casein protein meals affects interorgan amino acid metabolism and urea kinetics differently in pigs. , 1998, The Journal of nutrition.

[16]  J. Huneau,et al.  Energy nutrients modulate the splanchnic sequestration of dietary nitrogen in humans: a compartmental analysis. , 2001, American journal of physiology. Endocrinology and metabolism.

[17]  C Cobelli,et al.  Compartmental models: theory and practice using the SAAM II software system. , 1998, Advances in experimental medicine and biology.

[18]  H. Rabitz,et al.  Efficient input-output model representations , 1999 .

[19]  G. Biolo,et al.  Leucine and phenylalanine kinetics during mixed meal ingestion: a multiple tracer approach. , 1992, The American journal of physiology.

[20]  P. Pardalos,et al.  Handbook of global optimization , 1995 .

[21]  C Cobelli,et al.  Compartmental model of leucine kinetics in humans. , 1991, The American journal of physiology.

[22]  Giovanni Sparacino,et al.  Numerical non-identifiability regions of the minimal model of glucose kinetics: superiority of Bayesian estimation. , 2003, Mathematical biosciences.

[23]  A. Attie,et al.  The role of the LDL receptor in apolipoprotein B secretion. , 2000, The Journal of clinical investigation.

[24]  G. Biolo,et al.  Harry M. Vars Research Award. A new model to determine in vivo the relationship between amino acid transmembrane transport and protein kinetics in muscle. , 1992, JPEN. Journal of parenteral and enteral nutrition.

[25]  C Cobelli,et al.  A compartmental model of zinc metabolism in healthy women using oral and intravenous stable isotope tracers. , 1997, The American journal of clinical nutrition.

[26]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[27]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[28]  C. Gaudichon,et al.  Contribution of plasma proteins to splanchnic and total anabolic utilization of dietary nitrogen in humans. , 2003, American journal of physiology. Endocrinology and metabolism.

[29]  C Cobelli,et al.  Using what is accessible to measure that which is not: necessity of model of system. , 1998, Metabolism: clinical and experimental.

[30]  Stefano Tarantola,et al.  Sensitivity Analysis as an Ingredient of Modeling , 2000 .

[31]  C. Gaudichon,et al.  Approaches to quantifying protein metabolism in response to nutrient ingestion. , 2002, The Journal of nutrition.

[32]  R. Benamouzig,et al.  The gastro-ileal digestion of 15N-labelled pea nitrogen in adult humans , 1996, British Journal of Nutrition.

[33]  Feng-Sheng Wang,et al.  Evolutionary optimization with data collocation for reverse engineering of biological networks , 2005, Bioinform..

[34]  M. Auletta,et al.  Splanchnic and leg substrate exchange after ingestion of a natural mixed meal in humans. , 1999, Diabetes.

[35]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[36]  N. Akaike,et al.  Internal perfusion studies demonstrating GABA-induced chloride responses in frog primary afferent neurons. , 1984, The American journal of physiology.

[37]  A. Saltelli,et al.  Sensitivity Anaysis as an Ingredient of Modeling , 2000 .

[38]  D M Foster,et al.  Developing and testing integrated multicompartment models to describe a single-input multiple-output study using the SAAM II software system. , 1998, Advances in experimental medicine and biology.

[39]  J. Huneau,et al.  Compartmental modeling of postprandial dietary nitrogen distribution in humans. , 2000, American journal of physiology. Endocrinology and metabolism.

[40]  C. Cobelli,et al.  Global identifiability of linear compartmental models-a computer algebra algorithm , 1998, IEEE Transactions on Biomedical Engineering.

[41]  Douglas B. Kell,et al.  Non-linear optimization of biochemical pathways: applications to metabolic engineering and parameter estimation , 1998, Bioinform..

[42]  Claudio Cobelli,et al.  Global identifiability of nonlinear models of biological systems , 2001, IEEE Transactions on Biomedical Engineering.

[43]  B. Beaufrère,et al.  Influence of protein intake on whole body and splanchnic leucine kinetics in humans. , 1997, The American journal of physiology.

[44]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[45]  P. Tessari,et al.  A model of skeletal muscle leucine kinetics measured across the human forearm. , 1995, The American journal of physiology.

[46]  L Finkelstein,et al.  Validation of simple and complex models in physiology and medicine. , 1984, The American journal of physiology.