Physiology-based Mathematical Models for the Intensive Care Unit: Application to Mechanical Ventilation

This work takes us a step closer to realizing personalized medicine, complementing empirical and heuristic way in which clinicians typically work. This thesis presents mechanistic models of physiology. These models, given continuous signals from a patient, can be fine-tuned via parameter estimation methods so that the model's outputs match the patient's. We thus obtain a virtual patient mimicking the patient at hand. Therapeutic scenarios can then be applied and optimal diagnosis and therapy can thus be attained. As such, personalized medicine can then be achieved without resorting to costly genetics. In particular we have developed a novel comprehensive mathematical model of the cardiopulmonary system that includes cardiovascular circulation, respiratory mechanics, tissue and alveolar gas exchange, as well as short-term neural control. Validity of the model was proven by the excellent agreement with real patient data, under normo-physiological as well as hypercapnic and hypoxic conditions, taken from literature. As a concrete example, a submodel of the lung mechanics was fine-tuned using real patient data and personalized respiratory parameters (resistance, R_rs, and compliance, C_rs) were estimated continually. This allows us to compute the patient's effort (Work of Breathing), continuously and more importantly noninvasively. Finally, the use of Bayesian estimation techniques, which allow incorporation of population studies and prior information about model's parameters, was proposed in the contest of patient-specific physiological models. A Bayesian Maximum a Posteriori Probability (MAP) estimator was implemented and applied to a case-study of respiratory mechanics. Its superiority against the classical Least Squares method was proven in data-poor conditions using both simulated and real animal data. This thesis can serve as a platform for a plethora of applications for cardiopulmonary personalized medicine

Dynamic Assessment of Baroreflex Control of Heart Rate During Induction of Propofol Anesthesia Using a Point Process Method

In this article, we present a point process method to assess dynamic baroreflex sensitivity (BRS) by estimating the baroreflex gain as focal component of a simplified closed-loop model of the cardiovascular system. Specifically, an inverse Gaussian probability distribution is used to model the heartbeat interval, whereas the instantaneous mean is identified by linear and bilinear bivariate regressions on both the previous R−R intervals (RR) and blood pressure (BP) beat-to-beat measures. The instantaneous baroreflex gain is estimated as the feedback branch of the loop with a point-process filter, while the RRBP feedforward transfer function representing heart contractility and vasculature effects is simultaneously estimated by a recursive least-squares filter. These two closed-loop gains provide a direct assessment of baroreflex control of heart rate (HR). In addition, the dynamic coherence, cross bispectrum, and their power ratio can also be estimated. All statistical indices provide a valuable quantitative assessment of the interaction between heartbeat dynamics and hemodynamics. To illustrate the application, we have applied the proposed point process model to experimental recordings from 11 healthy subjects in order to monitor cardiovascular regulation under propofol anesthesia. We present quantitative results during transient periods, as well as statistical analyses on steady-state epochs before and after propofol administration. Our findings validate the ability of the algorithm to provide a reliable and fast-tracking assessment of BRS, and show a clear overall reduction in baroreflex gain from the baseline period to the start of propofol anesthesia, confirming that instantaneous evaluation of arterial baroreflex control of HR may yield important implications in clinical practice, particularly during anesthesia and in postoperative care.National Institutes of Health (U.S.) (Grant R01-HL084502)National Institutes of Health (U.S.) (Grant K25-NS05758)National Institutes of Health (U.S.) (Grant DP2- OD006454)National Institutes of Health (U.S.) (Grant T32NS048005)National Institutes of Health (U.S.) (Grant T32NS048005)National Institutes of Health (U.S.) (Grant R01-DA015644)Massachusetts General Hospital (Clinical Research Center, UL1 Grant RR025758

Automated Detection of Incomplete Exhalation as an Indirect Detection of Auto-PEEP on Mechanically Ventilated Adults

Auto-PEEP is auto positive end-expiratory pressure due to excessive amounts of alveolar gas produced by sustained recurrent incomplete exhalation. Incomplete exhalation occurs when the exhaled breath never reaches a flow rate of 0 L/min. The objective of this dissertation is to develop an automated detection system of auto-PEEP through incomplete exhalation as revealed by ventilator graphics for mechanically ventilated adults. Auto-PEEP can cause adverse effects if allowed to linger and if not quickly identified. An automated detection system will be instrumental in helping to quickly identify auto-PEEP. A computerized algorithm was developed to detect incomplete exhalation based on the following three parameters:1) Foi, was used to represent the value of the flow at the onset of inhalation, 2) ∆T, was used to represent the value of time difference between onset inhalation to the 0 L/min mark, and 3) slope threshold, a value set for the slope of change of flow over ∆T. Optimum parameters of the algorithm were achieved for Foi = -3 L/min, ∆T = 0.2 s, and slope threshold = 90 L-s/min. A novel data set was introduced to validate the algorithm, yielding no significant difference in true positive rates (t = 1.5, df = 12.402, p-value = 0.1408) and false positive rates (t = 1.9, df = 16.765, p-value = 0.0725) as outcomes for two-tailed t-tests comparing the novel and old data set. To determine the relationship between auto-PEEP and detection of sustained incomplete exhalation, a correlation of a linear model of the novel data set between auto-PEEP and the percentage of incomplete exhalation detection out of the existing breaths (index) was investigated. A linear model should interpret the index value that corresponds to significant auto-PEEP presence; unfortunately, no significant linear model was found between incomplete exhalation index and auto-PEEP (F1,62 = 1.67, p-value = 0.2010). However, there was a relationship between the intrinsic PEEP values and the incomplete exhalation index as functions of time. The automated detection algorithm produced by this work provides a non-invasive method of automatically detecting auto-PEEP

