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.

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)]

IL28B expression depends on a novel TT/-G polymorphism which improves HCV clearance prediction.

Approximately 3% of the world population is chronically infected with the hepatitis C virus (HCV), with potential development of cirrhosis and hepatocellular carcinoma. Despite the availability of new antiviral agents, treatment remains suboptimal. Genome-wide association studies (GWAS) identified rs12979860, a polymorphism nearby IL28B, as an important predictor of HCV clearance. We report the identification of a novel TT/-G polymorphism in the CpG region upstream of IL28B, which is a better predictor of HCV clearance than rs12979860. By using peripheral blood mononuclear cells (PBMCs) from individuals carrying different allelic combinations of the TT/-G and rs12979860 polymorphisms, we show that induction of IL28B and IFN-γ-inducible protein 10 (IP-10) mRNA relies on TT/-G, but not rs12979860, making TT/-G the only functional variant identified so far. This novel step in understanding the genetic regulation of IL28B may have important implications for clinical practice, as the use of TT/G genotyping instead of rs12979860 would improve patient management.

Prediction of clinical response after renal angioplasty: respective value of renal Doppler sonography and scintigraphy.

OBJECTIVE: The goal of our study was to compare Doppler sonography and renal scintigraphy as tools for predicting the therapeutic response in patients after undergoing renal angioplasty. SUBJECTS AND METHODS. Seventy-four hypertensive patients underwent clinical examination, Doppler sonography, and renal scintigraphy before and after receiving captopril in preparation for renal revascularization. The patients were evaluated for the status of hypertension 3 months after the procedure. The predictive values of the findings of clinical examination, Doppler sonography, renal scintigraphy, and angiography were assessed. RESULTS: For prediction of a favorable therapeutic outcome, abnormal results from renal scintigraphy before and after captopril administration had a sensitivity of 58% and specificity of 57%. Findings of Doppler sonography had a sensitivity of 68% and specificity of 50% before captopril administration and a sensitivity of 81% and specificity of 32% after captopril administration. Significant predictors of a cure or reduction of hypertension after revascularization were low unilateral (p = 0.014) and bilateral resistive (p = 0.016) indexes on Doppler sonography before (p = 0.009) and after (p = 0.028) captopril administration. On multivariate analysis, the best predictors were a unilateral resistive index of less than 0.65 (odds ratio [OR] = 3.7) after captopril administration and a kidney longer than 93 mm (OR = 7.8). The two best combined criteria to predict the favorable therapeutic outcome were a bilateral resistive index of less than 0.75 before captopril administration combined with a unilateral resistive index of less than 0.70 after captopril administration (sensitivity, 76%; specificity, 58%) or a bilateral resistive index of less than 0.75 before captopril administration and a kidney measuring longer than 90 mm (sensitivity, 81%; specificity, 50%). CONCLUSION: Measurements of kidney length and unilateral and bilateral resistive indexes before and after captopril administration were useful in predicting the outcome after renal angioplasty. Renal scintigraphy had no significant predictive value.

Response Prediction in Chronic Hepatitis C by Assessment of IP-10 and IL28B-Related Single Nucleotide Polymorphisms.

BACKGROUND: High baseline levels of IP-10 predict a slower first phase decline in HCV RNA and a poor outcome following interferon/ribavirin therapy in patients with chronic hepatitis C. Several recent studies report that single nucleotide polymorphisms (SNPs) adjacent to IL28B predict spontaneous resolution of HCV infection and outcome of treatment among HCV genotype 1 infected patients. METHODS AND FINDINGS: In the present study, we correlated the occurrence of variants at three such SNPs (rs12979860, rs12980275, and rs8099917) with pretreatment plasma IP-10 and HCV RNA throughout therapy within a phase III treatment trial (HCV-DITTO) involving 253 Caucasian patients. The favorable SNP variants (CC, AA, and TT, respectively) were associated with lower baseline IP-10 (P = 0.02, P = 0.01, P = 0.04) and were less common among HCV genotype 1 infected patients than genotype 2/3 (P<0.0001, P<0.0001, and P = 0.01). Patients carrying favorable SNP genotypes had higher baseline viral load than those carrying unfavorable variants (P = 0.0013, P = 0.029, P = 0.0004 respectively). Among HCV genotype 1 infected carriers of the favorable C, A, or T alleles, IP-10 below 150 pg/mL significantly predicted a more pronounced reduction of HCV RNA from day 0 to 4 (first phase decline), which translated into increased rates of RVR (62%, 53%, and 39%) and SVR (85%, 76%, and 75% respectively) among homozygous carriers with baseline IP-10 below 150 pg/mL. In multivariate analyses of genotype 1-infected patients, baseline IP-10 and C genotype at rs12979860 independently predicted the first phase viral decline and RVR, which in turn independently predicted SVR. CONCLUSIONS: Concomitant assessment of pretreatment IP-10 and IL28B-related SNPs augments the prediction of the first phase decline in HCV RNA, RVR, and final therapeutic outcome.

