Model-based optimal PEEP in mechanically ventilated ARDS patients in the intensive care unit

Pulmonary response prediction through personalized basis functions in a virtual patient model

BACKGROUND AND OBJECTIVE

Recruitment maneuvers with subsequent positive-end-expiratory-pressure (PEEP) have proven effective in recruiting lung volume and preventing alveoli collapse. However, determining a safe, effective, and patient-specific PEEP is not standardized, and this more optimal PEEP level evolves with patient condition, requiring personalised monitoring and care approaches to maintain optimal ventilation settings.

METHODS

This research examines 3 physiologically relevant basis function sets (exponential, parabolic, cumulative) to enable better prediction of elastance evolution for a virtual patient or digital twin model of MV lung mechanics, including novel elements to model and predict distension elastance. Prediction accuracy and robustness are validated against recruitment maneuver data from 18 volume-controlled ventilation (VCV) patients at 7 different baseline PEEP levels (0 to 12 cmH2O) and 14 pressure-controlled ventilation (PCV) patients at 4 different baseline PEEP levels (6 to 12 cmH2O), yielding 623 and 294 prediction cases, respectively. Predictions were made up to 12 cmH2O of added PEEP ahead, covering 6 × 2 cmH2O PEEP steps.

RESULTS

The 3 basis function sets yield median absolute peak inspiratory pressure (PIP) prediction error of 1.63 cmH2O for VCV patients, and median peak inspiratory volume (PIV) prediction error of 0.028 L for PCV patients. The exponential basis function set yields a better trade-off of overall performance across VCV and PCV prediction than parabolic and cumulative basis function sets from other studies. Comparing predicted and clinically measured distension prediction in VCV demonstrated consistent, robust high accuracy with R2 = 0.90-0.95.

CONCLUSIONS

The results demonstrate recruitment mechanics are best captured by an exponential basis function across different mechanical ventilation modes, matching physiological expectations, and accurately capture, for the first time, distension mechanics to within 5-10 % accuracy. Enabling the risk of lung injury to be predicted before changing ventilator settings. The overall outcomes significantly extend and more fully validate this digital twin or virtual mechanical ventilation patient model.

How Electronic Medical Record Integration Can Support More Efficient Critical Care Clinical Trials

Large volumes of data are collected on critically ill patients, and using data science to extract information from the electronic medical record (EMR) and to inform the design of clinical trials represents a new opportunity in critical care research. Using improved methods of phenotyping critical illnesses, subject identification and enrollment, and targeted treatment group assignment alongside newer trial designs such as adaptive platform trials can increase efficiency while lowering costs. Some tools such as the EMR to automate data collection are already in use. Refinement of data science approaches in critical illness research will allow for better clinical trials and, ultimately, improved patient outcomes.

Physiological trend analysis of a novel cardio-pulmonary model during a preload reduction manoeuvre

BACKGROUND AND OBJECTIVE

Mechanical ventilation causes adverse effects on the cardiovascular system. However, the exact nature of the effects on haemodynamic parameters is not fully understood. A recently developed cardio-vascular system model which incorporates cardio-pulmonary interactions is compared to the original 3-chamber cardiovascular model to investigate the exact effects of mechanical ventilation on haemodynamic parameters and to assess the trade-off of model complexity and model reliability between the 2 models.

METHODS

Both the cardio-pulmonary and three chamber models are used to identify cardiovascular system parameters from aortic pressure, left ventricular volume, airway flow and airway pressure measurements from 4 pigs during a preload reduction manoeuvre. Outputs and parameter estimations from both models are contrasted to assess the relative performance of each model and to further investigate the effects of mechanical ventilation on haemodynamic parameters.

RESULTS

Both models tracked measurements accurately as expected. There was no identifiable increase in error from the added complexity of the cardio-pulmonary model, with both models having a mean average error below 0.5% for all pigs. Identified left ventricle and vena cava elastances of the 3-chamber model was found to diverge exponentially with PEEP from identified left ventricle and vena cava elastances of the cardio-pulmonary model. The r² of the fit for each pig ranged from 0.888 to 0.998 for left ventricle elastance divergence and from 0.905 to 0.999 for vena cava elastance divergence. All other identified parameters showed no significant difference between models.

CONCLUSIONS

Despite the increase in model complexity, there was no loss in the cardio-pulmonary model's ability to accurately estimate haemodynamic parameters and reproduce system dynamics. Furthermore, the cardio-pulmonary model was able to demonstrate how mechanical ventilation affected parameter estimations as PEEP was increased. The 3-chamber model was shown to produce parameter estimations which diverged exponentially with PEEP, while the cardiopulmonary model estimations remained more stable, suggesting its ability to produce more physiologically accurate parameter estimations under higher PEEP conditions.

