Evaluation of model-based methods in estimating respiratory mechanics in the presence of variable patient effort

Virtual patient with temporal evolution for mechanical ventilation trial studies: A stochastic model approach

BACKGROUND AND OBJECTIVE

Healthcare datasets are plagued by issues of data scarcity and class imbalance. Clinically validated virtual patient (VP) models can provide accurate in-silico representations of real patients and thus a means for synthetic data generation in hospital critical care settings. This research presents a realistic, time-varying mechanically ventilated respiratory failure VP profile synthesised using a stochastic model.

METHODS

A stochastic model was developed using respiratory elastance (Ers) data from two clinical cohorts and averaged over 30-minute time intervals. The stochastic model was used to generate future Ers data based on current Ers values with added normally distributed random noise. Self-validation of the VPs was performed via Monte Carlo simulation and retrospective Ers profile fitting. A stochastic VP cohort of temporal Ers evolution was synthesised and then compared to an independent retrospective patient cohort data in a virtual trial across several measured patient responses, where similarity of profiles validates the realism of stochastic model generated VP profiles.

RESULTS

A total of 120,000 3-hour VPs for pressure control (PC) and volume control (VC) ventilation modes are generated using stochastic simulation. Optimisation of the stochastic simulation process yields an ideal noise percentage of 5-10% and simulation iteration of 200,000 iterations, allowing the simulation of a realistic and diverse set of Ers profiles. Results of self-validation show the retrospective Ers profiles were able to be recreated accurately with a mean squared error of only 0.099 [0.009-0.790]% for the PC cohort and 0.051 [0.030-0.126]% for the VC cohort. A virtual trial demonstrates the ability of the stochastic VP cohort to capture Ers trends within and beyond the retrospective patient cohort providing cohort-level validation.

CONCLUSION

VPs capable of temporal evolution demonstrate feasibility for use in designing, developing, and optimising bedside MV guidance protocols through in-silico simulation and validation. Overall, the temporal VPs developed using stochastic simulation alleviate the need for lengthy, resource intensive, high cost clinical trials, while facilitating statistically robust virtual trials, ultimately leading to improved patient care and outcomes in mechanical ventilation.

Computational physiological models for individualised mechanical ventilation: a systematic literature review focussing on quality, availability, and clinical readiness

BACKGROUND

Individualised optimisation of mechanical ventilation (MV) remains cumbersome in modern intensive care medicine. Computerised, model-based support systems could help in tailoring MV settings to the complex interactions between MV and the individual patient's pathophysiology. Therefore, we critically appraised the current literature on computational physiological models (CPMs) for individualised MV in the ICU with a focus on quality, availability, and clinical readiness.

METHODS

A systematic literature search was conducted on 13 February 2023 in MEDLINE ALL, Embase, Scopus and Web of Science to identify original research articles describing CPMs for individualised MV in the ICU. The modelled physiological phenomena, clinical applications, and level of readiness were extracted. The quality of model design reporting and validation was assessed based on American Society of Mechanical Engineers (ASME) standards.

RESULTS

Out of 6,333 unique publications, 149 publications were included. CPMs emerged since the 1970s with increasing levels of readiness. A total of 131 articles (88%) modelled lung mechanics, mainly for lung-protective ventilation. Gas exchange (n = 38, 26%) and gas homeostasis (n = 36, 24%) models had mainly applications in controlling oxygenation and ventilation. Respiratory muscle function models for diaphragm-protective ventilation emerged recently (n = 3, 2%). Three randomised controlled trials were initiated, applying the Beacon and CURE Soft models for gas exchange and PEEP optimisation. Overall, model design and quality were reported unsatisfactory in 93% and 21% of the articles, respectively.

CONCLUSION

CPMs are advancing towards clinical application as an explainable tool to optimise individualised MV. To promote clinical application, dedicated standards for quality assessment and model reporting are essential. Trial registration number PROSPERO- CRD42022301715 . Registered 05 February, 2022.

Mixed-integer quadratic programming approach for noninvasive estimation of respiratory effort profile during pressure support ventilation

Information about respiratory mechanics such as resistance, elastance, and muscular pressure is important to mitigate ventilator-induced lung injury. Particularly during pressure support ventilation, the available options to quantify breathing effort and calculate respiratory system mechanics are often invasive or complex. We herein propose a robust and flexible estimation of respiratory effort better than current methods. We developed a method for non-invasively estimating breathing effort using only flow and pressure signals. Mixed-integer quadratic programming (MIQP) was employed, and the binary variables were the switching moments of the respiratory effort waveform. Mathematical constraints, based on ventilation physiology, were set for some variables to restrict feasible solutions. Simulated and patient data were used to verify our method, and the results were compared to an established estimation methodology. Our algorithm successfully estimated the respiratory effort, resistance, and elastance of the respiratory system, resulting in more robust performance and faster solver times than a previously proposed algorithm that used quadratic programming (QP) techniques. In a numerical simulation benchmark, the worst-case errors for resistance and elastance were 25% and 23% for QP versus <0.1% and <0.1% for MIQP, whose solver times were 4.7 s and 0.5 s, respectively. This approach can estimate several breathing effort profiles and identify the respiratory system's mechanical properties in invasively ventilated critically ill patients.

