Journal of Healthcare Engineering
Volume 2020 (2020), Article ID 1414597, 11 pages
https://doi.org/10.1155/2020/1414597

Research Article

Dynamic Rule-Based Algorithm to Tune Insulin-on-Board Constraints for a Hybrid Artificial Pancreas System

Arthur Bertachi,1,2 Lyvia Biagi,1,2 Aleix Beneyto,¹ and Josep Vehí1,3

¹University of Girona, Girona, Spain
²Federal University of Technology–Paraná (UTFPR), Guarapuava, Brazil
³Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas, Madrid, Spain

Correspondence should be addressed to Josep Vehí; josep.vehi@udg.edu

Received 31 October 2018; Revised 16 April 2019; Accepted 14 May 2019; Published 11 January 2020

Academic Editor: Jiann-Shing Shieh

Copyright © 2020 Arthur Bertachi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The artificial pancreas (AP) is a system intended to control blood glucose levels through automated insulin infusion, reducing the burden of subjects with type 1 diabetes to manage their condition. To increase patients’ safety, some systems limit the allowed amount of insulin active in the body, known as insulin-on-board (IOB). The safety auxiliary feedback element (SAFE) layer has been designed previously to avoid overreaction of the controller and thus avoiding hypoglycemia. In this work, a new method, so-called “dynamic rule-based algorithm,” is presented in order to adjust the limits of IOB in real time. The algorithm is an extension of a previously designed method which aimed to adjust the limits of IOB for a meal with 60 grams of carbohydrates (CHO). The proposed method is intended to be applied on hybrid AP systems during 24 h operation. It has been designed by combining two different strategies to set IOB limits for different situations: (1) fasting periods and (2) postprandial periods, regardless of the size of the meal. The UVa/Padova simulator is considered to assess the performance of the method, considering challenging scenarios. In silico results showed that the method is able to reduce the time spent in hypoglycemic range, improving patients’ safety, which reveals the feasibility of the approach to be included in different control algorithms.

1. Introduction

Type 1 diabetes (T1D) is a chronic condition in which the pancreatic beta-cells either stop or reduce drastically the production of insulin. Insulin is a hormone whose function is to facilitate the glucose uptake from the bloodstream into the cells to be used or stored. Subjects with absence of insulin in the body face very high levels of blood glucose (BG) (hyperglycemia), which can lead to long-term micro- and macrovascular complications [1, 2]. Therefore, subjects living with T1D must inject insulin exogenously in order to regulate blood in a lifelong challenge [3], and intensive insulin therapy reduces the risk of long-term complications [4]. But maintaining blood glucose levels into near-normoglycemia is not a trivial task, and if insulin is overdosed, BG may fall to dangerously low levels (hypoglycemia), which can lead to serious hazards, such as diabetic coma or even death [5].

Over the last years, researchers have been working towards a closed-loop system to control BG automatically [6]. This system, known as artificial pancreas (AP), is usually composed of a continuous glucose monitor (CGM), a control algorithm, and continuous subcutaneous insulin infusion through a pump. Insulin-only AP systems consider only insulin infusion to control BG, and dual-hormone AP systems also consider glucagon infusion to elevate BG to reduce the risk of hypoglycemia. AP systems can also be classified into two different categories based on the degree of automation: hybrid closed-loop system, in which subjects are involved in the control loop and must announce meals or other disturbances to anticipate their effects, and fully closed-loop systems, where no actions are required from the patients [7, 8].

One of the main challenges of the AP system is achieving postprandial glucose control mainly because the insulin absorption through the subcutaneous tissue is slower than the appearance of glucose in the blood after a meal [9]. Due to this difference in dynamics of insulin action and carbohydrate (CHO) absorption, attempts to avoid hyperglycemic peaks are usually accompanied by hypoglycemic excursions [10]. Several approaches have been tested to overcome such issue [11–14], but while an ultrarapid insulin analogue is not available [15], postprandial control using subcutaneous route will continue to be a challenging situation for closed-loop systems. Diverse studies have included estimations of insulin concentration in the body to avoid excessive insulin stacking [16–19]; however, hypo- and hyperglycemia are still a hazard for AP systems, and novel approaches are still required.

Revert and colleagues introduced a safety auxiliary feedback element, so-called SAFE layer [20], to limit excessive insulin in the subcutaneous tissue, i.e., insulin-on-board (IOB). This layer is based on the sliding mode reference conditioning technique [21] and acts on the glucose reference signal when a specific constraint, related with the maximum IOB allowed (), is violated. Such technique has already been applied in different control schemes [22–25], but so far it is not clear the best methodology to tune . The selection of the constraint is critical in the design of the closed-loop system. This parameter regulates insulin infusion based on an estimation of the IOB. As higher , more insulin the controller will be allowed to deliver. Considering that once insulin is injected into the body, it cannot be removed, it will act naturally lowering BG levels.