Estimating the incidence of spontaneous breathing effort of mechanically ventilated patients using a non-linear auto regressive (NARX) model

BACKGROUND AND OBJECTIVE

Mechanical ventilation (MV) provides breathing support for acute respiratory distress syndrome (ARDS) patients in the intensive care unit, but is difficult to optimize. Too much, or too little of pressure or volume support can cause further ventilator-induced lung injury, increasing length of MV, cost and mortality. Patient-specific respiratory mechanics can help optimize MV settings. However, model-based estimation of respiratory mechanics is less accurate when patient exhibit un-modeled spontaneous breathing (SB) efforts on top of ventilator support. This study aims to estimate and quantify SB efforts by reconstructing the unaltered passive mechanics airway pressure using NARX model.

METHODS

Non-linear autoregressive (NARX) model is used to reconstruct missing airway pressure due to the presence of spontaneous breathing effort in mv patients. Then, the incidence of SB patients is estimated. The study uses a total of 10,000 breathing cycles collected from 10 ARDS patients from IIUM Hospital in Kuantan, Malaysia. In this study, there are 2 different ratios of training and validating methods. Firstly, the initial ratio used is 60:40 which indicates 600 breath cycles for training and remaining 400 breath cycles used for testing. Then, the ratio is varied using 70:30 ratio for training and testing data.

RESULTS AND DISCUSSION

The mean residual error between original airway pressure and reconstructed airway pressure is denoted as the magnitude of effort. The median and interquartile range of mean residual error for both ratio are 0.0557 [0.0230 - 0.0874] and 0.0534 [0.0219 - 0.0870] respectively for all patients. The results also show that Patient 2 has the highest percentage of SB incidence and Patient 10 with the lowest percentage of SB incidence which proved that NARX model is able to perform for both higher incidence of SB effort or when there is a lack of SB effort.

CONCLUSION

This model is able to produce the SB incidence rate based on 10% threshold. Hence, the proposed NARX model is potentially useful to estimate and identify patient-specific SB effort, which has the potential to further assist clinical decisions and optimize MV settings.

Quantifying neonatal patient effort using non-invasive model-based methods

Patient-specific spontaneous breathing effort (SB) is common in invasively mechanically ventilated (MV) adult patients, and especially common in preterm neonates who are not typically sedated. However, there is no proven, ethically feasible and non-invasive method to quantify SB effort in neonates, creating the potential for model-based measures. Lung mechanics and SB effort are segregated using a basis function model to identify passive lung mechanics, and an additional time-varying elastance model to identify patient-specific SB effort and asynchrony as negative and positive added elastances, respectively. Data from ten preterm neonates on standard MV care in the neonatal intensive care unit (NICU) are used to assess this model-based approach, using area under the curve (AUC) for positive (asynchrony) and negative (SB effort) time-varying elastance. Median [interquartile-range (IQR)] of passive pulmonary lung elastance was 3.82 [2.09-5.80] cmH₂O/ml. Median [IQR] AUC quantified SB effort was -0.32 [-0.43--0.12]cmH₂O/ml. AUC quantified asynchrony was 0.00 [0.00-0.01]cmH₂O/ml, and affected 28% of the 25,287 total breaths. This proof of concept model-based approach provides a non-invasive, computationally straightforward, and thus clinically feasible means to quantify patient-specific spontaneous breathing effort and asynchrony.

Over-distension prediction via hysteresis loop analysis and patient-specific basis functions in a virtual patient model

BACKGROUND AND OBJECTIVE

Recruitment maneuvers (RMs) with subsequent positive-end-expiratory-pressure (PEEP) have proven effective in recruiting lung volume and preventing alveolar collapse. However, a suboptimal PEEP could induce undesired injury in lungs by insufficient or excessive breath support. Thus, a predictive model for patient response under PEEP changes could improve clinical care and lower risks.

METHODS

This research adds novel elements to a virtual patient model to identify and predict patient-specific lung distension to optimise and personalise care. Model validity and accuracy are validated using data from 18 volume-controlled ventilation (VCV) patients at 7 different baseline PEEP levels (0-12cmH₂O), yielding 623 prediction cases. Predictions were made up to ΔPEEP = 12cmH₂O ahead covering 6x2cmH₂O PEEP steps.

RESULTS

Using the proposed lung distension model, 90% of absolute peak inspiratory pressure (PIP) prediction errors compared to clinical measurement are within 3.95cmH₂O, compared with 4.76cmH₂O without this distension term. Comparing model-predicted and clinically measured distension had high correlation increasing to R² = 0.93-0.95 if maximum ΔPEEP ≤ 6cmH₂O. Predicted dynamic functional residual capacity (V_frc) changes as PEEP rises yield 0.013L median prediction error for both prediction groups and overall R² of 0.84.

