AJP - Endo Journal of Applied Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Endocrinol Metab 293: E396-E409, 2007. First published April 24, 2007; doi:10.1152/ajpendo.00647.2006
0193-1849/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
293/1/E396    most recent
00647.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bertuzzi, A.
Right arrow Articles by Mingrone, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bertuzzi, A.
Right arrow Articles by Mingrone, G.

Insulin granule trafficking in beta-cells: mathematical model of glucose-induced insulin secretion

Alessandro Bertuzzi,1 Serenella Salinari,2 and Geltrude Mingrone3

1Institute of Systems Analysis and Computer Science, Consiglio Nazionale delle Ricerche, Rome; 2Department of Systems Analysis and Informatics, University of Rome "La Sapienza," Rome; and 3Institute of Internal Medicine, Catholic University School of Medicine, Rome, Italy

Submitted 27 November 2006 ; accepted in final form 16 April 2007


    ABSTRACT
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
A mathematical model that represents the dynamics of intracellular insulin granules in beta-cells is proposed. Granule translocation and exocytosis are controlled by signals assumed to be essentially related to ATP-to-ADP ratio and cytosolic Ca2+ concentration. The model provides an interpretation of the roles of the triggering and amplifying pathways of glucose-stimulated insulin secretion. Values of most of the model parameters were inferred from available experimental data. The numerical simulations represent a variety of experimental conditions, such as the stimulation by high K+ and by different time courses of extracellular glucose, and the predicted responses agree with published experimental data. Model capacity to represent data measured in a hyperglycemic clamp was also tested. Model parameter changes that may reflect alterations of beta-cell function present in type 2 diabetes are investigated, and the action of pharmacological agents that bind to sulfonylurea receptors is simulated.

beta-cell; insulin secretion; insulin granule dynamics; type 2 diabetes


IN THE PAST FEW YEARS, our knowledge of the events involved in the insulin secretion by pancreatic beta-cells has advanced considerably. New techniques emerged, such as the use of fluorescent proteins that can be targeted to secretory granules, the real-time imaging of granule trafficking in living cells, and the cell capacitance measurements that allow monitoring of the changes in the cell surface area that result from the fusion of granules with the cell membrane (41).

Since the early observations of the first and second phase of insulin secretion in response to a step in extracellular glucose (11, 13), pathways of stimulus-secretion coupling have been identified. In the triggering pathway, the glucose stimulus causes an increase in the ATP-to-ADP ratio. This leads to closure of ATP-dependent K+ (KATP) channels, cell membrane depolarization, and Ca2+ influx through voltage-sensitive Ca2+ channels. The consequent increase in cytosolic Ca2+ concentration ([Ca2+]i) stimulates exocytosis of insulin granules. Still not completely characterized are other signaling pathways, and several mediators may be involved in the KATP-independent Ca2+-dependent (amplifying) pathway (43, 21, 31, 30). ATP and Ca2+ cover a major role in the machinery that regulates granule translocation from the trans-Golgi network to the cell membrane, granule priming, fusion with the cell membrane, and release of contents in the extracellular space (41). Other important aspects of islet behavior are as follows: 1) the heterogeneity in the responsiveness to glucose of the beta-cell population (26) and 2) the presence of spontaneous sustained oscillations in the concentration of signaling molecules and in the insulin release, revealed both in isolated cells and islets, and in animal models as well as in human patients (19, 40).

Complex mathematical models of glucose-stimulated insulin secretion were proposed previously by Grodsky and coworkers (20, 37) (storage-limited model) and by Cerasi and coworkers (12, 36) (signal-limited model). Subsequently, models of insulin secretion were proposed that provide estimates of beta-cell function by using C-peptide data (49, 8, 33) and that are more easily utilizable in clinical practice. Complex models of the Hodgkin-Huxley type that represent the bursting behavior of the beta-cell electrical activity and the ATP and Ca2+ oscillations were also presented (46, 4, 38).

The mathematical model proposed in the present paper focuses on the dynamics of formation, translocation to cell membrane, and exocytosis of insulin granules in beta-cells. We have not attempted to model in detail the biochemical events that lead to changes in ATP concentration ([ATP]) and cytoplasmic Ca2+ concentration (or other signaling molecules) in response to stimulation by extracellular glucose or to the experimental manipulations that cause cell membrane depolarization. Instead, these processes are represented as time changes of rate coefficients regarded as "control" signals that regulate granule trafficking. The model provides a preliminary framework in which a number of experimental observations can be accommodated and analyzed in quantitative terms. The capacity of the model to analyze experimental data obtained in a hyperglycemic clamp was ascertained. In addition, to our knowledge, this is the first model that tries to represent the pharmacological action of drugs that stimulate insulin secretion.


    Glossary
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
State Variables (Equation Number Where Quantities First Appear)

I
Pool of proinsulin aggregates (Eq. 1)

V
Pool of granule membrane material (Eq. 1)

R
Reserve pool (Eq. 3)

D
Pool of docked granules (Eq. 4)

DIR
Pool of immediately releasable granules (Eq. 4)

F
Pool of granules fused with plasma membrane (Eq. 6)

{gamma}
Rate coefficient of granule externalization and priming, related to the ATP-to-ADP ratio (min–1) (Eq. 3)

{rho}
Rate coefficient of granule fusion with cell membrane, related to [Ca2+]i (min–1) (Eq. 5)

Additional Variables

C
Pool of unbound L-type Ca2+ channels

G
Extracellular glucose concentration (input variable) (mmol/l) (Eq. 7)

ISR
Insulin secretion rate (pmol/min) (Eq. 11)

IS
Insulin amount secreted in a time T during the first-phase response (pmol) (Eq. A4)

Gm
Measured plasma glucose concentration (mmol/l) (Eq. C1)

Ip
Measured plasma insulin concentration (pmol/l) (Eq. C1)

Ga
Estimated arterial glucose concentration (mmol/l) (Eq. C1)

Parameters in the Equations of Granule Dynamics

k
Rate constant of formation of proinsulin-containing granules in trans-Golgi network (min–1) (Eq. 1)

bI
Biosynthesis rate of proinsulin aggregates at a given glucose concentration (min–1) (Eq. 1)

{alpha}I
Rate constant of degradation of proinsulin aggregates (min–1) (Eq. 1)

bV
Rate of biosynthesis of granule membrane material (min–1) (Eq. 2)

{alpha}V
Rate constant of degradation of granule membrane material (min–1) (Eq. 2)

{tau}V
Time delay related to recycling of granule membrane material (min) (Eq. 2)

CT
Pool of total Ca2+ channels (Eq. 4)

k1+
Rate constant of association for the binding between granule and Ca2+ channel (min–1) (Eq. 4)

k1
Rate constant of dissociation for the binding between granule and Ca2+ channel (min–1) (Eq. 4)

{sigma}
Rate constant of insulin release from granules fused with cell membrane (min–1) (Eq. 2)

Parameters and Functions in the Equations Representing Stimulus-secretion Coupling

{eta}
Rate constant in the equation for {gamma} (min–1) (Eq. 7)

{gamma}b
Basal value of {gamma} (min–1) (Eq. 7)

{psi}
Oscillatory function that represents events inducing oscillations in {gamma} (min–1) (Eq. 7)

h{gamma}
Function of G representing the activatory action of glucose on {gamma} (min–1) (Eq. 7)

{tau}G
Time delay related to time required by glucose metabolism for activation of {gamma} (min) (Eq. 7)

G*
Glucose concentration threshold for the activation of {gamma} (mmol/l) (Eq. 8)

h
Maximal value of h{gamma} (min–1) (Eq. 8)

G
Glucose concentration over which h{gamma} remains constant and equal to h (mmol/l) (Eq. 8)

{zeta}
Rate constant in the equation for {rho} (min–1) (Eq. 9)

{rho}b
Basal value of {rho} (min–1) (Eq. 9)

h{rho}
Function of {gamma} representing the activatory action of {gamma} on {rho} (min–1) (Eq. 9)

k{rho}
Parameter representing the sensitivity of {rho} on the activatory action of {gamma} (Eq. 10)

Additional Parameters and Functions

I0
Insulin amount contained in a granule (amol) (Eq. 11)

Nc
Average number of beta-cells in an islet (Eq. 11)

Ni
Number of islets in pancreas (Eq. 11)

N
Total number of beta-cells in pancreas (Eq. 17)

f
Fraction of beta-cells responding to glucose stimulus (function of glucose concentration G) (Eq. 11)

fb
Basal value of the fraction f (Eq. A5)

Kf
Parameter in the function f(G) (mmol/l) (Eq. A5)

VG
Glucose distribution volume (liters) (Eq. C1)

VC
Volume of C-peptide accessible compartment (liters) (Eq. 17)

Q
Cardiac output (l/min) (Eq. C1)

SI
Insulin sensitivity index [l·min–1(pmol/l)–1] (Eq. C1)

Sd
Dynamic responsivity index (Eq. 17)

Ss
Static responsivity index (min–1) (Eq. 18)

a and b
Parameters in the approximate expression of the first-phase ISR (min–1) (Eq. A1)

xb
Denotes the basal value of the generic variable x

x
Denotes the steady-state value of the variable x


    DEVELOPMENT OF THE MATHEMATICAL MODEL
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
Granule Dynamics in a beta-Cell

Let us consider for the moment a single beta-cell. A schematic diagram of the formation and trafficking of insulin-containing granules is given in Fig. 1. Essentially, the model describes the kinetics of four intracellular pools of insulin granules: the reserve pool, the pool of docked granules, the pool of immediately releasable granules, and the pool of granules fused with cell membrane. In the following, the model of Fig. 1 will be illustrated and formulated as a mathematical model.