The Spanish Consortium on Artificial Pancreas and Diabetes Technology has been working over the last decade on the development of a new artificial pancreas system. In the first clinical trial, to evaluate the performance of the PD controller with the SAFE layer, an individualized constraint was designed to control postprandial BG levels after the consumption of a meal with 60 grams of CHO [26]. The closed-loop controller achieved better outcomes compared with the open-loop therapy, reducing significantly the time spent in hyperglycemia without increase the risk of hypoglycemia. However, a limitation of this study is that was tuned for this specific meal size. Therefore, novel approaches on how to select for meals of different sizes are required to cope with daily-life operation of an AP system. In this work, a novel approach to tune is presented. This new approach takes into account previous open-loop therapy to set for periods without meals and in case of the announcement of meals, can be raised to reduce hyperglycemia without leading to hypoglycemia in the late postprandial period. The proposal is evaluated in silico using the UVa/Padova simulator [27].

2. Materials and Methods

2.1. Control-Loop Scheme

In this section, the control scheme considered in this work is introduced. The control algorithm consists of two loops, as depicted by Figure 1. The inner loop is composed by a proportional derivative (PD) controller with an insulin feedback (IFB) loop. The outer loop contains the SAFE layer [20], which is inspired on the sliding mode reference conditioning technique [21]. This SAFE layer applies a discontinuous signal generated by a switching law when estimations of IOB () surpass a preset limit of IOB (). Later, the discontinuous signal is filtered generating a smooth reference signal () to be applied into the controller. Therefore, this safety layer is able to maintain IOB inside desired bounds.

Figure 1: Control scheme based on a PD controller with IFB and with the SAFE layer. PK, pharmacokinetic.

The control action produced by the PD controller is presented as follows:

where is the proportional gain, is the derivative time, is the glucose reference filtered after the action of the SAFE layer, and is the interstitial BG measurements provided by the CGM.

Then, the control action computed by the PD is augmented by two feed-forward signals: and . The signal is the insulin obtained from patients’ daily basal profile. The term is an impulse signal in case of the announcement of meals to compensate the disturbance caused by the ingestion of CHO.

The IFB algorithm [28] emulates the beta-cell physiology in healthy subjects, suppressing insulin secretion as plasma insulin concentration increases. The combination of the PD algorithm with both SAFE layer and IFB algorithm has already been investigated and shown to be more effective than when they are used separately [29]. The final control action signal provided to the insulin pump is

where γ is gain parameter and is the is the estimated deviation of plasma insulin from steady state conditions (basal levels).

The term in equation (1) is the conditioned reference due to the action of the SAFE layer to maintain IOB below . Since the IOB is inaccessible, an insulin absorption model is considered to estimate IOB [30], through the following equation:

where and are two compartments, is a constant related with the duration of insulin action (DIA), and is the estimation of IOB. The SAFE layer has a software-based nature and consists of two main elements: a switching block responsible to generate a discontinuous signal to maintain into the desired range and a first-order filter to smooth the discontinuous signal before being applied to the main controller.

Consider the sliding function defined by equation (4), the switching logic is defined as follows:

with W > 0 mg/dl.

Finally, the discontinuous signal is filtered by the following equation and generates smooth changes in the glucose reference signal:

Note that when , is greater than . In order to drive to the desired range, i.e., below , must be decreased. The addition of W in equation (6) generates greater than , diminishing the insulin suggestion provided by the main controller and thus reducing . When is below , no further action is provided by the outer loop, letting the controller work freely.

2.2. IOB Constraint Tuning

The selection of the constraint is a critical point in the design of the control system. In this work, a new tuning approach for is presented. In patients with T1D, insulin requirements vary during the day and also between days (intra- and interday variability). On traditional insulin pump therapy, physicians configure patients’ pump to deliver a steady flow of basal insulin to cope with intraday variability. In addition, insulin boluses are delivered when meals are informed by patients, to cover the disturbance caused on glycemic balance due to the ingestion of CHO [3].

The artificial pancreas under development by our research group has been evaluated clinically, where 20 T1D subjects (age 40.7 ± 10.4 years, T1D duration 22.2 ± 9.9 years, and A1c 7.8 ± 0.7%) used the CL system in front of a mixed meal containing 60 grams of CHO [26]. In this trial, an individualized tune was considered in order the improve postprandial glycemic control when compared with standard open-loop therapy. For this specific trial, where a single meal with 60 grams of CHO was consumed, was adjusted based on parameters taken from patients’ open-loop therapy and computed in an offline procedure. The procedure to compute was as follows: considering that patients were in basal levels of IOB, was computed as the estimation of IOB levels 90 minutes after the administration of an augmented bolus, by equation (3), to compensate this 60 grams meal. This augmented bolus was computed by adding to the standard bolus the amount of basal insulin that would have been delivered in the next hour in the case of being in open-loop therapy. Therefore, knowing in advance the size of the meal and all the parameters necessary to compute this augmented bolus, it was possible to compute in an offline fashion what would be the estimation of IOB levels 90 minutes after the meal bolus. This single value was applied as IOB limit during in this clinical trial [26].

