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 challenge of modeling nuclear receptor regulatory networks in mammalian cells.

Nuclear receptors are a major component of signal transduction in animals. They mediate the regulatory activities of many hormones, nutrients and metabolites on the homeostasis and physiology of cells and tissues. It is of high interest to model the corresponding regulatory networks. While molecular and cell biology studies of individual promoters have provided important mechanistic insight, a more complex picture is emerging from genome-wide studies. The regulatory circuitry of nuclear receptor regulated gene expression networks, and their response to cellular signaling, appear highly dynamic, and involve long as well as short range chromatin interactions. We review how progress in understanding the kinetics and regulation of cofactor recruitment, and the development of new genomic methods, provide opportunities but also a major challenge for modeling nuclear receptor mediated regulatory networks.

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.

Valeur prédictive de la surexpression/ amplification de Her2/neu pour un traitement ciblé du cancer du sein [Predictive value of Her2/neu expression/amplification for the targeted treatment of breast cancer]

Her2/neu is a tyrosine kinase receptor which stimulates cell growth. The receptor is overexpressed in about 20% of breast cancers. Her2/neu expression is an indicator of poor prognosis but also the target of the treatment of breast cancer using humanised anti-Her2/ neu antibodies. Only cancers overexpressing the protein will respond to this therapy, but which has significant (cardiac) side effects and is expensive. It is therefore important to test for the overexpression of the protein on breast cancer cells. This paper discusses how this can be done and ongoing research into new therapeutic options targeting the involved signaling pathways.

Testing the predictive adaptive response in a host-parasite system

1. Harsh environmental conditions experienced during development can reduce the performance of the same individuals in adulthood. However, the 'predictive adaptive response' hypothesis postulates that if individuals adapt their phenotype during development to the environments where they are likely to live in the future, individuals exposed to harsh conditions in early life perform better when encountering the same harsh conditions in adulthood compared to those never exposed to these conditions before. 2. Using the common vole (Microtus arvalis) as study organism, we tested how exposure to flea parasitism during the juvenile stage affects the physiology (haematocrit, resistance to oxidative stress, resting metabolism, spleen mass, and testosterone), morphology (body mass, testis mass) and motor performance (open field activity and swimming speed) of the same individuals when infested with fleas in adulthood. According to the 'predictive adaptive response' hypothesis, we predicted that voles parasitized at the adult stage would perform better if they had already been parasitized with fleas at the juvenile stage. 3. We found that voles exposed to fleas in adulthood had a higher metabolic rate if already exposed to fleas when juvenile, compared to voles free of fleas when juvenile and voles free of fleas in adulthood. Independently of juvenile parasitism, adult parasitism impaired adult haematocrit and motor performances. Independently of adult parasitism, juvenile parasitism slowed down crawling speed in adult female voles. 4. Our results suggest that juvenile parasitism has long-term effects that do not protect from the detrimental effects of adult parasitism. On the contrary, experiencing parasitism in early-life incurs additional costs upon adult parasitism measured in terms of higher energy expenditure, rather than inducing an adaptive shift in the developmental trajectory. 5. Hence, our study provides experimental evidence for long term costs of parasitism. We found no support for a predictive adaptive response in this host-parasite system.

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.

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.

Predictive factors of adrenal insufficiency in patients admitted to acute medical wards: a case control study.

BACKGROUND: Adrenal insufficiency is a rare and potentially lethal disease if untreated. Several clinical signs and biological markers are associated with glucocorticoid failure but the importance of these factors for diagnosing adrenal insufficiency is not known. In this study, we aimed to assess the prevalence of and the factors associated with adrenal insufficiency among patients admitted to an acute internal medicine ward. METHODS: Retrospective, case-control study including all patients with high-dose (250 μg) ACTH-stimulation tests for suspected adrenal insufficiency performed between 2008 and 2010 in an acute internal medicine ward (n = 281). Cortisol values <550 nmol/l upon ACTH-stimulation test were considered diagnostic for adrenal insufficiency. Area under the ROC curve (AROC), sensitivity, specificity, negative and positive predictive values for adrenal insufficiency were assessed for thirteen symptoms, signs and biological variables. RESULTS: 32 patients (11.4%) presented adrenal insufficiency; the others served as controls. Among all clinical and biological parameters studied, history of glucocorticoid withdrawal was the only independent factor significantly associated with patients with adrenal insufficiency (Odds Ratio: 6.71, 95% CI: 3.08 -14.62). Using a logistic regression, a model with four significant and independent variable was obtained, regrouping history of glucocorticoid withdrawal (OR 7.38, 95% CI [3.18 ; 17.11], p-value <0.001), nausea (OR 3.37, 95% CI [1.03 ; 11.00], p-value 0.044), eosinophilia (OR 17.6, 95% CI [1.02; 302.3], p-value 0.048) and hyperkalemia (OR 2.41, 95% CI [0.87; 6.69], p-value 0.092). The AROC (95% CI) was 0.75 (0.70; 0.80) for this model, with 6.3 (0.8 - 20.8) for sensitivity and 99.2 (97.1 - 99.9) for specificity. CONCLUSIONS: 11.4% of patients with suspected adrenal insufficient admitted to acute medical ward actually do present with adrenal insufficiency, defined by an abnormal response to high-dose (250 μg) ACTH-stimulation test. A history of glucocorticoid withdrawal was the strongest factor predicting the potential adrenal failure. The combination of a history of glucocorticoid withdrawal, nausea, eosinophilia and hyperkaliemia might be of interest to suspect adrenal insufficiency.

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.