CONCLUSIONS

Overall results demonstrate nonlinear distension mechanics are accurately captured in virtual lung mechanics patients for mechanical ventilation, for the first time. This result can minimise the risk of lung injury by predicting its potential occurrence of distension before changing ventilator settings. The overall outcomes significantly extend and more fully validate this virtual mechanical ventilation patient model.

Use of Machine Learning to Screen for Acute Respiratory Distress Syndrome Using Raw Ventilator Waveform Data

To develop and characterize a machine learning algorithm to discriminate acute respiratory distress syndrome from other causes of respiratory failure using only ventilator waveform data.

Design

Retrospective, observational cohort study.

Setting

Academic medical center ICU.

Patients

Adults admitted to the ICU requiring invasive mechanical ventilation, including 50 patients with acute respiratory distress syndrome and 50 patients with primary indications for mechanical ventilation other than hypoxemic respiratory failure.

Interventions

None.

Measurements and Main Results

Pressure and flow time series data from mechanical ventilation during the first 24-hours after meeting acute respiratory distress syndrome criteria (or first 24-hr of mechanical ventilation for non-acute respiratory distress syndrome patients) were processed to extract nine physiologic features. A random forest machine learning algorithm was trained to discriminate between the patients with and without acute respiratory distress syndrome. Model performance was assessed using the area under the receiver operating characteristic curve, sensitivity, specificity, positive predictive value, and negative predictive value. Analyses examined performance when the model was trained using data from the first 24 hours and tested using withheld data from either the first 24 hours (24/24 model) or 6 hours (24/6 model). Area under the receiver operating characteristic curve, sensitivity, specificity, positive predictive value, and negative predictive value were 0.88, 0.90, 0.71, 0.77, and 0.90 (24/24); and 0.89, 0.90, 0.75, 0.83, and 0.83 (24/6).

Conclusions

Use of machine learning and physiologic information derived from raw ventilator waveform data may enable acute respiratory distress syndrome screening at early time points after intubation. This approach, combined with traditional diagnostic criteria, could improve timely acute respiratory distress syndrome recognition and enable automated clinical decision support, especially in settings with limited availability of conventional diagnostic tests and electronic health records.

The Potential Role of Lung-Protective Ventilation in Preventing Postoperative Delirium in Elderly Patients Undergoing Prone Spinal Surgery: A Preliminary Study

Background

Postoperative delirium (POD) is a frequent complication in elderly patients, usually occurring within a few days after surgery. This study investigated the effect of lung-protective ventilation (LPV) on POD in elderly patients undergoing spinal surgery and the mechanism by which LPV suppresses POD.

Material/Methods

Seventy-one patients aged ≥65 years were randomized to receive LPV or conventional mechanical ventilation (MV), consisting of intermittent positive pressure ventilation following induction of anesthesia. The tidal volume in patients who received MV was 8 ml/kg predicted body weight (PBW), and the ventilation frequency was 12 times/min. The tidal volume in patients who received LPV was 6 ml/kg PBW, the positive end-expiratory pressure was 5 cmH₂O, and the ventilation frequency was 15 times/min, with a lung recruitment maneuver performed every 30 min. Blood samples were collected immediately before anesthesia induction (T₀), 10 min (T₁) and 60 min (T₂) after turning over, immediately after the operation (T₃), and 15 min after extubation (T₄) for blood gas analysis. Simultaneous cerebral oxygen saturation (rSO₂) and cerebral desaturation were recorded. Preoperative and postoperative serum concentrations of interleukin (IL)-6, IL-10 and glial fibrillary acidic protein (GFAP) were measured by ELISA. POD was assessed by nursing delirium screening score.

Results

Compared with the MV group, pH was lower and PaCO₂ higher in the LPV group at T₂. In addition PaO₂, SaO₂, and PaO₂/FiO₂ were higher at T₁, and T₄, and rSO₂ was higher at T₃, and T₄ in the LPV than in the MV group (P<0.05 each). Postoperative serum GFAP and IL-6 were lower and IL-10 higher in the LPV group. The incidences of cerebral desaturation and POD were significantly lower in the LPV group (P<0.05).

Conclusions

LPV may reduce POD in elderly patients undergoing spinal surgery by inhibiting inflammation and improving cerebral oxygen metabolism.

Parameter estimation in a minimal model of cardio-pulmonary interactions