Autonomic nervous system biomarkers from multi-modal and model-based signal processing in mental health and illness

Esta tesis se centra en técnicas de procesado multimodal y basado en modelos de señales para derivar parámetros fisiológicos, es decir, biomarcadores, relacionados con el sistema nervioso autónomo (ANS). El desarrollo de nuevos métodos para derivar biomarcadores de ANS no invasivos en la salud y la enfermedad mental ofrece la posibilidad de mejorar la evaluación del estrés y la monitorización de la depresión. Para este fin, el presente documento se estructura en tres partes principales. En la Parte I, se proporciona unaintroducción a la salud y la enfermedad mental (Cap. 1). Además, se presenta un marco teórico para investigar la etiología de los trastornos mentales y el papel del estrés en la enfermedad mental (Cap. 2). También se destaca la importancia de los biomarcadores no invasivos para la evaluación del ANS, prestando especial atención en la depresión clínica (Cap. 3, 4). En la Parte II, se proporciona el marco metodológico para derivar biomarcadores del ANS. Las técnicas de procesado de señales incluyen el análisis conjunto de la variabilidad del rítmo cardíaco (HRV) y la señal respiratoria (Cap. 6), técnicas novedosas para derivar la señal respiratoria del electrocardiograma (ECG) (Cap. 7) y un análisis robusto que se basa en modelar la forma de ondas del pulso del fotopletismograma (PPG) (Ch. 8). En la Parte III, los biomarcadores del ANS se evalúan en la quantificacióndel estrés (Cap. 9) y en la monitorización de la depresión (Ch. 10).Parte I: La salud mental no solo está relacionada con ese estado positivo de bienestar, en el que un individuo puede enfrentar a las situaciones estresantes de la vida, sino también con la ausencia de enfermedad mental. La enfermedad o trastorno mental se puede definir como un trastorno emocional, cognitivo o conductual que causa un deterioro funcional sustancial en una o más actividades importantes de la vida. Los trastornos mentales más comunes, que muchas veces coexisten, son la ansiedad y el trastorno depresivo mayor (MDD). La enfermedad mental tiene un impacto negativo en la calidad de vida, ya que se asocia con pérdidas considerables en la salud y el funcionamiento, y aumenta ignificativamente el riesgo de una persona de padecer enfermedades ardiovasculares.Un instigador común que subyace a la comorbilidad entre el MDD, la patologíacardiovascular y la ansiedad es el estrés mental. El estrés es común en nuestra vida de rítmo rapido e influye en nuestra salud mental. A corto plazo, ANS controla la respuesta cardiovascular a estímulos estresantes. La regulación de parámetros fisiológicos, como el rítmo cardíaco, la frecuencia respiratoria y la presión arterial, permite que el organismo responda a cambios repentinos en el entorno. Sin embargo, la adaptación fisiológica a un fenómeno ambiental que ocurre regularmente altera los sistemas biológicos involucrados en la respuesta al estrés. Las alteraciones neurobiológicas en el cerebro pueden alterar lafunción del ANS. La disfunción del ANS y los cambios cerebrales estructurales tienen un impacto negativo en los procesos cognitivos, emocionales y conductuales, lo que conduce al desarrollo de una enfermedad mental.Parte II: El desarrollo de métodos novedosos para derivar biomarcadores del ANS no invasivos ofrece la posibilidad de mejorar la evaluacón del estrés en individuos sanos y la disfunción del ANS en pacientes con MDD. El análisis conjunto de varias bioseñales (enfoquemultimodal) permite la cuantificación de interacciones entre sistemas biológicos asociados con ANS, mientras que el modelado de bioseãles y el análisis posterior de los parámetros del modelo (enfoque basado en modelos) permite la cuantificación robusta de cambios en mecanismos fisiológicos relacionados con el ANS. Un método novedoso, quetiene en cuenta los fenómenos de acoplo de fase y frecuencia entre la respiración y las señales de HRV para evaluar el acoplo cardiorrespiratorio no lineal cuadrático se propone en el Cap. 6.3. En el Cap. 7 se proponen nuevas técnicas paramejorar lamonitorización de la respiración. En el Cap. 8, para aumentar la robustez de algunas medidas morfológicas que reflejan cambios en el tonno arterial, se considera el modelado del pulso PPG como una onda principal superpuesta con varias ondas reflejadas.Parte III: Los biomarcadores del ANS se evalúan en la cuantificación de diferentes tipos de estrés, ya sea fisiológico o psicológico, en individuos sanos, y luego, en la monitorización de la depresión. En presencia de estrés mental (Cap. 9.1), inducido por tareas cognitivas, los sujetos sanos muestran un incremento en la frecuencia respiratoria y un mayor número de interacciones no lineales entre la respiración y la seãl de HRV. Esto podría estar asociado con una activación simpática, pero también con una respiración menos regular. En presencia de estrés hemodinámico (Cap. 9.2), inducido por un cambio postural, los sujetos sanos muestran una reducción en el acoplo cardiorrespiratoriono lineal cuadrático, que podría estar relacionado con una retracción vagal. En presencia de estrés térmico (Cap. 9.3), inducido por la exposición a emperaturas ambientales elevadas, los sujetos sanos muestran un aumento del equilibrio simpatovagal. Esto demuestra que los biomarcadores ANS son capaces de evaluar diferentes tipos de estrés y pueden explorarse más en el contexto de la monitorización de la depresión. En el Cap. 10, se evalúan las diferencias en la función del ANS entre elMDD y los sujetos sanos durante un protocolo de estrés mental, no solo con los valores brutos de los biomarcadores del ANS, sino también con los índices de reactividad autónoma, que reflejan la capacidad deun individuo para afrontar con una situación desafiante. Los resultados muestran que la depresión se asocia con un desequilibrio autonómico, que se caracteriza por una mayor actividad simpática y una reducción de la distensibilidad arterial. Los índices de reactividad autónoma cuantificados por cambios, entre etapas de estrés y de recuperación, en los sustitutos de la rigidez arterial, como la pérdida de amplitud de PPG en las ondas reflejadas, muestran el mejor rendimiento en términos de correlación con el grado de la depresión, con un coeficiente de correlación r = −0.5. La correlación negativa implicaque un mayor grado de depresión se asocia con una disminución de la reactividadautónoma. El poder discriminativo de los biomarcadores del ANS se aprecia también por su alto rendimiento diagnóstico para clasificar a los sujetos como MDD o sanos, con una precisión de 80.0%. Por lo tanto, se puede concluir que los biomarcadores del ANS pueden usarse para evaluar el estrés y que la distensibilidad arterial deteriorada podría constituir un biomarcador de salud mental útil en el seguimiento de la depresión.This dissertation is focused on multi-modal and model-based signal processing techniques for deriving physiological parameters, i.e. biomarkers, related to the autonomic nervous system (ANS). The development of novel approaches for deriving noninvasive ANS biomarkers in mental health and illness offers the possibility to improve the assessment of stress and the monitoring of depression. For this purpose, the present document is structured in three main parts. In Part I, an introduction to mental health and illness is provided (Ch. 1). Moreover, a theoretical framework for investigating the etiology of mental disorders and the role of stress in mental illness is presented (Ch. 2). The importance of noninvasive biomarkers for ANS assessment, paying particular attention in clinical depression, is also highlighted (Ch. 3, 4). In Part II, themethodological framework for deriving ANS biomarkers is provided. Signal processing techniques include the joint analysis of heart rate variability (HRV) and respiratory signals (Ch. 6), novel techniques for deriving the respiratory signal from electrocardiogram (ECG) (Ch. 7), and a robust photoplethysmogram(PPG)waveform analysis based on amodel-based approach (Ch. 8). In Part III, ANS biomarkers are evaluated in stress assessment (Ch. 9) and in the monitoring of depression (Ch. 10). Part I:Mental health is not only related to that positive state ofwell-being, inwhich an individual can cope with the normal stresses of life, but also to the absence of mental illness. Mental illness or disorder can be defined as an emotional, cognitive, or behavioural disturbance that causes substantial functional impairment in one or more major life activities. The most common mental disorders, which are often co-occurring, are anxiety and major depressive disorder (MDD). Mental illness has a negative impact on the quality of life, since it is associated with considerable losses in health and functioning, and increases significantly a person’s risk for cardiovascular diseases. A common instigator underlying the co-morbidity between MDD, cardiovascular pathology, and anxiety is mental stress. Stress is common in our fast-paced society and strongly influences our mental health. In the short term, ANS controls the cardiovascular response to stressful stimuli. Regulation of physiological parameters, such as heart rate, respiratory rate, and blood pressure, allows the organism to respond to sudden changes in the environment. However, physiological adaptation to a regularly occurring environmental phenomenon alters biological systems involved in stress response. Neurobiological alterations in the brain can disrupt the function of the ANS. ANS dysfunction and structural brain changes have a negative impact on cognitive, emotional, and behavioral processes, thereby leading to development of mental illness. Part II: The development of novel approaches for deriving noninvasive ANS biomarkers offers the possibility to improve the assessment of stress in healthy individuals and ANS dysfunction in MDD patients. Joint analysis of various biosignals (multi-modal approach) allows for the quantification of interactions among biological systems associated with ANS, while the modeling of biosignals and subsequent analysis of the model’s parameters (model-based approach) allows for the robust quantification of changes in physiological mechanisms related to the ANS. A novel method, which takes into account both phase and frequency locking phenomena between respiration and HRV signals, for assessing quadratic nonlinear cardiorespiratory coupling is proposed in Ch. 6.3. Novel techniques for improving the monitoring of respiration are proposed in Ch. 7. In Ch. 8, to increase the robustness for some morphological measurements reflecting arterial tone changes, the modeling of the PPG pulse as amain wave superposed with several reflected waves is considered. Part III: ANS biomarkers are evaluated in the assessment of different types of stress, either physiological or psychological, in healthy individuals, and then, in the monitoring of depression. In the presence of mental stress (Ch. 9.1), induced by cognitive tasks, healthy subjects show an increment in the respiratory rate and higher number of nonlinear interactions between respiration and HRV signal, which might be associated with a sympathetic activation, but also with a less regular breathing. In the presence of hemodynamic stress (Ch. 9.2), induced by a postural change, healthy subjects show a reduction in strength of the quadratic nonlinear cardiorespiratory coupling, whichmight be related to a vagal withdrawal. In the presence of heat stress (Ch. 9.3), induced by exposure to elevated environmental temperatures, healthy subjects show an increased sympathovagal balance. This demonstrates that ANS biomarkers are able to assess different types of stress and they can be further explored in the context of depression monitoring. In Ch. 10, differences in ANS function between MDD and healthy subjects during a mental stress protocol are assessed, not only with the raw values of ANS biomarkers but also with autonomic reactivity indices, which reflect the ability of an individual to copewith a challenging situation. Results show that depression is associated with autonomic imbalance, characterized by increased sympathetic activity and reduced arterial compliance. Autonomic reactivity indices quantified by changes, from stress to recovery, in arterial stiffness surrogates, such as the PPG amplitude loss in wave reflections, show the best performance in terms of correlation with depression severity, yielding to correlation coefficient r = −0.5. The negative correlation implies that a higher degree of depression is associated with a decreased autonomic reactivity. The discriminative power of ANS biomarkers is supported by their high diagnostic performance for classifying subjects as having MDD or not, yielding to accuracy of 80.0%. Therefore, it can be concluded that ANS biomarkers can be used for assessing stress and that impaired arterial compliance might constitute a biomarker of mental health useful in the monitoring of depression.<br /