Predicting mechanically ventilated patients future respiratory system elastance - A stochastic modelling approach

BACKGROUND AND OBJECTIVE

Respiratory mechanics of mechanically ventilated patients evolve significantly with time, disease state and mechanical ventilation (MV) treatment. Existing deterministic data prediction methods fail to comprehensively describe the multiple sources of heterogeneity of biological systems. This research presents two respiratory mechanics stochastic models with increased prediction accuracy and range, offering improved clinical utility in MV treatment.

METHODS

Two stochastic models (SM2 and SM3) were developed using retrospective patient respiratory elastance (E_rs) from two clinical cohorts which were averaged over time intervals of 10 and 30 min respectively. A stochastic model from a previous study (SM1) was used to benchmark performance. The stochastic models were clinically validated on an independent retrospective clinical cohort of 14 patients. Differences in predictive ability were evaluated using the difference in percentile lines and cumulative distribution density (CDD) curves.

RESULTS

Clinical validation shows all three models captured more than 98% (median) of future E_rs data within the 5th - 95th percentile range. Comparisons of stochastic model percentile lines reported a maximum mean absolute percentage difference of 5.2%. The absolute differences of CDD curves were less than 0.25 in the ranges of 5 < E_rs (cmH₂O/L) < 85, suggesting similar predictive capabilities within this clinically relevant E_rs range.

CONCLUSION

The new stochastic models significantly improve prediction, clinical utility, and thus feasibility for synchronisation with clinical interventions. Paired with other MV protocols, the stochastic models developed can potentially form part of decision support systems, providing guided, personalised, and safe MV treatment.

Virtual patient framework for the testing of mechanical ventilation airway pressure and flow settings protocol

BACKGROUND AND OBJECTIVE

Model-based and personalised decision support systems are emerging to guide mechanical ventilation (MV) treatment for respiratory failure patients. However, model-based treatments require resource-intensive clinical trials prior to implementation. This research presents a framework for generating virtual patients for testing model-based decision support, and direct use in MV treatment.

METHODS

The virtual MV patient framework consists of 3 stages: 1) Virtual patient generation, 2) Patient-level validation, and 3) Virtual clinical trials. The virtual patients are generated from retrospective MV patient data using a clinically validated respiratory mechanics model whose respiratory parameters (respiratory elastance and resistance) capture patient-specific pulmonary conditions and responses to MV care over time. Patient-level validation compares the predicted responses from the virtual patient to their retrospective results for clinically implemented MV settings and changes to care. Patient-level validated virtual patients create a platform to conduct virtual trials, where the safety of closed-loop model-based protocols can be evaluated.

RESULTS

This research creates and presents a virtual patient platform of 100 virtual patients generated from retrospective data. Patient-level validation reported median errors of 3.26% for volume-control and 6.80% for pressure-control ventilation mode. A virtual trial on a model-based protocol demonstrates the potential efficacy of using virtual patients for prospective evaluation and testing of the protocol.

CONCLUSION

The virtual patient framework shows the potential to safely and rapidly design, develop, and optimise new model-based MV decision support systems and protocols using clinically validated models and computer simulation, which could ultimately improve patient care and outcomes in MV.

Algebraic approximation of the distributed model for the pressure drop in the respiratory airways

The complexity of the human respiratory system causes that one of the main methods of analyzing the dynamic pulmonary phenomena and interpreting experimental results are simulations of its computational models. Among the most compound elements of these models, apart from the bronchial tree structure, is the phenomenon of flow limitation in flexible bronchi, which causes them to collapse with increasing flow, thus their properties, such as resistance, compliance and inertance, are highly nonlinear and time-varying. Commonly, this phenomenon is ignored, or a distributed model for the airway pressure drop is applied, simulated with a modified numerical solver of this differential equation (ODE). The disadvantages of this solution are the problems with taking into account the inherent singularity of the model and the long computation time due to iterative nature of the ODE procedure. The aim of the work was to derive an algebraic approximation of this distributed model, suitable for implementation in continuous dynamic models, to validate it by comparing the results of simulations with the respiratory system model including approximate and original (ODE solver) numerical procedures, as well as to evaluate possible acceleration of calculations. All simulations, including spontaneous breathing, mechanical ventilation with the optimal ventilatory waveform and forced expiration, proved that algebraic approximation yielded results negligibly differing from the ODE solution, and shortened the computation time by an order. The proposed approach is an attractive alternative in the case of computer implementations of pulmonary models, where simulations of flows and pressures in the complex respiratory system are of primary importance.

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.