Mechanical ventilation is a widely used breathing support for patients in intensive care. Its effects on the respiratory and cardiovascular systems are complex and difficult to predict. This work first presents a minimal mathematical model representing the mechanics of both systems and their interaction, in terms of flows, pressures and volumes. The aim of this model is to get insight on the two systems' status when mechanical ventilation settings, such as positive end-expiratory pressure, are changing. The parameters of the model represent cardiac elastances and vessel compliances and resistances. As a second step, these parameters are estimated from 16 experimental datasets. The data come from three pig experiments reproducing intensive care conditions, where a large range of positive end-expiratory pressures was imposed by the mechanical ventilator. The data used for parameter estimation is limited to information available in the intensive care unit, such as stroke volume, central venous pressure and systemic arterial pressure. The model is able to satisfactorily reproduce this experimental data, with mean relative errors ranging from 1 to 26%. The model also reproduces the dynamics of the cardio-vascular and respiratory systems, and their interaction. By looking at the estimated parameter values, one can quantitatively track how the two coupled systems mechanically react to changes in external conditions imposed by the ventilator. This work thus allows real-time, model-based management of ventilator settings.

Predictive Virtual Patient Modelling of Mechanical Ventilation: Impact of Recruitment Function

Mechanical ventilation is a life-support therapy for intensive care patients suffering from respiratory failure. To reduce the current rate of ventilator-induced lung injury requires ventilator settings that are patient-, time-, and disease-specific. A common lung protective strategy is to optimise the level of positive end-expiratory pressure (PEEP) through a recruitment manoeuvre to prevent alveolar collapse at the end of expiration and to improve gas exchange through recruitment of additional alveoli. However, this process can subject parts of the lung to excessively high pressures or volumes. This research significantly extends and more robustly validates a previously developed pulmonary mechanics model to predict lung mechanics throughout recruitment manoeuvres. In particular, the process of recruitment is more thoroughly investigated and the impact of the inclusion of expiratory data when estimating peak inspiratory pressure is assessed. Data from the McREM trial and CURE pilot trial were used to test model predictive capability and assumptions. For PEEP changes of up to and including 14 cmH₂O, the parabolic model was shown to improve peak inspiratory pressure prediction resulting in less than 10% absolute error in the CURE cohort and 16% in the McREM cohort. The parabolic model also better captured expiratory mechanics than the exponential model for both cohorts.

Biomedical engineer's guide to the clinical aspects of intensive care mechanical ventilation

Background

Mechanical ventilation is an essential therapy to support critically ill respiratory failure patients. Current standards of care consist of generalised approaches, such as the use of positive end expiratory pressure to inspired oxygen fraction (PEEP–FiO₂) tables, which fail to account for the inter- and intra-patient variability between and within patients. The benefits of higher or lower tidal volume, PEEP, and other settings are highly debated and no consensus has been reached. Moreover, clinicians implicitly account for patient-specific factors such as disease condition and progression as they manually titrate ventilator settings. Hence, care is highly variable and potentially often non-optimal. These conditions create a situation that could benefit greatly from an engineered approach. The overall goal is a review of ventilation that is accessible to both clinicians and engineers, to bridge the divide between the two fields and enable collaboration to improve patient care and outcomes. This review does not take the form of a typical systematic review. Instead, it defines the standard terminology and introduces key clinical and biomedical measurements before introducing the key clinical studies and their influence in clinical practice which in turn flows into the needs and requirements around how biomedical engineering research can play a role in improving care. Given the significant clinical research to date and its impact on this complex area of care, this review thus provides a tutorial introduction around the review of the state of the art relevant to a biomedical engineering perspective.

Discussion

This review presents the significant clinical aspects and variables of ventilation management, the potential risks associated with suboptimal ventilation management, and a review of the major recent attempts to improve ventilation in the context of these variables. The unique aspect of this review is a focus on these key elements relevant to engineering new approaches. In particular, the need for ventilation strategies which consider, and directly account for, the significant differences in patient condition, disease etiology, and progression within patients is demonstrated with the subsequent requirement for optimal ventilation strategies to titrate for patient- and time-specific conditions.

Conclusion

Engineered, protective lung strategies that can directly account for and manage inter- and intra-patient variability thus offer great potential to improve both individual care, as well as cohort clinical outcomes.

Next-generation, personalised, model-based critical care medicine: a state-of-the art review of in silico virtual patient models, methods, and cohorts, and how to validation them

