Development and optimisation of stochastic targeted (STAR) glycaemic control for pre-term infants in neonatal intensive care

Stochastic integrated model-based protocol for volume-controlled ventilation setting

BACKGROUND AND OBJECTIVE

Mechanical ventilation (MV) is the primary form of care for respiratory failure patients. MV settings are based on general clinical guidelines, intuition, and experience. This approach is not patient-specific and patients may thus experience suboptimal, potentially harmful MV care. This study presents the Stochastic integrated VENT (SiVENT) protocol which combines model-based approaches of the VENT protocol from previous works, with stochastic modelling to take the variation of patient respiratory elastance over time into consideration.

METHODS

A stochastic model of Ers is integrated into the VENT protocol from previous works to develop the SiVENT protocol, to account for both intra- and inter-patient variability. A cohort of 20 virtual MV patients based on retrospective patient data are used to validate the performance of this method for volume-controlled (VC) ventilation. A performance evaluation was conducted where the SiVENT and VENT protocols were implemented in 1080 instances each to compare the two protocols and evaluate the difference in reduction of possible MV settings achieved by each.

RESULTS

From an initial number of 189,000 possible MV setting combinations, the VENT protocol reduced this number to a median of 10,612, achieving a reduction of 94.4% across the cohort. With the integration of the stochastic model component, the SiVENT protocol reduced this number from 189,000 to a median of 9329, achieving a reduction of 95.1% across the cohort. The SiVENT protocol reduces the number of possible combinations provided to the user by more than 1000 combinations as compared to the VENT protocol.

CONCLUSIONS

Adding a stochastic model component into a model-based approach to selecting MV settings improves the ability of a decision support system to recommend patient-specific MV settings. It specifically considers inter- and intra-patient variability in respiratory elastance and eliminates potentially harmful settings based on clinically recommended pressure thresholds. Clinical input and local protocols can further reduce the number of safe setting combinations. The results for the SiVENT protocol justify further investigation of its prediction accuracy and clinical validation trials.

Insulin sensitivity in critically ill patients: are women more insulin resistant?

Background

Glycaemic control (GC) in intensive care unit is challenging due to significant inter- and intra-patient variability, leading to increased risk of hypoglycaemia. Recent work showed higher insulin resistance in female preterm neonates. This study aims to determine if there are differences in inter- and intra-patient metabolic variability between sexes in adults, to gain in insight into any differences in metabolic response to injury. Any significant difference would suggest GC and randomised trial design should consider sex differences to personalise care.

Methods

Insulin sensitivity (SI) levels and variability are identified from retrospective clinical data for men and women. Data are divided using 6-h blocks to capture metabolic evolution over time. In total, 91 male and 54 female patient GC episodes of minimum 24 h are analysed. Hypothesis testing is used to determine whether differences are significant (P < 0.05), and equivalence testing is used to assess whether these differences can be considered equivalent at a clinical level. Data are assessed for the raw cohort and in 100 Monte Carlo simulations analyses where the number of men and women are equal.

Results

Demographic data between females and males were all similar, including GC outcomes (safety from hypoglycaemia and high (> 50%) time in target band). Females had consistently significantly lower SI levels than males, and this difference was not clinically equivalent. However, metabolic variability between sexes was never significantly different and always clinically equivalent. Thus, inter-patient variability was significantly different between males and females, but intra-patient variability was equivalent.

Conclusion

Given equivalent intra-patient variability and significantly greater insulin resistance, females can receive the same benefit from safe, effective GC as males, but may require higher insulin doses to achieve the same glycaemia. Clinical trials should consider sex differences in protocol design and outcome analyses.

Safe and effective glycaemic control in premature infants: observational clinical results from the computerised STAR-GRYPHON protocol

OBJECTIVE

Previous studies examine clinical outcomes of insulin therapy in neonatal intensive care units (NICUs), without first developing safe and effective control protocols. This research quantifies the safety and performance of a computerised model-based control algorithmSTAR-GRYPHON (Stochastic TARgeted Glucose Regulation sYstem to Prevent Hyper- and hypO-glycaemia in Neonates).

DESIGN

Retrospective observational study of glycaemic control in very/extremely low birthweight infants treated with insulin from Christchurch Women's Hospital NICU between January 2013 and June 2017. Blood glucose (BG) outcomes and control performance is compared with retrospective data (n=22) and literature.

INTERVENTIONS

Insulin infusion doses were calculated from 3 to 4 hourly BG measurements using a computerised model-based control algorithm, STAR-GRYPHON.

MAIN OUTCOME MEASURES

Mean BG, time in targeted range and incidence of hypoglycaemia.

RESULTS

STAR-GRYPHON (n=35) had lower mean BG concentration (7.0mmol/L vs 7.9 mmol/L), higher %BG within the 4.0-8.0 mmol/L target range (71.1% vs 50.9%) and lower %BG <4.0 mmol/L (0.6% vs 2.1%). There were only 2 BG <2.6 mmol/L (over n=2, 5.5% of patients, 0.03% of all BG outcomes), one of which may be attributed to clinical error. These results show better control to target and lower incidence of hypoglycaemia than most literature results from intensive insulin therapy protocols or study groups in children and infants.