Due to the huge amount of insulin in meal boluses, violates , and a high frequency discontinuous signal is generated by the SAFE layer in order to return back to bounds. This action forces insulin delivery to zero for approximately 90 minutes, minimizing the effects of controller’s overcorrection. When (equation (4)), insulin infusion may be restored if the controller deems necessary.

However, a single value of may not be sufficient for 24-hour operation, especially due to the large intrapatient variability in T1D and to the different activities performed by subjects in their daily-life. The major problem observed in the strategy presented previously is that too high values of may cause that the SAFE layer be ineffective because it will act in very few conditions, e.g., only after a meal bolus. It is comprehensible having higher tuning during postprandial periods, especially because subjects tend to underestimate the CHO content in meals [31]. However, during late postprandial period, where the effects of meals have been covered either by the bolus or by the controller suggestions, such high values of may lead to excess of insulin in the body, increasing the risk of hypoglycemia.

In this work, a new method to tune to overcome the limitations presented by the former strategy is proposed. The new approach is called “dynamic rule-based” (DRB) algorithm and is intended to be used on hybrid artificial pancreas systems for 24-hour operation, where meals of any size are consumed. The proposed approach combines two different strategies to set for different situations: (1) fasting periods: where no big disturbance is expected, and the controller must deal mainly with intraday variability and (2) postprandial periods: where a substantial raise in BG levels is expected due to the consumption of a meal and insulin bolus may not be enough to compensate such disturbance. In summary, the DRB algorithm generates a time-varying based on patients’ basal insulin profile taken from open-loop therapy, and when a meal is announced, the algorithm evaluated in real-time if should be increased temporarily in order to reduce hyperglycemia. The following sections depict the algorithm, and a flowchart is also presented in Figure 2 to clarify the operation of the DRB algorithm.

Figure 2: Flowchart describing how the proposed dynamic rule-based algorithm works.

2.2.1. Fasting Periods

During fasting period (sleeping hours or daytime fasting periods), the absence of external disturbances makes the intraday variability the major challenge for the controller. The basal insulin profile from patients’ open-loop therapy is considered to create a baseline for :

where is the baseline for and is a gain that regulates the amplitude of , with . In case of , will not be allowed to be greater than it would have been during open-loop therapy. On the contrary, with , the insulin controller can suggest more insulin than what is programmed by the open-loop therapy. Thus, should be selected in order to protect patients from hypoglycemia but also allowing the control algorithm to suggest insulin above open-loop regimen when necessary. Thus, during fasting periods, is set to be equal to . In this work, the parameter is set to 1.3 for daytime period (06:00–23:00) and to 1.1 for night-time period (23:00–06:00).

Although this approach may be able to provide good glycemic control in front of intraday variability, patients in free-living conditions take CHO on several occasions during the day, requiring a greater amount of insulin for these periods. Therefore, the strategy of during fasting periods is combined with another one, intended for postprandial periods, presented in the following section.

2.2.2. Postprandial Periods

The tuning of during postprandial period is an extension to the method already validated clinically by Rossetti and colleagues [26]. Here, the method is generalized for meals with different amount of CHO and also includes a set of rules based on BG readings to modify . This set of rules has been designed to determine if it is necessary to increase , for cases which the bolus was not enough to drive BG to near-normoglycemia levels.

The method works as follows: in case of the announcement of a meal, an insulin bolus is delivered as a feed-forward action. This bolus is an augmented version of the standard bolus computation, by adding a portion of the future basal delivery according to the size of the meal, as shown in the following equation:

where is the content of CHO of the meal (in grams), is the insulin-to-CHO ratio, is the BG reference, and is the correction factor.

After the bolus, surpasses , and then is forced to zero, while is greater than , due to the action of the SAFE layer. The parameter (in minutes, equation (9)) is introduced to regulate the starting time after the meal from when may be increased. After minutes, BG readings start to be evaluated in order to check whether needs to be increased, aiming to drive BG below a selected target (). If in this moment, BG is greater than a threshold (), a new IOB limit is computed to control postprandial BG, based on equation (10). This value is maintained as , while BG is greater than . Finally, when BG returns to values below , returns to follow . Note that the parameters and can be adjusted intuitively by physicians. For a more aggressive postprandial control, these parameters should be decreased as follows:

The final tuning for is determined during a real-time procedure, based on the dynamic behavior of patients basal insulin profile and also in the set of rules activated after the announcement of meals, to increase during postprandial periods. Figure 2 depicts a flowchart of the proposed method to facilitate the understanding.

