Search tips
Search criteria 


Logo of bmengonBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleBioMedical Engineering OnLine
Biomed Eng Online. 2010; 9: 35.
Published online 2010 July 20. doi:  10.1186/1475-925X-9-35
PMCID: PMC2918623

In silico evaluation of a control system and algorithm for automated insulin infusion in the ICU setting



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.


Diabetes Mellitus affects 5.1 per cent of the world's adult population [1], with a prevalence of up to 12.4% among hospitalized patients [2].

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 [7].

The current standard of care for hospitalized diabetic patients (Figure (Figure1)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.

Figure 1
Current in-patient management of intravenous insulin infusion. Current in-patient management of intravenous insulin infusion(see text) Nurses have to measure capillary blood glucose every 1-6 hours, and after manually recording date, time and value, they ...

Thus, there is a need for automated protocols of intravenous insulin administration for hyperglycemic patients in the ICU [8,9].

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 Figures22 and and3.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).

Figure 2
System diagram of the "Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment. System diagram of the "Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment(see text). Glucose concentration measured by the subcutaneous sensor is ...
Figure 3
"Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment: Proposed architecture. Proposed architecture of the "Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment (see text). The nano amperimetric current produced by the subcutaneous ...

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 [10] 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 (on the use of CGMS in critically ill patients) [15], and of Steil [16], 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,

equation image

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 [16]. A white Gaussian noise to the sensor signal was added to our model, assuming a mean = 0.0 and a variance = 100 [17].

In their publication, Steil et. al. [16] 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 depicted a time constant (τ) of 5.77 minutes [16,18], which is thought to be primarily related to interstitial glucose equilibration.

α in equation 1 is equal to 0.23 nA × [mg × dl-1] -1 [19]. Using data provided by the manufacturer we included lower and upper limit saturations of 40 and 400 mg × dl-1, respectively.

I-C. Controller

The controller of the AIIS can be divided into three main parts: the insulin infusion pump controller, the glucose infusion pump controller, and the emergency infusion logic (Figures (Figures22 and and33).

The insulin infusion pump controller is based on the algorithm proposed by Marliss in 1977[20], in the form later published in 1981 by Broekhuyse [21], whose original equations were as follows,

equation image
equation image

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, G0[mg × dL-1] is the last available average glycemia, Gcrit 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 IRmin [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,

equation image

where the parameters K1 = 1 [mg × dL-1]-2 and K2 = 10 [nondimensional] [21], and the original method for calculating A [mg × dL-1] is as follows,

equation image

where Gi [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):

equation image

where G0 [mg × dL-1] is the last reading of glycemia, and G5 and G10 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 [21] and the glycemia and insulinemia profiles from a 69.5 Kg male type 1 diabetic patient being controlled by the artificial pancreas of Albisser [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 Figure4,4, but using equation (5) instead of equation (4).

Figure 4
"Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment: Simulation of mild positive challenge. "Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment: Simulation of mild positive challenge given by a standard meal containing 50 ...

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. [24].

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 [25]. The controller equation is given by equation 6,

equation image
equation image

where GR is the glucose infusion rate in ml × h-1, Gcrit is the critical serum glucose level in mg × dL-1 below which the infusion must begin, and G0 is the last 5-minute-averaged glycemia reading. KG is the controller gain. GR is always positive because the controller activates itself only when G0 falls below Gcrit. 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 (Gcrit) 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. IRmin≠ 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.

I-D Actuators

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 [26]. 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 (Figure33).

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) [29].

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 [27].

More recently, others [30] 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 tables11 and and2.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.

Table 1
Positive challenges.
Table 2
Negative challenge in 70-Kg type-1 diabetics.


I Positive Challenges

Table Table11 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 Table11 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) [33]. 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 Table1,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 Table11 and in Figure Figure5,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.

Figure 5
"Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment: Severe positive challenge. "Automatic Insulin Infusion System" (AIIS) for Inpatient Treatment: Severe positive challenge given by the sudden intravenous injection of 10 grams of glucose ...

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 (Table1).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 Table22 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. [34], 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 Table2,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) [37] and an improved version of the Medtronic Minimed CGMS (Guardian® RT) were introduced allowing real-time glucose readings [11]. In the year 2006, R. Hovorka [38] 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 [24] and subcutaneous [35,39] closed-loop insulin infusion systems. In fact, MPC for subcutaneous insulin infusion has already been tested by Schaller et. al. [40], 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 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 [41]. 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 [41], 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., [42] 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 [43]. 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 [43] 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 [45]. 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 (Figure6)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.

