Propagation velocity kinetics and repolarization alternans in a free-behaving sheep model of pacing-induced atrial fibrillation.

AIMS: Experimental models have reported conflicting results regarding the role of dispersion of repolarization in promoting atrial fibrillation (AF). Repolarization alternans, a beat-to-beat alternation in action potential duration, enhances dispersion of repolarization when propagation velocity is involved. METHODS AND RESULTS: In this work, original electrophysiological parameters were analysed to study AF susceptibility in a chronic sheep model of pacing-induced AF. Two pacemakers were implanted, each with a single right atrial lead. Right atrial depolarization and repolarization waves were documented at 2-week intervals. A significant and gradual decrease in the propagation velocity at all pacing rates and in the right atrial effective refractory period (ERP) was observed during the weeks of burst pacing before sustained AF developed when compared with baseline conditions. Right atrial repolarization alternans was observed, but because of the development of 2/1 atrioventricular block with far-field ventricular interference, its threshold could not be precisely measured. Non-sustained AF was not observed at baseline, but appeared during the electrical remodelling in association with a decrease in both ERP and propagation velocity. CONCLUSION: We report here on the feasibility of measuring ERP, atrial repolarization alternans, and propagation velocity kinetics and their potential in predicting susceptibility to AF in a free-behaving model of pacing-induced AF using the standard pacemaker technology.

Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility

Background: Germline genetic variation is associated with the differential expression of many human genes. The phenotypic effects of this type of variation may be important when considering susceptibility to common genetic diseases. Three regions at 8q24 have recently been identified to independently confer risk of prostate cancer. Variation at 8q24 has also recently been associated with risk of breast and colorectal cancer. However, none of the risk variants map at or relatively close to known genes, with c-MYC mapping a few hundred kilobases distally. Results: This study identifies cis-regulators of germline c-MYC expression in immortalized lymphocytes of HapMap individuals. Quantitative analysis of c-MYC expression in normal prostate tissues suggests an association between overexpression and variants in Region 1 of prostate cancer risk. Somatic c-MYC overexpression correlates with prostate cancer progression and more aggressive tumor forms, which was also a pathological variable associated with Region 1. Expression profiling analysis and modeling of transcriptional regulatory networks predicts a functional association between MYC and the prostate tumor suppressor KLF6. Analysis of MYC/Myc-driven cell transformation and tumorigenesis substantiates a model in which MYC overexpression promotes transformation by down-regulating KLF6. In this model, a feedback loop through E-cadherin down-regulation causes further transactivation of c-MYC.Conclusion: This study proposes that variation at putative 8q24 cis-regulator(s) of transcription can significantly alter germline c-MYC expression levels and, thus, contribute to prostate cancer susceptibility by down-regulating the prostate tumor suppressor KLF6 gene.

Personalized modeling for drug concentration prediction using Support Vector Machine

Building a personalized model to describe the drug concentration inside the human body for each patient is highly important to the clinical practice and demanding to the modeling tools. Instead of using traditional explicit methods, in this paper we propose a machine learning approach to describe the relation between the drug concentration and patients' features. Machine learning has been largely applied to analyze data in various domains, but it is still new to personalized medicine, especially dose individualization. We focus mainly on the prediction of the drug concentrations as well as the analysis of different features' influence. Models are built based on Support Vector Machine and the prediction results are compared with the traditional analytical models.

Student modeling based on fuzzy inference mechanisms

The paper presents a competence-based instructional design system and a way to provide a personalization of navigation in the course content. The navigation aid tool builds on the competence graph and the student model, which includes the elements of uncertainty in the assessment of students. An individualized navigation graph is constructed for each student, suggesting the competences the student is more prepared to study. We use fuzzy set theory for dealing with uncertainty. The marks of the assessment tests are transformed into linguistic terms and used for assigning values to linguistic variables. For each competence, the level of difficulty and the level of knowing its prerequisites are calculated based on the assessment marks. Using these linguistic variables and approximate reasoning (fuzzy IF-THEN rules), a crisp category is assigned to each competence regarding its level of recommendation.

Prediction of infant drug exposure through breastfeeding: population PK modeling and simulation of fluoxetine

