Abstract
Background:
The ability to simulate in silico experiments is crucial for fast and cost-effective preliminary studies prior to clinical trials. We present an in silico approach to the design of optimal pramlintide-to-insulin (P/I) ratios, using our computer simulator of the human metabolic system, with a population of virtual adult type 1 diabetes mellitus patients and with individual parameters modified to account for the dynamic effects of pramlintide.
Materials and Methods:
A model of pramlintide action on gastric emptying was built using data of 15 type 1 diabetes mellitus subjects studied twice with a standardized dual-tracer meal on placebo and pramlintide, which was incorporated in our type 1 diabetes simulator. Extensive in silico experiments on 100 virtual subjects were performed to optimize the co-administration of pramlintide and insulin prior to its submission to clinical trials; several P/I ratios were tested in terms of efficacy, in attenuating postprandial hyperglycemia, and in hypoglycemia safety.
Results:
In silico experiments estimated the optimal P/I ratio to be 9 μg of pramlintide per unit (U) of insulin. Additional simulations narrowing the investigated range indicated that P/I ratios of 8 and 10 μg/U would achieve similar performance. Moreover, simulation results suggested that in clinical trials, insulin boluses should be reduced by approximately 21% at a P/I ratio of 9 μg/U to account for the effects of pramlintide and avoid postprandial hypoglycemia.
Conclusions:
We can assert that a valid simulation model of pramlintide action was developed, leading to in silico estimation of optimal pramlintide:insulin co-administration ratio. Clinical trials will confirm (or adjust) this initial estimation.
Introduction
Typically, the in silico models of the human metabolism are compartmental (e.g., they represent the human body as a set of homogeneous compartments of specific concentrations and volumes linked by diffusion or rate-limited pathways). Classical examples include the Widmark 6 model of ethanol pharmacokinetics, which offers a straightforward interpretation with a constant ethanol clearance rate and the human body modeled as one compartment, or the more complex minimal model of glucose kinetics suggested by Bergman et al. 7 more than 30 years ago to assess insulin sensitivity in health and diabetes. Deeper understanding of the processes involved and novel measuring tools have allowed more precise measurement of glucose–insulin dynamics and the development of sophisticated nonlinear models, such as the vastly complex meal model, 8 which now serves as the basis for the most comprehensive computer simulator of the human metabolic system available to date. 9 –11
This simulator was developed with the specific goal to facilitate the preclinical testing of treatment strategies in diabetes. The simulation environment encompasses several metabolic subsystems, including gastrointestinal tract, renal function, hepatic glucose production, and others. 8,9 Furthermore, the simulation system was equipped with in silico images of the metabolic systems of 300 subjects with type 1 diabetes in three age groups: children, adolescents, and adults. Each in silico image is a set of 26 parameters uniquely identifying the metabolic system of a person. The distribution of these parameters across the general population provides the variability needed for defining the in silico population. 10 When the in silico model and the virtual population were built, the computer simulation of the observed biosystem became possible, leading to in silico trials involving virtual “subjects” rather than real people. Because such in silico trials can serve as cost-effective precursors guiding expensive and time-consuming clinical investigations by ruling out ineffective treatment approaches, in January 2008 this simulator received Food and Drug Administration (FDA) acceptance for the preclinical testing of control strategies in artificial pancreas studies. 12 The simulator was immediately put to its intended use: since 2008 several scientific publications have reported in silico results from the testing of various insulin treatment strategies, 8 –11 and several Investigational Device Exemptions were granted by the FDA for closed-loop clinical trials, using solely in silico experiments to test and validate their control systems. Two reviews 3,11 have presented the wide array of modeling approaches developed by us and others and their translation into the clinical practice. Thus, we were able to assert that (1) in silico modeling can produce credible preclinical results and (2) in silico testing yields results in a fraction of the time and resources required for animal trials.
Pramlintide is a synthetic analog of human amylin: a 37-amino-acid peptide hormone that is co-located with insulin in secretory granules and co-secreted by pancreatic β-cells in response to nutrient ingestion. Amylin has been shown to contribute to glucose control during the postprandial period by inhibiting glucagon secretion, slowing gastric emptying and reducing food intake via increased satiety, resulting in decreased appetite. 13 However, like insulin, amylin is deficient in individuals with diabetes requiring insulin therapy. 14 In these patients, maintaining postprandial glucose concentrations within the desired range remains challenging despite the introduction of rapid-acting insulin analogs and enhanced insulin delivery systems. Even so, it has been shown that pramlintide administration restores the effects of amylin and, in turn, lowers average postprandial blood glucose levels, substantially reduces postprandial blood glucose excursions, and decreases food intake, resulting in a reduction in mealtime insulin requirements. 15 –22 Pramlintide has been approved by the FDA, for use by patients with type 1 and type 2 diabetes who use mealtime insulin. The therapeutic effects and safety profile of pramlintide, when used with various rapid-acting insulin formulations, are well characterized. In pilot experiments, continuous subcutaneous co-administration of insulin and pramlintide appeared to be safe and well tolerated and reduced both fasting and postprandial glucose levels. 23 –25 However, an optimal pramlintide-to-insulin (P/I) ratio is yet to be determined.
In this article, in an effort to guide the development of a co-formulated pramlintide/insulin product, we present an in silico approach to the identification of optimal P/I ratios for use preprandially. We use our computer simulation environment with the population of virtual adults with type 1 diabetes and with individual parameters modified to account for the dynamic effects of pramlintide as observed in the clinic. 26,27
Materials and Methods
Dataset
A previously published dataset 27 was used to model the effect of pramlintide on glucose meal rate of appearance. The data included 15 subjects with type 1 diabetes (eight men and seven women, 37±2 years of age, body weight 76±3 kg). The study design included a standardized meal test containing 50 g of glucose enriched with 3 g of 6,6-dideuteroglucose, administered with and without 30 μg of pramlintide injected subcutaneously in the lower abdominal wall. Over the initial 90 min of the postprandial period, blood samples were taken at 15-min intervals and thereafter at 30-min intervals until completion of the experiment at 330 min. Concentrations of plasma glucose, tracer glucose, and insulin were measured. A detailed description of the subjects, protocol, and methods can be found in Woerle et al. 27
Modeling the effect of pramlintide
A model of oral glucose absorption (Fig. 1, upper panel) based on the previously introduced meal model 8 was used to quantify the effect of pramlintide on meal glucose rate of appearance (Ra). In brief, glucose transit through the stomach (Q sto) and intestine was modeled by two compartments, one for solid (Q sto1) and one for liquid (Q sto2) phase, whereas a single compartment was used to describe the gut (Q gut). k gri (min−1) is the rate of grinding, k empt(Q sto) (min−1) the rate constant of gastric emptying (a nonlinear function of Q sto), and k abs (min−1) the rate constant of intestinal absorption. Ra (mg/kg/min) is the rate of appearance of glucose in plasma. The parameter k gri was assumed equal to k max. 8 The model was numerically identified on placebo and pramlintide, yielding the parameter estimates presented in Figure 1, lower panel.

