Search

Sodium-dependent glucose transporter (SGLT) 2 is specifically expressed in the kidney, while SGLT1 is present in the kidneys and small intestine. SGLT2 inhibitors are a class of oral antidiabetic drugs that lower elevated plasma glucose levels by promoting the urinary excretion of excess glucose through the inhibition of renal glucose reuptake. The inhibition selectivity for SGLT2 over SGLT1 (SGLT2/1 selectivity) of marketed SGLT2 inhibitors is diverse, while SGLT2/1 selectivity of canagliflozin is relatively low. Although canagliflozin suppresses postprandial glucose levels, the degree of contribution for SGLT1 inhibition to this effect remains unproven. To analyze the effect of SGLT2 inhibitors on postprandial glucose level, we constructed a novel quantitative systems pharmacology (QSP) model, called human systemic glucose dynamics (HSGD) model, integrating intestinal absorption, metabolism, and renal reabsorption of glucose. This HSGD model reproduced the postprandial plasma glucose concentration-time profiles during a meal tolerance test under different clinical trial conditions. Simulations after canagliflozin administration showed a dose-dependent delay of time (Tmax,glc ) to reach maximum concentration of glucose (Cmax,glc ), and the delay of Tmax,glc disappeared when inhibition of SGLT1 was negated. In addition, contribution ratio of intestinal SGLT1 inhibition to the decrease in Cmax,glc was estimated to be 23%-28%, when 100 and 300 mg of canagliflozin are administered. This HSGD model enabled us to provide the partial contribution of intestinal SGLT1 inhibition to the improvement of postprandial hyperglycemia as well as to quantitatively describe the plasma glucose dynamics following SGLT2 inhibitors.

We present a general method for the analysis of the discharge trains of periodically forced noisy leaky integrate-and-fire neuron models. This approach relies on the iterations of a stochastic phase transition operator that generalizes the phase transition function used for the study of periodically forced deterministic oscillators to noisy systems. The kernel of this operator is defined in terms of the the first passage time probability density function of the Ornstein Uhlenbeck process through a suitable threshold. Numerically, it is computed as the solution of a singular integral equation. It is shown that, for the noisy system, quantities such as the phase distribution (cycle histogram), the interspike interval distribution, the autocorrelation function of the intervals, the autocorrelogram and the power spectrum density of the spike train, as well as the input-output cross-correlation and cross-spectral density can all be computed using the stochastic phase transition operator. A detailed description of the numerical implementation of the method, together with examples, is provided.

The force-velocity (F-V) relationship of filament sliding is traditionally used to define the inotropic condition of striated muscles. A simple circulation model combined with the Laplace heart was developed to get a deeper insight into the relationship between the F-V characteristics and the cardiac ventricular inotropy. The circulation model consists of a preload and an afterload compartments. The linear F-V relationship for filament sliding in the NL model (Negroni and Lascano 1996) was replaced by the exponential F-V relation observed by Piazzesi et al. (2002). We also modified the NL model to a hybrid model to benefit from the Ca(2+) cooperativity described by the Robinson model (Robinson et al. 2002). The model was validated by determining the diastolic ventricular pressure-volume relationship of the Laplace heart and the F-V relation of the new hybrid model. The computed parameters of the cardiac cycle agreed well with the physiological data. Computational results showed that the cross-bridge elongation (h in the NL model) temporally undershot the equilibrium h(c) during the ejection period and overshot it during the rapid refilling phase. Thereby the time course of ejection and refilling was retarded. In a simulation where the velocity of the mobile myosin head (dX/dt) was varied, the systolic peak pressure of the ventricle varied from a minimum value at dX/dt = 0 to a saturating value obtained with a constant h(c), providing in silico evidence for a functional impact of the cross-bridge sliding rate on the ventricular inotropy.