Figure 1
View larger version (5K):
[in this window]
[in a new window]

 
Fig. 1. Diagram of the kinetics of granule pools in the beta-cell. I, pool of proinsulin aggregates; V, pool of granule membrane material; R, reserve pool; D, pool of docked granules; C, pool of Ca2+ channels; DIR, pool of immediately releasable granules; F, pool of granules fused with plasma membrane. Other symbols explained in the text. The quantity {sigma}F represents the number of insulin granules that, in the unit time, release insulin in the extracellular space.

 
Let us denote by R(t) the number of insulin granules in the reserve pool at time t. Although it is still unclear how (pro)insulin aggregates are sorted and segregated into granule membranes (35), we assume that the rate of formation of insulin-containing granules in the trans-Golgi network may be expressed as kI(t)V(t), where k is a segregation rate constant, I(t) is the pool of "free" (i.e., not yet segregated into granule membranes) proinsulin aggregates, and V(t) is the pool of available free granule membrane material, not yet enclosing proinsulin, at time t. A conservation equation for I is written as follows:

Formula 1(1)
where the overdot indicates derivative with respect to time, {alpha}I is a degradation rate constant, and bI is a production rate that describes proinsulin biosynthesis and aggregation. Changes in proinsulin biosynthesis rate do not appear to be involved in the first and second phase of insulin secretion elicited by a step increase in ambient glucose concentration (42, 22). However, proinsulin synthesis has been reported to increase with glucose concentration (45), so bI is actually a function of the extracellular glucose concentration, G. The governing equation for V is

Formula 2(2)
where {alpha}V is a degradation rate constant and bV is the rate of production of the granule membranes. The last term in the right-hand side of Eq. 2 accounts, in simplified way, for the contribution to granule membrane formation resulting from the (complete) recycling of membrane material. This material becomes part of the plasma membrane at exocytosis and is successively removed and returned to the trans-Golgi network. We have denoted here by F the number of granules that are fused with cell membrane, by {sigma} the rate constant of the fusion process, and by {tau}V the constant time interval required for recycling. Incomplete recycling may easily be accounted for by means of a multiplicative coefficient smaller than 1.

To write the governing equation for the reserve pool, R, we assume that granules formed at the trans-Golgi network immediately enter the reserve pool and that granules leave this pool to become docked granules with a rate constant {gamma}. If granules are substantially stable, that is, they cannot be disaggregated or degraded, we can write:

Formula 3(3)

The mechanism by which granules are translocated from trans-Golgi to plasma membrane in beta-cells has not yet been completely elucidated, although it appears to involve unidirectional movement along microtubules by means of motor proteins (52). Experimental data from mouse and rat beta-cells show that only 1–2 granules/min per beta-cell are released at extracellular glucose concentration around 3 mmol/l, and around 20–30 granules/min per cell are released at the peak of the first phase of insulin secretion after a step increase in glucose (2, 7, 48). Thus granule transport must involve only a small fraction of the reserve pool that has been estimated to contain 9,000–13,000 granules (48). To keep the model simple, we include in the translocation also granule priming, so a unique rate coefficient, {gamma}, will represent the processes that produce docked granules that are competent for release. As discussed in the next section, {gamma} has to be considered a rate coefficient that is actually dependent on the glucose concentration G and represents a main rate-limiting step in the response to glucose.

According to a well-established view (2, 48), within the pool of docked and primed granules here denoted by D, a pool of immediately releasable granules, DIR, must be distinguished. The latter granules are likely to be tightly associated with L-type Ca2+ channels (2). Thus we assume (see Fig. 1) that the pool DIR is formed through the binding of docked granules with unbound L-type Ca2+ channels, C. For simplicity, we will also assume that the Ca2+ channels that become unbound at granule fusion are immediately available to bind new granules. Denoting by CT the constant pool of total Ca2+ channels (bound plus unbound), we have C = CT – DIR and we may write:

Formula 4(4)

Formula 5(5)
where k1+ and k1 are the rate constants of association and dissociation, respectively, for the binding between docked granule and Ca2+ channel, and {rho} is a rate coefficient that accounts for the factors that promote the fusion of granules with cell membrane. Actually, granule fusion is strictly related to changes in [Ca2+]i. As seen in the following sections, we consider {rho} a rate coefficient that is affected by changes in the ATP-to-ADP ratio or by experimental manipulations, such as the increase in extracellular K+ concentration ([K+]), that cause depolarization of the cell membrane and thus an increase in [Ca2+]i.

The rate of granule membrane fusion with cell membrane is given by the last term in the right-hand side of Eq. 5, i.e., {rho}DIR(t), so the equation for the pool of granules undergoing fusion with the cell membrane, F(t), is given by

Formula 6(6)
where {sigma} is the rate constant that regulates the duration of the fusion event. Although insulin cargo release follows the formation of the fusion pore with some delay (41), for simplicity, we lump the two events together by considering a unique rate constant. So {sigma}F(t) is the number of granules releasing insulin in the unit time at time t in a single beta-cell. When granules have completed the fusion process and released their content in the extracellular space, the material that forms granule membrane is recycled to the trans-Golgi network. This flow of internalized material, delayed by a time {tau}V, appears thus in Eq. 2.

Stimulation by Changes in Extracellular Glucose

The above model of granule dynamics must be complemented by equations that relate the glucose stimulus to the quantities that govern the machinery of granule trafficking. This has been done on the basis of the scheme in Fig. 2, modified from Fig. 1 in Ref. 21. The ATP-to-ADP ratio, and the other glucose-dependent factors that enhance granule supply from the reserve pool to the docked pool (21, 31, 30), are assumed to be globally represented by the phenomenological parameter {gamma} that appears in Eqs. 34. [Ca2+]i is represented by the rate constant {rho} in Eqs. 56. The fast oscillations in [Ca2+], related to changes in cell membrane potential, have been neglected in the present study. More comprehensive models of ATP and Ca2+ kinetics have been proposed (4, 38, 5).


Figure 2
View larger version (15K):
[in this window]
[in a new window]

 
Fig. 2. Schematic representation of the triggering and amplifying pathways that couple glucose stimulus to insulin secretion. Solid arrows, material fluxes; dotted arrows, control signals; +, stimulation; –, inhibition. The rate constants {gamma} and {rho} are the main signals that regulate granule trafficking and insulin secretion. Modified from Ref. 21.

 
Experimental evidence suggests that the glucose-induced increase in cytosolic [ATP] stimulates granule movement along microtubules by enhancing kinesin activity (52). Granules are recruited from the reserve pool to the plasma membrane and, in fact, sustained insulin release requires kinesin-dependent granule transport. Granule priming also appears to be accelerated by [ATP] increase (2, 41). Ca2+ contributes to granule translocation because kinesin motor activity is increased by the Ca2+-dependent dephosphorylation of kinesin heavy chain (16), and R-type Ca2+ channels are implicated in granule mobilization and priming (25). However, for simplicity, we will not include in the model this regulatory role of Ca2+ on granule transport. Moreover, it has been found that [ATP] undergoes oscillations, possibly originating in the glycolytic pathway, with a period of 5–10 min (19). By ignoring the actual molecular mechanisms underlying the glucose-induced supply of granules to the docked pool, we assume for {gamma} a first-order kinetics. In writing this equation, we account for: 1) the existence of a basal value of [ATP], and thus of {gamma}, 2) the possible presence of spontaneous oscillations, and 3) the ATP response elicited by the glucose stimulus. Thus we have:

Formula 7(7)
where {eta} is a rate constant, {gamma}b is the basal value at low glucose, and {psi} is an oscillatory forcing function that represents the events inducing [ATP] oscillations. Glucose activation is modeled by the function h{gamma}, and the time delay {tau}G is the time required by glucose metabolism (41, 48). Data on isolated beta-cells and cell clusters stimulated by glucose (26) show that cell response is maximal and almost constant for G >10–12 mmol/l and much lower below 6 mmol/l. So h{gamma} has been chosen equal to zero for G smaller than a threshold G*, linearly increasing for G in the glucose concentration range (G*, G), and then constant. We have:

Formula 8(8)
where h is the maximal value of h{gamma} and is reached at G = G.

The increase in the ATP-to-ADP ratio leads to the closure of KATP channels and increase in [Ca2+]i (see Fig. 2). So, the rate coefficient {rho} will depend on {gamma}. In writing the governing equation for {rho}, 1) we account for the existence of a basal value {rho}b, and 2) we assume linearity to model the action of [ATP] on [Ca2+]i. We obtain the following equation:

Formula 9(9)
where {zeta} is a rate constant, and the function h{rho}, which embodies the action of ATP on Ca2+, is simply defined as

Formula 10(10)
We remark that, for G below the threshold G* and {psi} = 0, the rate coefficients {gamma} and {rho} stay at their basal values, and the ISR may be thought to be determined by the "constitutive pathway" of insulin secretion. When G exceeds the threshold, h{gamma} is no longer equal to zero, and the ISR becomes determined by the "regulated pathway."

In summary, as shown in Fig. 2, extracellular glucose modulates the mitochondrial signal {gamma} according to Eqs. 78, and {gamma} regulates granule translocation and priming. In turn, {gamma} affects {rho}, and thus the fusion of DIR granules with cell membrane, according to Eqs. 910. So, in the present model, it is the parameter {gamma} that has a major role both in the triggering pathway and in the amplifying pathway (43, 21). Equations 710 represent in a very rough way the coupling of glucose metabolism to the signals that activate insulin secretion. In particular, we neglected the dependence of production and degradation rates of ATP on [Ca2+]i, as well as the kinetics of Ca2+ stores in the endoplasmic reticulum, which may cause oscillations in the concentrations with a phase shift between {gamma} and {rho} (1, 15).