Prediction of myocardial infarction size using the SYNTAX score in patients treated with primary percutaneous coronary intervention for acute ST- segment elevation myocardial infarction

Objectives The relevance of the SYNTAX score for the particular case of patients with acute ST- segment elevation myocardial infarction (STEMI) undergoing primary percutaneous coronary intervention (PPCI)  has previously only been studied in the setting of post hoc analysis of large prospective randomized clinical trials. A "real-life" population approach has never been explored before. The aim of this study was to evaluate the impact of the SYNTAX score for the prediction of the myocardial infarction size, estimated by the creatin-kinase (CK) peak value, using the SYNTAX score in patients treated with primary coronary intervention for acute ST-segment elevation myocardial infarction. Methods The primary endpoint of the study was myocardial infarction size as measured by the CK peak value. The SYNTAX score was calculated retrospectively in 253 consecutive patients with acute ST-segment elevation myocardial infarction (STEMI) undergoing primary percutaneous coronary intervention (PPCI) in a large tertiary referral center in Switzerland, between January 2009 and June 2010. Linear regression analysis was performed to compare myocardial infarction size with the SYNTAX score. This same endpoint was then stratified according to SYNTAX score tertiles: low <22 (n=178), intermediate [22-32] (n=60), and high >=33 (n=15). Results There were no significant differences in terms of clinical characteristics between the three groups. When stratified according to the SYNTAX score tertiles, average CK peak values of 1985 (low<22), 3336 (intermediate [22-32]) and 3684 (high>=33) were obtained with a p-value <0.0001. Bartlett's test for equal variances between the three groups was 9.999 (p-value <0.0067). A moderate Pearson product-moment correlation coefficient (r=0.4074) with a high statistical significance level (p-value <0.0001) was found. The coefficient of determination (R^2=0.1660) showed that approximately 17% of the variation of CK peak value (myocardial infarction size) could be explained by the SYNTAX score, i.e. by the coronary disease complexity. Conclusion In an all-comers population, the SYNTAX score is an additional tool in predicting myocardial infarction size in patients treated with primary percutaneous coronary intervention (PPCI). The stratification of patients in different risk groups according to SYNTAX enables to identify a high-risk population that may warrant particular patient care.

Implementation study of a clinical prediction score to rule out cardiogenic chest pain in community: Cluster randomized trial

1.1 Fundamentals Chest pain is a common complaint in primary care patients (1 to 3% of all consultations) (1) and its aetiology can be miscellaneous, from harmless to potentially life threatening conditions. In primary care practice, the most prevalent aetiologies are: chest wall syndrome (43%), coronary heart disease (12%) and anxiety (7%) (2). In up to 20% of cases, potentially serious conditions as cardiac, respiratory or neoplasic diseases underlie chest pain. In this context, a large number of laboratory tests are run (42%) and over 16% of patients are referred to a specialist or hospitalized (2).¦A cardiovascular origin to chest pain can threaten patient's life and investigations run to exclude a serious condition can be expensive and involve a large number of exams or referral to specialist -­‐ often without real clinical need. In emergency settings, up to 80% of chest pains in patients are due to cardiovascular events (3) and scoring methods have been developed to identify conditions such as coronary heart disease (HD) quickly and efficiently (4-­‐6). In primary care, a cardiovascular origin is present in only about 12% of patients with chest pain (2) and general practitioners (GPs) need to exclude as safely as possible a potential serious condition underlying chest pain. A simple clinical prediction rule (CPR) like those available in emergency settings may therefore help GPs and spare time and extra investigations in ruling out CHD in primary care patients. Such a tool may also help GPs reassure patients with more common origin to chest pain.

Markers of atherosclerosis and inflammation for prediction of coronary heart disease in older adults.

