Identification of Structured LTI MIMO State-Space Models

The identification of structured state-space model has been intensively studied for a long time but still has not been adequately addressed. The main challenge is that the involved estimation problem is a non-convex (or bilinear) optimization problem. This paper is devoted to developing an identification method which aims to find the global optimal solution under mild computational burden. Key to the developed identification algorithm is to transform a bilinear estimation to a rank constrained optimization problem and further a difference of convex programming (DCP) problem. The initial condition for the DCP problem is obtained by solving its convex part of the optimization problem which happens to be a nuclear norm regularized optimization problem. Since the nuclear norm regularized optimization is the closest convex form of the low-rank constrained estimation problem, the obtained initial condition is always of high quality which provides the DCP problem a good starting point. The DCP problem is then solved by the sequential convex programming method. Finally, numerical examples are included to show the effectiveness of the developed identification algorithm.Comment: Accepted to IEEE Conference on Decision and Control (CDC) 201

A mathematical model for breath gas analysis of volatile organic compounds with special emphasis on acetone

Recommended standardized procedures for determining exhaled lower respiratory nitric oxide and nasal nitric oxide have been developed by task forces of the European Respiratory Society and the American Thoracic Society. These recommendations have paved the way for the measurement of nitric oxide to become a diagnostic tool for specific clinical applications. It would be desirable to develop similar guidelines for the sampling of other trace gases in exhaled breath, especially volatile organic compounds (VOCs) which reflect ongoing metabolism. The concentrations of water-soluble, blood-borne substances in exhaled breath are influenced by: (i) breathing patterns affecting gas exchange in the conducting airways; (ii) the concentrations in the tracheo-bronchial lining fluid; (iii) the alveolar and systemic concentrations of the compound. The classical Farhi equation takes only the alveolar concentrations into account. Real-time measurements of acetone in end-tidal breath under an ergometer challenge show characteristics which cannot be explained within the Farhi setting. Here we develop a compartment model that reliably captures these profiles and is capable of relating breath to the systemic concentrations of acetone. By comparison with experimental data it is inferred that the major part of variability in breath acetone concentrations (e.g., in response to moderate exercise or altered breathing patterns) can be attributed to airway gas exchange, with minimal changes of the underlying blood and tissue concentrations. Moreover, it is deduced that measured end-tidal breath concentrations of acetone determined during resting conditions and free breathing will be rather poor indicators for endogenous levels. Particularly, the current formulation includes the classical Farhi and the Scheid series inhomogeneity model as special limiting cases.Comment: 38 page

Physiological modeling of isoprene dynamics in exhaled breath

Human breath contains a myriad of endogenous volatile organic compounds (VOCs) which are reflective of ongoing metabolic or physiological processes. While research into the diagnostic potential and general medical relevance of these trace gases is conducted on a considerable scale, little focus has been given so far to a sound analysis of the quantitative relationships between breath levels and the underlying systemic concentrations. This paper is devoted to a thorough modeling study of the end-tidal breath dynamics associated with isoprene, which serves as a paradigmatic example for the class of low-soluble, blood-borne VOCs. Real-time measurements of exhaled breath under an ergometer challenge reveal characteristic changes of isoprene output in response to variations in ventilation and perfusion. Here, a valid compartmental description of these profiles is developed. By comparison with experimental data it is inferred that the major part of breath isoprene variability during exercise conditions can be attributed to an increased fractional perfusion of potential storage and production sites, leading to higher levels of mixed venous blood concentrations at the onset of physical activity. In this context, various lines of supportive evidence for an extrahepatic tissue source of isoprene are presented. Our model is a first step towards new guidelines for the breath gas analysis of isoprene and is expected to aid further investigations regarding the exhalation, storage, transport and biotransformation processes associated with this important compound.Comment: 14 page

Physiological modeling of isoprene dynamics in exhaled breath

Human breath contains a myriad of endogenous volatile organic compounds (VOCs) which are reflective of ongoing metabolic or physiological processes. While research into the diagnostic potential and general medical relevance of these trace gases is conducted on a considerable scale, little focus has been given so far to a sound analysis of the quantitative relationships between breath levels and the underlying systemic concentrations. This paper is devoted to a thorough modeling study of the end-tidal breath dynamics associated with isoprene, which serves as a paradigmatic example for the class of low-soluble, blood-borne VOCs. Real-time measurements of exhaled breath under an ergometer challenge reveal characteristic changes of isoprene output in response to variations in ventilation and perfusion. Here, a valid compartmental description of these profiles is developed. By comparison with experimental data it is inferred that the major part of breath isoprene variability during exercise conditions can be attributed to an increased fractional perfusion of potential storage and production sites, leading to higher levels of mixed venous blood concentrations at the onset of physical activity. In this context, various lines of supportive evidence for an extrahepatic tissue source of isoprene are presented. Our model is a first step towards new guidelines for the breath gas analysis of isoprene and is expected to aid further investigations regarding the exhalation, storage, transport and biotransformation processes associated with this important compound.Comment: 14 page

