Purpose: The androgen receptor axis inhibitors (ARPI; e.g., enzalutamide, abiraterone acetate) are administered in daily practice for men with metastatic castration-resistant prostate cancer (mCRPC). However, not all patients respond, and mechanisms of both primary and acquired resistance remain largely unknown. Experimental Design: In the prospective trial MATCH-R (NCT02517892), 59 patients with mCRPC underwent whole-exome sequencing (WES) and/or RNA sequencing (RNA-seq) of samples collected before starting ARPI. Also, 18 patients with mCRPC underwent biopsy at time of resistance. The objectives were to identify genomic alterations associated with resistance to ARPIs as well as to describe clonal evolution. Associations of genomic and transcriptomic alterations with primary resistance were determined using Wilcoxon and Fisher exact tests. Results: WES analysis indicated that no single-gene genomic alterations were strongly associated with primary resistance. RNA-seq analysis showed that androgen receptor (AR) gene alterations and expression levels were similar between responders and nonresponders. RNA-based pathway analysis found that patients with primary resistance had a higher Hedgehog pathway score, a lower AR pathway score and a lower NOTCH pathway score than patients with a response. Subclonal evolution and acquisition of new alterations in AR-related genes or neuroendocrine differentiation are associated with acquired resistance. ARPIs do not induce significant changes in the tumor transcriptome of most patients; however, programs associated with cell proliferation are enriched in resistant samples. Conclusions: Low AR activity, activation of stemness programs, and Hedgehog pathway were associated with primary ARPIs’ resistance, whereas most acquired resistance was associated with subclonal evolution, AR-related events, and neuroendocrine differentiation. See related commentary by Slovin, p. 4323

Background: Current diabetes education methods are costly, time-consuming, and do not actively engage the patient. Here, we describe the development and verification of the physiological model for healthy subjects that forms the basis of the Eindhoven Diabetes Education Simulator (E-DES). E-DES shall provide diabetes patients with an individualized virtual practice environment incorporating the main factors that influence glycemic control: food, exercise, and medication. Method: The physiological model consists of 4 compartments for which the inflow and outflow of glucose and insulin are calculated using 6 nonlinear coupled differential equations and 14 parameters. These parameters are estimated on 12 sets of oral glucose tolerance test (OGTT) data (226 healthy subjects) obtained from literature. The resulting parameter set is verified on 8 separate literature OGTT data sets (229 subjects). The model is considered verified if 95% of the glucose data points lie within an acceptance range of ±20% of the corresponding model value. Results: All glucose data points of the verification data sets lie within the predefined acceptance range. Physiological processes represented in the model include insulin resistance and β-cell function. Adjusting the corresponding parameters allows to describe heterogeneity in the data and shows the capabilities of this model for individualization. Conclusion: We have verified the physiological model of the E-DES for healthy subjects. Heterogeneity of the data has successfully been modeled by adjusting the 4 parameters describing insulin resistance and β-cell function. Our model will form the basis of a simulator providing individualized education on glucose control.

Imbalanced immune responses are a prominent hallmark of cancer and autoimmunity. Myeloid cells can be overly sup-pressive, inhibiting protective immune responses or inactive not controlling autoreactive immune cells. Understanding the mechanisms that induce suppressive myeloid cells, such as myeloid-derived suppressor cells (MDSCs) and tolerogenic dendritic cells (TolDCs), can facilitate the development of immune-restoring therapeutic approaches. MDSCs are a major barrier for effective cancer immunotherapy by suppressing antitumor immune responses in cancer patients. TolDCs are administered to patients to promote immune tolerance with the intent to control autoimmune disease. Here, we investi-gated the development and suppressive/tolerogenic activity of human MDSCs and TolDCs to gain insight into signaling pathways that drive immunosuppression in these different myeloid subsets. Moreover, monocyte-derived MDSCs (M-MDSCs) generated in vitro were compared to M-MDSCs iso-lated from head-and-neck squamous cell carcinoma patients. PI3K-AKT signaling was identified as being crucial for the induction of human M-MDSCs. PI3K inhibition prevented the downregulation of HLA-DR and the upregulation of reactive oxygen species and MerTK. In addition, we show that the suppressive activity of dexamethasone-induced TolDCs is induced by beta-catenin-dependent Wnt signaling. The identifi- cation of PI3K-AKT and Wnt signal transduction pathways as respective inducers of the immunomodulatory capacity of M-MDSCs and TolDCs provides opportunities to overcome sup-pressive myeloid cells in cancer patients and optimize thera-peutic application of TolDCs. Lastly, the observed similarities between generated-and patient-derived M-MDSCs support the use of in vitro-generated M-MDSCs as powerful model to investigate the functionality of human MDSCs.

