Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Fat oxidation kinetics during exercise in obese individuals

Introduction An impaired ability to oxidize fat may be a factor in the obesity's aetiology (3). Moreover, the exercise intensity (Fatmax) eliciting the maximal fat oxidation rate (MFO) was lower in obese (O) compared with lean (L) individuals (4). However, difference in fat oxidation rate (FOR) during exercise between O and L remains equivocal and little is known about FORs during high intensities (>60% ) in O compared with L. This study aimed to characterize fat oxidation kinetics over a large range of intensities in L and O. Methods 12 healthy L [body mass index (BMI): 22.8±0.4] and 16 healthy O men (BMI: 38.9±1.4) performed submaximal incremental test (Incr) to determine whole-body fat oxidation kinetics using indirect calorimetry. After a 15-min resting period (Rest) and 10-min warm-up at 20% of maximal power output (MPO, determined by a maximal incremental test), the power output was increased by 7.5% MPO every 6-min until respiratory exchange ratio reached 1.0. Venous lactate and glucose and plasma concentration of epinephrine (E), norepinephrine (NE), insulin and non-esterified fatty acid (NEFA) were assessed at each step. A mathematical model (SIN) (1), including three variables (dilatation, symmetry, translation), was used to characterize fat oxidation (normalized by fat-free mass) kinetics and to determine Fatmax and MFO. Results FOR at Rest and MFO were not significantly different between groups (p≥0.1). FORs were similar from 20-60% (p≥0.1) and significantly lower from 65-85% in O than in L (p≤0.04). Fatmax was significantly lower in O than in L (46.5±2.5 vs 56.7±1.9 % respectively; p=0.005). Fat oxidation kinetics was characterized by similar translation (p=0.2), significantly lower dilatation (p=0.001) and tended to a left-shift symmetry in O compared with L (p=0.09). Plasma E, insulin and NEFA were significantly higher in L compared to O (p≤0.04). There were no significant differences in glucose, lactate and plasma NE between groups (p≥0.2). Conclusion The study showed that O presented a lower Fatmax and a lower reliance on fat oxidation at high, but not at moderate, intensities. This may be linked to a: i) higher levels of insulin and lower E concentrations in O, which may induce blunted lipolysis; ii) higher percentage of type II and a lower percentage of type I fibres (5), and iii) decreased mitochondrial content (2), which may reduce FORs at high intensities and Fatmax. These findings may have implications for an appropriate exercise intensity prescription for optimize fat oxidation in O. References 1. Cheneviere et al. Med Sci Sports Exerc. 2009 2. Holloway et al. Am J Clin Nutr. 2009 3. Kelley et al. Am J Physiol. 1999 4. Perez-Martin et al. Diabetes Metab. 2001 5. Tanner et al. Am J Physiol Endocrinol Metab. 2002

Retrovirus-induced target cell activation in the early phases of infection: the mouse mammary tumor virus model.

Mouse mammary tumor virus (MMTV) infects B lymphocytes and expresses a superantigen on the cell surface after integration of its reverse-transcribed genome. Superantigen-dependent B- and T-cell activation becomes detectable 2 to 3 days after infection. We show here that before this event, B cells undergo a polyclonal activation which does not involve massive proliferation. This first phase of B-cell activation is T cell independent. Moreover, during the first phase of activation, when only a small fraction of B cells is infected by MMTV(SW), viral DNA is detected only in activated B cells. Such a B-cell activation is also seen after injection of murine leukemia virus but not after injection of vaccinia virus, despite the very similar kinetics and intensity of the immune response. Since retroviruses require activated target cells to induce efficient infection, these data suggest that the early polyclonal retrovirus-induced target cell activation might play an important role in the establishment of retroviral infections.

Adaptive model selection using empirical complexities