Critical care, like many healthcare areas, is under a dual assault from significantly increasing demographic and economic pressures. Intensive care unit (ICU) patients are highly variable in response to treatment, and increasingly aging populations mean ICUs are under increasing demand and their cohorts are increasingly ill. Equally, patient expectations are growing, while the economic ability to deliver care to all is declining. Better, more productive care is thus the big challenge. One means to that end is personalised care designed to manage the significant inter- and intra-patient variability that makes the ICU patient difficult. Thus, moving from current "one size fits all" protocolised care to adaptive, model-based "one method fits all" personalised care could deliver the required step change in the quality, and simultaneously the productivity and cost, of care. Computer models of human physiology are a unique tool to personalise care, as they can couple clinical data with mathematical methods to create subject-specific models and virtual patients to design new, personalised and more optimal protocols, as well as to guide care in real-time. They rely on identifying time varying patient-specific parameters in the model that capture inter- and intra-patient variability, the difference between patients and the evolution of patient condition. Properly validated, virtual patients represent the real patients, and can be used in silico to test different protocols or interventions, or in real-time to guide care. Hence, the underlying models and methods create the foundation for next generation care, as well as a tool for safely and rapidly developing personalised treatment protocols over large virtual cohorts using virtual trials. This review examines the models and methods used to create virtual patients. Specifically, it presents the models types and structures used and the data required. It then covers how to validate the resulting virtual patients and trials, and how these virtual trials can help design and optimise clinical trial. Links between these models and higher order, more complex physiome models are also discussed. In each section, it explores the progress reported up to date, especially on core ICU therapies in glycemic, circulatory and mechanical ventilation management, where high cost and frequency of occurrence provide a significant opportunity for model-based methods to have measurable clinical and economic impact. The outcomes are readily generalised to other areas of medical care.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

Extrapolation of a non-linear autoregressive model of pulmonary mechanics

For patients with acute respiratory distress syndrome (ARDS), mechanical ventilation (MV) is an essential therapy in the intensive care unit (ICU). Suboptimal PEEP levels in MV can cause ventilator induced lung injury, which is associated with increased mortality, extended ICU stay, and high cost. The ability to predict the outcome of respiratory mechanics in response to changes in PEEP would thus provide a critical advantage in personalising and improving care. Testing the potentially dangerous high pressures would not be required to assess their impact. A nonlinear autoregressive (NARX) model was used to predict airway pressure in 19 data sets from 10 mechanically ventilated ARDS patients. Patient-specific NARX models were identified from pressure and flow data over one, two, three, or four adjacent PEEP levels in a recruitment manoeuvre. Extrapolation of NARX model elastance functions allowed prediction of patient responses to PEEP changes to higher or lower pressures. NARX model predictions were more successful than those using a well validated first order model (FOM). The most clinically important results were for extrapolation up one PEEP step of 2cmH₂O from the highest PEEP in the training data. When the NARX model was trained on one PEEP level, the mean RMS residual for the extrapolation PEEP level was 0.52 (90% CI: 0.47-0.57) cmH₂O, compared to 1.50 (90% CI: 1.38-1.62) cmH₂O for the FOM. When trained on four PEEP levels, the NARX result was 0.50 (90% CI: 0.42-0.58) cmH₂O, and was 1.95 (90% CI: 1.71-2.19) cmH₂O for the FOM. The results suggest that a full recruitment manoeuvre may not be required for the NARX model to obtain a useful estimate of the pressure waveform at higher PEEP levels. The methodology could thus allow clinicians to make informed decisions about ventilator PEEP settings while reducing the risk associated with high PEEP, and subsequent high peak airway pressures.

A Review on Human Respiratory Modeling

Input impedance of the respiratory system is measured by forced oscillation technique (FOT). Multiple prior studies have attempted to match the electromechanical models of the respiratory system to impedance data. Since the mechanical behavior of airways and the respiratory system as a whole are similar to an electrical circuit in a combination of series and parallel formats some theories were introduced according to this issue. It should be noted that, the number of elements used in these models might be less than those required due to the complexity of the pulmonary-chest wall anatomy. Various respiratory models have been proposed based on this idea in order to demonstrate and assess the different parts of respiratory system related to children and adults data. With regard to our knowledge, some of famous respiratory models in related to obstructive, restrictive diseases and also Acute Respiratory Distress Syndrome (ARDS) are reviewed in this article.

An estimation of mechanical stress on alveolar walls during repetitive alveolar reopening and closure