The variables “PP_state” and “flag_PP” included in Figure 2 were considered for implementation purposes. In the initialization of the system, both variables should be set to zero. “PP_state” indicates that a meal has been consumed, and that may be increased considering a set of rules to evaluated BG levels minutes after the announcement of the meal. The variable “flag_PP” guarantees that a single value of is computed for each meal. Figure 3 depicts the application of the method in one representative virtual patient during a simulation, with a meal containing 45 g of CHO at 07:30 (represented by the green triangle in Figure 3(a)). Note that 67.5 minutes after the meal (computed by equation (9)), the limit of IOB has been increased because BG levels were above . When was increased, the controller suggested more insulin in order to reduce the postprandial excursions since the insulin bolus was not enough for this specific meal. This period of time represents the first meal of a single patient during Scenario C, which is detailed in Section 2.3.3.

(a)

(b)

(c)

(d)

Figure 3: Operation of the dynamic rule-based algorithm during a postprandial period. Note that, after the meal (consumed at 07:30 and represented by the green triangle), the algorithm deemed necessary to increase because BG levels were above . (a) BG measurements in mg/dl. (b) Insulin delivery. (c) The conditioned reference signal . (d) Estimation of IOB levels and the constraint .

In summary, the DRB algorithm has been designed to make use of patients’ basal profile, which is an indicative of insulin requirements along the day, combined with a modified version of the approach already tested clinically with real patients, which achieved good results during postprandial control. In addition, BG measurements has been incorporated to set , aiming to track glucose back to regular values safely. All the relevant parameters used in the simulations are listed in Table 1.

Table 1: Parameters considered in this work for both inner and outer loops.

Parameter
Value
Unit

90

100

Γ
0.42

W
350

Τ
10

Λ
0.1

0.013

150

140

TDI, total daily insulin.

2.3. In Silico Evaluation

The proposed method is validated in the UVa/Padova simulator [27] on three challenging scenarios, including intrapatient variability in insulin sensitivity and in meal absorption rate [32].

Circadian variability has been included to simulate different requirements of insulin during the day and follows a sinusoidal variation. The parameters and , which are related with the insulin sensitivity are modified as follows:

where is the corresponding time-varying parameter; is the default individual parameter value ( or ), and is a uniformly distributed random number between 0 and 1.

Additionally, meal absorption rate and insulin absorption parameters (parameters , , , and ) assume different values (±30% around the standard value) after every single meal consumption. Further details about the aforementioned parameters are described elsewhere [27, 33].

2.3.1. Scenario A

This scenario is considered to compare the performance of the DRB algorithm () against the method already validated clinically, with a fixed value for (). In a 7-day scenario, the adult cohort consumed a single meal containing 60 grams of CHO per day, between 10:00 and 16:00, in order to assess the postprandial control under different intrapatient variability conditions. This amount of CHO has been selected to conduct a fairly comparison since the former strategy was clinically validated for this specific meal size.

2.3.2. Scenario B

This scenario is considered to compare the performance of the DRB algorithm () in front of meals with CHO content varying between 40 and 120 grams of CHO. In a 45-day scenario, the adult cohort consumed a single meal per day, between 08:00 and 19:00, in order to assess the postprandial control under different intrapatient variability conditions. It is considered just a single meal per day in order to avoid the accumulated effects of meals in the results.

2.3.3. Scenario C

A 14-day scenario is considered to assess the performance of the proposed method intended to mimic real-life operation of the AP systems. A total of three meals, with different amounts of CHO, are consumed per day at 7:30 (45 grams), 13:00 (90 grams), and 18:30 (50 grams). An error of on CHO counting has also been included to challenge the system. To apply meal absorption variability, the 10 mixed meal models from each patient were randomly assigned for each meal intake along with the simulation. A total of four different strategies are applied in this scenario: (1) the DRB algorithm to adjust (), (2) a fixed value of as used in Scenario A (), (3) is set to be equal to , without the rules considered for postprandial period (), and (4) the same insulin controller considered in other strategies, but without .

3. Results and Discussion

In this section, glycemic outcomes are presented for all the scenarios previously described. The performance of different methods to adjust is evaluated according the time spent into different glycemic ranges [34]. Results are computed based on CGM measurements. Individual metrics for the adult cohort from the simulator are computed, and then the results are presented as the median (25th–75th percentile) among the cohort. Additionally, the occurrence of hypoglycemic episodes (defined as at least 15 consecutive minutes with glucose below 70 mg/dl) is also analyzed.

3.1. Scenario A

Table 2 shows the metrics for Scenario A. These metrics assess the performance of and during postprandial period (i.e., 4 hours following the meal), once the last has been designed for such purpose. The metrics used for this evaluation are: average BG levels, percentage of time spent in different glycemic ranges, glycemic excursion (defined as the difference between the maximum BG level and the premeal BG level), and the total insulin delivered by the controller after the meal bolus.

Table 2: Population metrics for postprandial glycemic control in scenario A.

Mean glucose
Percentage of time spent in
Excursion (mg/dl)
Total basal (U)

(mg/dl)
70–140
70–180
>180
<70

141.48
46.58
96.13
3.87
0.00
62.00
3.56