Given $n$ independent replicates of a jointly distributed pair $(X,Y)\in {\cal R}^d \times {\cal R}$, we wish to select from a fixed sequence of model classes ${\cal F}_1, {\cal F}_2, \ldots$ a deterministic prediction rule $f: {\cal R}^d \to {\cal R}$ whose risk is small. We investigate the possibility of empirically assessingthe {\em complexity} of each model class, that is, the actual difficulty of the estimation problem within each class. The estimated complexities are in turn used to define an adaptive model selection procedure, which is based on complexity penalized empirical risk.The available data are divided into two parts. The first is used to form an empirical cover of each model class, and the second is used to select a candidate rule from each cover based on empirical risk. The covering radii are determined empirically to optimize a tight upper bound on the estimation error. An estimate is chosen from the list of candidates in order to minimize the sum of class complexity and empirical risk. A distinguishing feature of the approach is that the complexity of each model class is assessed empirically, based on the size of its empirical cover.Finite sample performance bounds are established for the estimates, and these bounds are applied to several non-parametric estimation problems. The estimates are shown to achieve a favorable tradeoff between approximation and estimation error, and to perform as well as if the distribution-dependent complexities of the model classes were known beforehand. In addition, it is shown that the estimate can be consistent,and even possess near optimal rates of convergence, when each model class has an infinite VC or pseudo dimension.For regression estimation with squared loss we modify our estimate to achieve a faster rate of convergence.

Atrial septal pacing for the termination of atrial fibrillation: study in a biophysical model of human atria.

AIMS: While successful termination by pacing of organized atrial tachycardias has been observed in patients, single site rapid pacing has not yet led to conclusive results for the termination of atrial fibrillation (AF). The purpose of this study was to evaluate a novel atrial septal pacing algorithm for the termination of AF in a biophysical model of the human atria. METHODS AND RESULTS: Sustained AF was generated in a model based on human magnetic resonance images and membrane kinetics. Rapid pacing was applied from the septal area following a dual-stage scheme: (i) rapid pacing for 10-30 s at pacing intervals 62-70% of AF cycle length (AFCL), (ii) slow pacing for 1.5 s at 180% AFCL, initiated by a single stimulus at 130% AFCL. Atrial fibrillation termination success rates were computed. A mean success rate for AF termination of 10.2% was obtained for rapid septal pacing only. The addition of the slow pacing phase increased this rate to 20.2%. At an optimal pacing cycle length (64% AFCL) up to 29% of AF termination was observed. CONCLUSION: The proposed septal pacing algorithm could suppress AF reentries in a more robust way than classical single site rapid pacing. Experimental studies are now needed to determine whether similar termination mechanisms and rates can be observed in animals or humans, and in which types of AF this pacing strategy might be most effective.

A wearable system for multi-segment foot kinetics measurement.

This study aims to design a wearable system for kinetics measurement of multi-segment foot joints in long-distance walking and to investigate its suitability for clinical evaluations. The wearable system consisted of inertial sensors (3D gyroscopes and 3D accelerometers) on toes, forefoot, hindfoot, and shank, and a plantar pressure insole. After calibration in a laboratory, 10 healthy elderly subjects and 12 patients with ankle osteoarthritis walked 50m twice wearing this system. Using inverse dynamics, 3D forces, moments, and power were calculated in the joint sections among toes, forefoot, hindfoot, and shank. Compared to those we previously estimated for a one-segment foot model, the sagittal and transverse moments and power in the ankle joint, as measured via multi-segment foot model, showed a normalized RMS difference of less than 11%, 14%, and 13%, respectively, for healthy subjects, and 13%, 15%, and 14%, for patients. Similar to our previous study, the coronal moments were not analyzed. Maxima-minima values of anterior-posterior and vertical force, sagittal moment, and power in shank-hindfoot and hindfoot-forefoot joints were significantly different between patients and healthy subjects. Except for power, the inter-subject repeatability of these parameters was CMC>0.90 for healthy subjects and CMC>0.70 for patients. Repeatability of these parameters was lower for the forefoot-toes joint. The proposed measurement system estimated multi-segment foot joints kinetics with acceptable repeatability but showed difference, compared to those previously estimated for the one-segment foot model. These parameters also could distinguish patients from healthy subjects. Thus, this system is suggested for outcome evaluations of foot treatments.

Interplay between T cell receptor binding kinetics and the level of cognate peptide presented by major histocompatibility complexes governs CD8+ T cell responsiveness.