Figure 6
A commercially available subcutaneous glucose sensor. Five views (A, B, C, D and E). (A) A commercially available subcutaneous glucose sensor. (B) Exploded-view of a generic glucose sensor structure. (C) The tip of the glucose sensor with its three layers. ...

List of Abbreviations used

α: (greek letter alpha) Sensitivity (nA × [mg × dl-1] -1) [Eq. 1]; A, Differential Parameter (mg × dL-1) [Eq. 4 and 5]; AIIS: Automatic Insulin Infusion System (for inpatient treatment); B: Blood glucose level (mg × dL-1) at which half-maximum infusion rate is chosen to occur [Eq. 2a]; CGMS: Continuous Glucose Monitoring System; DF: Differential Factor (mg × dL-1) [Eq. 2a and 3]; FDA: Food and Drug Administration (USA); G: 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]; GFR: Glomerular Filtration Rate (mL × min-1), a measure of renal function; GLUT-4: Glucose Transporter - 4 (present in muscular and adipose tissues); GR: 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; GR': Infusion rate in mL × h-1 as determined by the glucose infusion controller [Eq. 6b]; ICU: Intensive Care Unit; IR: Insulin Infusion Rate (mIU × min-1) [Eq.2]; IR': Insulin infusion rate (mL × h-1) determined by the insulin infusion controller; K1: = 1 [mg × dL-1]-2 , a parameter [Eq. 3] [21]; K2: = 10 [adimensional], a parameter [Eq.3] [21]; KG: Controller gain ([mL × h-1] × [mg × dL-1]-1) [Eq. 6a]; KVO: "Keep Vein Open". Minimal infusion rate (= 1.0 mL × h-1) needed to keep venous access from being occluded by thrombosis; M: Maximal possible insulin infusion rate (mIU × min-1) [Eq. 2]; NHGB: Net Hepatic Glucose Balance (mg × min-1); SGM: Subcutaneous Glucose Monitor; S: Slope (dL × mg-1) of the hyperbolic tangent curve [Eq. 2]; SP: Peripheral sensitivity to insulin (adimensional; 0.0 ≤ SP ≤ 1.0); Sh: Hepatic sensitivity to insulin (adimensional; 0.0 ≤ Sh ≤ 1.0); τ: (greek letter tau). Time constant (minutes) [Eq. 1]; TD: Transport lag (minutes) [Eq. 1]

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

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.

Author's information

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)


The authors are grateful to Mrs Carolina Torres for her excellent Secretarial Assistance.