The beta-Cell Population

According to Eq. 6, the insulin amount released in the unit time by a single beta-cell in the extracellular space is given by I0{sigma}F, where I0 is the insulin amount contained in a granule, if complete emptying of granules that undergo fusion is assumed. Actually, by targeting fluorescent probes to granule membrane or granule lumen (50), it has been found that membrane fusion is not always associated with insulin release. Thus only a fraction (15–40%) of granule insulin amount, with this fraction being possibly modulated by the stimulus intensity, may be released.

To compute the insulin amount released in the unit time by the whole beta-cell population (ISR), we must account for both the size of the population (total number of beta-cells) and the fraction of cells that actually release insulin in the given experimental situation. Denoting by Nc the average number of beta-cells per islet and by Ni the number of islets, the total number N of beta-cells may be written as N = Nc x Ni. In the experiments with isolated islets, if a single islet is considered, we have Ni = 1.

When the intact pancreas in vivo is considered, the number of active beta-cells is regulated by several factors. The beta-cell number is known to undergo time changes in relation to the balance among cell replication, cell death by apoptosis, and neogenesis by differentiation from ductal epithelium. Regulation of these processes appears to be altered in obese and type 2 diabetic subjects (9, 14). Estimates of the total number N of pancreatic beta-cells have been reported for mouse (1.06 x 106) and rat (2.76 x 106) (6). In studies on human pancreata (14), the number of islet equivalents (IEq, islet with average size of 150 µm diameter) per kilogram of body weight was found to be 6,640 IEq/kg (decreased to 2,641 IEq/kg in type 2 diabetes).

Another important factor is the heterogeneous responsiveness of the cell population to glucose (neglecting other secretagogues). The fraction of responsive islet cells and clusters increases with the extracellular glucose (26), and cell subpopulations with different intrinsic sensitivity to glucose are likely to exist. In addition, because glucose diffuses in the islet volume and is consumed by beta-cells, even in isolated islets, the cells will be exposed to different glucose concentrations both in the basal state and during changes in glucose concentration (3). In spite of the rich vascularization of pancreas, extracellular glucose concentration will decrease as the distance from vessels increases, and a wide range of glucose concentrations is likely to be found. Moreover, it is known that insulin itself, acting as a paracrine signal, sensitizes beta-cells to glucose stimuli and may thus participate in cell recruitment (47).

To keep the model as simple as possible, we have adopted a rough description of beta-cell recruitment by glucose stimulus, by assuming that only a fraction, f, of the total cell population responds to glucose. This fraction will increase as the glucose concentration G increases, with a time delay for recruitment equal to {tau}G, so we have f(t) = f[G(t{tau}G)]. Under the further assumption that beta-cells respond to stimulation in a synchronous manner, which may be realistic on a not-too-fast time scale thanks to cell-to-cell electrical and chemical interactions (38), we may finally write for the insulin secretion rate of the whole cell population at time t the following expression:

Formula 11(11)

Given the time course of glucose concentration G, Eqs. 111 predict the time evolution of the control rate coefficients {gamma} and {rho}, the evolution of the quantities that describe granule dynamics (I, V, R, D, DIR, F), and the evolution of the recruited cell fraction f and hence of the ISR. The following parameters may be considered constant, independent of the degree of stimulation of the beta-cell: k, {alpha}I, {alpha}V,{tau}V, CT, k1+, k1, {sigma}, {tau}G, {eta},{gamma}b, G*, G, h, {zeta}, {rho}b, k{rho}, and I0 (parameters of granule dynamics and signaling model); Nc and Ni (parameters of the beta-cell population). By contrast, other quantities in the model, namely bI, bV, {psi}, and f, may be affected by changes in the concentration of glucose (and other secretagogues) or as a consequence of changes in the experimental conditions. An equation for f as a function of G is given in APPENDIX A, and bV will be taken time and glucose independent.

Steady-state Conditions

To represent the basal condition (glucose concentration Gb), or a situation in which glucose concentration is kept constant for a long time at a basal or suprabasal value, we may consider the steady-state solution of the model in which bI, {gamma}, {rho}, and f take values corresponding to the assigned value of G, {psi} = 0, and the pools I, V, R, D, DIR, and F are constant. These constant values are obtained by setting to zero the time derivatives in model equations. We observe that, in the steady state, we have the following chain of equalities:

Formula 12(12)
Taking Eq. 12 into account, from Eqs. 13 we see that I, V, and R are given in the steady state by the following expressions:

Formula 13(13)
Thus we have

Formula 14(14)
Using Eqs. 46, D, DIR, and F can be conveniently expressed in terms of the steady-state value of R:

Formula 15(15)
To guarantee the positivity of D in the steady state, the constraint {rho}CT{gamma}R > 0 has to be satisfied.

For the ISR in the steady state, Eq. 14 gives

Formula 15
showing that, as intuitive, the rate of insulin secretion in the long time is essentially determined by the proinsulin biosynthesis rate, bI. Because bI and f depend on G, the above equation is the "dose-response relationship" predicted by the model when glucose levels (and all other quantities) are kept constant for a long time. As we will see in the following, this relationship may not reflect the ISR data measured from an early part of the second phase of insulin secretion.

Assessment of Model Parameter Values

For some of the parameters of the model of Eqs. 111, we may use experimental measurements or estimates, e.g., I0 = 1.6 amol in rat granules (41) and Nc = 1,000–2,000 cells for an islet (2). The granule fusion event plus the insulin cargo release has been reported to require a time of the order of a few seconds (41), so we may take an estimate for {sigma} in the range 6–30 min–1. Moreover, CT = 500 was reported as an upper estimate for the number of Ca2+ channels (2). We will take {tau}G = 1 min (2).

Other indications about the values of model parameters can be obtained in view of the available experimental data concerning the basal state in isolated islets. At the glucose concentration of 3 mmol/l, only 1–2 granules·min–1·cell–1 are released in the mouse or rat (48). We will assume that, at the basal state, {sigma}F (number of granules released per minute by a beta-cell) and the other quantities equal to {sigma}F that appear in Eq. 12 are equal to 1. With this choice, assuming Nc = 103 active cells/islet and f = 1 (that is, 100% of islet cells release insulin), the ISR of a single islet predicted by Eq. 11 at the basal state is equal to 1.6 x 10–3 pmol/min (or 9.28 pg/min), which is in the range of the values experimentally observed (22, 44, 47).

To determine the baseline values of model parameters, we proceed by taking {sigma}F = 1. The pool of immediately releasable granules, DIR, is reported to contain ~50 granules in the basal state (2, 7, 48). By taking DIR = 50, from Eq. 12 we estimate that the value of {rho} in the basal state, denoted as {rho}b, is {rho}b = 0.02 min–1. A plausible estimate is that the reserve pool, R, contains 10,000 granules, and that 1,000 granules are docked with the plasma membrane: thus we take D + DIR =1,000. The chosen R value leads to the estimate for {gamma} in the basal state (denoted by {gamma}b) as {gamma}b =10–4 min–1. To proceed further, we take as plausible the value of 10 for both I and V at the basal state (it may be seen that this assumption does not affect the general pattern of the first and second phase of insulin secretion). So kIV = 1 (from Eq. 12) implies k = 10–2 min–1. Also, from Eq. 13 we have bIkbV/(kbV + {alpha}I{alpha}V) = 1. The quantity kbV/(kbV + {alpha}I{alpha}V) ranges from 0 to 1 and increases with bV, so we may look at it as a factor that takes into account the variability in the vesicle membrane production rate related, for instance, to variability in lipid metabolism within the beta-cell. Taking for this factor the value 0.25 at the basal state, we deduce for the basal bI, denoted as bI,b, the estimate bI,b = 4 min–1 and, using Eq. 1 written in the basal state, we also obtain {alpha}I = 0.3 min–1. Let us now assume that the degradation rate constant for V is larger than the degradation rate constant for I, e.g., {alpha}V = 2{alpha}I, and use a similar reasoning as above. This leads to the following estimates: {alpha}V = 0.6 min–1 and bV,b = 6 min–1.

Concerning {tau}V, k1+, k1, and CT, we have tentatively assumed, on the basis of simulation results, the following baseline values: {tau}V = 5 min, k1+ = 1.447 x 10–5 min–1, k1 = 0.10375 min–1, CT = 500. With these parameter values, the (basal) steady-state values of model variables are thus as follows: I = 10, V = 10, R = 10,000, D = 950, DIR = 50, F = 3.33 x 10–2, and ISR = 9.28 pg·min–1·islet–1 (f = 1). The values of the other parameters were chosen as discussed in APPENDIX A.


    RESULTS
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
Stimulation by High K+

Data of islet response to stimulation by high extracellular [K+] in low glucose and in the presence of diazoxide have been given by various authors (see Ref. 21). Henquin et al. (22) report time courses in mouse islets of [Ca2+]i and of the insulin secretion rate in response to a step increase in extracellular K+ and to a sequence of impulsive K+ increments with pulses lasting 6 min and interpulse intervals of 6 min. The continuous increase in K+ caused a sustained [Ca2+]i rise, whereas the insulin secretion rate exhibited a large initial peak followed by a slow return to the initial value. The impulsive increments in K+ caused an intermittent rise of [Ca2+]i and a sequence of decreasing peaks of the insulin secretion rate. The pattern of the response, but not its basal value, was influenced by the glucose level (0 or 3 mmol/l).