Through a rational design approach, we generated a panel of HLA-A*0201/NY-ESO-1(157-165)-specific T cell receptors (TCR) with increasing affinities of up to 150-fold from the wild-type TCR. Using these TCR variants which extend just beyond the natural affinity range, along with an extreme supraphysiologic one having 1400-fold enhanced affinity, and a low-binding one, we sought to determine the effect of TCR binding properties along with cognate peptide concentration on CD8(+) T cell responsiveness. Major histocompatibility complexes (MHC) expressed on the surface of various antigen presenting cells were peptide-pulsed and used to stimulate human CD8(+) T cells expressing the different TCR via lentiviral transduction. At intermediate peptide concentration we measured maximum cytokine/chemokine secretion, cytotoxicity, and Ca(2+) flux for CD8(+) T cells expressing TCR within a dissociation constant (K(D)) range of ∼1-5 μM. Under these same conditions there was a gradual attenuation in activity for supraphysiologic affinity TCR with K(D) < ∼1 μM, irrespective of CD8 co-engagement and of half-life (t(1/2) = ln 2/k(off)) values. With increased peptide concentration, however, the activity levels of CD8(+) T cells expressing supraphysiologic affinity TCR were gradually restored. Together our data support the productive hit rate model of T cell activation arguing that it is not the absolute number of TCR/pMHC complexes formed at equilibrium, but rather their productive turnover, that controls levels of biological activity. Our findings have important implications for various immunotherapies under development such as adoptive cell transfer of TCR-engineered CD8(+) T cells, as well as for peptide vaccination strategies.

Fat oxidation kinetics during exercise in obese individuals

Introduction An impaired ability to oxidize fat may be a factor in the obesity's aetiology (3). Moreover, the exercise intensity (Fatmax) eliciting the maximal fat oxidation rate (MFO) was lower in obese (O) compared with lean (L) individuals (4). However, difference in fat oxidation rate (FOR) during exercise between O and L remains equivocal and little is known about FORs during high intensities (>60% ) in O compared with L. This study aimed to characterize fat oxidation kinetics over a large range of intensities in L and O. Methods 12 healthy L [body mass index (BMI): 22.8±0.4] and 16 healthy O men (BMI: 38.9±1.4) performed submaximal incremental test (Incr) to determine whole-body fat oxidation kinetics using indirect calorimetry. After a 15-min resting period (Rest) and 10-min warm-up at 20% of maximal power output (MPO, determined by a maximal incremental test), the power output was increased by 7.5% MPO every 6-min until respiratory exchange ratio reached 1.0. Venous lactate and glucose and plasma concentration of epinephrine (E), norepinephrine (NE), insulin and non-esterified fatty acid (NEFA) were assessed at each step. A mathematical model (SIN) (1), including three variables (dilatation, symmetry, translation), was used to characterize fat oxidation (normalized by fat-free mass) kinetics and to determine Fatmax and MFO. Results FOR at Rest and MFO were not significantly different between groups (p≥0.1). FORs were similar from 20-60% (p≥0.1) and significantly lower from 65-85% in O than in L (p≤0.04). Fatmax was significantly lower in O than in L (46.5±2.5 vs 56.7±1.9 % respectively; p=0.005). Fat oxidation kinetics was characterized by similar translation (p=0.2), significantly lower dilatation (p=0.001) and tended to a left-shift symmetry in O compared with L (p=0.09). Plasma E, insulin and NEFA were significantly higher in L compared to O (p≤0.04). There were no significant differences in glucose, lactate and plasma NE between groups (p≥0.2). Conclusion The study showed that O presented a lower Fatmax and a lower reliance on fat oxidation at high, but not at moderate, intensities. This may be linked to a: i) higher levels of insulin and lower E concentrations in O, which may induce blunted lipolysis; ii) higher percentage of type II and a lower percentage of type I fibres (5), and iii) decreased mitochondrial content (2), which may reduce FORs at high intensities and Fatmax. These findings may have implications for an appropriate exercise intensity prescription for optimize fat oxidation in O. References 1. Cheneviere et al. Med Sci Sports Exerc. 2009 2. Holloway et al. Am J Clin Nutr. 2009 3. Kelley et al. Am J Physiol. 1999 4. Perez-Martin et al. Diabetes Metab. 2001 5. Tanner et al. Am J Physiol Endocrinol Metab. 2002