We are currently in the middle of a major shift in biomedical research: unprecedented and rapidly growing amounts of data may be obtained today, from in vitro, in vivo and clinical studies, at molecular, physiological and clinical levels. To make use of these large-scale, multi-level datasets, corresponding multi-level mathematical models are needed, i.e. models that simultaneously capture multiple layers of the biological, physiological and disease-level organization (also referred to as quantitative systems pharmacology—QSP—models). However, today's multi-level models are not yet embedded in end-usage applications, neither in drug research and development nor in the clinic. Given the expectations and claims made historically, this seemingly slow adoption may seem surprising. Therefore, we herein consider a specific example—type 2 diabetes—and critically review the current status and identify key remaining steps for these models to become mainstream in the future. This overview reveals how, today, we may use models to ask scientific questions concerning, e.g., the cellular origin of insulin resistance, and how this translates to the whole-body level and short-term meal responses. However, before these multi-level models can become truly useful, they need to be linked with the capabilities of other important existing models, in order to make them ‘personalized’ (e.g. specific to certain patient phenotypes) and capable of describing long-term disease progression. To be useful in drug development, it is also critical that the developed models and their underlying data and assumptions are easily accessible. For clinical end-usage, in addition, model links to decision-support systems combined with the engagement of other disciplines are needed to create user-friendly and cost-efficient software packages.

The Metabolic Syndrome (MetS) is a complex, multifactorial disorder that develops slowly over time presenting itself with large differences among MetS patients. We applied a systems biology approach to describe and predict the onset and progressive development of MetS, in a study that combined in vivo and in silico models. A new data-driven, physiological model (MINGLeD: Model INtegrating Glucose and Lipid Dynamics) was developed, describing glucose, lipid and cholesterol metabolism. Since classic kinetic models cannot describe slowly progressing disorders, a simulation method (ADAPT) was used to describe longitudinal dynamics and to predict metabolic concentrations and fluxes. This approach yielded a novel model that can describe long-term MetS development and progression. This model was integrated with longitudinal in vivo data that was obtained from male APOE*3-Leiden.CETP mice fed a high-fat, high-cholesterol diet for three months and that developed MetS as reflected by classical symptoms including obesity and glucose intolerance. Two distinct subgroups were identified: those who developed dyslipidemia, and those who did not. The combination of MINGLeD with ADAPT could correctly predict both phenotypes, without making any prior assumptions about changes in kinetic rates or metabolic regulation. Modeling and flux trajectory analysis revealed that differences in liver fluxes and dietary cholesterol absorption could explain this occurrence of the two different phenotypes. In individual mice with dyslipidemia dietary cholesterol absorption and hepatic turnover of metabolites, including lipid fluxes, were higher compared to those without dyslipidemia. Predicted differences were also observed in gene expression data, and consistent with the emergence of insulin resistance and hepatic steatosis, two well-known MetS co-morbidities. Whereas MINGLeD specifically models the metabolic derangements underlying MetS, the simulation method ADAPT is generic and can be applied to other diseases where dynamic modeling and longitudinal data are available. Obesity is becoming a growing health problem, with the risk of running into a disease state called the Metabolic Syndrome (MetS). It is difficult to study MetS in humans as it develops over a long period of time while many processes and organs are involved. Moreover, each patient may present itself with different characteristics. We combined data from animal experiments with a computer model and computer simulations to study the development of MetS. The new model correctly predicts which changes in the underlying metabolic processes could explain the MetS symptoms. Unexpectedly, two different subgroups were identified: those with high cholesterol and high triglycerides, and those without. The computer model found that in those who develop lipid abnormalities, both dietary cholesterol absorption and hepatic liver fluxes were higher.