Structural identifiability analysis and reparameterisation (parameter reduction) of a cardiovascular feedback model

Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination

BACKGROUND

Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error (χ 2 ) of LV myocardial strain, strain rate, and cavity volume.

RESULTS

A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients (χ 2  < 1.6), but minimum parameter reproducibility was poor (ICC min  = 0.01). Iterative reduction yielded a reproducible (ICC min  = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs (χ 2  < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05).

CONCLUSIONS

By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.

Classical structural identifiability methodology applied to low-dimensional dynamic systems in receptor theory

Mathematical modelling has become a key tool in pharmacological analysis, towards understanding dynamics of cell signalling and quantifying ligand-receptor interactions. Ordinary differential equation (ODE) models in receptor theory may be used to parameterise such interactions using timecourse data, but attention needs to be paid to the theoretical identifiability of the parameters of interest. Identifiability analysis is an often overlooked step in many bio-modelling works. In this paper we introduce structural identifiability analysis (SIA) to the field of receptor theory by applying three classical SIA methods (transfer function, Taylor Series and similarity transformation) to ligand-receptor binding models of biological importance (single ligand and Motulsky-Mahan competition binding at monomers, and a recently presented model of a single ligand binding at receptor dimers). New results are obtained which indicate the identifiable parameters for a single timecourse for Motulsky-Mahan binding and dimerised receptor binding. Importantly, we further consider combinations of experiments which may be performed to overcome issues of non-identifiability, to ensure the practical applicability of the work. The three SIA methods are demonstrated through a tutorial-style approach, using detailed calculations, which show the methods to be tractable for the low-dimensional ODE models.

A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability

There is an inherent tension in Quantitative Systems Pharmacology (QSP) between the need to incorporate mathematical descriptions of complex physiology and drug targets with the necessity of developing robust, predictive and well-constrained models. In addition to this, there is no “gold standard” for model development and assessment in QSP. Moreover, there can be confusion over terminology such as model and parameter identifiability; complex and simple models; virtual populations; and other concepts, which leads to potential miscommunication and misapplication of methodologies within modeling communities, both the QSP community and related disciplines. This perspective article highlights the pros and cons of using simple (often identifiable) vs. complex (more physiologically detailed but often non-identifiable) models, as well as aspects of parameter identifiability, sensitivity and inference methodologies for model development and analysis. The paper distills the central themes of the issue of identifiability and optimal model size and discusses open challenges.

A tutorial on model informed approaches to cardiovascular safety with focus on cardiac repolarisation

Drugs can affect the cardiovascular (CV) system either as an intended treatment or as an unwanted side effect. In both cases, drug-induced cardiotoxicities such as arrhythmia and unfavourable hemodynamic effects can occur, and be described using mathematical models; such a model informed approach can provide valuable information during drug development and can aid decision-making. However, in order to develop informative models, it is vital to understand CV physiology. The aims of this tutorial are to present (1) key background biological and medical aspects of the CV system, (2) CV electrophysiology, (3) CV safety concepts, (4) practical aspects of development of CV models and (5) regulatory expectations with a focus on using model informed and quantitative approaches to support nonclinical and clinical drug development. In addition, we share several case studies to provide practical information on project strategy (planning, key questions, assumptions setting, and experimental design) and mathematical models development that support decision-making during drug discovery and development.

Physiologically Based Pharmacokinetic Model Qualification and Reporting Procedures for Regulatory Submissions: A Consortium Perspective

This work provides a perspective on the qualification and verification of physiologically based pharmacokinetic (PBPK) platforms/models intended for regulatory submission based on the collective experience of the Simcyp Consortium members. Examples of regulatory submission of PBPK analyses across various intended applications are presented and discussed. European Medicines Agency (EMA) and US Food and Drug Administration (FDA) recent draft guidelines regarding PBPK analyses and reporting are encouraging, and to advance the use and acceptability of PBPK analyses, more clarity and flexibility are warranted.

Deterministic identifiability of population pharmacokinetic and pharmacokinetic-pharmacodynamic models

Identifiability is an important component of pharmacokinetic-pharmacodynamic (PKPD) model development. Structural identifiability is concerned with the uniqueness of the model parameters for a set of perfect input-output data and deterministic identifiability with the precision of parameter estimation given imperfect input-output data. We introduce two subcategories of deterministic identifiability, external and internal, and consider factors that distinguish between these forms. We define external deterministic identifiability as a function of externally controllable variables, i.e., the design, and internal deterministic identifiability as a function of the model and its parameter values. The concepts are explored using three common PK and PKPD models, and verified for their precision for the selected set of parameter values under optimal design.

Good Practices in Model-Informed Drug Discovery and Development: Practice, Application, and Documentation

This document was developed to enable greater consistency in the practice, application, and documentation of Model-Informed Drug Discovery and Development (MID3) across the pharmaceutical industry. A collection of "good practice" recommendations are assembled here in order to minimize the heterogeneity in both the quality and content of MID3 implementation and documentation. The three major objectives of this white paper are to: i) inform company decision makers how the strategic integration of MID3 can benefit R&D efficiency; ii) provide MID3 analysts with sufficient material to enhance the planning, rigor, and consistency of the application of MID3; and iii) provide regulatory authorities with substrate to develop MID3 related and/or MID3 enabled guidelines.

Modeling and Simulation Approaches for Cardiovascular Function and Their Role in Safety Assessment

Systems pharmacology modeling and pharmacokinetic-pharmacodynamic (PK/PD) analysis of drug-induced effects on cardiovascular (CV) function plays a crucial role in understanding the safety risk of new drugs. The aim of this review is to outline the current modeling and simulation (M&S) approaches to describe and translate drug-induced CV effects, with an emphasis on how this impacts drug safety assessment. Current limitations are highlighted and recommendations are made for future effort in this vital area of drug research.