Alveolar overdistension and mechanical stresses generated by repetitive opening and closing of small airways and alveoli have been widely recognized as two primary mechanistic factors that may contribute to the development of ventilator-induced lung injury. A long-duration exposure of alveolar epithelial cells to even small, shear stresses could lead to the changes in cytoskeleton and the production of inflammatory mediators. In this paper, we have made an attempt to estimate in situ the magnitudes of mechanical stresses exerted on the alveolar walls during repetitive alveolar reopening by using a tape-peeling model of McEwan and Taylor (35). To this end, we first speculate the possible ranges of capillary number ( Ca) ≡ μU/ γ (a dimensionless combination of surface tension γ, fluid viscosity μ, and alveolar opening velocity U) during in vivo alveolar opening. Subsequent calculations show that increasing respiratory rate or inflation rate serves to increase the values of mechanical stresses. For a normal lung, the predicted maximum shear stresses are <15 dyn/cm2 at all respiratory rates, whereas for a lung with elevated surface tension or viscosity, the maximum shear stress will notably increase, even at a slow respiratory rate. Similarly, the increased pressure gradients in the case of elevated surface or viscosity may lead to a pressure drop >300 dyn/cm2 across a cell, possibly inducing epithelial hydraulic cracks. In addition, we have conceived of a geometrical model of alveolar opening to make a prediction of the positive end-expiratory pressure (PEEP) required to splint open a collapsed alveolus, which as shown by our results, covers a wide range of pressures, from several centimeters to dozens of centimeters of water, strongly depending on the underlying pulmonary conditions. The establishment of adequate regional ventilation-to-perfusion ratios may prevent recruited alveoli from reabsorption atelectasis and accordingly, reduce the required levels of PEEP. The present study and several recent animal experiments likewise suggest that a lung-protective ventilation strategy should not only include small tidal volume and plateau pressure limitations but also consider such cofactors as ventilation frequency and inflation rate.

Evaluation of lung recruitment maneuvers in acute respiratory distress syndrome using computer simulation

Introduction

Direct comparison of the relative efficacy of different recruitment maneuvers (RMs) for patients with acute respiratory distress syndrome (ARDS) via clinical trials is difficult, due to the heterogeneity of patient populations and disease states, as well as a variety of practical issues. There is also significant uncertainty regarding the minimum values of positive end-expiratory pressure (PEEP) required to ensure maintenance of effective lung recruitment using RMs. We used patient-specific computational simulation to analyze how three different RMs act to improve physiological responses, and investigate how different levels of PEEP contribute to maintaining effective lung recruitment.

Methods

We conducted experiments on five ‘virtual’ ARDS patients using a computational simulator that reproduces static and dynamic features of a multivariable clinical dataset on the responses of individual ARDS patients to a range of ventilator inputs. Three recruitment maneuvers (sustained inflation (SI), maximal recruitment strategy (MRS) followed by a titrated PEEP, and prolonged recruitment maneuver (PRM)) were implemented and evaluated for a range of different pressure settings.

Results

All maneuvers demonstrated improvements in gas exchange, but the extent and duration of improvement varied significantly, as did the observed mechanism of operation. Maintaining adequate post-RM levels of PEEP was seen to be crucial in avoiding cliff-edge type re-collapse of alveolar units for all maneuvers. For all five patients, the MRS exhibited the most prolonged improvement in oxygenation, and we found that a PEEP setting of 35 cm H₂O with a fixed driving pressure of 15 cm H₂O (above PEEP) was sufficient to achieve 95% recruitment. Subsequently, we found that PEEP titrated to a value of 16 cm H₂O was able to maintain 95% recruitment in all five patients.

Conclusions

There appears to be significant scope for reducing the peak levels of PEEP originally specified in the MRS and hence to avoid exposing the lung to unnecessarily high pressures. More generally, our study highlights the huge potential of computer simulation to assist in evaluating the efficacy of different recruitment maneuvers, in understanding their modes of operation, in optimizing RMs for individual patients, and in supporting clinicians in the rational design of improved treatment strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13054-014-0723-6) contains supplementary material, which is available to authorized users.

A patient-specific airway branching model for mechanically ventilated patients

BACKGROUND

Respiratory mechanics models have the potential to guide mechanical ventilation. Airway branching models (ABMs) were developed from classical fluid mechanics models but do not provide accurate models of in vivo behaviour. Hence, the ABM was improved to include patient-specific parameters and better model observed behaviour (ABMps).

METHODS

The airway pressure drop of the ABMps was compared with the well-accepted dynostatic algorithm (DSA) in patients diagnosed with acute respiratory distress syndrome (ARDS). A scaling factor (α) was used to equate the area under the pressure curve (AUC) from the ABMps to the AUC of the DSA and was linked to patient state.

RESULTS

The ABMps recorded a median α value of 0.58 (IQR: 0.54-0.63; range: 0.45-0.66) for these ARDS patients. Significantly lower α values were found for individuals with chronic obstructive pulmonary disease (P < 0.001).

CONCLUSION

The ABMps model allows the estimation of airway pressure drop at each bronchial generation with patient-specific physiological measurements and can be generated from data measured at the bedside. The distribution of patient-specific α values indicates that the overall ABM can be readily improved to better match observed data and capture patient condition.

Model-based setting of inspiratory pressure and respiratory rate in pressure-controlled ventilation