Kinetics of K release from soils of Brazilian coffee regions: effect of organic acids

Kinetic studies on soil potassium release can contribute to a better understanding of K availability to plants. This study was conducted to evaluate K release rates from the whole soil, clay, silt, and sand fractions of B-horizon samples of a basalt-derived Oxisol and a sienite-derived Ultisol, both representative soils from coffee regions of Minas Gerais State, Brazil. Potassium was extracted from each fraction after eight different shaking time periods (0-665 h) with either 0.001 mol L-1 citrate or oxalate at a 1:10 solid:solution ratio. First-order, Elovich, zero-order, and parabolic diffusion equations were used to parameterize the time dependence of K release. For the Oxisol, the first-order equation fitted best to the experimental data of K release, with similar rates for all fractions and independent of the presence of citrate or oxalate in the extractant solution. For all studied Ultisol fractions, in which K release rates increased when extractions were performed with citrate solution, the Elovich model described K release kinetics most adequately. The highest potassium release rate of the Ultisol silt fraction was probably due to the transference of "non-exchangeable" K to the extractant solution, whereas in the Oxisol exchangeable potassium represented the main K source in all studied fractions.

A coherent neurobiological framework for functional neuroimaging provided by a model integrating compartmentalized energy metabolism.

Functional neuroimaging has undergone spectacular developments in recent years. Paradoxically, its neurobiological bases have remained elusive, resulting in an intense debate around the cellular mechanisms taking place upon activation that could contribute to the signals measured. Taking advantage of a modeling approach, we propose here a coherent neurobiological framework that not only explains several in vitro and in vivo observations but also provides a physiological basis to interpret imaging signals. First, based on a model of compartmentalized energy metabolism, we show that complex kinetics of NADH changes observed in vitro can be accounted for by distinct metabolic responses in two cell populations reminiscent of neurons and astrocytes. Second, extended application of the model to an in vivo situation allowed us to reproduce the evolution of intraparenchymal oxygen levels upon activation as measured experimentally without substantially altering the initial parameter values. Finally, applying the same model to functional neuroimaging in humans, we were able to determine that the early negative component of the blood oxygenation level-dependent response recorded with functional MRI, known as the initial dip, critically depends on the oxidative response of neurons, whereas the late aspects of the signal correspond to a combination of responses from cell types with two distinct metabolic profiles that could be neurons and astrocytes. In summary, our results, obtained with such a modeling approach, support the concept that both neuronal and glial metabolic responses form essential components of neuroimaging signals.

Fat oxidation kinetics : effect of exercise

ABSTRACT Fat oxidation kinetics: effect of exercise. During graded exercise, absolute whole body fat oxidation rates increase from low to moderate intensities, and then markedly decline at high intensities, implying an exercise intensity (Fatmax) at which the fat oxidation rate is maximal (MFO). The main aim of the present work was to examine the effect of exercise on whole body fat oxidation kinetics. For this purpose, a sinusoidal mathematical model (SIN) has been developped in the first study to provide an accurate description of the shape of fat oxidation kinetics during graded exercise, represented as a function of exercise intensity, and to determine Fatmax and MFO. The SIN model incorporates three independent variables (i.e., dilatation, symmetry, and translation) that correspond to main expected modulations of the basic fat oxidation curve because of factors such as mode of exercise or training status. The results of study 1 showed that the SIN model was a valuable tool to determine Fatmax and MFO, and to precisely characterize and quantify the different shape of fat oxidation kinetics through its three variables. The effectiveness of the SIN model to detect differences in fat oxidation kinetics induced by a specific factor was then confirmed in the second study, which quantitatively described and compared fat oxidation kinetics in two different popular modes of exercise: running and cycling. It was found that the mean fat oxidation kinetics during running was characterized by a greater dilatation and a rightward asymmetry compared with the symmetric parabolic curve in cycling. In the two subsequent studies, the effect of a prior endurance exercise of different intensities and durations on whole body fat oxidation kinetics was examined. Study 3 determined the impact of a 1-h continuous exercise bout at an exercise intensity corresponding to Fatmax on fat oxidation kinetics during a subsequent graded test, while study 4 investigated the effect of an exercise leading to a more pronounced muscle glycogen depletion. The results of these two latter studies showed that fat oxidation rates, MFO, and Fatmax were enhanced following endurance exercise, but were increased to a greater extent with a more severe mucle glycogen depletion, inducing therefore modifications in the postexercise fat oxidation kinetics (i.e., greater dilatation and rightward asymmetry). In perspective, further studies have been suggested 1) to assess physiological meaning of the three independent variables of the SIN model; and 2) to compare the effect of two different training programs on fat oxidation kinetics in obese subjects.