(136.4–153.0)
(24.7–56.5)
(93.8–100.0)
(0.0–6.3)
(0.0–0.0)
(55.9–67.3)
(3.0–4.0)

141.40
47.77
96.43
3.57
0.00
63.36
3.78

(136.7–150.8)
(25.6–53.6)
(94.0–100.0)
(0.0–6.0)
(0.0–0.0)
(53.6–68.7)
(3.1–4.1)

Note that both strategies achieved similar results. Although obtained numerical results slightly superior than the proposed method, no significant difference was observed in any of the metrics. The new strategy (), besides being applicable for meals of any size, achieved equivalent outcomes when compared with the former strategy, which was designed for meals with 60 grams of CHO and achieved good results in real patients.

However, analyzing the whole simulation period (7-days), the new strategy was able to eliminate the occurrence of hypoglycemic events, while five events were observed for . This fact reinforces the hypothesis that modifying over the day increases patients safety.

3.2. Scenario B

Table 3 shows the metrics during postprandial period in Scenario B for the two methods applied in this Scenario B.

Table 3: Population metrics for postprandial glycemic control in scenario B.

Mean glucose
Percentage of time spent in
Excursion (mg/dl)
Total basal (U)

(mg/dl)
70–140
70–180
>180
<70

146.79
41.52
85.86
14.14
0.00
71.79
3.19

(144.7–154.7)
(35.4–48.5)
(77.7–90.2)
(9.8–20.8)
(0.0–0.0)
(66.1–74.7)
(2.1–3.6)

148.96
38.10
84.38
15.63
0.00
75.43
3.43

(147.0–156.2)
(25.3–47.3)
(77.1–86.3)
(13.7–22.9)
(0.0–0.0)
(69.0–79.4)
(2.6–4.1)

Different from Scenario A, in Scenario B, meals from different sizes were considered to evaluate the postprandial performance of the controller, which is more realistic. Results showed a slight superiority of the against the in all the glycemic metrics analyzed. Notice that the percentage of time spent in tight glycemic range increased from 38.1% to 41.52%, while the percentage of time spent in hyperglycemia reduces at the same time, from 15.63% to 14.14%. This analysis shows that the previous tuning of may not be sufficient for varied size of meals.

3.3. Scenario C

Table 4 shows the results for the 14-day scenario comparing the system with the DRB algorithm () against the other strategies. It displays the mean glucose, percentage of time spent in different glycemic ranges and the number of hypoglycemic events for the entire cohort.

Table 4: Population metrics comparing the performance of the system with the dynamic rule-based algorithm () against other strategies, in scenario C.

Without

Daytime (06:00–23:00)

    Mean glucose
140.30 (134.2–143.5)
135.57 (133.0–140.3)
142.72 (134.7–145.3)
116.92 (114.0–121.6)

    % of time spent in

    70–140
56.58 (49.3–66.3)
60.17 (56.1–65.4)
53.43 (46.3–64.6)
71.53 (67.4–79.4)

    70–180
92.10 (85.5–95.3)
92.70 (85.9–96.6)
90.91 (85.2–95.0)
90.88 (83.5–93.7)

    >180
7.90 (4.7–14.5)
7.30 (3.4–14.1)
9.09 (5.0–14.8)
5.25 (3.0–8.4)

    <70
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)
5.27 (0.4–9.9)

    <54
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)
1.61 (0.0–6.3)

    # Hypoglycemic events
3
4
4
112

Night-time (06:00–23:00)

    Mean glucose
114.61 (112.1–118.8)
106.51 (104.0–107.2)
114.81 (112.9–126.1)
102.37 (100.1–103.9)

    % of time spent in

    70–140
98.09 (94.3–99.1)
99.36 (96.0–100.0)
97.75 (87.2–99.2)
98.04 (96.0–99.3)

    70–180
100.00 (100.0–100.0)
100.00 (100.0–100.0)
100.00 (100.0–100.0)
98.13 (97.1–99.4)

    >180
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)

    <70
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)
1.87 (0.6–2.9)

    <54
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.00 (0.0–0.0)
0.13 (0.0–0.9)

    # Hypoglycemic events
0
4
0
11

The results display the solid performance of the strategies using the SAFE layer when compared with the insulin controller without IOB limitation, mainly to avoid hypoglycemia. The arm with the strategy achieved the lowest amount of hypoglycemic events when compared with the other methods. The strategy is the one with the lowest mean glucose values during daytime, but the fixed value of is not enough to avoid hypoglycemic events, especially during night-time. On the contrary, the arm considered only the adjustment for presented in Section 2.2.1. Such approach did not consider the rule-based algorithm to increase during postprandial periods, and this is reflected in the slightly worse outcomes for this strategy when compared with , during daytime. The arm achieved lower mean glucose values and spent less time in the hyperglycemic range, when compared with the arm.