This type of stimulation is modeled here by assuming that the time course of {rho} reflects the time course of [Ca2+]i reported in Ref. 22, whereas {gamma} retains its basal value {gamma}b. Thus, starting from the basal value {rho}b, the time course of {rho} in the case of continuous stimulation has been obtained as the response of a linear first-order system with rate constant {zeta} to a step function (the step increase in K+). A component increasing linearly with slope s{rho} has been added to mimic more closely the measured [Ca2+]i in the time interval 0–50 min. So, if the step increase in K+ occurs at t = t1, we have {rho}(t) = {rho}b for t < t1 and

Formula 16(16)
see the plot in Fig. 3A. The values of Formula 16 and s{rho} were chosen on the basis of simulation results, for a K+ increase from 4.8 to 30 mmol/l (22).


Figure 3
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 3. Simulated islet response to K+ stimulation. A (continuous stimulation) and C (intermittent stimulation): time course of {rho} (dotted line), D (dashed-dotted line), and DIR (solid line). B and D: insulin secretion rate (ISR) time course with baseline parameters (solid line) and alternative parameters (dotted line). Baseline parameters: k = 10–2 min–1, bI,b = 4 min–1, {alpha}I = 0.3 min–1, {alpha}V = 0.6 min–1, bV,b = 6 min–1, {gamma}b = 10–4 min–1, {rho}b = 0.02 min–1, Formula 1 = 1 min–1, s{rho} = 4 x 10–3 min–2, {zeta} = 1 min–1, {tau}V = 5 min, k

Formula 1= 1.447 x 10–5 min–1, kFormula 1= 0.10375 min–1, CT = 500, {sigma} = 30 min–1, I0 = 1.6 amol, Nc = 1,000, and f = 1. Alternative values: kFormula 1= 5.788 x 10–5 min–1, kFormula 1= 0.255 min–1, CT = 300, and other parameters as in baseline.

 
The function {rho}(t) is in this case the input signal for the model of Eqs. 16. The model response to this stimulation (taking 1,000 beta-cells in an islet and f = 1) is shown in Fig. 3, to be compared with data reported in Ref. 22. Figure 3A gives the time course of pools D and DIR. The latter pool is seen to be rapidly depleted upon stimulation. The response of I and V (not reported) shows that these pools are affected after the time lag, {tau}V, required for the recycling from the plasma membrane to the pool V, which is in fact transiently increased. The pool R remains virtually unchanged. The predicted ISR is shown in Fig. 3B (solid line). The number of granules released per beta-cell, given by the integral of {sigma}F over the time, is equal to 67 granules in the first 5 min of the response. According to this model, this type of stimulation does not affect the steady-state values of I, V, R, F, and ISR, which are recovered in the long time: only the steady-state values of D and DIR are decreased (see Eqs. 13 and 15). If a Ca2+-mediated enhancement of microtubule-based granule transport (16) were included in the model, high K+ stimulation would produce a component of granule translocation to plasma membrane with a decrease in R.

For the intermittent stimulation, Fig. 3C shows the time course of {rho} (a sequence of 4 pulses lasting 6 min with values of Formula 16, {zeta}, and s{rho} being the same as for the continuous stimulation). The predicted response is illustrated in Fig. 3, C and D. Note in Fig. 3C the refilling of pool DIR during the interpulse interval. When CT, kFormula 16, and kFormula 16 were set to alternative values that produce the same amount of granules in the various pools at the steady state, the pattern of the response was qualitatively similar. The ISR is reported in Fig. 3, B and D, by the dotted line (83 granules released in the first 5 min). The simulated ISR of the islet in Fig. 3D reproduces the declining amplitude of the peaks experimentally observed at zero glucose (22). The data show a slight broadening of the peaks, possibly related to the incomplete synchrony of beta-cells in an islet.

Stimulation by Extracellular Glucose

As seen in the previous section, insulin secretion induced by K+ rise can be simulated in the model by an increase in the rate coefficient {rho} and shows a first-phase response only. In that case, the ISR pattern is monophasic, and, after the peak, the secretion rate declines toward the basal value. The secretory response to a step increase in extracellular glucose is, instead, biphasic. The second phase of glucose-stimulated insulin secretion had been initially considered to be sustained by proinsulin biosynthesis. Subsequently, this interpretation was refuted (18, 42). Confirmed by model simulations (data not shown), an increase in the rate of proinsulin biosynthesis bI, concomitant with the increase in {rho}, does not produce a secondary response if {gamma} is kept at {gamma}b. As previously discussed, a critical role is played by the increase in the rate coefficient {gamma} modulated by extracellular glucose, because {gamma} intervenes in both the triggering and the amplifying pathway (21, 43).

To compute the model response to glucose stimuli, the functions h{gamma}(G) and f(G) have to be specified. To this aim, we rely on data of insulin secretion by the perfused rat pancreas, reported by Grodsky (20). In particular, we used the data of insulin amount released during the first phase of the response to different values of a step change in glucose concentration (IS, see Fig. 4). An approximate expression of the first-phase ISR and of the insulin amount secreted, as predicted by the present model, is derived in APPENDIX A. APPENDIX A also gives the form chosen for the function f (Eq. A5) and describes the procedure followed to estimate the parameters of h{gamma}(G) and f(G) from Grodsky's data. Figure 4 shows the IS data (20) with the fitting curve of Eq. A4 (parameter values reported in APPENDIX A). The number of beta-cells used in the simulations in this and the following section is equal to the estimate reported for the rat.