Gene expression profiling of human glioblastomas for the identification of predictive factors and characterisation of molecular mechanism implicated in response to therapy and outcome

ABSTRACT Poor outcome for glioblastoma patients is largely due to resistance to chemoradiation therapy. While epigenetic inactivation of MGMT mediated DNA repair is highly predictive for benefit from the alkylating agent therapy Temozolomide, additional mechanisms for resistance associated with molecular alterations exist. Furthermore, new concepts in cancer suggest that resistance to treatment may be linked to cancer stem cells that escape therapy and act as source for tumour recurrence. We determined gene expression signatures associated with outcome in glioblastoma patients enrolled in a phase II and phase III clinical trial establishing the new combination therapy of radiation plus concomitant and adjuvant Temozolomide. Correlating stable gene clusters emerging from unsupervised analysis with survival of 42 treated patients identified a number of biological processes associated with outcome. Most prominent, a gene cluster dominated by HOX genes and comprising PROM1, was associated with resistance. PROM1 encodes CD133, a marker for a subpopulation of tumour cells enriched for glioblastoma stem- like cells. The core of this correlated HOX cluster was comprised in the top genes of a "self-renewal signature" defined in a mouse model for MLL-AF9 initiated leukaemia. The association of the HOX gene cluster with tumour resistance was confirmed in two external data sets of 146 malignant glioma As additional resistance factors we identified over-expression of the epidermal growth factor receptor gene, EGFR, while increased gene expression related to biological features of tumour host interaction, including markers for tumour vascular and cell adhesion, and innate immune response, were associated with better outcome. The "self-renewal" signature associated with resistance to the new combination chemoradiation therapy provides first clinical evidence that glioma stem like cells may implicated in resistance in a uniformly treated cohort of glioblastoma patients. This study underlines the need to target the tumour stem cell compartment, and provides some testable hypothesis for biological mechanisms relevant for malignant behaviour of glioblastoma that may be targeted in new treatment approaches. Résumé Le glioblastome, tumeur cérébrale primaire maligne la plus fréquente, est connue pour son mauvais pronostique. Des avancées chimiothérapeutiques récentes avec des agents alkylants comme le témozolomide (TMZ), ont permis une amélioration notable dans la survie de certains patients. Les bénéficiaires ont la caractéristique commune de présenter une particularité génétique, la methylation du MGMT (methylguanine methyltransferase). Néanmoins, d'autres mécanismes de résistance en fonction des aberrations moléculaires existent. Nous avons établi les profils d'expressions génétiques des patients traités par irradiation et TMZ dans des études cliniques de phase II et III. En combinant des méthodes non-supervisées et supervisées, de l'étude de la cohorte des patients traités nous avons découvert des groupes de gènes associés à la survie. Un ensemble de gènes contenant les gènes Hox semble lié au mécanisme de résistance au traitement. Récemment, les gènes Hox ont été décrits comme faisant partie d"une signature d'autorenouvellement (self-renewal) des cellules souches cancéreuses de la leucémie. L'autorenouvellement est un processus grâce auquel les cellules souches se maintiennent tout au long de la vie. Cette association à la résistance est confirmée dans deux autres études indépendantes. Un autre facteur de résistance au traitement est la surexpression du gène EGFR. D'autre part, deux groupes de gènes associés à la relation entre hôte-tumeur tels que les marqueurs des vaisseaux tumoraux et de la réponse immunitaire innée s'avèrent avoir un effet positif sur la survie des patients traités. La découverte de la signature d'autorenouvellement comme facteur de résistance à la nouvelle chimio-radiothérapie offre une preuve clinique que les cellules souches cancéreuses sont impliquées dans la résistance au traitement. If est donc logique de penser que le traitement ciblé contre des cellules souches cancéreuses va dans l'avenir permettre des thérapies anticancéreuses plus performantes.

Modeling the forensic two-trace problem with Bayesian networks

The forensic two-trace problem is a perplexing inference problem introduced by Evett (J Forensic Sci Soc 27:375-381, 1987). Different possible ways of wording the competing pair of propositions (i.e., one proposition advanced by the prosecution and one proposition advanced by the defence) led to different quantifications of the value of the evidence (Meester and Sjerps in Biometrics 59:727-732, 2003). Here, we re-examine this scenario with the aim of clarifying the interrelationships that exist between the different solutions, and in this way, produce a global vision of the problem. We propose to investigate the different expressions for evaluating the value of the evidence by using a graphical approach, i.e. Bayesian networks, to model the rationale behind each of the proposed solutions and the assumptions made on the unknown parameters in this problem.