Grants: financed in part by DIPUC Grant # 003/2003


  • Imprimerie L, Vanmelle SA, editor. Diabetes Atlas, Global estimates of Diabetes. International Diabetes Federation. 2. Gent/Mariakerke, Belgium; 2003. The prevalence of diabetes; pp. 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. doi: 10.2337/diacare.27.2.553. [PubMed] [Cross Ref]
  • 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. doi: 10.1136/hrt.2003.025460. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1097/01.CCM.0000189741.44071.25. [PubMed] [Cross Ref]
  • 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. doi: 10.1056/NEJMoa011300. [PubMed] [Cross Ref]
  • Bloomgarden ZT. Inpatient diabetes control: rationale. Diabetes Care. 2004;27:2074–2080. doi: 10.2337/diacare.27.8.2074. [PubMed] [Cross Ref]
  • 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. doi: 10.1067/mtc.2003.181. [PubMed] [Cross Ref]
  • 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. doi: 10.1345/aph.1E308. [PubMed] [Cross Ref]
  • Wild C, Narath M. Evaluating and planning ICUs: methods and approaches to differentiate between need and demand. Health Policy. 2005;71:289–301. doi: 10.1016/j.healthpol.2003.12.020. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/S0168-8227(99)00113-8. [PubMed] [Cross Ref]
  • Warren J, McLaughlin D. Medtronic announces FDA approval of real-time continuous glucose monitoring system. News Release, Medtronic, Inc; 2006. [cited June 5th, 2006]
  • Howorka K. In: Functional Insulin Treatment, Principles, teaching approach and practice. 2. Howorka K, editor. Springer-Verlag, Berlin Heidelberg; 1996. Initial algorithms, "K" and the blood glucose target value; pp. 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. doi: 10.2337/diacare.26.3.582. [PubMed] [Cross Ref]
  • 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. doi: 10.2337/diacare.25.5.889. [PubMed] [Cross Ref]
  • 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. doi: 10.1109/TITB.2003.808509. [PubMed] [Cross Ref]
  • Steil G, Pantelon A, Rebrin K. Closed-loop insulin delivery - the path to physiological glucose control. Adv Drug Deliver Rev. 2004;56:125–144. doi: 10.1016/j.addr.2003.08.011. [PubMed] [Cross Ref]
  • 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. doi: 10.1089/152091503763816436. [PubMed] [Cross Ref]
  • 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. [PubMed]
  • 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. doi: 10.2337/diabetes.26.7.663. [PubMed] [Cross Ref]
  • 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. doi: 10.1109/TBME.1981.324661. [PubMed] [Cross Ref]
  • Albisser A, Leibel B, Ewart T, Davidovac Z, Botz C, Zingg W. An artificial endocrine pancreas. Diabetes. 1974;23:389–396. [PubMed]
  • 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. [PubMed]
  • Parker RS, Doyle-III FJ, Peppas NA. The Intravenous route to blood glucose control. IEEE Eng Med Biol. 2001;20:65–73. doi: 10.1109/51.897829. [PubMed] [Cross Ref]
  • 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. doi: 10.1111/j.1525-1594.1979.tb01061.x. [PubMed] [Cross Ref]
  • IMED Gemini Infusion Systems. Alaris Medical Systems. 2000. Accessed March 30th, 2009.
  • Lehmann E, Deutsch T. A physiological model of glucose-insulin interaction in type 1 diabetes mellitus. J Biomed Eng. 1992;14:235–242. doi: 10.1016/0141-5425(92)90058-S. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/1350-4533(94)90038-8. [PubMed] [Cross Ref]
  • 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. doi: 10.2337/diabetes.27.10.1027. [PubMed] [Cross Ref]
  • Dalla Man C, Rizza RA, Cobelli C. Meal simulation model of the glucose-insulin system. IEEE T Bio-Med Eng. 2007;54:1740–1749. doi: 10.1109/TBME.2007.893506. [PubMed] [Cross Ref]
  • 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. doi: 10.2337/diabetes.52.1.133. [PubMed] [Cross Ref]
  • Sherwin RS, Sacca L, Eigler N, Hendler R, Fradkin J, Shamoon H, Walesky M, Felig P. In: Carbohydrate metabolism. Cobelli C, Bergman RN, editor. Chapter 7. John Wiley & Sons Ltd., Chichester, UK; 1981. Influence of glucoregulatory hormones and blood glucose per se on the regulation of hepatic glucose production; pp. 127–152.
  • Marshall SM, Walker M, Alberti KGMM. In: International Textbook of Diabetes Mellitus, KGMM. 2. Alberti, Zimmet P, DeFronzo RA, Keen H, editor. John Wiley & Sons Chichester, West Sussex, England; 1997. Diabetic ketoacidosis and hyperglycaemic non-ketotic coma; pp. 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. doi: 10.2337/diabetes.32.4.331. [PubMed] [Cross Ref]
  • 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. doi: 10.1056/NEJM199309303291401. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/0168-8227(95)01064-K. [PubMed] [Cross Ref]
  • Maran A, Poscia A. Continuous subcutaneous glucose monitoring: The GlucoDay® system. Diabetes Nutr Metab. 2002;15:429–433. [PubMed]
  • Hovorka R. Continuous glucose monitoring and closed-loop systems. Diabetic Med. 2006;23:1–12. doi: 10.1111/j.1464-5491.2005.01672.x. [PubMed] [Cross Ref]
  • Steil G, Rebrin K, Mastrototaro J. Metabolic modeling and the closed-loop insulin delivery problem. Diab Res Clin Pract. 2006;74:S183–S186. doi: 10.1016/S0168-8227(06)70028-6. [PubMed] [Cross Ref]
  • 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. doi: 10.1111/j.1464-5491.2006.01695.x. [PubMed] [Cross Ref]
  • 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. doi: 10.2337/dc07-2401. [PubMed] [Cross Ref]
  • 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. (doi: 10.1186/cc3981) [PMC free article] [PubMed] [Cross Ref]
  • 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. [PMC free article] [PubMed]
  • 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. (doi: 10.1186/cc6868) [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.3748/wjg.15.4105. [PMC free article] [PubMed] [Cross Ref]

Articles from BioMedical Engineering OnLine are provided here courtesy of BioMed Central