Search

Background and objective: Structural identifiability is a concept that considers whether the structure of a model together with a set of input-output relations uniquely determines the model parameters. In the mathematical modelling of biological systems, structural identifiability is an important concept since biological interpretations are typically made from the parameter estimates. For a system defined by ordinary differential equations, several methods have been developed to analyse whether the model is structurally identifiable or otherwise. Another well-used modelling framework, which is particularly useful when the experimental data are sparsely sampled and the population variance is of interest, is mixed-effects modelling. However, established identifiability analysis techniques for ordinary differential equations are not directly applicable to such models. Methods: In this paper, we present and apply three different methods that can be used to study structural identifiability in mixed-effects models. The first method, called the repeated measurement approach, is based on applying a set of previously established statistical theorems. The second method, called the augmented system approach, is based on augmenting the mixed-effects model to an extended state-space form. The third method, called the Laplace transform mixed-effects extension, is based on considering the moment invariants of the systems transfer function as functions of random variables. Results: To illustrate, compare and contrast the application of the three methods, they are applied to a set of mixed-effects models. Conclusions: Three structural identifiability analysis methods applicable to mixed-effects models have been presented in this paper.

The concept of structural identifiability for state-space models is expanded to cover mixed-effects state-space models. Two methods applicable for the analytical study of the structural identifiability of mixed-effects models are presented. The two methods are based on previously established techniques for non-mixed-effects models; namely the Taylor series expansion and the input-output form approach. By generating an exhaustive summary, and by assuming an infinite number of subjects, functions of random variables can be derived which in turn determine the distribution of the system's observation function(s). By considering the uniqueness of the analytical statistical moments of the derived functions of the random variables, the structural identifiability of the corresponding mixed-effects model can be determined. The two methods are applied to a set of examples of mixed-effects models to illustrate how they work in practice.

This study presents a dose-response-time (DRT) analysis based on a large preclinical biomarker dataset on the interaction between nicotinic acid (NiAc) and free fatty acids (FFA). Data were collected from studies that examined different rates, routes, and modes of NiAc provocations on the FFA time course. All information regarding the exposure to NiAc was excluded in order to demonstrate the utility of a DRT model. Special emphasis was placed on the selection process of the biophase model. An inhibitory I-max-model, driven by the biophase amount, acted on the turnover rate of FFA. A second generation NiAc/FFA model, which encompasses integral (slow buildup of tolerance - an extension of the previously used NiAc/FFA turnover models) and moderator (rapid and oscillatory) feedback control, was simultaneously fitted to all time courses in normal rats. The integral feedback control managed to capture an observed 90% adaptation (i.e., almost a full return to baseline) when 10 days constant-rate infusion protocols of NiAc were used. The half-life of the adaptation process had a 90% prediction interval between 3.5-12 in the present population. The pharmacodynamic parameter estimates were highly consistent when compared to an exposure-driven analysis, partly validating the DRT modelling approach and suggesting the potential of DRT analysis in areas where exposure data are not attainable. Finally, new numerical algorithms, which rely on sensitivity equations to robustly and efficiently compute the gradients in the parameter optimization, were successfully used for the mixed-effects approach in the parameter estimation.