In silico evaluation of a control system and algorithm for automated insulin infusion in the ICU setting
BioMedical Engineering OnLine volume 9, Article number: 35 (2010)
It is known that tight control of glucose in the Intensive Care Unit reduces morbidity and mortality not only in diabetic patients but also in those non-diabetics who become transiently hyperglycemic. Taking advantage of a recently marketed subcutaneous glucose sensor we designed an Automatic Insulin Infusion System (AIIS) for inpatient treatment, and tested its stability under simulated clinical conditions.
The system included: reference glucose, glucose sensor, insulin and glucose infusion controllers and emergency infusion logic. We carried out computer simulations using Matlab/Simulink®, in both common and worst-case conditions.
The system was capable of controlling glucose levels without entering in a phase of catastrophic instability, even under severe simulated challenges. Care was taken to include in all simulations the 5-10 minute delay of the subcutaneous glucose signal when compared to the real-time serum glucose signal, a well-known characteristic of all subcutaneous glucose sensors.
When tested in-Silico, a commercially available subcutaneous glucose sensor allowed the stable functioning of a proportional-derivative Automatic Insulin Infusion System, which was able to maintain glucose within acceptable limits when using a well-established glucose response model simulating a patient. Testing of the system in vivo using animal models is now warranted.
It has been shown that "tight control of blood glucose" (i.e. between 80 and 110 mg × dL-1 = 4.4-6.1 mmol × L-1) in critically ill patients reduces dramatically the mortality from 8.0% to 4.6% in Intensive Care Unit (ICU) patients. Additionally, it reduces bloodstream infections by 46%, renal failure by 41%, transfusions by 50% and polyneuropathy by 44%. Interestingly, these results apply to those individuals who -being diabetic or not- have blood glucose on admission > 110 mg × dL-1 [6.1 mmol·L-1] (up to 76% of ICU patients) [3–6]. Tight control of blood glucose is also useful in Coronary Care Units, where 33% of patients are diabetics and 33% have glucose intolerance. As a token of the importance of the subject, in Portland, Oregon, USA, tight control of blood glucose with intravenous insulin infusion reduced wound infection by 50%, hospital stay by 56% and mortality by 67% in diabetic patients who underwent open heart surgery .
The current standard of care for hospitalized diabetic patients (Figure 1) requires the measurement of capillary blood glucose at least 4 times a day, and by the use of a simple algorithm ("slide scale"), adjustment of the insulin infusion rate, usually by means of an infusion pump. In cases in which "tight control" is needed, capillary blood sampling has to be done every hour, increasing personnel workload. Hypoglycemia (i.e., serum glucose < 60 mg/dl) is a frequent side effect, and must be corrected by intravenous glucose.
In this context, we studied the behavior of an automated intravenous insulin infusion system (AIIS) designed by the authors, that uses an FDA-approved subcutaneous glucose sensor [10, 11]. This sensor has been shown to be useful for outpatient care by enabling individuals with type-1 diabetes to manually adjust their subcutaneous insulin doses according to the glucose levels displayed by the device [11, 12].
Our aim was to use Matlab/Simulink® to develop a computer program for the AIIS, and to challenge this system repeatedly in a fashion that in some cases was similar to the worst clinical conditions possible to find in patients with diabetes in the ICU. What we considered a success was the AIIS being able to control blood glucose without becoming instable even in the worst conceivable circumstances.
Our concerns about potential instability of the AIIS originated in a characteristic of the subcutaneous glucose sensor, i.e. its 5 to 10-min. delay when compared to changes in blood glucose. We hypothesized that this delay, by altering the phase of feedback loop of the AIIS, could result in catastrophic instability. However, the results of the simulations showed that the system behaves with a high level of stability under a wide range of clinical conditions likely to be encountered in the ICU setting.
I. Architecture of the AIIS
Our closed-loop control system is shown in Figures 2 and 3. The glycemia of the patient acquired by the glucose sensor [13, 14] is compared with a desired (reference) blood glucose level and the difference between these two values (the error) is processed by the controller, generating the signals for the actuators (insulin and glucose [dextrose in water] pumps).
I-A. Reference Glycemia
Normal fasting glycemia in adults varies between 70 and 99 mg × dL-1. In the case of an hospitalized patient, clinicians tend to keep the glycemia ( G ) below 120 mg × dL-1, but away from hypoglycemia. As symptoms of hypoglycemia start at G < 60 mg × dL-1, we set the reference level at 100 mg × dL-1, with an acceptable range of 90 to 120 mg × dL-1, near the "tight control" of blood glucose level (80-110 mg × dL-1) and safely away from hypoglycemia.
I-B Glucose Sensor
We choose the Medtronic Minimed (California, USA) Continuous Glucose Monitoring System (CGMS), a subcutaneous needle sensor capable of providing readings representative of blood glucose concentrations every 10 seconds, that are averaged every 5 minutes, totalizing 288 glucose values per day, to be either stored or displayed [11, 13, 14]. Its operation is based on the oxidation of glucose in the presence of oxygen and the enzyme glucose oxidase, generating an electrical current that is proportional to interstitial and blood glucose levels.
The sensor has been shown to be of value to improve outpatient control by deferred physician-supervised analysis. Its version with glucose reading display  allows patients to manually modify their own subcutaneous insulin dose (when using syringes) or infusion rate (when using a wearable insulin pump).
Based on the work of Chee et.al. (on the use of CGMS in critically ill patients) , and of Steil et.al. , we have developed a mathematical model of the CGMS sensor. For a sensor signal S(t) reacting to glucose level G(t) with sensitivity α (nano Ampères per mg/dL), the model equation becomes,
In equation 1, τ is the time constant that defines an immediate proportional response, and operationally represents the time elapsed for the signal from the sensor to reach 63% of the equilibrium when compared to the intravenous glucose concentration. TD is a transport lag . A white Gaussian noise to the sensor signal was added to our model, assuming a mean = 0.0 and a variance = 100 .
In their publication, Steil et. al.  included a graph of the simultaneous measurement of plasma and subcutaneous glucose, showing no transport lag (i.e., TD = 0.0) [16, 18]. Of the two glucose sensors presently being developed by Medtronic Minimed, the subcutaneous glucose sensor has been shown to have a time constant ranging from 2 to 12 min. The graph by Steil et.al. depicted a time constant (τ) of 5.77 minutes [16, 18], which is thought to be primarily related to interstitial glucose equilibration.
The insulin infusion pump controller is based on the algorithm proposed by Marliss et.al. in 1977, in the form later published in 1981 by Broekhuyse et.al. , whose original equations were as follows,
In equation 2a, IR is the insulin infusion rate effectively being pumped [mU × min-1], M is the maximum infusion rate possible [mU × min-1], S [dL × mg-1] is the slope of the hyperbolic tangent curve, and B [mg × dL-1] is the blood glucose level at which half maximum infusion rate is chosen to occur (as determined from the tanh function curve).
In equation 2b (see "Emergency Infusion Logic", below), IR' is the infusion rate in ml × hr-1 as determined by the insulin infusion controller, G 0 [mg × dL-1] is the last available average glycemia, G crit is the critical serum glucose level in mg × dL-1 below which intravenous glucose infusion must begin, (see also "glucose infusion controller" and equation 6, below), and IR min [mU × min-1] is the minimum rate of insulin infusion.
In equation 2a, DF is the Differential Factor [mg × dL-1], which is computed as,
where the parameters K 1 = 1 [mg × dL-1]-2 and K 2 = 10 [nondimensional] , and the original method for calculating A [mg × dL-1] is as follows,
where G i [mg × dL-1] is the average glycemia determined i minutes ago.
Since the CMGS System provides one measurement of glycemia every 5 minutes rather that at every minute, we propose the following modification of equation (4):
where G 0 [mg × dL-1] is the last reading of glycemia, and G 5 and G 10 are the previous readings [mg × dL-1] five and ten minutes earlier respectively. The coefficients accompanying the glycemia readings were obtained using as a reference the experimental data published by Broekhuyse et.al.  and the glycemia and insulinemia profiles from a 69.5 Kg male type 1 diabetic patient being controlled by the artificial pancreas of Albisser et.al. [22, 23] which used the algorithms corresponding to the aforementioned equations (2),(3) and (4). Needing blood samples for serum glucose every minute, Albisser's artificial pancreas was able to function for a limited period of time. However, because it allowed a fair glycemic control, we considered adequate to emulate its profile, but using the three glycemia readings obtained every 5 minutes as shown in Equation (5). With Albisser's data and the least square approach to find the optimum parameters, we reproduced the glycemic control as shown in Figure 4, but using equation (5) instead of equation (4).
From the point of view of the control theory, considering the current glycemia in addition to the differences with previous values in equations (2), (3) and (5), allowed us to build a proportional-derivative control model, as reviewed by Parker et. al. .
The Glucose Infusion Pump Controller: This subsystem was included in our AIIS as an emergency measure to avoid hypoglycemia. In normal circumstances, the insulin infusion alone should be enough to keep an adequate glycemic control. This emergency subsystem is activated when serum glucose falls below 60 mg × dL-1 (hypoglycemia). It has been determined that for the glucose infusion in these circumstances, proportional control (P) provides adequate performance . The controller equation is given by equation 6,
where GR is the glucose infusion rate in ml × h-1, G crit is the critical serum glucose level in mg × dL-1 below which the infusion must begin, and G 0 is the last 5-minute-averaged glycemia reading. K G is the controller gain. GR is always positive because the controller activates itself only when G 0 falls below G crit . The definition of GR' is given in the next paragraph, below.
Emergency Infusion Logic: It has been included in our system to both shut down the insulin infusion to a minimum and activate the glucose infusion whenever a critically low glycemia (G crit ) is detected. Once the glycemia returns to a value that is above the critical level, insulin infusion returns back to normal operation and glucose infusion is disabled. The activation or deactivation of each of the infusion pumps were summarized in equation 2b for insulin (see above) and in equation 6b for glucose, where GR' is the infusion rate in ml × hr-1 as determined by the glucose infusion controller, and GR is the infusion rate effectively leaving the glucose pump.
The insulin infusion has a minimum value that is different from zero (i.e. IR min ≠ 0; see Equation 2, above) for two reasons. First, considering that the insulin and glucose infusions take place through a catheter inserted in a vein, it is essential to ensure a minimum permanent flow to avoid the occlusion of the catheter (KVO, "keep vein open" feature). Second, in diabetic patients lacking endogenous production of insulin (Type-1 diabetes mellitus), it is essential to ensure a minimal serum concentration of insulin in order to avoid a massive lipolysis with ensuing ketoacidosis.
As described above, the AIIS has two infusion mechanisms regulating blood glucose concentration: insulin infusion and glucose infusion using the IMED Gemini PC-1 infusion pump produced by Alaris Medical Systems (San Diego, CA, U.S.A.). This pump fulfills the specifications of our system in terms of range of infusion rate, computer communication capabilities and accuracy. The ideal operation of the device for the simulations considered a minimum infusion rate of 1.0 ml × h-1, a maximum rate of 500 ml × h-1, and rate steps of 0.1 ml × h-1 . Although the IMED Gemini PC-2 model could have been more convenient due to its double channel capacity that allows the infusion of two solutions simultaneously and independently, for simplicity we used two separate IMED Gemini PC-1 pumps. (Figure 3).
I-E Diabetic Patient
In order to carry out the computer simulations aimed to test the behavior of our system, we used a mathematical model of human glucose physiology based on the work by Lehmann, Deutsch & Hermanyi [27, 28], which assumes a patient completely lacking endogenous insulin secretion. It contains a single compartment of extracellular glucose, where it enters via both intestinal absorption and hepatic glucose production. Similarly, the gastric emptying rate that provides the glucose flux from the stomach to the small intestine is controlled by a process defined as a function of the carbohydrate content of the ingested food. The glucose is removed from the extracellular space by insulin-independent utilization in the central nervous system and red blood cells, and also by the insulin-dependent utilization by the liver and peripheral tissues (muscle and fat). Additionally, glucose excretion is considered to take place through the kidneys when the glycemia exceeds renal glucose threshold.
The model accounts separately for the glucose entering the peripheral tissues and the glucose entering/leaving the liver, an approach that allows to assign a range of values to patient-specific parameters, such as hepatic insulin sensitivity [26, 27]. In this context, the model summarizes the glucose processing by the liver as the net hepatic glucose balance (NHGB), which is the algebraic sum of the glucose produced by both neoglucogenesis and glycogen breakdown, minus the glucose removed by glycogen synthesis. The NHGB is a function of both the glucose and the insulin plasma concentrations (see below) .
In fact, the Lehman & Deutsch model is based on four equations. The first one is a differential equation describing the change in plasma insulin concentration (dI/dt). The second is also a differential equation that now describes the build-up and deactivation of the active insulin pool (dIa/dt). The third and fourth equations describe the rate of insulin absorption from the subcutaneous tissue. As our AIIS injects insulin directly into the bloodstream, we did not use these last two equations. On the other hand, the overall rate of peripheral plus insulin-dependent glucose utilization (Gout , mg × Kg-1 × h-1) depends on an equation whose two main variables are the current glucose level (G, mg × dL-1) and the insulin concentration (I*eq , mU.I. × L-1). In the same equation, the most important constant is the slope of peripheral glucose utilization versus insulin line (c, 0.045 mg × min-1 × mIU-1 × L).
When the patient eats carbohydrates, the Lehman & Deutsch model allows for calculation of the rate of change of glucose in the gut (Ggut , mg × dL-1), which depends on the glucose load and the rate of gastric emptying. This Ggut in turn, allows the continuous calculation of the systemic appearance of glucose via absorption from the gut (Gin , mg × dL-1 × h-1).
The NHGB (mg × h-1), which can be either positive or negative, depends on effective plasma insulin, that is, plasma insulin (I*eq , mU.I. × L-1), times the patient-specific hepatic insulin sensitivity (Sh), an adimensional number that has a normalized value between 0 and 1. Published tables [27, 29] allow for the calculation of NHGB for any effective insulin level.
Finally, the change of glycemia with time (dG/dt, mg × dL-1 × h-1 ) is given by a differential equation depending on Gin(t), plus NHGB(t), minus Gout(t), and minus Gren(t). Where Gren is the rate of renal excretion of glucose (mg × h-1).
The main advantage of the Lehman & Deutsch model resides in its explicit anatomical compartments, taking into account the function of the body organs involved. For the same reason, it is easy to introduce modifications in the model, as it was necessary with our own simulations. The Lehman & Deutsch model was validated clinically in 1994 by their own authors .
More recently, others  have tried to improve the understanding of the glucose-insulin interactions by means of an approach based in 4 unit processes, i.e., liver, gastrointestinal tract, muscle & adipose tissue and beta cell. In our view however, this interesting new model underestimates the endogenous glucose production [31, 32].
II.- Simulations of the Behavior of the AIIS in Hospitalized Patients
To test the behavior of the AIIS in normal and in the worst possible scenarios, we designed the series of cases summarized in tables 1 and 2. Simulations were carried out in Matlab/Simulink® v6.5 from Mathworks Inc., using a Pentium IV 1,800 MHz computer, equipped with 256 MB RAM. Under these conditions, each 24-hour simulation took about 15 seconds to complete.
I Positive Challenges
Table 1 shows the positive challenges that we expected to be controlled by the AIIS, i.e. all perturbations tending to increase glycemia. We assumed that a single perturbation occurred at the hour 10:00 within the 24-hour simulation span.
I-A Mild Challenge
The mild challenge simulates a patient eating normally an average meal, containing 50 grams of complex carbohydrates (starch), the digestion of which allows the glucose to enter the bloodstream in a slow fashion. As shown in Table 1 and in Figure 4-A, in the case of a 70-Kg diabetic with average hepatic sensitivity (Sh = 0.5), the glycemia rises temporarily to a level close to 150 mg × dL-1, comparable to the transient elevation of serum glucose after a meal in a normal adult. As soon as the AIIS detects this rise in glycemia, the insulin infusion rate rapidly increases in concordance with the derivative component of the control algorithm (Figure 4-B). The same algorithm makes insulin infusion rate to decrease sharply as soon as the first derivative of the glycemia reaches the positive-to-negative critical point, without falling below 1 ml × h-1, thus avoiding the risk of ketoacidosis (a complication that might arise if plasma insulin concentration was allowed to fall below 0.5 μU × ml-1). In fact, as shown in Figure 4-C, plasma insulin never falls below the lower limit of the normal fasting range (6-26 μU × ml-1) . Furthermore, the aforementioned minimal infusion rate of 1 ml/hr corresponds to the parameter "K.V.O." of the IMED Gemini PC-1 infusion pump. The lack of an undershoot of glycemia after this mild challenge precluded the glucose infusion to start.
I-B Moderate Challenge
A moderate positive demand was given by the ingestion of 50 grams of glucose in the form of 250 grams of grapes. As seen in the 14th row of Table 1, in the case of a 70-Kg diabetic with average hepatic sensitivity (Sh = 0.5) the glycemia rose rapidly up to a value close to 170 mg × dL-1, and a peak duration of 30 min. The hyperglycemia was then controlled by the AIIS by means of a sharp and sustained increase in the insulin infusion rate, that corrected glycemia after a TS of 60 minutes. As with the mild challenge, only the insulin infusion was required whereas the glucose pump remained inactive.
I-C Severe Challenge
For the severe positive challenge we selected the most abrupt method of administering glucose to a patient: injecting 20 ml of 50% dextrose-in water into the subclavian vein in 30 seconds. As seen both in the 23d row of Table 1 and in Figure 5, in the case of a 70-Kg diabetic with average hepatic sensitivity (Sh = 0.5), a very steep rise in glycemia ensued after glucose injection (Figure 5-A). The fast rise in the insulin infusion rate that was ordered by the AIIS controller (Figure 5-B) was able to keep the hyperglycemic overshoot within tolerable limits, i.e., 180 mg × dL-1 for 15 minutes. For some minutes after the end of the sharp insulin response, serum insulin remained high (Figure 5-C), resulting in a reduction in serum glucose to levels below 90 mg × dL-1, but without reaching a true undershoot value thanks to (a) a 45-minute period of minimum insulin infusion(KVO), and (b) the activation of the glucose pump. After a TS of 120 minutes, serum glucose returned to the 90-120 mg × dL-1 range.
I-D Inter-Individual Variation In The Performance of the AIIS
In the year 1994, when the Lehman & Deutsch model was validated clinically, two parameters, Sh y Sp, were introduced in order to reduce the differences between observed and predicted blood glucose. Sh y Sp are, respectively, hepatic and peripheral sensitivities to insulin. Both are adimensional numbers ranging from 0.0 to 1.0.
In order to explore the inter-individual variation in the behavior of the AIIS, we challenged our own model for patients with body weight above and below 70 kg, and hepatic sensitivity (Sh ) above and below 0.5 (Table 1). For the mild positive challenge, the overshoot varied between 143 and 168 mg × dL-1, lasting for 22-102 min., with settling times between 42 and 144 min. For the moderate positive challenges, the overshoot varied between 162 and 183 mg × dL-1, lasting for 30-86 min., with settling times between 60 and 108 min.
For the severe positive challenges, these figures were: overshoot of 170-201 mg × dL-1 , lasting 29-41 min., with a settling time of 30-120 min. There were no "bounce back" undershoots.
II Negative Challenges
Table 2 shows the negative challenge (a perturbation that tends to decrease glycemia) that we expected to be controlled by the AIIS. We assumed that this occurred at the hour 10:00 as a single perturbation within the 24-hour simulation span. The strong negative challenge was the intravenous injection of 10 units of regular insulin. To put this in mathematical terms, we applied the two-compartment model for intravenous insulin bolus described by Kobayashi et al. , using the parameters shown in their Appendix A. After injection, serum insulin concentration rose to 3,300 μU × ml-1 , making the glycemia to fall below 60 mg × dL-1. However, once more the controller of the AIIS responded instantly, both by almost stopping insulin infusion (down to "KVO") and aggressively activating glucose infusion. As a result, the undershoot (glycemia < 60 mg × dL-1) was corrected after 5 min., with a Ts of 90 minutes for it to return to the 90-120 mg × dL-1 range.
In Table 2, in the same 70-Kg patient having hepatic sensitivities (Sh ) of 0.7 and 0.3, the negative challenge produced glycemias of 10 and 31 mg × dL-1 respectively, that lasted 2 and 7 min., with settling times of 22 and 45 min.
A quarter of a century ago, the introduction of capillary blood glucose monitors allowed Type-1 diabetics to intermittently measure their glycemia and adjust their subcutaneous insulin dosage 3 to 4 times a day using empirical algorithms. This so-called "Intensified Insulin Therapy", currently the standard of care for outpatients and an example of "open-loop glucose control", has resulted in a significant reduction in chronic microvascular complications in diabetic patients [35, 36]. Continuous glucose monitoring became possible in 1999 by means of the Continuous Glucose Monitoring System (Medtronic Minimed CGMS, Northridge, CA, USA) [10, 19]. In the last 5 years, the subcutaneous glucose monitoring systems using microdialysis (GlucoDay, A. Menarini Diagnostics, Italy)  and an improved version of the Medtronic Minimed CGMS (Guardian® RT) were introduced allowing real-time glucose readings . In the year 2006, R. Hovorka  summarized 15 clinical studies of closed-loop control in Type 1 diabetics carried out between 1960 and 2004, of which at least nine used subcutaneous sensors. These studies were oriented to outpatient Type-1 diabetics; with insulin being administered either subcutaneously or intraperitoneally according to empirical, semi-automatic and automatic algorithms.
The reader might argue why we have concentrated ourselves on intravenous insulin infusion in hospitalized patients (IVII-H), instead of modeling a system of continuous subcutaneous insulin infusion (CSII-H).
Although CSII has been used in outpatients as well as in hospitalized patients, its use is not advisable in patients requiring Intensive & Intermediate Care, due to the impredictable absorption of insulin from the subcutaneous tissue due to vasoconstriction. The Intravenous Insulin Infusion, on the contrary, allows for 100 per cent bioavailability of insulin even in patients with decreased perfusion of subcutaneous tissues.
In the methods section, we mentioned the use of the Lehman & Deutsch [26, 27] model, in our view a valid choice, considering that some more recent approaches to glucose metabolism tend to underestimate endogenous glucose production at zero insulin [30–32].
As it was mentioned in the introduction, our aim was to explore the closed-loop glucose control under conditions observed in ICU patients. The AIIS included a subcutaneous sensor with the characteristics of the Guardian® RT (Medtronic Minimed) and two infusion pumps for the intravenous administration of regular insulin or glucose. Our "Automatic Insulin Infusion System" proved to be stable even when subjected to the worst conceivable challenges in the clinical setting.
Our main concern when modeling our AIIS was the possibility of catastrophic instability due to the 5-10 min. delay of the subcutaneous glucose sensor. From the point of view of the control theory, it is well known that an important source of instability is precisely the delay in the feedback subsystem. The fact that the AIIS did not become instable even with the worst worst-case challenges, encourages our team to take the next step, i.e., to build a prototype for animal testing.
With respect to the negative challenge, we observed that, the higher the insulin sensitivity, the deepest the glycemia undershoot was, being compensated by a shorter hypoglycemia, and shorter settling times.
In the future, other protocols, perhaps a "Model Predictive Control" (MPC) might be used for inpatient automatic insulin infusion systems, as it has been already proposed for both intravenous  and subcutaneous [35, 39] closed-loop insulin infusion systems. In fact, MPC for subcutaneous insulin infusion has already been tested by Schaller et. al. , who used either intravenous or subcutaneous (simulated by a 30-minute delay) glucose measurements. As we have used only intravenous insulin infusion in our in-Silico model of the AIIS, the MPC system of Schaller et.al. will need to be adapted for intravenous insulin when used in hospitalized patients.
We are aware that the subcutaneous glucose sensors that constitute the essential components of all continuous glucose monitoring devices do not perfectly reflect changes in blood glucose. In fact, during euglycemia, the mean absolute relative differences (MARD's) have been proven by others to be, respectively, 15.2, 21.2, 15.3 and 15.6 per cent, for the four currently available CGMs, i.e., Guardian (Medtronic, Northridge, CA, USA), Dexcom (DexCom, San Diego, CA, USA), Navigator (Abbott Diabetes Care, Alameda, CA, USA), and Glucoday (A. Menarini Diagnostics, Florence, Italy). Of these, the first three are subcutaneous needle-type, and the last one uses a microdyalysis system . However the needle-type system that we decided to include in our model (the same as the one used in Guardian®) has an accuracy that is not significantly different from the others , and constitute now the best available choice.
Comparison of the performance of the AIIS to other existing solutions can be done at three different levels:
In the first level, the averaged results of 1583 ICU patients in 12 "expert-based control studies" reviewed by Moijering S., et.al.  were compared with the results of the AIIS. These studies used sliding scale or titration-based protocols developed on the basis of previous clinical experience . Their reported mean blood glucose of 7.92 ± 1.32 mmol × L-1 (142.6 ± 23.7 mg × dL-1) was obtained with capillary sampling every 1.44 - 3.0 hours. On the other hand, our proposed AIIS was aimed to a lower blood glucose of 100 mg × dL-1 (90 to 120 mg × dL-1), was not dependent on repeated blood sampling, and had a more rigorous criteria for hypoglycemia (i.e., < 60 mg × dL-1).
At the second comparison level are the "model-based methods". In the year 2008, Chase et.al.  published the results of their SPRINT (Specialized Relative INsulin and Nutrition Tables), that regulated insulin administration by being fed data (both current and during the previous 1-2 hours) from 3 sources, i.e., (a) capillary glucose (every 90 minutes); (b) intravenous boluses of insulin over 15-30 seconds, and (c) enteral nutrition by nasogastric tube (not exceeding 90 mL × hour-1). The SPRINT has the merit of having been clinically tested [43, 44], albeit with an average of 16 capillary samples per day. However, the SPRINT has been designed for tube-fed patients only.
The third comparison level is the "artificial pancreas", of which there is only one in current clinical use - the Japanese STG-22. This device is composed of three main parts: the sensor system, a computer and an operation system capable of delivering intravenous solutions of insulin or glucose. The sensor system obtains blood samples every 5 minutes by means of an indwelling double-lumen catheter, that conducts the blood towards a glucose-oxidase sensor. The glucose data is fed into the computer, programmed with a proportional-derivative algorithm whose output controls the delivery of insulin or glucose to the patient. Since the year 1983, the STG-22 has been used in more that 14,000 patients, including clinical and experimental applications . The single and most transcendental difference with our AIIS is the glucose sensing system, as our device, instead of relying on repeated blood sampling, uses a much less invasive subcutaneous glucose sensor.
The fact that after severely challenging the AIIS it did not become instable, even in the context of using a subcutaneous glucose sensor with a 5 to 10 min. delay, means to us that a prototype can be constructed for animal testing, under the expectations of using it in the ICU setting in order to achieve tight control of blood glucose without frequent blood sampling.
Indeed, we believe that the closed-loop control, with subcutaneous glucose sensing (Figure 6) and intravenous insulin administration constitutes a promising approach for developing new devices oriented to the tight control of glycemia in hospitalized patients, particularly in the intensive care unit. This approach has two advantages. First, it provides a new tool to reduce mortality and morbidity in the ICU. Last, but not least, this system shall be tested in an environment with highly trained personnel and 24-hour surveillance, so the experience thus gained might one day extend its benefits to the outpatient setting, i.e., to millions of insulin-requiring diabetics worldwide.
JO is a BscEE (2006) and MscEE (2008)
MG is a Full Professor of Electrical Engineering, with a PhD in Electrical Engineering (University of Arizona, U.S.A. 1989)
GB is Associate Professor of Medicine (Department of Respiratory Diseases), Respiratory Diseases Specialist (The Ohio State University College of Medicine, USA, 1993), and has a PhD in Medical Biochemistry (The Ohio State University Graduate School, USA, 1993)
PO is Associate Professor of Medicine (Department of Diabetes, Nutrition & Metabolism), Endocrinologist (The Ohio State University College of Medicine, USA, 1993) and has a MSc in Biomedical Engineering (The Ohio State University Graduate School, USA, 1993)
(greek letter alpha) Sensitivity (nA × [mg × dl-1]-1) [Eq. 1]
Differential Parameter (mg × dL-1) [Eq. 4 and 5]
Automatic Insulin Infusion System (for inpatient treatment)
Blood glucose level (mg × dL-1) at which half-maximum infusion rate is chosen to occur [Eq. 2a]
Continuous Glucose Monitoring System
Differential Factor (mg × dL-1) [Eq. 2a and 3]
Food and Drug Administration (USA)
Glycemia (blood glucose concentration, mg × dL-1) [Eq. 1]
- G0 :
Last 5-minute averaged glycemia (mg × dL-1) [Eq. 5 and 6]
- Gcrit :
Critical glycemia (mg × dL-1) below which the glucose infusion controller activates itself [Eq.2b and 6b]
- Gi :
Glycemia (mg × dL-1) determined "i" minutes ago [Eq. 4 and 5]
Glomerular Filtration Rate (mL × min-1), a measure of renal function
Glucose Transporter - 4 (present in muscular and adipose tissues)
Glucose infusion rate (mL × h-1) [Eqs. 6a and 6b]. Note that it is in milliliters per hour (and not in mL × min-1) as the commercially available pumps are graduated in this way
Infusion rate in mL × h-1 as determined by the glucose infusion controller [Eq. 6b]
Intensive Care Unit
Insulin Infusion Rate (mIU × min-1) [Eq.2]
Insulin infusion rate (mL × h-1) determined by the insulin infusion controller
- K1 :
= 1 [mg × dL-1]-2 , a parameter [Eq. 3] 
- K2 :
= 10 [adimensional], a parameter [Eq.3] 
- KG :
Controller gain ([mL × h-1] × [mg × dL-1]-1) [Eq. 6a]
"Keep Vein Open". Minimal infusion rate (= 1.0 mL × h-1) needed to keep venous access from being occluded by thrombosis
Maximal possible insulin infusion rate (mIU × min-1) [Eq. 2]
Net Hepatic Glucose Balance (mg × min-1)
Subcutaneous Glucose Monitor
Slope (dL × mg-1) of the hyperbolic tangent curve [Eq. 2]
Peripheral sensitivity to insulin (adimensional 0.0 ≤ SP ≤ 1.0)
Hepatic sensitivity to insulin (adimensiona l0.0 ≤ Sh ≤ 1.0)
(greek letter tau). Time constant (minutes) [Eq. 1]
- TD :
Transport lag (minutes) [Eq. 1]
Imprimerie L, Vanmelle SA (Eds): The prevalence of diabetes In Diabetes Atlas, Global estimates of Diabetes. International Diabetes Federation. 2nd edition. Gent/Mariakerke, Belgium; 2003:28–30.
Clement S, Braithewaite SS, Magee MF, Ahmann A, Smith EP, Schafer RG, Hirsch IB: Management of diabetes and hyperglycaemia in hospitals. Diabetes Care 2004, 27: 553–591. 10.2337/diacare.27.2.553
Hanania G, Cambou J-P, Guéret P, Vaur L, Blanchard D, Lablanche J-M, Boutalbi Y, Humbert R, Clerson P, Genès N, Danchin N: Management and in-hospital outcome of patients with acute myocardial infarction admitted to intensive care units at the turn of the century: results from the French nationwide USIC 2000 registry. Heart 2004, 90: 1404–1410. 10.1136/hrt.2003.025460
Whitcomb BW, Pradhan EK, Pittas AG, Roghmann RC, Perencevich RN: Impact of admission hyperglycaemia on hospital mortality in various intensive care unit populations. Crit Care Med 2005, 33: 2772–2777. 10.1097/01.CCM.0000189741.44071.25
Van Den Berghe G, Wouters P, Weekers F, Verwaest C, Bruyninckx F, Schetz M, Vlasselaers D, Ferdinande P, Lauwers P, Bouillon R: Intensive insulin therapy in the critically ill patients. N Eng J Med 2001, 345: 1359–1367. 10.1056/NEJMoa011300
Bloomgarden ZT: Inpatient diabetes control: rationale. Diabetes Care 2004, 27: 2074–2080. 10.2337/diacare.27.8.2074
Furnary AP, Gao G, Grunkemeier GL, Wu Y, Zerr KJ, Bookin SO, Floten HS, Starr A: Continuous insulin infusion reduces mortality in patients with diabetes undergoing coronary artery bypass grafting. J Thorac Cardiovasc Surg 2003, 125: 1007–1021. 10.1067/mtc.2003.181
Conner TM, Flesner-Gurley KR, Barner JC: Hyperglycaemia in the hospital setting: the case for improved control among non-diabetics. Ann Pharmacother 2005, 39: 492–501. 10.1345/aph.1E308
Wild C, Narath M: Evaluating and planning ICUs: methods and approaches to differentiate between need and demand. Health Policy 2005, 71: 289–301. 10.1016/j.healthpol.2003.12.020
Bode BW, Gross TM, Thornton KR, Mastrototaro JJ: Continuous glucose monitoring used to adjust diabetes therapy improves glycosilated hemoglobin: a pilot study. Diab Res Clin Pract 1999, 46: 183–190. 10.1016/S0168-8227(99)00113-8
Warren J, McLaughlin D: Medtronic announces FDA approval of real-time continuous glucose monitoring system. News Release, Medtronic, Inc; 2006. [http://Medtronic.com/Newsroom/NewsReleaseDetails.do?itemId1123700847661&format=pdf&lang]
Howorka K: Initial algorithms, "K" and the blood glucose target value. In Functional Insulin Treatment, Principles, teaching approach and practice. 2nd edition. Edited by: Howorka K. Springer-Verlag, Berlin Heidelberg; 1996:45–49. Section 4.1
Guerci B, Floriot M, Böhme P, Durain D, Benichou M, Jellimann S, Drouin P: Clinical performance of CGMS in Type-1 diabetic patients treated by continuous subcutaneous insulin infusion using insulin analogs. Diabetes Care 2003, 26: 582–589. 10.2337/diacare.26.3.582
Monsod TP, Flanagan DE, Rife F, Saenz R, Caprio S, Sherwin RS, Tamborlane WV: Do sensor glucose levels accurately predict plasma glucose concentrations during hypoglycemia and hyperinsulinemia? Diabetes Care 2002, 25: 889–893. 10.2337/diacare.25.5.889
Chee F, Fernando T, Van Heerden P: Closed-loop glucose control in critically ill patients using Continuous Glucose Monitoring System (CGMS) in real time. IEEE Trans Inf Technol Biomed 2003, 7: 43–53. 10.1109/TITB.2003.808509
Steil G, Pantelon A, Rebrin K: Closed-loop insulin delivery - the path to physiological glucose control. Adv Drug Deliver Rev 2004, 56: 125–144. 10.1016/j.addr.2003.08.011
Gross T, Mastrototaro J, Shin J: Improved performance of the Medtronic MiniMed CGMS using redesigned glucose sensors and updated calibration software. Diabetes 2003, 52(Suppl 1):A95.
Steil GM, Rebrin K, Mastrototaro J, Bernaba B, Saad MF: Determination of plasma glucose during rapid glucose excursions with a subcutaneous glucose sensor. Diabetes Technol Ther 2003, 5: 27–31. 10.1089/152091503763816436
Rebrin K, Steil GM, Van Antwerp W, Mastrototaro J: Subcutaneous glucose predicts plasma glucose independent of insulin: Implications for continuous monitoring. Am J Physiol 1999, 277: E561-E571.
Marliss EB, Murray FT, Stokes EF, Zinman B, Nakhooda AF, Denoga A, Leibel BS, Albisser AM: Normalization of glycaemia in diabetics during meals with insulin and glucagon delivery by the artificial pancreas. Diabetes 1977, 26: 663–672. 10.2337/diabetes.26.7.663
Broekhuyse H, Nelson J, Zinman B, Albisser A: Comparision of algorithms for the closed-loop control of blood glucose using the artificial beta cell. IEEE T Bio-Med Eng 1981, 28: 678–687. 10.1109/TBME.1981.324661
Albisser A, Leibel B, Ewart T, Davidovac Z, Botz C, Zingg W: An artificial endocrine pancreas. Diabetes 1974, 23: 389–396.
Albisser A, Leibel B, Ewart T, Davidovac Z, Botz C, Zingg W, Schipper H, Gander R: Clinical control of diabetes by the artificial pancreas. Diabetes 1977, 23: 397–404.
Parker RS, Doyle-III FJ, Peppas NA: The Intravenous route to blood glucose control. IEEE Eng Med Biol 2001, 20: 65–73. 10.1109/51.897829
Yamasaki Y, Shichiri M, Kawamori R, Goriya Y, Sasai T, Morishima T, Nomura M, Tohdo R, Abe H: Counterregulatory system in an artificial endocrine pancreas. Glucose infusion algorithm. Artif Organs 1979, 3: 265–270. 10.1111/j.1525-1594.1979.tb01061.x
IMED Gemini Infusion Systems: Alaris Medical Systems. 2000. [http://www.alarismed.com/na/products/pdfs/Geminisingle.pdf]
Lehmann E, Deutsch T: A physiological model of glucose-insulin interaction in type 1 diabetes mellitus. J Biomed Eng 1992, 14: 235–242. 10.1016/0141-5425(92)90058-S
Lehmann E, Hermanyi I, Deutsch T: Retrospective validation of a physiological model of glucose-insulin interaction in diabetes mellitus. Med Eng Phys 1994, 16: 193–202. 10.1016/1350-4533(94)90038-8
Guyton JR, Foster RO, Soeldner JS, Tan MH, Kahn CB, Konez L, Gleason RE: A model of glucose-insulin homeostasis in man that incorporates the heterogeneous fast pool theory of pancreatic insulin release. Diabetes 1978, 27: 1027–1042. 10.2337/diabetes.27.10.1027
Dalla Man C, Rizza RA, Cobelli C: Meal simulation model of the glucose-insulin system. IEEE T Bio-Med Eng 2007, 54: 1740–1749. 10.1109/TBME.2007.893506
Boden G, Cheung P, Homko C: Effects of acute insulin excess and deficiency on gluconeogenesis and glycogenolysis in type 1 diabetes. Diabetes 2003, 52: 133–137. 10.2337/diabetes.52.1.133
Sherwin RS, Sacca L, Eigler N, Hendler R, Fradkin J, Shamoon H, Walesky M, Felig P: Influence of glucoregulatory hormones and blood glucose per se on the regulation of hepatic glucose production. In Carbohydrate metabolism. Volume Chapter 7. Edited by: Cobelli C, Bergman RN. John Wiley & Sons Ltd., Chichester, UK; 1981:127–152.
Marshall SM, Walker M, Alberti KGMM: Diabetic ketoacidosis and hyperglycaemic non-ketotic coma. In International Textbook of Diabetes Mellitus, KGMM. 2nd edition. Edited by: Alberti, Zimmet P, DeFronzo RA, Keen H. John Wiley & Sons Chichester, West Sussex, England; 1997:1215–1229.
Kobayashi T, Sawano S, Itoh T, Kosaka K, Kasuya Y: The pharmacokinetics of insulin after continuous subcutaneous infusion or bolus subcutaneous injection in diabetic patients. Diabetes 1983, 32: 331–336. 10.2337/diabetes.32.4.331
The Diabetes Control and Complications Trial Research Group: The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus. N Eng J Med 1993, 329: 977–986. 10.1056/NEJM199309303291401
Ohkubo Y, Kishikawa H, Araki E, Myyata T, Isami S, Motoyoshy S, Kojima Y, Furuyoshi N, Shichiri M: Intensive insulin therapy prevents the progression of diabetic microvascular complications in Japanese patients with non-insulin-dependent diabetes mellitus: A randomized prospective 6-year study. Diab Res Clin Pract 1995, 28: 103–117. 10.1016/0168-8227(95)01064-K
Maran A, Poscia A: Continuous subcutaneous glucose monitoring: The GlucoDay® system. Diabetes Nutr Metab 2002, 15: 429–433.
Hovorka R: Continuous glucose monitoring and closed-loop systems. Diabetic Med 2006, 23: 1–12. 10.1111/j.1464-5491.2005.01672.x
Steil G, Rebrin K, Mastrototaro J: Metabolic modeling and the closed-loop insulin delivery problem. Diab Res Clin Pract 2006, 74: S183-S186. 10.1016/S0168-8227(06)70028-6
Schaller HC, Schaupp L, Bodenlenz M, Wilinska ME, Chassin LJ, Wach P, Vering T, Hovorka R, Pieber TR: On-line adaptive algorithm with glucose prediction capacity for subcutaneous closed-loop control of glucose: Evaluation under fasting conditions in patients with type 1 diabetes. Diabetic Med 2006, 23: 90–93. 10.1111/j.1464-5491.2006.01695.x
Kovatchev B, Heinemann L, Anderson S, Clarke W: Comparison of the numerical and clinical accuracy of four continuous glucose monitors. Diabetes Care 2008, 31: 1160–1164. 10.2337/dc07-2401
Meijering S, Corstjens A, Tullecken JE, Meertens JHJM, Zijlstra JG, Ligtemberg JJM: Towards a feasible algorithm for tight glycemic control in critically ill patients: a systematic review of the literature. Critical Care 2001, 10: R19. (doi: 10.1186/cc3981) 10.1186/cc3981
Chase JG, LeCompte A, Shaw GM, Blakemore A, Wong J, Lin J, Hann CE: A bencmark data set for model-based glycemic control in critical care. J Diabetes Sci Technol 2008, 2(4):584–594.
Chase JG, Shaw G, LeCompte A, Lonergan T, Willacy M, Wong XW, Lin J, Lee D, Hann C: Implementation and evaluation of the SPRINT protocol for tight glycemic control in critically ill patients: A clinical practice change. Critical Care 2008, 12: R49. (doi: 10.1186/cc6868) 10.1186/cc6868
Nishida K, Shimoda S, Ichinose K, Araki E, Motoaki S: What is endocrine artificial pancreas? - Mechanism and history. World J Gastroenterol 2009, 15(33):4105–4110. 10.3748/wjg.15.4105
The authors are grateful to Mrs Carolina Torres for her excellent Secretarial Assistance.
Grants: financed in part by DIPUC Grant # 003/2003
The authors declare that they have no competing interests.
JO made his MSc Thesis on Electronic Engineering on this subject. MG was adviser to JO in his thesis, contributed to write and reviewed the manuscript. GB contributed to write and critically reviewed the manuscript. PO was adviser to JO in his thesis, contributed to write the manuscript and drew the illustrations. All authors have read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Ortiz, J.L., Guarini, M.W., Borzone, G.R. et al. In silico evaluation of a control system and algorithm for automated insulin infusion in the ICU setting. BioMed Eng OnLine 9, 35 (2010). https://doi.org/10.1186/1475-925X-9-35
- Insulin Infusion
- Model Predictive Control
- Glucose Infusion
- Continuous Subcutaneous Insulin Infusion
- Intensive Care Unit Setting