BACKGROUND: Risks of significant infant drug exposurethrough breastmilk are poorly defined for many drugs, and largescalepopulation data are lacking. We used population pharmacokinetics(PK) modeling to predict fluoxetine exposure levels ofinfants via mother's milk in a simulated population of 1000 motherinfantpairs.METHODS: Using our original data on fluoxetine PK of 25breastfeeding women, a population PK model was developed withNONMEM and parameters, including milk concentrations, wereestimated. An exponential distribution model was used to account forindividual variation. Simulation random and distribution-constrainedassignment of doses, dosing time, feeding intervals and milk volumewas conducted to generate 1000 mother-infant pairs with characteristicssuch as the steady-state serum concentrations (Css) and infantdose relative to the maternal weight-adjusted dose (relative infantdose: RID). Full bioavailability and a conservative point estimate of1-month-old infant CYP2D6 activity to be 20% of the adult value(adjusted by weigth) according to a recent study, were assumed forinfant Css calculations.RESULTS: A linear 2-compartment model was selected as thebest model. Derived parameters, including milk-to-plasma ratios(mean: 0.66; SD: 0.34; range, 0 - 1.1) were consistent with the valuesreported in the literature. The estimated RID was below 10% in >95%of infants. The model predicted median infant-mother Css ratio was0.096 (range 0.035 - 0.25); literature reported mean was 0.07 (range0-0.59). Moreover, the predicted incidence of infant-mother Css ratioof >0.2 was less than 1%.CONCLUSION: Our in silico model prediction is consistent withclinical observations, suggesting that substantial systemic fluoxetineexposure in infants through human milk is rare, but further analysisshould include active metabolites. Our approach may be valid forother drugs. [supported by CIHR and Swiss National Science Foundation(SNSF)]

The impact of fibrosis and steatosis on early viral kinetics in HCV genotype 1-infected patients treated with Peg-IFN-alfa-2a and ribavirin.

Summary.  Hepatitis C viral (HCV) kinetics after initiation of interferon-based therapy provide valuable insights for understanding virus pathogenesis, evaluating treatment antiviral effectiveness and predicting treatment outcome. Adverse effects of liver fibrosis and steatosis on sustained virological response have been frequently reported, yet their impacts on the early viral kinetics remain unclear. In this study, associations between histology status and early viral kinetics were assessed in 149 HCV genotype 1-infected patients treated with pegylated interferon alfa-2a and ribavirin (DITTO trial). In multivariate analyses adjusted for critical factors such as IL28B genotype and baseline viral load, presence of significant fibrosis (Ishak stage > 2) was found to independently reduce the odds of achieving an initial reduction (calculated from day 0 to day 4) in HCV RNA of ≥2 logIU/mL (adjusted OR 0.03, P = 0.004) but was not associated with the second-phase slope of viral decline (calculated from day 8 to day 29). On the contrary, presence of liver steatosis was an independent risk factor for not having a rapid second-phase slope, that is, ≥0.3 logIU/mL/week (adjusted OR 0.22, P = 0.012) but was not associated with the first-phase decline. Viral kinetic modelling theory suggests that significant fibrosis primarily impairs the treatment antiviral effectiveness in blocking viral production by infected cells, whereas the presence of steatosis is associated with a lower net loss of infected cells. Further studies will be necessary to identify the biological mechanisms underlain by these findings.

Aerodynamic drag modeling of alpine skiers performing giant slalom turns.

PURPOSE: Aerodynamic drag plays an important role in performance for athletes practicing sports that involve high-velocity motions. In giant slalom, the skier is continuously changing his/her body posture, and this affects the energy dissipated in aerodynamic drag. It is therefore important to quantify this energy to understand the dynamic behavior of the skier. The aims of this study were to model the aerodynamic drag of alpine skiers in giant slalom simulated conditions and to apply these models in a field experiment to estimate energy dissipated through aerodynamic drag. METHODS: The aerodynamic characteristics of 15 recreational male and female skiers were measured in a wind tunnel while holding nine different skiing-specific postures. The drag and the frontal area were recorded simultaneously for each posture. Four generalized and two individualized models of the drag coefficient were built, using different sets of parameters. These models were subsequently applied in a field study designed to compare the aerodynamic energy losses between a dynamic and a compact skiing technique. RESULTS: The generalized models estimated aerodynamic drag with an accuracy of between 11.00% and 14.28%, and the individualized models estimated aerodynamic drag with an accuracy between 4.52% and 5.30%. The individualized model used for the field study showed that using a dynamic technique led to 10% more aerodynamic drag energy loss than using a compact technique. DISCUSSION: The individualized models were capable of discriminating different techniques performed by advanced skiers and seemed more accurate than the generalized models. The models presented here offer a simple yet accurate method to estimate the aerodynamic drag acting upon alpine skiers while rapidly moving through the range of positions typical to turning technique.

Oxygen uptake kinetics and middle distance swimming performance.