CONCLUSIONS

Model-based protocols can safely and effectively control BG in very premature infants and should be used in future studies to determine the effect of insulin therapy on clinical outcomes.

Model-based glycaemic control: methodology and initial results from neonatal intensive care

Very/extremely premature infants often experience glycaemic dysregulation, resulting in abnormally elevated (hyperglycaemia) or low (hypoglycaemia) blood glucose (BG) concentrations, due to prematurity, stress, and illness. STAR-GRYPHON is a computerised protocol that utilises a model-based insulin sensitivity parameter to directly tailor therapy for individual patients and their changing conditions, unlike other common insulin protocols in this cohort. From January 2013 to January 2015, 13 patients totalling 16 hyperglycaemic control episodes received insulin under STAR-GRYPHON. A significant improvement in control was achieved in comparison to a retrospective cohort, with a 26% absolute improvement in BG within the targeted range and no hypoglycaemia. This improvement was obtained predominantly due to the reduction of hyperglycaemia (%BG>10.0 mmol/l: 5.6 vs. 17.7%, p<0.001), and lowering of the median per-patient BG [6.9 (6.1-7.9) vs. 7.8 (6.6-9.1) mmol/l, p<0.001, Mann-Witney U test]. While cohort-wide control results show good control overall, there is high intra-patient variability in BG behaviour, resulting in overly conservative treatments for some patients. Patient insulin sensitivity differs between and within patients over time, with some patients having stable insulin sensitivity, while others change rapidly. These results demonstrate the trade-off between safety and performance in a highly variable and fragile cohort.

Untangling glycaemia and mortality in critical care

Background

Hyperglycaemia is associated with adverse outcomes in the intensive care unit, and initial studies suggested outcome benefits of glycaemic control (GC). However, subsequent studies often failed to replicate these results, and they were often unable to achieve consistent, safe control, raising questions about the benefit or harm of GC as well as the nature of the association of glycaemia with mortality and clinical outcomes. In this study, we evaluated if non-survivors are harder to control than survivors and determined if glycaemic outcome is a function of patient condition and eventual outcome or of the glycaemic control provided.

Methods

Clinically validated, model-based, hour-to-hour insulin sensitivity (SI) and its hour-to-hour variability (%ΔSI) were identified over the first 72 h of therapy in 145 patients (119 survivors, 26 non-survivors). In hypothesis testing, we compared distributions of SI and %ΔSI in 6-hourly blocks for survivors and non-survivors. In equivalence testing, we assessed if differences in these distributions, based on blood glucose measurement error, were clinically significant.

Results

SI level was never equivalent between survivors and non-survivors (95% CI of percentage difference in medians outside ±12%). Non-survivors had higher SI, ranging from 9% to 47% higher overall in 6-h blocks, and this difference became statistically significant as glycaemic control progressed. %ΔSI was equivalent between survivors and non-survivors for all 6-hourly blocks (95% CI of difference in medians within ±12%) and decreased in general over time as glycaemic control progressed.

Conclusions

Whereas non-survivors had higher SI levels, variability was equivalent to that of survivors over the first 72 h. These results indicate survivors and non-survivors are equally controllable, given an effective glycaemic control protocol, suggesting that glycaemia level and variability, and thus the association between glycaemia and outcome, are essentially determined by the control provided rather than by underlying patient or metabolic condition.

Electronic supplementary material

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

Insulin kinetics and the Neonatal Intensive Care Insulin-Nutrition-Glucose (NICING) model

BACKGROUND

Models of human glucose-insulin physiology have been developed for a range of uses, with similarly different levels of complexity and accuracy. STAR (Stochastic Targeted) is a model-based approach to glycaemic control. Elevated blood glucose concentrations (hyperglycaemia) are a common complication of stress and prematurity in very premature infants, and have been associated with worsened outcomes and higher mortality. This research identifies and validates the model parameters for model-based glycaemic control in neonatal intensive care.

METHODS

C-peptide, plasma insulin, and BG from a cohort of 41 extremely pre-term (median age 27.2 [26.2-28.7] weeks) and very low birth weight infants (median birth weight 839 [735-1000] g) are used alongside C-peptide kinetic models to identify model parameters associated with insulin kinetics in the NICING (Neonatal Intensive Care Insulin-Nutrition-Glucose) model. A literature analysis is used to determine models of kidney clearance and body fluid compartment volumes. The full, final NICING model is validated by fitting the model to a cohort of 160 glucose, insulin, and nutrition data records from extremely premature infants from two different NICUs (neonatal intensive care units).

RESULTS

Six model parameters related to insulin kinetics were identified. The resulting NICING model is more physiologically descriptive than prior model iterations, including clearance pathways of insulin via the liver and kidney, rather than a lumped parameter. In addition, insulin diffusion between plasma and interstitial spaces is evaluated, with differences in distribution volume taken into consideration for each of these spaces. The NICING model was shown to fit clinical data well, with a low model fit error similar to that of previous model iterations.