Mechanical ventilation carries the risk of ventilator-induced-lung-injury (VILI). To minimize the risk of VILI, ventilator settings should be adapted to the individual patient properties. Mathematical models of respiratory mechanics are able to capture the individual physiological condition and can be used to derive personalized ventilator settings. This paper presents model-based calculations of inspiration pressure (pI), inspiration and expiration time (tI, tE) in pressure-controlled ventilation (PCV) and a retrospective evaluation of its results in a group of mechanically ventilated patients. Incorporating the identified first order model of respiratory mechanics in the basic equation of alveolar ventilation yielded a nonlinear relation between ventilation parameters during PCV. Given this patient-specific relation, optimized settings in terms of minimal pI and adequate tE can be obtained. We then retrospectively analyzed data from 16 ICU patients with mixed pathologies, whose ventilation had been previously optimized by ICU physicians with the goal of minimization of inspiration pressure, and compared the algorithm's 'optimized' settings to the settings that had been chosen by the physicians. The presented algorithm visualizes the patient-specific relations between inspiration pressure and inspiration time. The algorithm's calculated results highly correlate to the physician's ventilation settings with r = 0.975 for the inspiration pressure, and r = 0.902 for the inspiration time. The nonlinear patient-specific relations of ventilation parameters become transparent and support the determination of individualized ventilator settings according to therapeutic goals. Thus, the algorithm is feasible for a variety of ventilated ICU patients and has the potential of improving lung-protective ventilation by minimizing inspiratory pressures and by helping to avoid the build-up of clinically significant intrinsic positive end-expiratory pressure.

Analysis of different model-based approaches for estimating dFRC for real-time application

BACKGROUND

Acute Respiratory Distress Syndrome (ARDS) is characterized by inflammation, filling of the lung with fluid and the collapse of lung units. Mechanical ventilation (MV) is used to treat ARDS using positive end expiratory pressure (PEEP) to recruit and retain lung units, thus increasing pulmonary volume and dynamic functional residual capacity (dFRC) at the end of expiration. However, simple, non-invasive methods to estimate dFRC do not exist.

METHODS

Four model-based methods for estimating dFRC are compared based on their performance on two separate clinical data cohorts. The methods are derived from either stress-strain theory or a single compartment lung model, and use commonly controlled or measured parameters (lung compliance, plateau airway pressure, pressure-volume (PV) data). Population constants are determined for the stress-strain approach, which is implemented using data at both single and multiple PEEP levels. Estimated values are compared to clinically measured values to assess the reliability of each method for each cohort individually and combined.

RESULTS

The stress-strain multiple breath (at multiple PEEP levels) method produced an overall correlation coefficient R2 = 0.966. The stress-strain single breath method produced R2 = 0.530. The single compartment single breath method produced R2 = 0.415. A combined method at single and multiple PEEP levels produced R2 = 0.963.

CONCLUSIONS

The results suggest that model-based, single breath and non-invasive approaches to estimating dFRC may be viable in a clinical scenario, ensuring no interruption to MV. The models provide a means of estimating dFRC at any PEEP level. However, model limitations and large estimation errors limit the use of the methods at very low PEEP.

Optimization of mechanical ventilator settings for pulmonary disease states

The selection of mechanical ventilator settings that ensure adequate oxygenation and carbon dioxide clearance while minimizing the risk of ventilator-associated lung injury (VALI) is a significant challenge for intensive-care clinicians. Current guidelines are largely based on previous experience combined with recommendations from a limited number of in vivo studies whose data are typically more applicable to populations than to individuals suffering from particular diseases of the lung. By combining validated computational models of pulmonary pathophysiology with global optimization algorithms, we generate in silico experiments to examine current practice and uncover optimal combinations of ventilator settings for individual patient and disease states. Formulating the problem as a multiobjective, multivariable constrained optimization problem, we compute settings of tidal volume, ventilation rate, inspiratory/expiratory ratio, positive end-expiratory pressure and inspired fraction of oxygen that optimally manage the tradeoffs between ensuring adequate oxygenation and carbon dioxide clearance and minimizing the risk of VALI for different pulmonary disease scenarios.

Physiological relevance and performance of a minimal lung model: an experimental study in healthy and acute respiratory distress syndrome model piglets

Background

Mechanical ventilation (MV) is the primary form of support for acute respiratory distress syndrome (ARDS) patients. However, intra- and inter- patient-variability reduce the efficacy of general protocols. Model-based approaches to guide MV can be patient-specific. A physiological relevant minimal model and its patient-specific performance are tested to see if it meets this objective above.

Methods