Figure 4
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 4. Insulin amount secreted in the first-phase response to different values of a step glucose stimulation (IS). Filled squares, data (perfused rat pancreas, means ± SE) replotted from Ref. 20; continuous line, model fitting (parameter values reported in APPENDIX A).

 
Figure 5 shows examples of the model response to changes in glucose concentration. Figure 5A shows the simulated step change of G from 1 to 16.7 mmol/l (300 mg/dl in Grodsky's data) with the corresponding time course of D and DIR. Figure 5B reports the corresponding ISR (solid line), together with the ISR predicted with step in G equal to 150 and 500 mg/dl. In the 16.7 mmol/l glucose step, 570 granules per beta-cell are released in 50 min, to be compared with the value of 680 granules reported in Ref. 48 for rat pancreatic islets that exhibit a large second phase. The ISR predicted by the model exhibits both the first and the second phase of insulin secretion. However, the amplifying effect can be "deleted" in the model by letting {gamma} and {rho} evolve according to Eqs. 710 but imposing that {gamma} is kept to its basal value {gamma}b in the Eqs. 34 that govern granule translocation. In this case, the response only shows a first phase, followed by a decline in the ISR (Fig. 5B, dashed-dotted line). The difference between the solid and dashed-dotted line in Fig. 5B may thus be related to the amplifying action of glucose on insulin secretion. Figure 5, C and D, shows the response to two other patterns of glucose stimulation. The resensitization effect resulting from a previous prolonged stimulation is shown in Fig. 5C, and the model prediction agrees fairly well with experimental data (20). The ISR patterns obtained when G has a fast or a slow ramp rise are shown in Fig. 5D, where the delay in the response is because of the time required to reach G*. As also found experimentally (20, 37), the ramp stimulus greatly reduces or virtually abolishes the first phase of insulin secretion.


Figure 5
View larger version (16K):
[in this window]
[in a new window]

 
Fig. 5. Pancreatic response to step increase in glucose concentration. A: time course of G (from 1 to 16.7 mmol/l, dashed line), D (dotted line), and DIR (solid line). B: ISR time course for a step in G to 300 mg/dl (solid line), 500 mg/dl (dotted line), and 150 mg/dl (dashed line). Dashed-dotted line: time course of the ISR when the G step does not induce granule translocation and priming. C: ISR time course when glucose stimulation (16.7 mmol/l) is interrupted in a rest period from 60 to 65 min. D: ISR response to glucose rise as a fast ramp (slope 50 mg·dl–1·min–1, solid line) and a slow ramp (slope 5 mg·dl–1·min–1, dotted line). Baseline parameters (rat pancreas): {eta} = {zeta} = 4 min–1, {tau}G = 1 min, G* = 4.58 mmol/l, G = 10 mmol/l, h = 3.93 x 10–3 min–1, k{rho} = 350, fb = 0.05, Kf = 3.43 mmol/l, N = 2.76 x 106. Other parameters are same as the baseline parameters of Fig. 3.

 
A point to stress is that neither the increase in the proinsulin biosynthesis rate bI, which has been observed to follow glucose increase (45), nor the recycling of granule membrane material appreciably affect the second phase of the response in the short time. A large increase in bI, concomitant with the increase in G, does not modify appreciably the ISR response in a time interval of 1–2 h after the glucose step. By contrast, as expected from Eq. 14, a rise in bI modifies the ISR in the long term. Figure 6A reports the response to a step increase in G computed up to 12 h. After a phase of ~3 h in which the secondary response increases, the ISR tends slowly to the initial value with a marked depletion of the pool R. If, upon glucose stimulation, the biosynthesis rate is increased, the ISR tends to a larger value with a diminished depletion of the pool R, which is now refilled at a greater rate by proinsulin biosynthesis (see Eq. 13). The bI increase with G will be maintained in the simulations that follow. A large value of {tau}V, which slows down granule membrane recycling, tends to delay this effect.


Figure 6
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 6. Glucose stimulation. A: response to a step in G from 1 to 16.7 mmol/l, given at t = 10 min and lasting 12 h, bI kept at bI,b (ISR, solid line; R, dashed line), and bI increased to 8 x bI,b (ISR, dotted line; R, dashed-dotted line). B: response to a G step from 1 to 16.7 mmol/l, given at t = 1 h in the presence of the endogenous oscillation {psi}(t) in {gamma}. ISR, solid line; {rho}, dotted line. After t= 4 h, {rho} is set equal to 2 min–1 to simulate extracellular K+ rise. Other parameters are as in A (bI increased to 8 x bI,b).

 
The simulations of Fig. 6A show that the ISR measured at 1–3 h after the glucose step, which may be assumed to represent the stationary response and reported in a dose-response curve ISR vs. G, might be different from the actual stationary response of the ISR. In fact, because the long response time of the second phase related to the small value of {gamma} (see APPENDIX B), the stationary state may be attained only after a very long time.

Figure 6B shows the effect on {rho} and ISR of a sinusoidal oscillation possibly arising in the glycolytic pathway, represented by the function {psi} in Eq. 7. Because the oscillations of beta-cell signaling molecules appear to be enhanced by extracellular glucose (19), we tentatively assumed the amplitude of {psi} equal to (1/2)h{gamma}(G) (G changed stepwise from 1 to 16.7 mmol/l). The oscillation period was set to 10 min. When G exceeds the threshold G*, the oscillation becomes clearly visible in both {rho} and ISR. Depending on the amplitude, frequency, and site (oscillations may also be present in the fraction f of active beta-cells), the oscillatory component may represent a substantial or a major part of the ISR (40). In Fig. 6B, we also simulate the effect of the experimental manipulation consisting in the increase of extracellular K+, in the presence of diazoxide, when glucose is at a stimulatory concentration (19). As in the preceding section, K+ rise is modeled by an increase in {rho} (we set {rho} = 2 at t = 4 h). When {rho} is kept at a high value, the second-phase response is preserved since {gamma} remains at a high level and oscillates, but the oscillations in the ISR disappear. This behavior is consistent with the observations (19).

Simulation of Derangements of Physiological Mechanisms and Effect of Drugs

As seen in APPENDIX A, the product k{rho}h, but not the values of the two factors, was estimated from Grodsky's data. Figure 7A shows that, when h is increased and k{rho} decreased with k{rho}h kept constant, the first-phase response to a step increase in G is unchanged, whereas the second phase has a larger amplitude and a more prompt decay. Pool R is more completely depleted as h increases. This behavior agrees with the analysis in APPENDIX B and may be seen as the result of an increase in the sensitivity of second-phase ISR with respect to glucose. The pattern of the response at large h, where the second phase may largely exceed the peak of the first phase, is more reminiscent of the ISR observed in the rat than of that observed in the mouse (22). For a given G step, a larger h causes a larger rise in {gamma}, and this agrees with the idea that the amplifying action may be stronger in the rat than in the mouse (22). Figure 7A also shows the ISR obtained with the alternative values of CT, k1+, and k1, given in the legend of Fig. 3. Because the product k1+CT is larger with the alternative parameters than with baseline parameters, the second-phase ISR has a faster increase, as predicted by the initial value of the derivative (see APPENDIX B).


Figure 7
View larger version (16K):
[in this window]
[in a new window]

 
Fig. 7. ISR responses after glucose step from 1 to 16.7 mmol/l (bI increased to 8 x bI,b). A: h changed with k{rho}h kept constant and {psi} = 0. Solid line, baseline parameters; dotted line, h x 5; dashed line, h/5; dashed-dotted line, baseline parameters with alternative values of CT, k1+, and k1. B: solid line, ISR response to the step in G given at 1.5 h and {psi}(t) as in Fig. 6B, with h/2 min–1, k{rho}/5 min–1 and {sigma} = 3 min–1. After t = 4 h, to simulate the effect of a drug that blocks ATP-dependent K+ (KATP) channels, k{rho} is exponentially increased from 1/5 of the baseline value to the baseline value (time constant of the exponential, 10 min). Dotted-line plot, ISR response to the step in G given at 1 h, with bV = 1.5 bV,b and the other parameters as in Fig. 6B.

 
In the simulations of Fig. 7A, the parameter changes modified only the second phase of the response. Figure 7B shows the effect of changes that affect both the first and the second phase. The dotted-line plot reports the ISR response predicted with bV = 1.5 bV,b, which may represent an increased synthesis rate of granule membranes. The basal values of granule pools, as well as the basal ISR and the whole ISR response, are increased. In contrast, the solid-line plot was obtained by assuming a decreased effectiveness of stimulus-secretion coupling. A decreased h may simulate an impaired effectiveness of processes that regulate ATP metabolism, such as a loss of function of glucokinase, and decrease in k{rho} simulates impairment of Ca2+ transport in the cell related, for instance, to a defective closure of KATP channels. Moreover, a smaller {sigma} may represent impairment of processes involved in the exocytotic machinery (41). With h, k{rho}, and {sigma} decreased with respect to the baseline values (Fig. 7B, solid-line plot up to t = 4 h), the glucose-stimulated ISR almost completely lacks the first phase, and both the second phase and the oscillation amplitude are markedly reduced. These parameter changes may reflect alterations in beta-cell function present in type 2 diabetes (41).

In the solid-line plot of Fig. 7B, we also simulate, within the extremely simplified representation of the stimulus-secretion coupling adopted in the present model, the action of pharmacological agents that bind to sulfonylurea receptors (27, 24) and cause closure of KATP channels, leading to increased conductance of voltage-sensitive Ca2+ channels and increased Ca2+ influx. An increase in the model parameter k{rho} translates into an increase in {rho}, that is, in the [Ca2+]i. Thus, the drug action may be simulated by an increase in k{rho}. The solid-line plot of Fig. 7B shows, up to t = 4 h, the response to a glucose step of a beta-cell with impaired stimulus-secretion coupling. At t = 4 h, to simulate the drug effect, k{rho} is exponentially increased to the baseline value. This change causes a rise in {rho} and thus a transient increase in the ISR, as it would occur with a rise in extracellular K+. The amplitude of the oscillations is restored.

In general, the effect of antidiabetic drugs should guarantee that the signal {rho} is large enough to allow exocytosis of the docked granules that are fed by the reserve pool at rate {gamma}R, as remarked in Ref. 21. If the quantity {rho}CT{gamma}R in Eq. 15 for D becomes small (and positive), the pool D may in fact become exceedingly large.

C-peptide Secretion During Hyperglycemic Clamp in Humans

As a preliminary validation of the present model, we simulated the C-peptide response to a hyperglycemic clamp in healthy humans. Previously obtained data of glucose, insulin, and C-peptide were used (32). The C-peptide secretion rate was computed according to the model of Eqs. 111, and the whole body C-peptide kinetics were represented by the usual two-compartment model (17, 39) with the standard kinetic parameters computed for each subject (51).

To determine the time course of glucose concentration that actually stimulates the in vivo pancreatic response, we estimated the time course of arterial glucose concentration, Ga, as shown in APPENDIX C. In the insulin secretion model, the number N of beta-cells in human pancreas was set according to the estimate in Ref. 14. A parameter identification was not attempted; instead, a few parameters were modified with respect to the baseline values previously reported. In particular, we used the alternative values for k1+, k1, and CT, and we set h = 1.38 x 10–3 min–1, a value smaller than that used in the simulations of rat pancreatic response. Other parameters were chosen in the following ranges: {rho}b = 0.01–0.02 min–1, k{rho} was two to four times the value in the rat, and bV = 6–8 min–1. We found that C-peptide data were adequately approximated by assuming f constant at a low value (f = 0.06–0.07). This may be related to a very slow recruitment of the beta-cells during the course of the experiment (fasting subjects) or to inhibitory signals acting on the islets in vivo. To obtain the observed value of C-peptide preceding the clamp, we added to the ISR a "constitutive" insulin secretion rate given by I0{sigma}FfN with {sigma}F = 1 and f = 0.05.

Figure 8 shows, for one of the subjects, the experimental data of glucose and C-peptide together with the time course of Ga, reconstructed according to Eq. C1. Figure 8B. reports the time courses of the ISR and of the C-peptide concentration as predicted by the model when the glucose stimulus in Eq. 7 is Ga (solid lines) or Gm (dotted lines). The initial more rapid increase of Ga with respect to Gm may be a factor that accounts for the occurrence of an appreciable first-phase ISR and the observed rapid increase of the C-peptide concentration. No significant differences in the C-peptide concentration calculated by the model with Ga or Gm were observed, likely because of the long time constant of plasma C-peptide kinetics (51).


Figure 8
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 8. C-peptide response to hyperglycemic clamp in human subjects (data from Ref. 32). A: measured glucose concentration (triangles) with the spline approximation (dotted line) and time course of reconstructed arterial glucose concentration (solid line). B: measured C-peptide (CP) concentration (filled squares) with the time course predicted by the model and the predicted ISR (solid lines). The dotted-line plots represent the model response when the glucose stimulus is the measured glucose concentration, Gm.

 
Interpretation of beta-Cell Responsivity Indexes

Using the approximate expressions of the first- and second-phase ISR given in APPENDIXES A and B, we may derive the expressions for the dynamic and the static responsivity index of the beta-cell, denoted as {Phi}d and {Phi}s in Refs 8, 49, and 10, that are predicted by the present model for a step increase in glucose concentration.

According to Ref. 10, we define a "dynamic responsivity index," Sd, as the ratio

Formula 16
where IS is the insulin amount secreted during the first-phase response, given by Eq. A4, and VC is the volume of the C-peptide-accessible compartment. Note that G* has a meaning similar to the threshold parameter h in Ref. 10. The index Sd is dimensionless, like the index {Phi}d. As seen from Eq. A4, the index Sd depends on G, so we consider the case of a small G step (a different expression may be derived for large G). By approximating exp(–aT) in Eq. A3 as 1 – aT, we get Z(G) = Formula 16DIR,bT. Furthermore, neglecting {rho}b in Formula 16, we obtain from Eq. A4 the equation for Sd:

Formula 17(17)
The above expression shows that Sd essentially depends on the size of the pool of immediately releasable granules and on the product of the [ATP] sensitivity to glucose, h/(G – G*), times the [Ca2+]i sensitivity to [ATP], k{rho}. These parameters are defined in Eqs. 8 and 10.

The "static responsivity index," in the spirit of the index {Phi}s in Ref. 10, takes into account the dynamics of the second-phase ISR. Thus we have obtained a dynamic index by computing the derivative at t = 0 of the second-phase ISR and dividing by the rate constant, kFormula 17 CT, of ISR increase. From Eqs. B3 and B4, we get for this ratio the quantity I0Formula 17RbfN. Proceeding as for the dynamic index, the static responsivity index Ss (with dimension of time–1) has the expression

Formula 18(18)
that depends on the [ATP] sensitivity to glucose and the size of the reserve pool. As noted in APPENDIX B, Eq. 18 is obtained under the assumption of large G (free [Ca2+]i at its maximum level).

With the values of model parameters used in the present study to simulate the response to the hyperglycemic clamp in a normal human subject, we obtain Sd {approx} 670 x 10–9 and Ss {approx} 26 x 10–9 min–1. These values are in the range of the {Phi}d and {Phi}s values estimated in Refs. 8, 49, and 10 under different experimental conditions and by means of a different model.


    DISCUSSION
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
In this study, a model of insulin granule trafficking in pancreatic beta-cells is proposed. Essentially, the model describes the formation of insulin granules and the dynamics of four intracellular granule pools: the reserve pool, the pool of docked granules, the pool of immediately releasable granules, and the pool of granules fused with cell membrane. Granule translocation from the trans-Golgi network to the cell membrane and exocytosis of granules fused with the cell membrane are regulated by two rate constants, {gamma} and {rho}, assumed to be related to the ATP-to-ADP ratio and other mitochondrial signals, and to the free [Ca2+]i, respectively. As anticipated by Grodsky (20), the "labile" insulin may be grossly identified with the insulin contained in the DIR pool. The "stable" insulin is the R pool, with the D pool, that feeds the DIR pool as insulin is secreted, being in an intermediate position. The model links in a quantitative, although simplistic, way the dynamics of granule release by the single beta-cell, with the recruitment of a subpopulation of glucose-responsive cells within the whole beta-cell population.

The different roles of the triggering pathway and of the amplifying pathway in the glucose-stimulated insulin secretion are modeled according to Ref. 21. The rate coefficient {gamma}, which represents events leading from glucose influx into the cell to ATP production, has an essential role in promoting the granule externalization and priming that supports the second-phase ISR. This rate coefficient, in particular the parameters {eta} and h in its governing equation, is a rate-limiting step in the glucose-stimulated insulin secretion. The representative parameters of the part of the triggering pathway from the ATP-to-ADP ratio to [Ca2+]i are {zeta} and k{rho}, which thus represent another rate-limiting step in the insulin secretion (48). We stress that the model does not include the regulatory role of Ca2+ in granule translocation and thus in the second-phase insulin release (16, 25), nor Ca2+ implication in ATP metabolism (1, 15). Moreover, the events on the fast time scale of the rapid fluctuations in beta-cell membrane potential (38, 46) are not accounted for.

The pattern of the pancreatic response to glucose is shown in the present study by the expression of the ISR in the first and the second phase (see APPENDIXES A and B), by the simulation tests, and by the responsivity indexes of Eqs. 1718. The first phase mainly depends on the insulin amount in the DIR pool, and on the product of [ATP] sensitivity to glucose, h/(G – G*), times the [Ca2+]i sensitivity to [ATP], k{rho}. Accordingly, the dynamic responsivity index Sd accounts for the fusion of DIR granules (release of labile insulin), essentially regulated by the increase in [Ca2+]i. Notably, agents that produce membrane depolarization, e.g., through closure of KATP channels, or potentiate the triggering pathway at some point, such as gastric inhibitory peptide (GIP) and glucagon-like peptide (GLP)-1, may result in an increased Sd. This action is, however, transient if it is not accompanied by adequate increase in the proinsulin synthesis rate and provision of granules competent for release. In contrast, the static responsivity index Ss mainly accounts for granule provision, which is regulated by the size of the reserve pool and the increase in the ATP-to-ADP ratio. The rate constant {gamma} is a synthetic descriptor of the signals that regulate these events, and its increase, together with the increase in bI, is able to produce a durable release of insulin without degranulation of the beta-cell.

The secretory response is also determined by the number, fN, of responsive beta-cells. Given the total number N of beta-cells in pancreas, possibly subject to slow but considerable changes in pathological conditions (14), the fraction f of recruited cells will change with the time in response to a variety of stimuli (glucose, insulin, GIP, GLP-1, etc.). The complex phenomena and interactions among the cells distributed in pancreas volume [for instance, the diffusion and consumption of glucose in the islets with consequent heterogeneity of glucose concentration (3) and electrical and chemical interactions (38)] are not truly represented in the present model. These phenomena are described by a dependence of the fraction f of active cells on G, at most with a time delay. It may be speculated that cell recruitment has a slower dynamics, also related to the concentration of insulin and other substances in the interstitial fluid that surrounds the beta-cells. In vivo, because of inhibitory signals, the maximal value attained by f is likely to be substantially smaller than the unity, and was estimated to be of the order of 6–7%.

Although the present model does not allow to discriminate among the various pathways involved in the pathophysiology of type 2 diabetes, it might provide some insight into the coupling between energy changes and insulin secretion. Therefore, should the mitochondrial derangement described in type 2 diabetes for skeletal muscle (29) occur also in beta-cell, then the parameters regulating the kinetics of the signals {gamma} and {rho} would be altered with a subsequent reduction in insulin secretion. The problem also arises of understanding how the various events involved in the trafficking of intracellular granules are coordinated and how the physiological amount of granules in the various pools is maintained over time. There may be control mechanisms devoted to "sensing" the size of granule pools to avoid granule depletion or excessive crowding and devoted to check the status of the network of microtubules that vehicle granule translocation.

The model proposed allows a good fitting of experimental data obtained during a hyperglycemic clamp, suggesting that it may be applied to in vivo experimental conditions where glucose is stably increased over time. However, the experimental pattern of the ISR and the value of the responsivity indexes, determined for instance in an oral glucose tolerance test (8, 10), will also be determined by factors other than glucose, such as the signals arising in the enteroinsular axis. The present model can represent the effect of these signals through changes in the equations for the control variables {gamma} and {rho}. For instance, the incretin effect may be translated into an increase in k{rho}. The k{rho} increase may also be related to an increase in the Ca2+ efficacy on exocytosis.

A similar consideration is also true for the effect of pharmacological agents on insulin secretion. The mechanism of action of different drugs that stimulate insulin secretion can be modeled, at least partially, by taking into account the drug effect at different levels, for instance on the KATP channels or on the metabolic energy level as it happens for acetylcarnitine. According to the observations in Ref. 21, it is important that drugs stimulating insulin secretion act not only at the level of the amplifying pathway but also of the triggering pathway; otherwise, the final effect would be that of an imbalance in the equilibrium between these pathways. Therefore, tailor-made drugs in diabetes should be designed to obtain the maximum effect.

A full description of these events cannot be given by the present model because the representation of the signaling pathways involved in the insulin secretion is still oversimplified. For instance, the role of long-chain acyl-CoAs and protein kinases (43, 53) is not taken into account. A more detailed description of the kinetics of signaling molecules should improve the capability of the model to investigate the events that lead to beta-cell dysfunction and predict the effects of agents stimulating insulin secretion.


    APPENDIX A: ANALYSIS OF FIRST-PHASE RESPONSE
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
We derive here an approximate expression for the insulin amount released in the first-phase response following a step in glucose concentration. The parameters appearing in this expression will be estimated from data reported by Grodsky (20) (see Fig. 4). We assume, for the sake of simplicity, that the dynamics of {gamma} and {rho} are much faster than granule dynamics (large {eta} and {zeta} values). Although the measured [ATP] and free [Ca2+] in beta-cell cytosol after a glucose step show rise times of a few minutes (1), glucose appears to provoke the formation of subplasma membrane ATP microdomains, where ATP rapidly responds to a glucose increase with a half-time = 20 s (28). The state of nearby KATP and L-type Ca2+ channels is likely to rapidly change, leading to fast Ca2+ influx at sites close to the granules of the DIR pool. The transient increase in [Ca2+]i at these sites, and thus in {rho}, may be completed in 1 min. We will take as baseline values {eta} = {zeta} = 4 min–1, so the transient of {gamma} and {rho} will require ~1 min. Simulations of a recently proposed model for mitochondrial ATP production (5) show that mitochondrial ATP and [Ca2+] rapidly respond to a step increase in fructose bisphosphate, with a transient that is completed in <1 min.

Following the step change in glucose concentration to a suprabasal value G, and after the time delay {tau}G, {gamma} and {rho} will rapidly attain their steady-state values, denoted as Formula 18 and Formula 18, respectively. From Eqs. 7 [with {psi}(t) = 0] and 9, we have Formula 18 = {gamma}b + h{gamma}(G) and Formula 18 = {rho}b + k{rho}h{gamma}(G). Let us take the time origin so that {gamma}(t) = Formula 18 and {rho}(t) = Formula 18 for t > 0. Assuming that pool D does not change appreciably during the first-phase response, we solve Eq. 5 with D(t) held constant at its initial (basal) value Db, {rho} equal to Formula 18, and initial condition DIR(0) = DIR,b. The pool DIR will be depleted according to

Formula A1(A1)
where a = k1+ Db +k1 + Formula A1, with Formula A1 being a function of G, and b = k1+DbCT. Then we solve Eq. 6 for F, with initial condition at the basal value Fb, {rho} = Formula A1, and DIR(t) given by Eq. A1, obtaining

Formula A2(A2)

To compare the prediction of the present model with the data of insulin amount released by the perfused rat pancreas during the first-phase response, we must compute the integral of the ISR over a time interval, say of length T, that just contains the first ISR peak. Recalling Eq. 11, we compute the total number of insulin granules released by a single beta-cell in the interval T, that is, the integral of {sigma}F(t) from t = 0 to t = T. Assuming a << {sigma} and Fb small, this integral, denoted as Z(G), is given by

Formula A3(A3)
Thus the insulin amount, IS(G), released in the first-phase response to a glucose step G can be written as

Formula A4(A4)
where N here is the number of beta-cells in rat pancreas. The function f was chosen of the form

Formula A5(A5)
where fb is the fraction of responding cells at glucose concentration smaller than the threshold G* (for simplicity assumed equal to that in Eq. 8), and Kf is a constant related to the effectiveness of recruitment. According to Eq. A5, {approx}100% of beta-cells may be recruited for a large enough G.

Eq. A4 was used to fit Grodsky's data. Taking the baseline values of model parameters, we have still to estimate the quantities G*, G, the product k{rho}h, fb, and Kf. To simplify the estimation problem, we set G = 10 mmol/l according to data in Ref. 26 and fb = 0.05. Least-squares fitting, used to estimate the remaining parameters, provided the following values: G* = 4.58 mmol/l, k{rho}h = 1.38 min–1, and Kf = 3.43 mmol/l (see Fig. 4). By model simulation we chose, as baseline parameters, h = 3.93 x 10–3 min–1 and k{rho} = 350. Finally, we note that a simple expression for the first-phase ISR can no longer be obtained if the dynamics of {gamma} and {rho} are not much faster than granule dynamics. Model simulations show that the peak ISR becomes smaller and is delayed in time when {eta} and {zeta} decrease, but the peak values in Fig. 5B can be recovered by increasing the product k{rho}h. If, for instance, {eta} and {zeta} are decreased from the baseline value of 4 min–1 to {eta} = {zeta} = 1 min–1, the baseline value of k{rho}h must be multiplied by 2.5.


    APPENDIX B: ANALYSIS OF SECOND-PHASE RESPONSE
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
By a similar analysis as in APPENDIX A, an approximate expression for the second-phase ISR can be obtained as follows. Let us assume that, during this phase, the pool R decays exponentially from the basal value Rb to a final value Formula A5, where Formula A5 is expressed as in Eq. 13 with the final values of {gamma} and bI,Formula A5 and Formula A51, respectively, and the other parameters unchanged. The rate constant of the exponential is Formula A5, so we have (initial time t = 0)

Formula 1(B1)
Using the quasistationary approximation Formula 1IR = 0 and summing the equations for D and DIR, we may write

Formula 2(B2)
For Formula 2 >> kFormula 2 D + kFormula 2(valid when G is large and the concentration of free Ca2+ is at its maximum level), the differential equation in B2 may be written as a linear differential equation for D that can be solved using Eq. B1 and with D(0) = Db, obtaining

Formula 3(B3)
By the further quasistationary approximation Formula 3 = 0, we finally get

Formula 4(B4)
that may be used with Eq. B3 to compute the approximate expression of the second-phase ISR. The derivative of {sigma}F(t) at t= 0 is given by kFormula 4 CT (Formula 4RbkFormula 4 CTDb).

If kFormula 4 CT = 7.24 x 10–3 min–1, and Formula 4 ≤ {gamma}b + h = 3.93 x 10–3 min–1 (baseline parameters), the descending front of second-phase ISR has a large time constant (1/Formula 4 > 4 h). The time constant becomes smaller if h is larger, as shown in Fig. 7A. With the baseline parameters, the negative term in the above derivative may be neglected for G >10 mmol/l.


    APPENDIX C: ESTIMATION OF ARTERIAL GLUCOSE CONCENTRATION
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 
An estimate of arterial glucose concentration was obtained by using the measured plasma glucose concentration, Gm, and by computing the net tissue glucose uptake from the insulin sensitivity index and the plasma insulin concentration. Neglecting non-insulin-dependent uptake, the following equation can be written:

Formula 1(C1)
where VG is the glucose distribution volume, Q is the cardiac output, SI is the insulin sensitivity, and Ip and Ip,b denote the current and the basal insulin concentration in plasma. The time course of Ga was thus obtained by computing the derivative Formula 1m from a smooth approximation of Gm. The parameters were assessed for each subject as follows: VG was obtained by using the estimate of 0.16 l/kg body wt (23), and cardiac output was given by Q = 0.2 x kg body wt0.75 (Brody's formula); SI was available from Ref. 32; Ip,b was computed as the mean Ip value in the last 30 min preceding the clamp. An improved estimate of Ga might be obtained by means of a more complex circulatory model (34).


    FOOTNOTES
 

Address for reprint requests and other correspondence: G. Mingrone, Istituto di Medicina Interna, Università Cattolica del Sacro Cuore, Largo A. Gemelli, 8, 00168 Rome, Italy (e-mail: gmingrone{at}rm.unicatt.it)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 Glossary
 DEVELOPMENT OF THE MATHEMATICAL...
 RESULTS
 DISCUSSION
 APPENDIX A: ANALYSIS OF...
 APPENDIX B: ANALYSIS OF...
 APPENDIX C: ESTIMATION OF...
 REFERENCES
 

  1. Ainscow EK, Rutter GA. Glucose-stimulated oscillations in free cytosolic ATP concentration imaged in single islet beta-cells. Evidence for a Ca2+-dependent mechanism. Diabetes 51: S162–S170, 2002.[Abstract/Free Full Text]
  2. Barg S, Eliasson L, Renström E, Rorsman P. A subset of 50 secretory granules in close contact with L-type Ca2+ channels accounts for the first-phase insulin secretion in mouse beta-cells. Diabetes 51: S74–S82, 2002.[Abstract/Free Full Text]
  3. Bertram R, Pernarowsky M. Glucose diffusion in pancreatic islets of Langerhans. Biophys J 74: 1722–1731, 1998.[Web of Science][Medline]
  4. Bertram R, Sherman A. A calcium-based phantom bursting model for pancreatic islets. Bull Math Biol 66: 1313–1344, 2004.[CrossRef][Web of Science][Medline]
  5. Bertram R, Pedersen MG, Luciani DS, Sherman A. A simplified model for mitochondrial ATP production. J Theor Biol 243: 575–586, 2006.[CrossRef][Web of Science][Medline]
  6. Bock T, Svenstrup K, Pakkenberg B, Buschard K. Unbiased estimation of total beta-cell number and mean beta-cell volume in rodent pancreas. Acta Pathol Microbiol Immunol Scand 107: 791–799, 1999.
  7. Bratanova-Tochkova TK, Cheng H, Daniel S, Gunawardana S, Liu YJ, Mulvaney-Musa J, Schermerhorn T, Straub SG, Yajima H, Sharp GWG. Triggering and augmentation mechanisms, granule pools, and biphasic insulin secretion. Diabetes 51: S83–S90, 2002.[Abstract/Free Full Text]
  8. Breda E, Cavaghan MK, Toffolo G, Polonsky KS, Cobelli C. Oral glucose tolerance test minimal model indexes of beta-cell function and insulin sensitivity. Diabetes 50: 150–158, 2001.[Abstract/Free Full Text]
  9. Butler AE, Janson J, Bonner-Weir S, Ritzel R, Rizza RA, Butler PC. beta-Cell deficit and increased beta-cell apoptosis in humans with type 2 diabetes. Diabetes 52: 102–110, 2003.[Abstract/Free Full Text]
  10. Campioni M, Toffolo G, Shuster LT, Service FJ, Rizza RA, Cobelli C. Incretin effect potentiates beta-cell responsivity to glucose as well as to its rate of change: OGTT and matched intravenous study. Am J Physiol Endocrinol Metab 292: E54–E60, 2007.[Abstract/Free Full Text]
  11. Cerasi E, Luft R. The plasma insulin response to glucose infusion in healthy subjects and in diabetes mellitus. Acta Endocrinol 55: 278–304, 1967.[Medline]
  12. Cerasi E, Fick G, Rudemo M. A mathematical model for the glucose induced insulin release in man. Eur J Clin Invest 4: 267–278, 1974.[Web of Science][Medline]
  13. Curry DL, Bennett LL, Grodsky GM. Dynamics of insulin secretion by the perfused rat pancreas. Endocrinology 83: 572–584, 1968.[Abstract/Free Full Text]
  14. Deng S, Vatamaniuk M, Huang X, Doliba N, Lian MM, Frank A, Velidedeoglu E, Desai NM, Koeberlein B, Wolf B, Barker CF, Naji A, Matschinsky FM, Markmann JF. Structural and functional abnormalities in the islets isolated from type 2 diabetic subjects. Diabetes 53: 624–632, 2004.[Abstract/Free Full Text]
  15. Detimary P, Gilon P, Henquin JC. Interplay between cytoplasmic Ca2+ and the ATP/ADP ratio: a feedback control mechanism in mouse pancreatic islets. Biochem J 333: 269–274, 1998.[Web of Science][Medline]
  16. Donelan MJ, Morfini G, Julyan R, Sommers S, Hays L, Kajio H, Briaud I, Easom RA, Molkentini JD, Brady ST, Rhodes CJ. Ca2+-dependent dephosphorylation of kinesin heavy chain on beta-granules in pancreatic beta-cells. J Biol Chem 277: 24232–24242, 2002.[Abstract/Free Full Text]
  17. Eaton RP, Allen RC, Schade DS, Erickson KM, Standefer J. Prehepatic insulin production in man: kinetic analysis using peripheral connecting peptide behavior. J Clin Endocrinol Metab 51: 520–528, 1980.[Abstract/Free Full Text]
  18. Garcia-Barrado MJ, Ravier MA, Rolland JF, Gilon P, Nenquin M, Henquin JC. Inhibition of protein synthesis sequentially impairs distinct steps of stimulus-secretion coupling in pancreatic beta-cells. Endocrinology 142: 299–307, 2001.[Abstract/Free Full Text]
  19. Gilon P, Ravier MA, Jonas JC, Henquin JC. Control mechanisms of the oscillations of insulin secretion in vitro and in vivo. Diabetes 51: S144–S151, 2002.[Abstract/Free Full Text]
  20. Grodsky GM. A threshold distribution hypothesis for packet storage of insulin and its mathematical model. J Clin Invest 51: 2047–2059, 1972.[Web of Science][Medline]
  21. Henquin JC. Triggering and amplifying pathways of regulation of insulin secretion by glucose. Diabetes 49: 1751–1760, 2000.[Abstract]
  22. Henquin JC, Ishiyama N, Nenquin M, Ravier MA, Jonas JC. Signals and pools underlying biphasic insulin secretion. Diabetes 51: S60–S67, 2002.[Abstract/Free Full Text]
  23. Hovorka R, Shojaee-Moradie F, Carroll PV, Chassin LJ, Gowrie IJ, Jackson NC, Tudor RS, Umpleby AM, Jones RH. Partitioning glucose distribution/transport, disposal, and endogenous production during IVGTT. Am J Physiol Endocrinol Metab 282: E992–E1007, 2002.[Abstract/Free Full Text]
  24. Hu S. Interaction of nateglinide with KATP channel in beta-cells underlies its unique insulinotropic action. Eur J Pharmacol 442: 163–171, 2002.[CrossRef][Web of Science][Medline]
  25. Jing X, Li DQ, Olofsson CS, Salehi A, Surve VV, Caballero J, Ivarsson R, Lundquist I, Pereverzev A, Schneider T, Rorsman P, Renström E. CaV2.3 calcium channels control second-phase insulin release. J Clin Invest 115: 146–154, 2005.[CrossRef][Web of Science][Medline]
  26. Jonkers FC, Henquin JC. Measurements of cytoplasmic Ca2+ in islet cell clusters show that glucose rapidly recruits beta-cells and gradually increases the individual cell response. Diabetes 50: 540–550, 2001.[Abstract/Free Full Text]
  27. Juhl CB, Pørksen N, Hollingdal M, Sturis J, Pincus S, Veldhuis JD, Dejgaard A, Schmitz O. Repaglinide acutely amplifies pulsatile insulin secretion by augmentation of burst mass with no effect on burst frequency. Diab Care 23: 675–681, 2000.[Abstract/Free Full Text]
  28. Kennedy HJ, Pouli AE, Ainscow EK, Jouaville LS, Rizzuto R, Rutter GA. Glucose generates sub-plasma membrane ATP microdomains in single islet beta-cells. J Biol Chem 274: 13281–13291, 1999.[Abstract/Free Full Text]
  29. Lowell BB, Shulman GI. Mitochondrial dysfunction and type 2 diabetes. Science 307: 384–387, 2005.[Abstract/Free Full Text]
  30. MacDonald MJ, Fahien LA, Brown LJ, Hasan NM, Buss JD, Kendrick MA. Perspective: emerging evidence for signaling roles of mitochondrial anaplerotic products in insulin secretion. Am J Physiol Endocrinol Metab 288: E1–E15, 2005.[Abstract/Free Full Text]
  31. Maechler P, Wollheim CB. Mitochondrial signals in glucose-stimulated insulin secretion in the beta cell. J Physiol 529: 49–56, 2000.[Abstract/Free Full Text]
  32. Manco M, Bertuzzi A, Salinari S, Scarfone A, Calvani M, Greco AV, Mingrone G. The ingestion of saturated fatty acids and triacylglycerols acutely affects insulin secretion and insulin sensitivity in human subjects. Br J Nutr 92: 895–903, 2004.[CrossRef][Web of Science][Medline]
  33. Mari A, Camastra S, Toschi E, Giancaterini A, Gastaldelli A, Mingrone G, Ferrannini E. A model for glucose control of insulin secretion in 24 h of free living. Diabetes 50: S164–S168, 2001.[Free Full Text]
  34. Mari A, Stojanovska L, Proietto J, Thorburn AW. A circulatory model for calculating non-steady-state glucose fluxes. Validation and comparison with compartmental models. Comput Meth Prog Biomed 71: 269–281, 2003.[CrossRef]
  35. Molinete M, Irminger JC, Tooze SA, Halban PA. Trafficking/sorting and granule biogenesis in the beta-cell. Cell Dev Biol 11: 243–251, 2000.[CrossRef]
  36. Nesher R, Cerasi E. Modeling phasic insulin release. Immediate and time-dependent effects of glucose. Diabetes 51: S53–S59, 2002.[Abstract/Free Full Text]
  37. O'Connor MDL, Landahl H, Grodsky GM. Comparison of storage- and signal-limited models of pancreatic insulin secretion. Am J Physiol Regul Integr Comp Physiol 238: R378–R389, 1980.[Abstract/Free Full Text]
  38. Pedersen MG, Bertram R, Sherman A. Intra- and inter-islet synchronization of metabolically driven insulin secretion. Biophys J 89: 107–119, 2005.[CrossRef][Web of Science][Medline]
  39. Polonsky KS, Licinio-Paixao J, Given BD, Pugh W, Rue P, Galloway J, Karrison T, Frank B. Use of biosynthetic human C-peptide in the measurement of insulin secretion rates in normal volunteers and type I diabetic patients. J Clin Invest 77: 98–105, 1986.[Web of Science][Medline]
  40. Pørksen N, Grøfte T, Greisen J, Mengel A, Juhl C, Veldhuis JD, Schmitz O, Rössle M, Vilstrup H. Human insulin release processes measured by intraportal sampling. Am J Physiol Endocrinol Metab 282: E695–E702, 2002.[Abstract/Free Full Text]
  41. Rorsman P, Renström E. Insulin granule dynamics in pancreatic beta cells. Diabetologia 46: 1029–1045, 2003.[CrossRef][Web of Science][Medline]
  42. Sando H, Borg J, Steiner DF. Studies on the secretion of newly synthesized proinsulin and insulin from isolated rat islets of Langerhans. J Clin Invest 51: 1476–1485, 1972.[Web of Science][Medline]
  43. Sato Y, Henquin JC. The K+-ATP channel-independent pathway of regulation of insulin secretion by glucose. In search of the underlying mechanism. Diabetes 47: 1713–1721, 1998.[Abstract]
  44. Sato Y, Anello M, Henquin JC. Glucose regulation of insulin secretion independent of the opening or closure of adenosine triphosphate-sensitive K+ channels in beta-cells. Endocrinology 140: 2252–2257, 1999.[Abstract/Free Full Text]
  45. Schuit FC, In 't Veld PA, Pipeleers DG. Glucose stimulates proinsulin biosynthesis by a dose-dependent recruitment of pancreatic beta cells. Proc Nat Acad Sci USA 85: 3865–3869, 1988.[Abstract/Free Full Text]
  46. Sherman A. Contributions of modeling to understanding of stimulus-secretion coupling in pancreatic beta-cells. Am J Physiol Endocrinol Metab 271: E362–E372, 1996.[Abstract/Free Full Text]
  47. Srivastava S, Goren HJ. Insulin constitutively secreted by beta-cells is necessary for glucose-stimulated insulin secretion. Diabetes 52: 2049–2056, 2003.[Abstract/Free Full Text]
  48. Straub SG, Sharp GWG. Hypothesis: one rate-limiting step controls the magnitude of both phases of glucose-stimulated insulin secretion. Am J Physiol Cell Physiol 287: C565–C571, 2004.[Abstract/Free Full Text]
  49. Toffolo G, Breda E, Cavaghan MK, Ehrmann DA, Polonsky KS, Cobelli C. Quantitative indexes of beta-cell function during graded up&down glucose infusion from C-peptide minimal models. Am J Physiol Endocrinol Metab 280: E2–E10, 2001.[Abstract/Free Full Text]
  50. Tsuboi T, Rutter GA. Multiple forms of "kiss-and-run" exocytosis revealed by evanescent wave microscopy. Curr Biol 13: 563–567, 2003.[CrossRef][Web of Science][Medline]
  51. Van Cauter E, Mestrez F, Sturis J, Polonsky KS. Estimation of insulin secretion rates from C-peptide levels. Comparison of individual and standard kinetic parameters for C-peptide clearance. Diabetes 41: 368–377, 1992.[Abstract]
  52. Varadi A, Ainscow EK, Allan VJ, Rutter GA. Involvement of conventional kinesin in glucose-stimulated secretory granule movements and exocytosis in clonal pancreatic beta-cells. J Cell Sci 115: 4177–4189, 2002.[Abstract/Free Full Text]
  53. Yaney GC, Corkey BE. Fatty acid metabolism and insulin secretion in pancreatic beta cells. Diabetologia 46: 1297–1312, 2003.[CrossRef][Web of Science][Medline]



This article has been cited by other articles:


Home page
Phil Trans R Soc AHome page
M. G. Pedersen, A. Corradin, G. M Toffolo, and C. Cobelli
A subcellular model of glucose-stimulated pancreatic insulin secretion
Phil Trans R Soc A, October 13, 2008; 366(1880): 3525 - 3543.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
293/1/E396    most recent
00647.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bertuzzi, A.
Right arrow Articles by Mingrone, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bertuzzi, A.
Right arrow Articles by Mingrone, G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2007 by the American Physiological Society.