Also during daytime, it can be observed that achieved higher mean glucose values when compared with the insulin controller without , considerably above the . Although these values are higher than in healthy people, they are still acceptable considering the recommendations of the American Diabetes Association (ADA) [35], which is A1C <7% for nonpregnant adults. The median GMI (glucose management indicator) achieved by the proposed method during daytime is 6.7%, which gives the approximate A1C levels based on CGM measurements [36]. During night-time, was able to lower mean glucose, without any hypoglycemic event. Considering 24-hour period for the entire scenario, the proposed method achieved a median GMI of 6.5%, within the limits recommendable by the ADA.

The control-variability grid-analysis (CVGA) [37], in Figure 4, allows a graphical visualization regarding the glycemic variability within an observational period of 24 hours. Note that the approach achieved the highest percentage of points falling either in Zone A or Zone B, with 97.85%. Additionally, Figure 5 shows the dynamics of BG and insulin delivery for a single day, comparing with the arm without IOB limits, to better illustrate the influence of the in the control action.

(a)

(b)

(c)

(d)

Figure 4: Control-variability grid-analysis for different arms evaluated in scenario C: (a) , (b) , (c) , and (d) without .

(a)

(b)

Figure 5: Comparison of the dynamics of BG and insulin delivery between the controller with the DRB algorithm and without for a single patient in scenario C: (a) BG readings; (b) insulin delivery.

Observing the results in Table 4 and Figure 4 it is possible to note that the action of the SAFE layer avoids the overreaction of the insulin controller due to the rise of glucose, caused by the meals. Although the results obtained by the DRB algorithm are only slightly better when compared with the other strategies which included the SAFE layer, it was observed a reduction on the occurrence of hypoglycemic events without leading to excessive hyperglycemia. All the three hypoglycemic events observed in the arm were caused by the meal bolus and not by any insulin suggested by the controller after the meal. Therefore, the hypoglycemic episodes were very likely caused by overestimation of CHO content to compute the bolus.

The starting point of the proposed method was another strategy which has already been extensively tested both in silico and clinically, making the task of achieving better results even more difficult. Nevertheless, it has been possible to improve slightly the performance of the AP system in this in silico study, being able to apply different limits of IOB according to the CHO content of meals and by applying a lower during night-time. Any improvement in postprandial glycemic control can reflect in a reduction on the risks associated with long-term complications, given its correlation with A1C levels [38], and the avoidance of nocturnal hypoglycemia is a major concern in T1D treatment. In addition, it also allows modification on the parameters of the proposed algorithm to be performed intuitively by physicians, if they deem necessary to further improvements on glycemic control. However, clinical trials involving real patients must be conducted to assess the performance of the proposed algorithm under real-life operation, in which patients may forget to announce meals, and the performance of the DRB algorithm may not be suitable to avoid hyperglycemia. Furthermore, the DRB algorithm may not be applicable in fully automatic AP systems, in which patients do not need to announce meals. Finally, the adjustments of is a major task for AP systems which considers the SAFE layer, and the proper adjustment of this constraint also plays an important role even during physical activity [39, 40], by reducing the amount of injected insulin during and after exercise.

4. Conclusions

The dynamic rule-based algorithm proposed in this work has been designed to tune in real-time limiting of IOB to safely control BG levels. This algorithm is intended for 24-hour operation, which includes postprandial and fasting periods. During postprandial periods, it allows the increase of IOB limits when more insulin is required to have a more aggressive controller but yet safe. The strategy has been evaluated in silico under challenging conditions and achieved satisfactory performance, with emphasis on the reduction of hypoglycemic events during nocturnal period and without excessive hyperglycemia during postprandial period for meals with different CHO contents. Although, in this paper, the proposed strategy has been applied and evaluated in a PD controller, such approach could also be used by other algorithms, since the SAFE layer can be added to control algorithms of any nature.

Abbreviations

ADA: | American Diabetes Association
AP: | Artificial pancreas
BG: | Blood glucose
CGM: | Continuous glucose monitor
CHO: | Carbohydrates
DIA: | Duration of insulin action
DRB: | Dynamic rule-based
IFB: | Insulin feedback
IOB: | Insulin-on-board
PD: | Proportional derivative
PK: | Pharmacokinetic
SAFE: | Safe auxiliary feedback element
T1D: | Type 1 diabetes.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work has been partially funded by the Spanish Government (grant no. DPI2016-78831-C2-2-R), National Counsel of Technological and Scientific Development, CNPq-Brazil (grant nos. 207688/2014-1 and 202050/2015-7), and Formación de Professorado Universitario (grant no. FPU0244 2015).

References

M. Lind, A.-M. Svensson, M. Kosiborod et al., “Glycemic control and excess mortality in type 1 diabetes,” New England Journal of Medicine, vol. 371, no. 21, pp. 1972–1982, 2014.

F. Zaccardi, D. R. Webb, T. Yates, and M. J. Davies, “Pathophysiology of type 1 and type 2 diabetes mellitus: a 90-year perspective,” Postgraduate Medical Journal, vol. 92, no. 1084, pp. 63–69, 2016.