Healthy anesthetized piglets weighing 24.0 kg [IQR: 21.0-29.6] underwent a step-wise PEEP increase manoeuvre from 5cmH₂O to 20cmH₂O. They were ventilated under volume control using Engström Care Station (Datex, General Electric, Finland), with pressure, flow and volume profiles recorded. ARDS was then induced using oleic acid. The data were analyzed with a Minimal Model that identifies patient-specific mean threshold opening and closing pressure (TOP and TCP), and standard deviation (SD) of these TOP and TCP distributions. The trial and use of data were approved by the Ethics Committee of the Medical Faculty of the University of Liege, Belgium.

Results and discussions

3 of the 9 healthy piglets developed ARDS, and these data sets were included in this study. Model fitting error during inflation and deflation, in healthy or ARDS state is less than 5.0% across all subjects, indicating that the model captures the fundamental lung mechanics during PEEP increase. Mean TOP was 42.4cmH₂O [IQR: 38.2-44.6] at PEEP = 5cmH₂O and decreased with PEEP to 25.0cmH₂O [IQR: 21.5-27.1] at PEEP = 20cmH₂O. In contrast, TCP sees a reverse trend, increasing from 10.2cmH₂O [IQR: 9.0-10.4] to 19.5cmH₂O [IQR: 19.0-19.7]. Mean TOP increased from average 21.2-37.4cmH₂O to 30.4-55.2cmH₂O between healthy and ARDS subjects, reflecting the higher pressure required to recruit collapsed alveoli. Mean TCP was effectively unchanged.

Conclusion

The minimal model is capable of capturing physiologically relevant TOP, TCP and SD of both healthy and ARDS lungs. The model is able to track disease progression and the response to treatment.

Structural identifiability and practical applicability of an alveolar recruitment model for ARDS patients

Patient-specific mathematical models of respiratory mechanics can offer substantial insight into patient state and pulmonary dynamics that are not directly measurable. Thus, they offer significant potential to evaluate and guide patient-specific lung protective ventilator strategies for acute respiratory distress syndrome (ARDS) patients. To assure bedside applicability, the model must be computationally efficient and identifiable from the limited available data, while also capturing dominant dynamics and trends observed in ARDS patients. In this study, an existing static recruitment model is enhanced by considering alveolar distension and implemented in a novel time-continuous dynamic respiratory mechanics model. The model was tested for structural identifiability and a hierarchical gradient descent approach was used to fit the model to low-flow test responses of 12 ARDS patients. Finally, a comprehensive practical identifiability analysis was performed to evaluate the impact of data quality on the model parameters. Identified parameter values were physiologically plausible and very accurately reproduced the measured pressure responses. Structural identifiability of the model was proven, but practical identifiability analysis of the results showed a lack of convexity on the error surface indicating that successful parameter identification is currently not assured in all test sets. Overall, the model presented is physiologically and clinically relevant, captures ARDS dynamics, and uses clinically descriptive parameters. The patient-specific models show the ability to capture pulmonary dynamics directly relevant to patient condition and clinical guidance. These characteristics currently cannot be directly measured or established without such a validated model.

Model-based PEEP optimisation in mechanical ventilation

BACKGROUND

Acute Respiratory Distress Syndrome (ARDS) patients require mechanical ventilation (MV) for breathing support. Patient-specific PEEP is encouraged for treating different patients but there is no well established method in optimal PEEP selection.

METHODS

A study of 10 patients diagnosed with ALI/ARDS whom underwent recruitment manoeuvre is carried out. Airway pressure and flow data are used to identify patient-specific constant lung elastance (E lung) and time-variant dynamic lung elastance (E drs) at each PEEP level (increments of 5 cm H2O), for a single compartment linear lung model using integral-based methods. Optimal PEEP is estimated using E lung versus PEEP, Edrs-Pressure curve and E drs Area at minimum elastance (maximum compliance) and the inflection of the curves (diminishing return). Results are compared to clinically selected PEEP values. The trials and use of the data were approved by the New Zealand South Island Regional Ethics Committee.

RESULTS

Median absolute percentage fitting error to the data when estimating time-variant E drs is 0.9% (IQR = 0.5-2.4) and 5.6% [IQR: 1.8-11.3] when estimating constant E lung. Both E lung and E drs decrease with PEEP to a minimum, before rising, and indicating potential over-inflation. Median E drs over all patients across all PEEP values was 32.2 cmH2O/l [IQR: 26.1-46.6], reflecting the heterogeneity of ALI/ARDS patients, and their response to PEEP, that complicates standard approaches to PEEP selection. All E drs-Pressure curves have a clear inflection point before minimum E drs, making PEEP selection straightforward. Model-based selected PEEP using the proposed metrics were higher than clinically selected values in 7/10 cases.

CONCLUSION

Continuous monitoring of the patient-specific E lung and E drs and minimally invasive PEEP titration provide a unique, patient-specific and physiologically relevant metric to optimize PEEP selection with minimal disruption of MV therapy.