OBJECTIVE: The aim of this study was to determine whether V˙O(2) kinetics and specifically, the time constant of transitions from rest to heavy (τ(p)H) and severe (τ(p)S) exercise intensities, are related to middle distance swimming performance. DESIGN: Fourteen highly trained male swimmers (mean ± SD: 20.5 ± 3.0 yr; 75.4 ± 12.4 kg; 1.80 ± 0.07 m) performed an discontinuous incremental test, as well as square wave transitions for heavy and severe swimming intensities, to determine V˙O(2) kinetics parameters using two exponential functions. METHODS: All the tests involved front-crawl swimming with breath-by-breath analysis using the Aquatrainer swimming snorkel. Endurance performance was recorded as the time taken to complete a 400 m freestyle swim within an official competition (T400), one month from the date of the other tests. RESULTS: T400 (Mean ± SD) (251.4 ± 12.4 s) was significantly correlated with τ(p)H (15.8 ± 4.8s; r=0.62; p=0.02) and τ(p)S (15.8 ± 4.7s; r=0.61; p=0.02). The best single predictor of 400 m freestyle time, out of the variables that were assessed, was the velocity at V˙O(2max)vV˙O(2max), which accounted for 80% of the variation in performance between swimmers. However, τ(p)H and V˙O(2max) were also found to influence the prediction of T400 when they were included in a regression model that involved respiratory parameters only. CONCLUSIONS: Faster kinetics during the primary phase of the V˙O(2) response is associated with better performance during middle-distance swimming. However, vV˙O(2max) appears to be a better predictor of T400.

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.

New Trends in Public Facility Location Modeling

The past four decades have witnessed an explosive growth in the field of networkbased facility location modeling. This is not at all surprising since location policy is one of the most profitable areas of applied systems analysis in regional science and ample theoretical and applied challenges are offered. Location-allocation models seek the location of facilities and/or services (e.g., schools, hospitals, and warehouses) so as to optimize one or several objectives generally related to the efficiency of the system or to the allocation of resources. This paper concerns the location of facilities or services in discrete space or networks, that are related to the public sector, such as emergency services (ambulances, fire stations, and police units), school systems and postal facilities. The paper is structured as follows: first, we will focus on public facility location models that use some type of coverage criterion, with special emphasis in emergency services. The second section will examine models based on the P-Median problem and some of the issues faced by planners when implementing this formulation in real world locational decisions. Finally, the last section will examine new trends in public sector facility location modeling.

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

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

GeoKernels: Modeling of Spatial Data on GeoManifolds

This paper presents a review of methodology for semi-supervised modeling with kernel methods, when the manifold assumption is guaranteed to be satisfied. It concerns environmental data modeling on natural manifolds, such as complex topographies of the mountainous regions, where environmental processes are highly influenced by the relief. These relations, possibly regionalized and nonlinear, can be modeled from data with machine learning using the digital elevation models in semi-supervised kernel methods. The range of the tools and methodological issues discussed in the study includes feature selection and semisupervised Support Vector algorithms. The real case study devoted to data-driven modeling of meteorological fields illustrates the discussed approach.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Making sense of AMPA receptor trafficking by modeling molecular mechanisms of synaptic plasticity.

Synaptic plasticity involves a complex molecular machinery with various protein interactions but it is not yet clear how its components give rise to the different aspects of synaptic plasticity. Here we ask whether it is possible to mathematically model synaptic plasticity by making use of known substances only. We present a model of a multistable biochemical reaction system and use it to simulate the plasticity of synaptic transmission in long-term potentiation (LTP) or long-term depression (LTD) after repeated excitation of the synapse. According to our model, we can distinguish between two phases: first, a "viscosity" phase after the first excitation, the effects of which like the activation of NMDA receptors and CaMKII fade out in the absence of further excitations. Second, a "plasticity" phase actuated by an identical subsequent excitation that follows after a short time interval and causes the temporarily altered concentrations of AMPA subunits in the postsynaptic membrane to be stabilized. We show that positive feedback is the crucial element in the core chemical reaction, i.e. the activation of the short-tail AMPA subunit by NEM-sensitive factor, which allows generating multiple stable equilibria. Three stable equilibria are related to LTP, LTD and a third unfixed state called ACTIVE. Our mathematical approach shows that modeling synaptic multistability is possible by making use of known substances like NMDA and AMPA receptors, NEM-sensitive factor, glutamate, CaMKII and brain-derived neurotrophic factor. Furthermore, we could show that the heteromeric combination of short- and long-tail AMPA receptor subunits fulfills the function of a memory tag.