Drug-eluting beads loaded with antiangiogenic agents for chemoembolization: in vitro sunitinib loading and release and in vivo pharmacokinetics in an animal model.

PURPOSE: The combination of embolic beads with a multitargeted tyrosine kinase inhibitor that inhibits tumor vessel growth is suggested as an alternative and improvement to the current standard doxorubicin-eluting beads for use in transarterial chemoembolization. This study demonstrates the in vitro loading and release kinetics of sunitinib using commercially available embolization microspheres and evaluates the in vitro biologic efficacy on cell cultures and the resulting in vivo pharmacokinetics profiles in an animal model. MATERIALS AND METHODS: DC Bead microspheres, 70-150 µm and 100-300 µm (Biocompatibles Ltd., Farnham, United Kingdom), were loaded by immersion in sunitinib solution. Drug release was measured in saline in a USP-approved flow-through apparatus and quantified by spectrophotometry. Activity after release was confirmed in cell culture. For pharmacokinetics and in vivo toxicity evaluation, New Zealand white rabbits received sunitinib either by intraarterial injection of 100-300 µm sized beads or per os. Plasma and liver tissue drug concentrations were assessed by liquid chromatography-tandem mass spectroscopy. RESULTS: Sunitinib loading on beads was close to complete and homogeneous. A total release of 80% in saline was measured, with similar fast-release profiles for both sphere sizes. After embolization, drug plasma levels remained below the therapeutic threshold (< 50 ng/mL), but high concentrations at 6 hours (14.9 µg/g) and 24 hours (3.4 µg/g) were found in the liver tissue. CONCLUSIONS: DC Bead microspheres of two sizes were efficiently loaded with sunitinib and displayed a fast and almost complete release in saline. High liver drug concentrations and low systemic levels indicated the potential of sunitinib-eluting beads for use in embolization.

Fat oxidation kinetics during exercise in obese individuals.

Introduction An impaired ability to oxidize fat may be a factor in the obesity's aetiology (3). Moreover, the exercise intensity (Fatmax) eliciting the maximal fat oxidation rate (MFO) was lower in obese (O) compared with lean (L) individuals (4). However, difference in fat oxidation rate (FOR) during exercise between O and L remains equivocal and little is known about FORs during high intensities (>60% ) in O compared with L. This study aimed to characterize fat oxidation kinetics over a large range of intensities in L and O. Methods 12 healthy L [body mass index (BMI): 22.8±0.4] and 16 healthy O men (BMI: 38.9±1.4) performed submaximal incremental test (Incr) to determine whole-body fat oxidation kinetics using indirect calorimetry. After a 15-min resting period (Rest) and 10-min warm-up at 20% of maximal power output (MPO, determined by a maximal incremental test), the power output was increased by 7.5% MPO every 6-min until respiratory exchange ratio reached 1.0. Venous lactate and glucose and plasma concentration of epinephrine (E), norepinephrine (NE), insulin and non-esterified fatty acid (NEFA) were assessed at each step. A mathematical model (SIN) (1), including three variables (dilatation, symmetry, translation), was used to characterize fat oxidation (normalized by fat-free mass) kinetics and to determine Fatmax and MFO. Results FOR at Rest and MFO were not significantly different between groups (p≥0.1). FORs were similar from 20-60% (p≥0.1) and significantly lower from 65-85% in O than in L (p≤0.04). Fatmax was significantly lower in O than in L (46.5±2.5 vs 56.7±1.9 % respectively; p=0.005). Fat oxidation kinetics was characterized by similar translation (p=0.2), significantly lower dilatation (p=0.001) and tended to a left-shift symmetry in O compared with L (p=0.09). Plasma E, insulin and NEFA were significantly higher in L compared to O (p≤0.04). There were no significant differences in glucose, lactate and plasma NE between groups (p≥0.2). Conclusion The study showed that O presented a lower Fatmax and a lower reliance on fat oxidation at high, but not at moderate, intensities. This may be linked to a: i) higher levels of insulin and lower E concentrations in O, which may induce blunted lipolysis; ii) higher percentage of type II and a lower percentage of type I fibres (5), and iii) decreased mitochondrial content (2), which may reduce FORs at high intensities and Fatmax. These findings may have implications for an appropriate exercise intensity prescription for optimize fat oxidation in O. References 1. Cheneviere et al. Med Sci Sports Exerc. 2009 2. Holloway et al. Am J Clin Nutr. 2009 3. Kelley et al. Am J Physiol. 1999 4. Perez-Martin et al. Diabetes Metab. 2001 5. Tanner et al. Am J Physiol Endocrinol Metab. 2002