J. Walsh, R. Roberts, C. Varma, and T. Bailey, Using Insulin: Everything You Need for Success with Insulin, Torrey Pines Press, San Diego, CA, USA, 2003.

D. M. Nathan and DCCT/EDIC Research Group, “The diabetes control and complications trial/epidemiology of diabetes interventions and complications study at 30 years: Overview,” Diabetes Care, vol. 37, no. 1, pp. 9–16, 2014.

A. J. Graveling and B. M. Frier, “The risks of nocturnal hypoglycaemia in insulin-treated diabetes,” Diabetes Research and Clinical Practice, vol. 133, pp. 30–39, 2017.

E. Bekiari, K. Kitsios, H. Thabit et al., “Artificial pancreas treatment for outpatients with type 1 diabetes: systematic review and meta-analysis,” BMJ, vol. 361, p. k1310, 2018.

H. Thabit and R. Hovorka, “Coming of age: the artificial pancreas for type 1 diabetes,” Diabetologia, vol. 59, no. 9, pp. 1795–1805, 2016.

F. J. Doyle, L. M. Huyett, J. B. Lee, H. C. Zisser, and E. Dassau, “Closed-loop artificial pancreas systems: engineering the algorithms,” Diabetes Care, vol. 37, no. 5, pp. 1191–1197, 2014.

V. Gingras, N. Taleb, A. Roy-Fleming, L. Legault, and R. Rabasa-Lhoret, “The challenges of achieving postprandial glucose control using closed-loop systems in patients with type 1 diabetes,” Diabetes, Obesity and Metabolism, vol. 20, no. 2, pp. 245–256, 2018.

G. C. Goodwin, A. M. Medioli, D. S. Carrasco, B. R. King, and Y. Fu, “A fundamental control limitation for linear positive systems with application to type 1 diabetes treatment,” Automatica, vol. 55, pp. 73–77, 2015.

Y. M. Luijf, A. C. van Bon, J. B. Hoekstra, and J. H. DeVries, “Premeal injection of rapid-acting insulin reduces postprandial glycemic excursions in type 1 diabetes,” Diabetes Care, vol. 33, no. 10, pp. 2152–2155, 2010.

P. Rossetti, F. J. Ampudia-Blasco, A. Laguna et al., “Evaluation of a novel continuous glucose monitoring-based method for mealtime insulin dosing-the iBolus-in subjects with type 1 diabetes using continuous subcutaneous insulin infusion therapy: a randomized controlled trial,” Diabetes Technology & Therapeutics, vol. 14, no. 11, pp. 1043–1052, 2012.

M. Reddy, P. Pesl, M. Xenou et al., “Clinical safety and feasibility of the advanced bolus calculator for type 1 diabetes based on case-based reasoning: a 6-week nonrandomized single-arm pilot study,” Diabetes Technology & Therapeutics, vol. 18, no. 8, pp. 487–493, 2016.

K. Turksoy, I. Hajizadeh, S. Samadi et al., “Real-time insulin bolusing for unannounced meals with artificial pancreas,” Control Engineering Practice, vol. 59, pp. 159–164, 2017.

M. M. Al-Tabakha, “Future prospect of insulin inhalation for diabetic patients: the case of afrezza versus exubera,” Journal of Controlled Release, vol. 215, pp. 25–38, 2015.

C. Ellingsen, E. Dassau, H. Zisser et al., “Safety constraints in an artificial pancreatic β cell: an implementation of model predictive control with insulin on board,” Journal of Diabetes Science and Technology, vol. 3, no. 3, pp. 536–544, 2009.

I. Hajizadeh, M. Rashid, S. Samadi et al., “Adaptive and personalized plasma insulin concentration estimation for artificial pancreas systems,” Journal of Diabetes Science and Technology, vol. 12, no. 3, pp. 639–649, 2018.

C. Toffanin, H. Zisser, F. J. Doyle, and E. Dassau, “Dynamic insulin on board: incorporation of circadian insulin sensitivity variation,” Journal of Diabetes Science and Technology, vol. 7, no. 4, pp. 928–940, 2013.

I. Hajizadeh, M. Rashid, K. Turksoy et al., “Plasma insulin estimation in people with type 1 diabetes mellitus,” Industrial & Engineering Chemistry Research, vol. 56, no. 35, pp. 9846–9857, 2017.

A. Revert, F. Garelli, J. Pico et al., “Safety auxiliary feedback element for the artificial pancreas in type 1 diabetes,” IEEE Transactions on Biomedical Engineering, vol. 60, no. 8, pp. 2113–2122, 2013.

F. Garelli, R. J. Mantz, and H. D. Battista, Advanced Control for Constrained Processes and Systems, The Institution of Engineering and Technology, London, UK, 2011.

F. León-Vargas, F. Garelli, H. De Battista, and J. Vehí, “Postprandial response improvement via safety layer in closed-loop blood glucose controllers,” Biomedical Signal Processing and Control, vol. 16, pp. 80–87, 2015.