Although both inflammatory and atherosclerosis markers have been associated with coronary heart disease (CHD) risk, data directly comparing their predictive value are limited. The authors compared the value of 2 atherosclerosis markers (ankle-arm index (AAI) and aortic pulse wave velocity (aPWV)) and 3 inflammatory markers (C-reactive protein (CRP), interleukin-6 (IL-6), and tumor necrosis factor-alpha (TNF-alpha)) in predicting CHD events. Among 2,191 adults aged 70-79 years at baseline (1997-1998) from the Health, Aging, and Body Composition Study cohort, the authors examined adjudicated incident myocardial infarction or CHD death ("hard" events) and "hard" events plus hospitalization for angina or coronary revascularization (total CHD events). During 8 years of follow-up between 1997-1998 and June 2007, 351 participants developed total CHD events (197 "hard" events). IL-6 (highest quartile vs. lowest: hazard ratio = 1.82, 95% confidence interval: 1.33, 2.49; P-trend < 0.001) and AAI (AAI </= 0.9 vs. AAI 1.01-1.30: hazard ratio = 1.57, 95% confidence interval: 1.14, 2.18) predicted CHD events above traditional risk factors and modestly improved global measures of predictive accuracy. CRP, TNF-alpha, and aPWV had weaker associations. IL-6 and AAI accurately reclassified 6.6% and 3.3% of participants, respectively (P's </= 0.05). Results were similar for "hard" CHD, with higher reclassification rates for AAI. IL-6 and AAI are associated with future CHD events beyond traditional risk factors and modestly improve risk prediction in older adults.

The renal hemodynamic effects of ibuprofen in the newborn rabbit.

In early childhood, nonsteroidal anti-inflammatory drugs are mainly used to either prevent or treat premature labor of the mother and patent ductus arteriosus of the newborn infant. The most frequently used prostaglandin-synthesis inhibitor is indomethacin. Fetuses exposed to indomethacin in utero have been born with renal developmental defects, and in both the unborn child and the term and premature newborn this drug may compromise renal glomerular function. The latter has in the past also been observed when i.v. indomethacin or i.v. acetylsalicylic acid (aspirin) were administered to newborn rabbits. The present experiments were designed to evaluate whether ibuprofen has less renal side effects than indomethacin, as claimed. Three groups of anesthetized, ventilated, normoxemic neonatal rabbits were infused with increasing doses of ibuprofen (0.02, 0.2, 2.0 mg/kg body weight) and the following renal parameters were measured: urine volume, urinary sodium excretion, GFR, and renal plasma flow. Renal blood flow, filtration fraction, and the renal vascular resistance were calculated according to standard formulae. Intravenous ibuprofen caused a dose-dependent, significant reduction in urine volume, GFR, and renal blood flow with a fall in filtration fraction in the animals receiving the highest dose of ibuprofen (2 mg/kg body weight). There was a very steep rise in renal vascular resistance. Urinary sodium excretion decreased. These experiments in neonatal rabbits clearly show that acute i.v. doses of ibuprofen also have significant renal hemodynamic and functional side effects, not less than seen previously with indomethacin.

Robust response transformations based on optimal prediction

Nonlinear regression problems can often be reduced to linearity by transforming the response variable (e.g., using the Box-Cox family of transformations). The classic estimates of the parameter defining the transformation as well as of the regression coefficients are based on the maximum likelihood criterion, assuming homoscedastic normal errors for the transformed response. These estimates are nonrobust in the presence of outliers and can be inconsistent when the errors are nonnormal or heteroscedastic. This article proposes new robust estimates that are consistent and asymptotically normal for any unimodal and homoscedastic error distribution. For this purpose, a robust version of conditional expectation is introduced for which the prediction mean squared error is replaced with an M scale. This concept is then used to develop a nonparametric criterion to estimate the transformation parameter as well as the regression coefficients. A finite sample estimate of this criterion based on a robust version of smearing is also proposed. Monte Carlo experiments show that the new estimates compare favorably with respect to the available competitors.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Prediction of postnatal outcome in fetuses with pulmonary hypoplasia by fetal-MRI: Comparison between lung volumetry and pulmonary signal intensity values

Purpose: Pulmonary hypoplasia is a determinant parameter for extra-uterine life. In the last years, MRI appears as a complement to US in order to evaluate the degree of pulmonary hypoplasia in foetuses with congenital anomalies, by using different methods - fetal lung volumetry (FLV), lung-to-liver signal intensity ratio (LLSIR)-. But until now, information about the correlation between the MRI prediction and the real postnatal outcome is limited. Methods and materials: We retrospectively reviewed the fetal MRI performed at our Institution in the last 8 years and selected the cases with suspicion of fetal pulmonary hypoplasia (n = 30). The pulmonary volumetry data of these foetuses were collected and the lung-to-liver signal intensity ratio (LLSIR) measures performed. These data were compared with those obtained from a control group of 25 foetuses considered as normal at MRI. The data of the study group were also correlated with the autopsy records or the post-natal clinical information of the patients. Results: As expected, the control group showed higher FLV and LLSIR values than the problem group at all gestational ages. Higher values of FLV and LLSIR were associated with a better post-natal outcome. Sensitivity, specificity, positive and negative predictive values and accuracy for the relative LLSIR and the relative FLV showed no significant differences. Conclusion: Our data show that not only the FLV but also the relative LLSIR inform about the degree of fetal lung development. This information may help to predict the fetal outcome and to evaluate the need for neonatal intensive care.