Inhibition of glutathione efflux in the perfused rat liver and isolated hepatocytes by organic anions and bilirubin. Kinetics, sidedness, and molecular forms.

Using isolated, in situ, single-pass perfused rat livers, incubations of freshly isolated hepatocytes, and sinusoidal membrane-enriched vesicles, we and others have shown the saturability of transport (efflux) of hepatic glutathione (GSH). These observations have implicated a carrier mechanism. Our present studies were designed to provide further evidence in support of a carrier mechanism for hepatic GSH efflux by demonstrating competition by liver-specific ligands which are taken up by hepatocytes. Perfusing livers with different substances, we found that: (a) sulfobromophthalein-GSH (BSP-GSH) had a dose-dependent and fully reversible inhibitory effect on GSH efflux, while GSH alone did not have any effect; (b) taurocholate had no inhibitory effect; (c) all of the organic anions studied, i.e., BSP, rose bengal, indocyanine green, and unconjugated bilirubin (UCB), manifested potent, dose-dependent inhibitory effects, with absence of toxic effects and complete reversibility of inhibition in the case of UCB. The inhibitory effects of UCB could be overcome partially by raising (CoCl2-induced) hepatic GSH concentration. Because of the physiological importance of UCB, we conducted a detailed study of its inhibitory kinetics in the isolated hepatocyte model in the range of circulating concentrations of UCB. Studies with Cl- -free media, to inhibit the uptake of UCB by hepatocytes, showed that the inhibition of GSH efflux by UCB is apparently from inside the cell. This point was confirmed by showing that the inhibition is overcome only when bilirubin-loaded cells are cleared of bilirubin (incubation with 5% bovine serum albumin). Using Gunn rat hepatocytes and purified bilirubin mono- and diglucuronides, we found that both UCB and glucuronide forms of bilirubin inhibit GSH efflux in a dose-dependent manner. We conclude that the organic anions, although taken up by a mechanism independent of GSH, may competitively inhibit the carrier for GSH efflux from inside the hepatocyte.

Chemoconvection patterns in the Methylene-blue-glucose system: Weakly nonlinear analysis

The oxidation of solutions of glucose with methylene-blue as a catalyst in basic media can induce hydrodynamic overturning instabilities, termed chemoconvection in recognition of their similarity to convective instabilities. The phenomenon is due to gluconic acid, the marginally dense product of the reaction, which gradually builds an unstable density profile. Experiments indicate that dominant pattern wavenumbers initially increase before gradually decreasing or can even oscillate for long times. Here, we perform a weakly nonlinear analysis for an established model of the system with simple kinetics, and show that the resulting amplitude equation is analogous to that obtained in convection with insulating walls. We show that the amplitude description predicts that dominant pattern wavenumbers should decrease in the long term, but does not reproduce the aforementioned increasing wavenumber behavior in the initial stages of pattern development. We hypothesize that this is due to horizontally homogeneous steady states not being attained before pattern onset. We show that the behavior can be explained using a combination of pseudo-steady-state linear and steady-state weakly nonlinear theories. The results obtained are in qualitative agreement with the analysis of experiments.