R. Hu and C. Li, “An improved PID algorithm based on insulin-on-board estimate for blood glucose control with type 1 diabetes,” Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 281589, 8 pages, 2015.

R. Sánchez-Peña, P. Colmegna, F. Garelli et al., “Artificial pancreas: clinical study in Latin America without premeal insulin boluses,” Journal of Diabetes Science and Technology, vol. 12, no. 5, pp. 914–925, 2018.

F. León-Vargas, F. Garelli, H. De Battista, and J. Vehí, “Postprandial blood glucose control using a hybrid adaptive PD controller with insulin-on-board limitation,” Biomedical Signal Processing and Control, vol. 8, no. 6, pp. 724–732, 2013.

P. Rossetti, C. Quirós, V. Moscardó et al., “Closed-loop control of postprandial glycemia using an insulin-on-board limitation through continuous action on glucose target,” Diabetes Technology & Therapeutics, vol. 19, no. 6, pp. 355–362, 2017.

C. D. Man, F. Micheletto, D. Lv, M. Breton, B. Kovatchev, and C. Cobelli, “The UVA/PADOVA type 1 diabetes simulator,” Journal of Diabetes Science and Technology, vol. 8, no. 1, pp. 26–34, 2014.

G. M. Steil, C. C. Palerm, N. Kurtz et al., “The effect of insulin feedback on closed loop glucose control,” The Journal of Clinical Endocrinology & Metabolism, vol. 96, no. 5, pp. 1402–1408, 2011.

I. Sala-Mira, J. L. Díez, and J. Bondia, “Insulin limitation in the artificial pancreas by sliding mode reference conditioning and insulin feedback: an in silico comparison,” IFAC-PapersOnLine, vol. 50, no. 1, pp. 7743–7748, 2017.

M. E. Wilinska, L. J. Chassin, H. C. Schaller, L. Schaupp, T. R. Pieber, and R. Hovorka, “Insulin kinetics in type-1 diabetes: continuous and bolus delivery of rapid acting insulin,” IEEE Transactions on Biomedical Engineering, vol. 52, no. 1, pp. 3–12, 2005.

A. S. Brazeau, H. Mircescu, K. Desjardins et al., “Carbohydrate counting accuracy and blood glucose variability in adults with type 1 diabetes,” Diabetes Research and Clinical Practice, vol. 99, no. 1, pp. 19–23, 2013.

P. Herrero, J. Bondia, O. Adewuyi et al., “Enhancing automatic closed-loop glucose control in type 1 diabetes with an adaptive meal bolus calculator-in silico evaluation under intra-day variability,” Computer Methods and Programs in Biomedicine, vol. 146, pp. 125–131, 2017.

C. Dalla Man, R. A. Rizza, and C. Cobelli, “Meal simulation model of the glucose-insulin system,” IEEE Transactions on Biomedical Engineering, vol. 54, no. 10, pp. 1740–1749, 2007.

G. Agiostratidou, H. Anhalt, D. Ball et al., “Standardizing clinically meaningful outcome measures beyond HbA1cfor type 1 diabetes: a consensus report of the american association of clinical endocrinologists, the american association of diabetes educators, the american diabetes association, the endocrine society, JDRF international, the leona m. and harry b. helmsley charitable trust, the pediatric endocrine society, and the t1d exchange,” Diabetes Care, vol. 40, no. 12, pp. 1622–1630, 2017.

American Diabetes Association, “Glycemic targets: standards of medical care in diabetes—2018,” Diabetes Care, vol. 41, pp. S55–S64, 2018.

R. M. Bergenstal, R. W. Beck, K. L. Close et al., “Glucose management indicator (GMI): a new term for estimating A1C from continuous glucose monitoring,” Diabetes Care, vol. 41, no. 11, pp. 2275–2280, 2018.

L. Magni, D. M. Raimondo, C. D. Man et al., “Evaluating the efficacy of closed-loop glucose regulation via control-variability grid analysis,” Journal of Diabetes Science and Technology, vol. 2, no. 4, pp. 630–635, 2008.

E. B. Ketema and K. T. Kibret, “Correlation of fasting and postprandial plasma glucose with HbA1c in assessing glycemic control; systematic review and meta-analysis,” Archives of Public Health, vol. 73, no. 1, 2015.

A. Bertachi, A. Beneyto, C. M. Ramkissoon, and J. Vehí, “Assessment of mitigation methods to reduce the risk of hypoglycemia for announced exercise in a uni-hormonal artificial pancreas,” Diabetes Technology & Therapeutics, vol. 20, no. 4, pp. 285–295, 2018.

C. M. Ramkissoon, A. Bertachi, A. Beneyto, J. Bondia, and J. Vehi, “Detection and control of unannounced exercise in the artificial pancreas without additional physiological signals,” IEEE Journal of Biomedical and Health Informatics, p. 1, 2019, https://ieeexplore.ieee.org/document.