Compartmental model identification based on an empirical Bayesian approach: The case of thiamine kinetics in rats

Compartmental models are a very popular tool for the analysis of experiments in living systems. There are three main aspects that have to be taken into account: the degree of detail of the model, its a priori identifiability and thea posteriori (numerical) identifiability. In some cases, where standard approaches are adopted, the models can be eithera priori ora posteriori unidentifiable. The paper proposes model identification within a Bayesian framework, to solvea posteriori unidentifiability problems. In particular, a stochastic simulation algorithm is proposed to perform a Bayesian identification of compartmental models, and an empirical Bayesian technique is proposed to propagate information among multiple experiments. The power of this methodology was demonstrated by evaluating the kinetics of thiamine under several experimental conditions. The complexity of the existing model (nine parameters) and limited experimental data (8/12 for each model) causeda posteriori identifiability problems when standard approaches were adopted. The application of the methodology identifies all 28 models (four tissues under seven different conditions).

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

[2]  J. Cooper,et al.  The role of thiamine in nervous tissue , 1979, Neurochemical Research.

[3]  C. Cobelli,et al.  Parameter and structural identifiability concepts and ambiguities: a critical review and analysis. , 1980, The American journal of physiology.

[4]  J. DiStefano,et al.  Identifiability of Model Parameter , 1985 .

[5]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[6]  Ewart R. Carson,et al.  The mathematical modeling of metabolic and endocrine systems : model formulation, identification, and validation , 1983 .

[7]  Carlo Reggiani,et al.  Thiamine content and turnover rates of some rat nervous regions, using labeled thiamine as a tracer , 1980, Brain Research.

[8]  A. Rukhin Bayes and Empirical Bayes Methods for Data Analysis , 1997 .

[9]  R. Bellazzi,et al.  In Vivo Study of the Kinetics of Thiamine and its Phosphoesters in the Deafferented Rat Cerebellum , 1997, Metabolic Brain Disease.

[10]  Giuseppe De Nicolao,et al.  Dynamic Probabilistic Networks for Modelling and Identifying Dynamic Systems: a MCMC Approach , 1997, Intell. Data Anal..

[11]  Adrian F. M. Smith,et al.  Bayesian Analysis of Linear and Non‐Linear Population Models by Using the Gibbs Sampler , 1994 .

[12]  G. Rindi,et al.  Nervous tissue thiamine metabolism in vivo. I. Transport of thiamine and thiamine monophosphate from plasma to different brain regions of the rat , 1984, Brain Research.

[13]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[14]  G. Rindi,et al.  Nervous tissue thiamine metabolism in vivo. III. Influence of ethanol intake on the dynamics of thiamine and its phosphoesters in different brain regions and sciatic nerve of the rat , 1987, Brain Research.

[15]  G. Rindi,et al.  Effects of phenytoin on the in vivo kinetics of thiamine and its phosphoesters in rat nervous tissues , 1993, Brain Research.

[16]  L. Bettendorff Thiamine in excitable tissues: Reflections on a non-cofactor role , 1994, Metabolic Brain Disease.

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

[18]  Gelfand,et al.  AD-A 254 769 BAYESIAN ANALYSIS OF LINEAR AND NONLINEAR POPULATION MODELS USING THE GIBBS SAMPLER , 2022 .

[19]  G. Rindi,et al.  Nervous tissue thiamine metabolism in vivo. II. Thiamine and its phosphoesters dynamics in different brain regions and sciatic nerve of the rat , 1984, Brain Research.

[20]  K. R. Godfrey,et al.  Chapter 1 – IDENTIFIABILITY OF MODEL PARAMETERS , 1987 .

[21]  G. Rindi,et al.  Thiamine monophosphate, a normal constituent of rat plasma. , 1968, The Journal of nutrition.