An investigation into the effects of commencing haemodialysis in the critically ill

&lt;b&gt;Introduction:&lt;/b&gt; We have aimed to describe haemodynamic changes when haemodialysis is instituted in the critically ill. 3 hypotheses are tested: 1)The initial session is associated with cardiovascular instability, 2)The initial session is associated with more cardiovascular instability compared to subsequent sessions, and 3)Looking at unstable sessions alone, there will be a greater proportion of potentially harmful changes in the initial sessions compared to subsequent ones. &lt;b&gt;Methods:&lt;/b&gt; Data was collected for 209 patients, identifying 1605 dialysis sessions. Analysis was performed on hourly records, classifying sessions as stable/unstable by a cutoff of &gt;+/-20% change in baseline physiology (HR/MAP). Data from 3 hours prior, and 4 hours after dialysis was included, and average and minimum values derived. 3 time comparisons were made (pre-HD:during, during HD:post, pre-HD:post). Initial sessions were analysed separately from subsequent sessions to derive 2 groups. If a session was identified as being unstable, then the nature of instability was examined by recording whether changes crossed defined physiological ranges. The changes seen in unstable sessions could be described as to their effects: being harmful/potentially harmful, or beneficial/potentially beneficial. &lt;b&gt;Results:&lt;/b&gt; Discarding incomplete data, 181 initial and 1382 subsequent sessions were analysed. A session was deemed to be stable if there was no significant change (&gt;+/-20%) in the time-averaged or minimum MAP/HR across time comparisons. By this definition 85/181 initial sessions were unstable (47%, 95% CI SEM 39.8-54.2). Therefore Hypothesis 1 is accepted. This compares to 44% of subsequent sessions (95% CI 41.1-46.3). Comparing these proportions and their respective CI gives a 95% CI for the standard error of the difference of -4% to 10%. Therefore Hypothesis 2 is rejected. In initial sessions there were 92/1020 harmful changes. This gives a proportion of 9.0% (95% CI SEM 7.4-10.9). In the subsequent sessions there were 712/7248 harmful changes. This gives a proportion of 9.8% (95% CI SEM 9.1-10.5). Comparing the two unpaired proportions gives a difference of -0.08% with a 95% CI of the SE of the difference of -2.5 to +1.2. Hypothesis 3 is rejected. Fisher’s exact test gives a result of p=0.68, reinforcing the lack of significant variance. &lt;b&gt;Conclusions:&lt;/b&gt; Our results reject the claims that using haemodialysis is an inherently unstable choice of therapy. Although proportionally more of the initial sessions are classed as unstable, the majority of MAP and HR changes are beneficial in nature