Reconstructing asynchrony for mechanical ventilation using a hysteresis loop virtual patient model

BACKGROUND

Patient-specific lung mechanics during mechanical ventilation (MV) can be identified from measured waveforms of fully ventilated, sedated patients. However, asynchrony due to spontaneous breathing (SB) effort can be common, altering these waveforms and reducing the accuracy of identified, model-based, and patient-specific lung mechanics.

METHODS

Changes in patient-specific lung elastance over a pressure-volume (PV) loop, identified using hysteresis loop analysis (HLA), are used to detect the occurrence of asynchrony and identify its type and pattern. The identified HLA parameters are then combined with a nonlinear mechanics hysteresis loop model (HLM) to extract and reconstruct ventilated waveforms unaffected by asynchronous breaths. Asynchrony magnitude can then be quantified using an energy-dissipation metric, Easyn, comparing PV loop area between model-reconstructed and original, altered asynchronous breathing cycles. Performance is evaluated using both test-lung experimental data with a known ground truth and clinical data from four patients with varying levels of asynchrony.

RESULTS

Root mean square errors for reconstructed PV loops are within 5% for test-lung experimental data, and 10% for over 90% of clinical data. Easyn clearly matches known asynchrony magnitude for experimental data with RMS errors < 4.1%. Clinical data performance shows 57% breaths having Easyn > 50% for Patient 1 and 13% for Patient 2. Patient 3 only presents 20% breaths with Easyn > 10%. Patient 4 has Easyn = 0 for 96% breaths showing accuracy in a case without asynchrony.

CONCLUSIONS

Experimental test-lung validation demonstrates the method's reconstruction accuracy and generality in controlled scenarios. Clinical validation matches direct observations of asynchrony in incidence and quantifies magnitude, including cases without asynchrony, validating its robustness and potential efficacy as a clinical real-time asynchrony monitoring tool.

Protocol conception for safe selection of mechanical ventilation settings for respiratory failure Patients

BACKGROUND AND OBJECTIVE

Mechanical ventilation is the primary form of care provided to respiratory failure patients. Limited guidelines and conflicting results from major clinical trials means selection of mechanical ventilation settings relies heavily on clinician experience and intuition. Determining optimal mechanical ventilation settings is therefore difficult, where non-optimal mechanical ventilation can be deleterious. To overcome these difficulties, this research proposes a model-based method to manage the wide range of possible mechanical ventilation settings, while also considering patient-specific conditions and responses.

METHODS

This study shows the design and development of the "VENT" protocol, which integrates the single compartment linear lung model with clinical recommendations from landmark studies, to aid clinical decision-making in selecting mechanical ventilation settings. Using retrospective breath data from a cohort of 24 patients, 3,566 and 2,447 clinically implemented VC and PC settings were extracted respectively. Using this data, a VENT protocol application case study and clinical comparison is performed, and the prediction accuracy of the VENT protocol is validated against actual measured outcomes of pressure and volume.

RESULTS

The study shows the VENT protocols' potential use in narrowing an overwhelming number of possible mechanical ventilation setting combinations by up to 99.9%. The comparison with retrospective clinical data showed that only 33% and 45% of clinician settings were approved by the VENT protocol. The unapproved settings were mainly due to exceeding clinical recommended settings. When utilising the single compartment model in the VENT protocol for forecasting peak pressures and tidal volumes, median [IQR] prediction error values of 0.75 [0.31 - 1.83] cmH₂O and 0.55 [0.19 - 1.20] mL/kg were obtained.

CONCLUSIONS

Comparing the proposed protocol with retrospective clinically implemented settings shows the protocol can prevent harmful mechanical ventilation setting combinations for which clinicians would be otherwise unaware. The VENT protocol warrants a more detailed clinical study to validate its potential usefulness in a clinical setting.

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.

Stochastic Modelling of Respiratory System Elastance for Mechanically Ventilated Respiratory Failure Patients

While lung protective mechanical ventilation (MV) guidelines have been developed to avoid ventilator-induced lung injury (VILI), a one-size-fits-all approach cannot benefit every individual patient. Hence, there is significant need for the ability to provide patient-specific MV settings to ensure safety, and optimise patient care. Model-based approaches enable patient-specific care by identifying time-varying patient-specific parameters, such as respiratory elastance, Ers, to capture inter- and intra-patient variability. However, patient-specific parameters evolve with time, as a function of disease progression and patient condition, making predicting their future values crucial for recommending patient-specific MV settings. This study employs stochastic modelling to predict future Ers values using retrospective patient data to develop and validate a model indicating future intra-patient variability of Ers. Cross validation results show stochastic modelling can predict future elastance ranges with 92.59 and 68.56% of predicted values within the 5-95% and the 25-75% range, respectively. This range can be used to ensure patients receive adequate minute ventilation should elastance rise and minimise the risk of VILI should elastance fall. The results show the potential for model-based protocols using stochastic model prediction of future Ers values to provide safe and patient-specific MV. These results warrant further investigation to validate its clinical utility.