CONCLUSIONS

Insulin kinetic parameters have been identified, and the NICING model is presented for glycaemic control neonatal intensive care. The resulting NICING model is more complex and physiologically relevant, with no loss in bedside-identifiability or ability to capture and predict metabolic dynamics.

Brain mass estimation by head circumference and body mass methods in neonatal glycaemic modelling and control

INTRODUCTION

Hyperglycaemia is a common complication of stress and prematurity in extremely low-birth-weight infants. Model-based insulin therapy protocols have the ability to safely improve glycaemic control for this group. Estimating non-insulin-mediated brain glucose uptake by the central nervous system in these models is typically done using population-based body weight models, which may not be ideal.

METHOD

A head circumference-based model that separately treats small-for-gestational-age (SGA) and appropriate-for-gestational-age (AGA) infants is compared to a body weight model in a retrospective analysis of 48 patients with a median birth weight of 750g and median gestational age of 25 weeks. Estimated brain mass, model-based insulin sensitivity (SI) profiles, and projected glycaemic control outcomes are investigated. SGA infants (5) are also analyzed as a separate cohort.

RESULTS

Across the entire cohort, estimated brain mass deviated by a median 10% between models, with a per-patient median difference in SI of 3.5%. For the SGA group, brain mass deviation was 42%, and per-patient SI deviation 13.7%. In virtual trials, 87-93% of recommended insulin rates were equal or slightly reduced (Δ<0.16mU/h) under the head circumference method, while glycaemic control outcomes showed little change.

CONCLUSION

The results suggest that body weight methods are not as accurate as head circumference methods. Head circumference-based estimates may offer improved modelling accuracy and a small reduction in insulin administration, particularly for SGA infants.

On the problem of patient-specific endogenous glucose production in neonates on stochastic targeted glycemic control

BACKGROUND

Both stress and prematurity can induce hyperglycemia in the neonatal intensive care unit, which, in turn, is associated with worsened outcomes. Endogenous glucose production (EGP) is the formation of glucose by the body from substrates and contributes to blood glucose (BG) levels. Due to the inherent fragility of the extremely low birth weight (ELBW) neonates, true fasting EGP cannot be explicitly determined, introducing uncertainty into glycemic models that rely on quantifying glucose sources. Stochastic targeting, or STAR, is one such glycemic control framework.

METHODS

A literature review was carried out to gather metabolic and EGP values on preterm infants with a gestational age (GA) &lt;32 weeks and a birth weight (BW) &lt;2 kg. The data were analyzed for EGP trends with BW, GA, BG, plasma insulin, and glucose infusion (GI) rates. Trends were modeled and compared with a literature-derived range of population constant EGP models using clinically validated virtual trials on retrospective clinical data.

RESULTS

No clear relationship was found for EGP and BW, GA, or plasma insulin. Some evidence of suppression of EGP with increasing GI or BG was seen. Virtual trial results showed that population-constant EGP models fit clinical data best and gave tighter control performance to a target band in virtual trials.

CONCLUSIONS

Variation in EGP cannot easily be quantified, and EGP is sufficiently modeled as a population constant in the neonatal intensive care insulin-nutrition-glucose model. Analysis of the clinical data and fitting error suggests that ELBW hyperglycemic preterm neonates have unsuppressed EGP in the higher range than that seen in literature.

Nasogastric aspiration as an indicator for feed absorption in model-based glycemic control in neonatal intensive care

BACKGROUND

STAR (stochastic targeted) is a glycemic control model-based framework for critically ill neonates that has shown benefits in reducing hypoglycemia and hyperglycemia. STAR uses a stochastic matrix method to forecast future changes in a patient's insulin sensitivity and then applies this result to a physiological model to select an optimal insulin treatment. Nasogastric aspiration may be used as an indicator to suggest periods of care when enteral feed absorption is compromised, improving the performance of glycemic control. An analysis has been carried out to investigate the effect of poorly absorbed feeds on glycemic control.

METHOD

Clinical data were collected from eight patients on insulin therapy and enteral feed, which included large or significantly milky aspirates. Patients had a median gestational age of 25 weeks and postnatal age of 5.5 days. Virtual patients were created using the NICING model, and insulin sensitivity (SI) profiles were fit. Alternative feed profiles were generated whereby enteral feed absorption was redistributed with time to account for poor feed absorption. The effect of poor feed absorption, as indicated by aspirates, is investigated.

RESULTS

The average percentage change of SI 4 h before a significant aspirate was 1.16%, and 1.49% in the 4 h following the aspirate. No distinct relationship was found between the fractional change in SI and the volume of the aspirate. Accounting for aspirates had a clinically negligible impact on glycemic control in virtual trials.

CONCLUSION

Accounting for aspirates by manipulating enteral feed profiles had a minimal influence on both modeling and controlling glycemia in neonates. The impact of this method is clinically insignificant, suggesting that a population constant for the rate of glucose absorption in the gut adequately models feed absorption within the STAR framework.