Risk Scores for Predicting Mortality in Flail Chest

The objective of this thesis was to develop two risk scores which could predict the individual risk of in-hospital mortality for patients with flail chest using data from the Ontario Trauma Registry. The first study describes the univariate analyses conducted to identify mortality predictors. The second study details the logistic regression analysis that generated a risk score. Finally, the third study describes the decision tree analysis that produced the second risk score. The two risk scores were then compared. In summary, these three studies show that a minority of flail chest patients are currently undergoing operative repair and that a risk score may be a useful adjunct for surgeons to determine the individual risk of in-hospital mortality in patients requiring operative repair for flail chest

Understanding Peripheral Blood Pressure Signals: A Statistical Learning Approach

Proper estimation of body fluid status for human or animal subjects has always been a challenging problem. Accurate and timely estimate of body fluid can prevent life threatening conditions under trauma and severe dehydration. The main objective of this research is the estimation, classification and detection of dehydration in human and animal subjects using peripheral blood pressure (PBP) signals. Peripheral venous pressure (PVP) and peripheral arterial pressure (PAP) signals have been investigated in this research. Both PVP and PAP signals are PBP signals. A dataset of PVP signals was collected using standard peripheral intravenous catheters from human subjects suffering from hypertrophic pyloric stenosis. Using this dataset, we successfully classified dehydrated subjects from hydrated subjects using regularized logistic regression on frequency domain data of the PVP signals. During the data acquisition process, the PVP signals was corrupted by noise and blood clot. So, we developed an unsupervised anomaly detection algorithm for PVP signals using hidden Markov model and Kalman filter. This anomaly detection algorithm removed the human bias in data-preprocessing. Another dataset of PAP and PVP signals was collected from pigs under anesthesia using the Millar catheter. We proposed a integral pulse frequency modulation (IPFM) based signal model for both PAP and PVP signals. The proposed model-synthesized signal is highly correlated with the experimental data. The model-synthesized signals also performs similar to experimental signals under classification tasks. We also examine the model estimated parameters both qualitatively and quantitatively. This model can also quantify the effect of respiratory rate on heart rate variability. Increasing doses of anesthesia has similar effect of getting hydrated from dehydration