Model-based estimation of negative inspiratory driving pressure in patients receiving invasive NAVA mechanical ventilation

BACKGROUND

Optimisation of mechanical ventilation (MV) and weaning requires insight into underlying patient breathing effort. Current identifiable models effectively describe lung mechanics, such as elastance (E) and resistance (R) at the bedside in sedated patients, but are less effective when spontaneous breathing is present. This research derives and regularises a single compartment model to identify patient-specific inspiratory effort.

METHODS

Constrained second-order b-spline basis functions (knot width 0.05 s) are used to describe negative inspiratory drive (Pp, cmH₂O) as a function of time. Breath-breath Pp are identified with single E and R values over inspiration and expiration from n = 20 breaths for N = 22 patients on NAVA ventilation. Pp is compared to measured electrical activity of the diaphragm (Eadi) and published results.

RESULTS

Average per-patient root-mean-squared model fit error was (median [interquartile range, IQR]) 0.9 [0.6-1.3] cmH₂O, and average per-patient median Pp was -3.9 [-4.5- -3.0] cmH₂O, with range -7.9 - -1.9 cmH₂O. Per-patient E and R were 16.4 [13.6-21.8] cmH₂O/L and 9.2 [6.4-13.1] cmH₂O.s/L, respectively. Most patients showed an inspiratory volume threshold beyond which Pp started to return to baseline, and Pp at peak Eadi (end-inspiration) was often strongly correlated with peak Eadi (R²=0.25-0.86). Similarly, average transpulmonary pressure was consistent breath-breath in most patients, despite differences in peak Eadi and thus peak airway pressure.

CONCLUSIONS

The model-based inspiratory effort aligns with electrical muscle activity and published studies showing neuro-muscular decoupling as a function of pressure and/or volume. Consistency in coupling/dynamics were patient-specific. Quantification of patient and ventilator work of breathing contributions may aid optimisation of MV modes and weaning.

Patient asynchrony modelling during controlled mechanical ventilation therapy

BACKGROUND AND OBJECTIVE

Mechanical ventilation therapy of respiratory failure patients can be guided by monitoring patient-specific respiratory mechanics. However, the patient's spontaneous breathing effort during controlled ventilation changes airway pressure waveform and thus affects the model-based identification of patient-specific respiratory mechanics parameters. This study develops a model to estimate respiratory mechanics in the presence of patient effort.

METHODS

Gaussian effort model (GEM) is a derivative of the single-compartment model with basis function. GEM model uses a linear combination of basis functions to model the nonlinear pressure waveform of spontaneous breathing patients. The GEM model estimates respiratory mechanics such as Elastance and Resistance along with the magnitudes of basis functions, which accounts for patient inspiratory effort.

RESULTS AND DISCUSSION

The GEM model was tested using both simulated data and a retrospective observational clinical trial patient data. GEM model fitting to the original airway pressure waveform is better than any existing models when reverse triggering asynchrony is present. The fitting error of GEM model was less than 10% for both simulated data and clinical trial patient data.

CONCLUSION

GEM can capture the respiratory mechanics in the presence of patient effect in volume control ventilation mode and also can be used to assess patient-ventilator interaction. This model determines basis functions magnitudes, which can be used to simulate any waveform of patient effort pressure for future studies. The estimation of parameter identification GEM model can further be improved by constraining the parameters within a physiologically plausible range during least-square nonlinear regression.

Respiratory muscular response to obstructive maneuvers in non-invasively ventilated healthy subjects

The hypothesis of this study is that muscular activity measured through surface electromyography (sEMG) is useful to estimate the work of breathing (WOB) and respiratory mechanics. Thirty-two healthy volunteers were non-invasively ventilated, and an airflow resistor was attached to the airway circuit. sEMG signals from diaphragm, intercostal and sternocleidomastoid muscles were processed and compared with WOB changes. The airway resistance was increased from a median of 9.58 to 22.51 cmH₂O/L/s adding a resistance of 20 cmH₂O/L/s, achieving the lower compliance too. The respiratory mechanics changes implied linear increases in WOB, with Pearson correlation of 88.43% respect to changes in resistance. Muscles increased their activity in agreement with changes of WOB, being higher the increment in diaphragm followed by sternocleidomastoid. The non-invasively monitored respiratory muscles activity allowed evaluating the changes in WOB when it depends on addition of obstructive loads, confirming that it could be used to improve the available respiratory mechanics and WOB monitoring tools.