Comparative Assessment of Epidemiological Models for Analyzing and Forecasting Infectious Disease Outbreaks

Mathematical modeling offers a quantitative framework for analyzing mechanisms underlying infectious disease transmission and explaining patterns in epidemiological data. Models are also commonly applied in outbreak investigations for assessing intervention and control strategies and generating epidemic forecasts in real time. However, successful application of mathematical models depends on the ability to reliably estimate key transmission and severity parameters, which are critical for guiding public health interventions. Overall, the three studies presented provide a thorough guide for assessing and utilizing mathematical models for describing infectious disease outbreak trends. In the first study, we describe the process for analyzing identifiability of parameters of interest in mechanistic disease transmission models. In the second study, we expand this idea to simple phenomenological models and explore the idea of overdispersion in the data and how to determine an appropriate error structure within the analyses. In the third study, we use previously validated phenomenological models to generate short-term forecasts of the ongoing COVID-19 pandemic. During infectious disease epidemics, public health authorities rely on modeling results to inform intervention decisions and resource allocation. Therefore, we highlight the importance of interpreting modeling results with caution, particularly regarding theoretical aspects of mathematical models and parameter estimation methods. Further, results from modeling studies should be presented with quantified uncertainty and interpreted in terms of the assumptions and limitations of the model, methods, and data used. The methodology presented in this dissertation provides a thorough guide for conducting model-based inferences and presenting the uncertainty associated with parameter estimation results

Lost in modelling and simulation?

Over the past few decades, physiologically-based pharmacokinetic modelling (PBPK) has been anticipated to be a powerful tool to improve the productivity of drug discovery and development. However, recently, multiple systematic evaluation studies independently suggested that the predictive power of current oral absorption (OA) PBPK models needs significant improvement. There is some disagreement between the industry and regulators about the credibility of OA PBPK modelling. Recently, the editorial board of AMDET&DMPK has announced the policy for the articles related to PBPK modelling (Modelling and simulation ethics). In this feature article, the background of this policy is explained: (1) Requirements for scientific writing of PBPK modelling, (2) Scientific literacy for PBPK modelling, and (3) Middle-out approaches. PBPK models are a useful tool if used correctly. This article will hopefully help advance the science of OA PBPK models.

Inverse Problems in data-driven multi-scale Systems Medicine: application to cancer physiology

Systems Medicine is an interdisciplinary framework involving reciprocal feedback between clinical investigation and mathematical modeling/analysis. Its aim is to improve the understanding of complex diseases by integrating knowledge and data across multiple levels of biological organization. This Thesis focuses on three inverse problems, arising from three kinds of data and related to cancer physiology, at different scales: tissues, cells, molecules. The general assumption of this piece of research is that cancer is associated toa path ological glucose consumption and, in fact, its functional behavior can be assessed by nuclear medicine experiments using [18F]-fluorodeoxyglucose (FDG) as a radioactive tracer mimicking the glucose properties. At tissue-scale, this Thesis considers the Positron Emission Tomography (PET) imaging technique, and deals with two distinct issues within compartmental analysis. First, this Thesis presents a compartmental approach, referred to as reference tissue model, for the estimation of FDG kinetics inside cancer tissues when the arterial blood input of the system is unknown. Then, this Thesis proposes an efficient and reliable method for recovering the compartmental kinetic parameters for each PET image pixel in the context of parametric imaging, exploiting information on the tissue physiology. Standard models in compartmental analysis assume that phosphorylation and dephosphorylation of FDG occur in the same intracellular cytosolic volume. Advances in cell biochemistry have shown that the appropriate location of dephosphorylation is the endoplasmic reticulum (ER). Therefore, at cell-scale, this Thesis formalizes a biochemically-driven compartmental model accounting for the specific role played by the ER, and applies it to the analysis of in vitro experiments on FDG uptake by cancer cell cultures obtained with a LigandTracer (LT) device. Finally, at molecule-scale, this Thesis provides a preliminary mathematical investigation of a chemical reaction network (CRN), represented by a huge Molecular Interaction Map (MIM), describing the biochemical interactions occurring between signaling proteins in specific pathways within a cancer cell. The main issue addressed in this case is the network parameterization problem, i.e. how to determine the reaction rate coefficients from protein concentration data

Cardiac cell modelling: Observations from the heart of the cardiac physiome project

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field