Cardiopulmonary responses to maximal aerobic exercise in patients with cystic fibrosis (article)

This is the final version. Available from Public Library of Science via the DOI in this record.The dataset associated with this article is located in ORE at: https://doi.org/10.24378/exe.1105Cystic fibrosis (CF) is a debilitating chronic condition, which requires complex and expensive disease management. Exercise has now been recognised as a critical factor in improving health and quality of life in patients with CF. Hence, cardiopulmonary exercise testing (CPET) is used to determine aerobic fitness of young patients as part of the clinical management of CF. However, at present there is a lack of conclusive evidence for one limiting system of aerobic fitness for CF patients at individual patient level. Here, we perform detailed data analysis that allows us to identify important systems-level factors that affect aerobic fitness. We use patients’ data and principal component analysis to confirm the dependence of CPET performance on variables associated with ventilation and metabolic rates of oxygen consumption. We find that the time at which participants cross the gas exchange threshold (GET) is well correlated with their overall performance. Furthermore, we propose a predictive modelling framework that captures the relationship between ventilatory dynamics, lung capacity and function and performance in CPET within a group of children and adolescents with CF. Specifically, we show that using Gaussian processes (GP) we can predict GET at the individual patient level with reasonable accuracy given the small sample size of the available group of patients. We conclude by presenting an example and future perspectives for improving and extending the proposed framework. The modelling and analysis have the potential to pave the way to designing personalised exercise programmes that are tailored to specific individual needs relative to patient’s treatment therapies.Wellcome TrustEngineering and Physical Sciences Research Counci