The relative variation of gastrointestinal tract parameters due to pramlintide was calculated as:
where
Assuming normality, joint parameter distribution was derived from mean and the covariance of Δpij . To achieve intersubject variability, we used this parameter distribution to randomly generate n=100 distinct sets of parameters (i.e., n=100 in silico “adults” with type 1 diabetes responding to pramlintide). Thus, the effect of pramlintide on Ra was incorporated in the simulator via modification of the parameters of the virtual subjects corresponding to the empirically determined effects of administration of pramlintide. 27
In silico experiments
The simulation scenario included a meal with a carbohydrate content of 50 g given at 8:00 h to n=100 in silico adults with type 1 diabetes. Concurrent with the meal, placebo and several different P/I ratios were administered: P/I=3, 6, 8, 9, 10, and 12 μg of pramlintide/unit (U) of insulin (μg/U). On placebo, the virtual subjects received a premeal insulin bolus based on each individual's insulin-to-carbohydrate ratio (ICR). With pramlintide, two in silico experiments were performed.
Experiment 1
The virtual subjects received a premeal insulin bolus identical to the bolus used on placebo, without adjustment for the effects of pramlintide. As seen in the next section this resulted in postprandial hypoglycemia, which necessitated Experiment 2.
Experiment 2
The virtual subjects received premeal insulin bolus lowered to account for the effects of pramlintide. The adjustment of the insulin dose (i.e., of each individual's ICR) was done iteratively for each subject, until postprandial hypoglycemia due to insulin overdose was avoided. Then the same P/I ratios were administered again.
Data analysis
Efficacy (attenuation of postprandial hyperglycemia) and safety (reduction of hypoglycemia) of different P/I ratios have been quantitatively evaluated using control variability-grid analysis (CVGA).
28
In addition, to select the best P/I candidates, a CVGA Zone Quotient index was calculated using the formula:
where A, B, C, D, and E are the percentages of subjects in the CVGA Zones A–E. 28
CVGA provides a simultaneous visual and numerical assessment of the overall quality of glycemic regulation in the entire population. In particular, CVGA zones represent accurate control (Zone A), benign deviation into hyper- or hypoglycemia (Zone B), overcorrection of hyper- or hypoglycemia (Zone C), failure to deal with hyper- or hypoglycemia (Zone D), and erroneous control (Zone E). 28
Results
Prediction of glucose Ra
Figure 2 compares the experimentally observed Ra of glucose on placebo and pramlintide (upper panel) and the model-predicted Ra on placebo and pramlintide (lower panel), indicating good agreement between experimental data and model prediction, and therefore model validity.

A comparison of the rate parameters k min (gastric emptying) and k abs (intestinal absorption) off and on pramlintide indicated a significantly slower gastric emptying and intestinal absorption rates on pramlintide (P=0.003 and P=0.015, respectively) (Fig. 1, lower panel). Comparing parameters b and c shows a significantly increased quantity of glucose retained by the stomach on pramlintide (P=0.0002 and P=0.0001), respectively.
Experiment 1: no adjustment of insulin for pramlintide effect
Figure 3, upper panel, presents the average simulated profiles of plasma glucose concentration and glucose Ra for placebo and all tested P/I ratios. The traces indicate that a P/I ratio of 3 was not better than insulin alone (placebo)—it delayed slightly the postprandial blood glucose peak but did not attenuate its magnitude. For P/I ratios of 6, 8, 9, 10, and 12 μg/U, hypoglycemic events occurred in 9%, 10%, 11%, 12%, and 15% of the in silico subjects; thus subjects' ICRs had to be adjusted to compensate for the effects of pramlintide.

Experiment 2: individual ICR adjusted for pramlintide effects
Figure 3, lower panel, presents average simulated profiles of plasma glucose concentrations and meal Ra, after ICR adjustment. As seen in Figure 3, lower-left panel, hypoglycemia, typically occurring 2 h after premeal insulin, was significantly attenuated, which allowed further evaluation of the relative effectiveness of the various P/I ratios in Experiment 2. We found that: • A P/I ratio of 3 was no more effective than placebo, with 62% versus 71% in CVGA Zone A. • P/I ratios of 6, 8, 9, 10, and 12 μg/U resulted in significant improvements of postprandial glucose control: 81%, 89%, 90%, 89%, and 88% of subjects in CVGA Zone A, respectively, and no hypoglycemia, as shown in Figure 4.

Control variability-grid analysis for placebo and pramlintide-to-insulin (P/I) ratios of 6, 8, 9, 10, and 12 μg/unit. Maximum and minimum blood glucose (BG) units are in mg/dL.
During Experiment 2, the ICR ratio was increased on average by a mean of 26±2%, which resulted in an average reduction in insulin dosing of 20±2% for P/I ratios of 6–12. In Figure 5, upper panel, the total average dose expressed in micrograms of pramlintide administered with the meal and the average percentage reduction in insulin bolus needed to minimize hypoglycemia at different P/I ratios are reported for the different P/I ratios.

The CVGA Zone Quotient is shown in Figure 5, lower panel, for placebo, and P/I ratios ranging from 3 to 12 μg/U. The Zone Quotient peaked at a P/I ratio of 9 μg/U, indicating that this would be an optimal mixture of pramlintide and insulin in terms of safety and efficacy.
Discussion
Despite advances in insulin formulations and improvements in insulin delivery and glucose monitoring devices, as well as robust data demonstrating the importance of good glycemic control, multiple recent clinical studies have shown that the majority of patients with type 1 diabetes are not achieving A1C goals with insulin alone; this holds true for observational studies, randomized controlled trials, and registries. 29 –33 Given that we are now over 20 years past the publication of the landmark results from the Diabetes Control and Complications Trial, 34 these data are a powerful indicator that insulin alone is not adequate to achieve the desired level of control in the general population. Major obstacles to achieving this goal include the fear of hypoglycemia and weight gain, which limit the willingness of patients and/or their providers to use more insulin. Furthermore, for approximately 20% of patients who do achieve goal, they do so by maintaining a constant awareness of the blood glucose levels, analyzing the impact of what they have done recently, what they plan to eat, and what they anticipate doing during the next 12-h period. Patients who are not able to practice this level of diligence, for whatever reason, are unlikely to succeed. 29
The difficulty in controlling glucose concentration with insulin alone is somewhat expected because glucose regulation in health is achieved by complex feedback mechanisms involving several other hormones (e.g., counterregulatory hormones and amylin). Amylin, a second peptide hormone secreted with insulin by pancreatic β-cells, suppresses glucagon secretion during the postprandial period, 16 slows gastric emptying, 35 which decreases the rate of appearance of ingested nutrients in the peripheral circulation, and reduces food intake by a centrally mediated mechanism. 19 These effects, which work to improve glucose control during the postprandial period, resulting in a reduction in the dose of short-acting mealtime insulin required, are restored by preprandial administration of pramlintide, an analog of human amylin. Pramlintide was approved by the FDA as an adjunct to insulin to improve glycemic control in patients who use mealtime insulin. However, in its current formulation, pramlintide must be administered as a separate injection, increasing the patient's injection burden, and can only be given as a fixed dose regardless the insulin dose. Because amylin and insulin are secreted in a fixed molar ratio by the healthy human pancreas, 14 the potential to co-formulate pramlintide with short-acting insulin is attractive as it may come closer to mimicking normal physiology and it would reduce the number of injections required. A fixed-ratio product could also potentially be administered by a subcutaneous infusion similar to today's insulin pumps.
Designing a clinical study that would identify an appropriate fixed ratio for pramlintide and insulin would be prohibitively expensive because countless dose ratios would need to be tested. A cost-efficient alternative is in silico experimentation—an increasing trend in drug design, typically including bioinformatics techniques that allow rapid selection of drug candidates based on high-throughput search in massive databases containing molecules, molecular fragments, or molecular structures. 1 In this article the meaning of in silico drug optimization is different: a model of the human metabolic system is used to study the co-dynamics of a combination of compounds with separately known kinetic properties. In that sense, the approach presented here is new and only possible because of recent advances in the in silico modeling of insulin–glucose dynamics in diabetes. 3,10 The availability of data from sophisticated tracer studies 26,27 allowed the reconstruction of the effects of pramlintide on the glucose rate of appearance following standardized meals. In turn, this allowed for dynamic modeling of the action of pramlintide on the gastrointestinal function and for the formal description of its dynamics via differential equations. A key element of this description is the assembly of an in silico population—a comprehensive set of virtual individuals with widespread characteristics (individual parameters) approximating the variability of pramlintide action observed in vivo.
Thus, deriving from previous technology developments and recent studies of pramlintide dynamics, we have updated our simulator of type 1 diabetes to incorporate the possibility of “administering” pramlintide to our virtual “subjects.” Using this technology, extensive in silico experiments have been performed to optimize the co-administration of pramlintide and insulin prior to its submission to clinical trials. The optimal P/I ratio was estimated to be 9 μg of pramlintide/U of insulin. Additional simulations narrowing the investigated range indicated that P/I ratios of 8 and 10 μg/U would achieve similar performance.
During the simulations, a well-known clinical observation was reproduced: to avoid postprandial hypoglycemia, insulin doses needed to be reduced by approximately 21% at the P/I ratio of 9 μg/U, which is consistent with the reduction of insulin during the titration of pramlintide in clinical practice. 21,22 This in silico effect was not preprogrammed—it was a result from the simulation experiments, further confirming the validity of our model.
Although modeling can certainly provide insight into individualization of the P/I ratio and pramlintide doses, the purpose of this study was to suggest a generally optimal ratio, usable in a future co-formulation of the two compounds to be delivered as a single mixture. We can therefore conclude that a valid simulation model of the action of pramlintide was developed, leading to in silico estimation of optimal pramlintide–insulin co-administration ratio. Clinical trials are needed to confirm (or adjust) this initial estimation.
Footnotes
Acknowledgments
This work was supported by grant 15-2011-45 (“Insulin-Pramlintide Combination for the Treatment of Type I Diabetes”) from the JDRF to Amylin.
Author Disclosure Statement
O.K., E.C., and K.H are employees of Bristol-Myers Squibb. F.M., C.D.M., B.K., and C.C. hold a pending patent on the simulation model. J.S. declares no competing financial interests exist.