Understanding Peripheral Blood Pressure Signals: A Statistical Learning Approach

Proper estimation of body fluid status for human or animal subjects has always been a challenging problem. Accurate and timely estimate of body fluid can prevent life threatening conditions under trauma and severe dehydration. The main objective of this research is the estimation, classification and detection of dehydration in human and animal subjects using peripheral blood pressure (PBP) signals. Peripheral venous pressure (PVP) and peripheral arterial pressure (PAP) signals have been investigated in this research. Both PVP and PAP signals are PBP signals. A dataset of PVP signals was collected using standard peripheral intravenous catheters from human subjects suffering from hypertrophic pyloric stenosis. Using this dataset, we successfully classified dehydrated subjects from hydrated subjects using regularized logistic regression on frequency domain data of the PVP signals. During the data acquisition process, the PVP signals was corrupted by noise and blood clot. So, we developed an unsupervised anomaly detection algorithm for PVP signals using hidden Markov model and Kalman filter. This anomaly detection algorithm removed the human bias in data-preprocessing. Another dataset of PAP and PVP signals was collected from pigs under anesthesia using the Millar catheter. We proposed a integral pulse frequency modulation (IPFM) based signal model for both PAP and PVP signals. The proposed model-synthesized signal is highly correlated with the experimental data. The model-synthesized signals also performs similar to experimental signals under classification tasks. We also examine the model estimated parameters both qualitatively and quantitatively. This model can also quantify the effect of respiratory rate on heart rate variability. Increasing doses of anesthesia has similar effect of getting hydrated from dehydration