Generating Virtual Patients by Multivariate and Discrete Re-Sampling Techniques

Does Sample Size, Sampling Strategy, or Handling of Concentrations Below the Lower Limit of Quantification Matter When Externally Evaluating Population Pharmacokinetic Models?

BACKGROUND AND OBJECTIVES

Precision dosing requires selecting the appropriate population pharmacokinetic model, which can be assessed through external evaluations (EEs). The lack of understanding of how different study design factors influence EE study outcomes makes it challenging to select the most suitable model for clinical use. This study aimed to evaluate the impact of sample size, sampling strategy, and handling of concentrations below the lower limit of quantification (BLQ) on the outcomes of EE for four population pharmacokinetic models using vancomycin and tobramycin as examples.

METHODS

Three virtual patient populations undergoing vancomycin or tobramycin therapy were simulated with varying sample size and sampling scenarios. The three approaches used to handle BLQ data were to (1) discard them, (2) impute them as LLOQ/2, or (3) use a likelihood-based approach. EEs were performed with NONMEM and R.

RESULTS

Sample size did not have an important impact on the EE results for a given scenario. Increasing the number of samples per patient did not improve predictive performance for two out of the three evaluated models. Evaluating a model developed with rich sampling did not result in better performance than those developed with regular therapeutic drug monitoring. A likelihood-based method to handle BLQ samples impacted the outcomes of the EE with lower bias for predicted troughs.

CONCLUSIONS

This study suggests that a large sample size may not be necessary for an EE study, and models selected based on TDM may be more generalizable. The study highlights the need for guidelines for EE of population pharmacokinetic models for clinical use.

Quantitative systems pharmacology modeling of macrophage-targeted therapy combined with PD-L1 inhibition in advanced NSCLC

Immune checkpoint inhibitors remained the standard-of-care treatment for advanced non-small cell lung cancer (NSCLC) for the past decade. In unselected patients, anti-PD-(L)1 monotherapy achieved an overall response rate of about 20%. In this analysis, we developed a pharmacokinetic and pharmacodynamic module for our previously calibrated quantitative systems pharmacology model (QSP) to simulate the effectiveness of macrophage-targeted therapies in combination with PD-L1 inhibition in advanced NSCLC. By conducting in silico clinical trials, the model confirmed that anti-CD47 treatment is not an optimal option of second- and later-line treatment for advanced NSCLC resistant to PD-(L)1 blockade. Furthermore, the model predicted that inhibition of macrophage recruitment, such as using CCR2 inhibitors, can potentially improve tumor size reduction when combined with anti-PD-(L)1 therapy, especially in patients who are likely to respond to anti-PD-(L)1 monotherapy and those with a high level of tumor-associated macrophages. Here, we demonstrate the application of the QSP platform on predicting the effectiveness of novel drug combinations involving immune checkpoint inhibitors based on preclinical or early-stage clinical trial data.

Envisioning the Future of Personalized Medicine: Role and Realities of Digital Twins

Digital twins have emerged as a groundbreaking concept in personalized medicine, offering immense potential to transform health care delivery and improve patient outcomes. It is important to highlight the impact of digital twins on personalized medicine across the understanding of patient health, risk assessment, clinical trials and drug development, and patient monitoring. By mirroring individual health profiles, digital twins offer unparalleled insights into patient-specific conditions, enabling more accurate risk assessments and tailored interventions. However, their application extends beyond clinical benefits, prompting significant ethical debates over data privacy, consent, and potential biases in health care. The rapid evolution of this technology necessitates a careful balancing act between innovation and ethical responsibility. As the field of personalized medicine continues to evolve, digital twins hold tremendous promise in transforming health care delivery and revolutionizing patient care. While challenges exist, the continued development and integration of digital twins hold the potential to revolutionize personalized medicine, ushering in an era of tailored treatments and improved patient well-being. Digital twins can assist in recognizing trends and indicators that might signal the presence of diseases or forecast the likelihood of developing specific medical conditions, along with the progression of such diseases. Nevertheless, the use of human digital twins gives rise to ethical dilemmas related to informed consent, data ownership, and the potential for discrimination based on health profiles. There is a critical need for robust guidelines and regulations to navigate these challenges, ensuring that the pursuit of advanced health care solutions does not compromise patient rights and well-being. This viewpoint aims to ignite a comprehensive dialogue on the responsible integration of digital twins in medicine, advocating for a future where technology serves as a cornerstone for personalized, ethical, and effective patient care.

Virtual Patient Simulation Using Copula Modeling

Virtual patient simulation is increasingly performed to support model-based optimization of clinical trial designs or individualized dosing strategies. Quantitative pharmacological models typically incorporate individual-level patient characteristics, or covariates, which enable the generation of virtual patient cohorts. The individual-level patient characteristics, or covariates, used as input for such simulations should accurately reflect the values seen in real patient populations. Current methods often make unrealistic assumptions about the correlation between patient's covariates or require direct access to actual data sets with individual-level patient data, which may often be limited by data sharing limitations. We propose and evaluate the use of copulas to address current shortcomings in simulation of patient-associated covariates for virtual patient simulations for model-based dose and trial optimization in clinical pharmacology. Copulas are multivariate distribution functions that can capture joint distributions, including the correlation, of covariate sets. We compare the performance of copulas to alternative simulation strategies, and we demonstrate their utility in several case studies. Our work demonstrates that copulas can reproduce realistic patient characteristics, both in terms of individual covariates and the dependence structure between different covariates, outperforming alternative methods, in particular when aiming to reproduce high-dimensional covariate sets. In conclusion, copulas represent a versatile and generalizable approach for virtual patient simulation which preserve relationships between covariates, and offer an open science strategy to facilitate re-use of patient data sets.

Bayesian Inference-Based Gaussian Mixture Models With Optimal Components Estimation Towards Large-Scale Synthetic Data Generation for In Silico Clinical Trials

Goal: To develop a computationally efficient and unbiased synthetic data generator for large-scale in silico clinical trials (CTs). Methods: We propose the BGMM-OCE, an extension of the conventional BGMM (Bayesian Gaussian Mixture Models) algorithm to provide unbiased estimations regarding the optimal number of Gaussian components and yield high-quality, large-scale synthetic data at reduced computational complexity. Spectral clustering with efficient eigenvalue decomposition is applied to estimate the hyperparameters of the generator. A case study is conducted to compare the performance of BGMM-OCE against four straightforward synthetic data generators for in silico CTs in hypertrophic cardiomyopathy (HCM). Results: The BGMM-OCE generated 30000 virtual patient profiles having the lowest coefficient-of-variation (0.046), inter- and intra-correlation differences (0.017, and 0.016, respectively) with the real ones in reduced execution time. Conclusions: BGMM-OCE overcomes the lack of population size in HCM which obscures the development of targeted therapies and robust risk stratification models.

Feasibility and Assessment of a Machine Learning-Based Predictive Model of Outcome After Lumbar Decompression Surgery

STUDY DESIGN

Retrospective study at a unique center.

OBJECTIVE

The aim of this study is twofold, to develop a virtual patients model for lumbar decompression surgery and to evaluate the precision of an artificial neural network (ANN) model designed to accurately predict the clinical outcomes of lumbar decompression surgery.

METHODS

We performed a retrospective study of complete Electronic Health Records (EHR) to identify potential unfavorable criteria for spine surgery (predictors). A cohort of synthetics EHR was created to classify patients by surgical success (green zone) or partial failure (orange zone) using an Artificial Neural Network which screens all the available predictors.

RESULTS

In the actual cohort, we included 60 patients, with complete EHR allowing efficient analysis, 26 patients were in the orange zone (43.4%) and 34 were in the green zone (56.6%). The average positive criteria amount for actual patients was 8.62 for the green zone (SD+/- 3.09) and 10.92 for the orange zone (SD 3.38). The classifier (a neural network) was trained using 10,000 virtual patients and 2000 virtual patients were used for test purposes. The 12,000 virtual patients were generated from the 60 EHR, of which half were in the green zone and half in the orange zone. The model showed an accuracy of 72% and a ROC score of 0.78. The sensitivity was 0.885 and the specificity 0.59.

CONCLUSION

Our method can be used to predict a favorable patient to have lumbar decompression surgery. However, there is still a need to further develop its ability to analyze patients in the "failure of treatment" zone to offer precise management of patient health before spinal surgery.

The PHARMACOM-EPI Framework for Integrating Pharmacometric Modelling Into Pharmacoepidemiological Research Using Real-World Data: Application to Assess Death Associated With Valproate

In pharmacoepidemiology, it is usually expected that the observed association should be directly or indirectly related to the pharmacological effects of the drug/s under investigation. Pharmacological effects are, in turn, strongly connected to the pharmacokinetic and pharmacodynamic properties of a drug, which can be characterized and investigated using pharmacometric models. Recently, the use of pharmacometrics has been proposed to provide pharmacological substantiation of pharmacoepidemiological findings derived from real-world data. However, validated frameworks suggesting how to combine these two disciplines for the aforementioned purpose are missing. Therefore, we propose PHARMACOM-EPI, a framework that provides a structured approach on how to identify, characterize, and apply pharmacometric models with practical details on how to choose software, format dataset, handle missing covariates/dosing data, how to perform the external evaluation of pharmacometric models in real-world data, and how to provide pharmacological substantiation of pharmacoepidemiological findings. PHARMACOM-EPI was tested in a proof-of-concept study to pharmacologically substantiate death associated with valproate use in the Danish population aged ≥ 65 years. Pharmacological substantiation of death during a follow-up period of 1 year showed that in all individuals who died (n = 169) individual predictions were within the subtherapeutic range compared with 52.8% of those who did not die (n = 1,084). Of individuals who died, 66.3% (n = 112) had a cause of death possibly related to valproate and 33.7% (n = 57) with well-defined cause of death unlikely related to valproate. This proof-of-concept study showed that PHARMACOM-EPI was able to provide pharmacological substantiation for death associated with valproate use in the study population.

Constructing a representative in-silico population for paediatric simulations: Application to HIV-positive African children

AIMS

Simulations are an essential tool for investigating scenarios in pharmacokinetics-pharmacodynamics. The models used during simulation often include the effect of highly correlated covariates such as weight, height and sex, and for children also age, which complicates the construction of an in silico population. For this reason, a suitable and representative patient population is crucial for the simulations to produce meaningful results. For simulation in paediatric patients, international growth charts from the World Health Organization (WHO) and the Centers for Disease Control and Prevention (CDC) provide a reference, but these may not always be representative for specific populations, such as malnourished children with HIV or acutely unwell children.

METHODS

We present a workflow to construct a virtual paediatric patient population using WHO and CDC growth charts, suggest piecewise linear functions to adjust the median of the growth charts by sex and age, and suggest visual diagnostics to compare with the target population. We applied this workflow in a population of 1206 HIV-positive African children, consisting of 19 742 observations with weight ranging from 3.8 to 79.7 kg, height from 55.5 to 180 cm, and an age between 0.40 and 18 years.

RESULTS

Before adjustment, the WHO and CDC charts produced weights and heights higher compared to the observed data. After applying our methodology, we could simulate weight, height, sex and age combinations in good agreement with the observed data.

CONCLUSION

The methodology presented here is flexible and may be applied to other scenarios where WHO and CDC growth standards might not be appropriate. In addition we provide R scripts and a large ready-to-use paediatric population.

A computational pipeline for data augmentation towards the improvement of disease classification and risk stratification models: A case study in two clinical domains

Virtual population generation is an emerging field in data science with numerous applications in healthcare towards the augmentation of clinical research databases with significant lack of population size. However, the impact of data augmentation on the development of AI (artificial intelligence) models to address clinical unmet needs has not yet been investigated. In this work, we assess whether the aggregation of real with virtual patient data can improve the performance of the existing risk stratification and disease classification models in two rare clinical domains, namely the primary Sjögren's Syndrome (pSS) and the hypertrophic cardiomyopathy (HCM), for the first time in the literature. To do so, multivariate approaches, such as, the multivariate normal distribution (MVND), and straightforward ones, such as, the Bayesian networks, the artificial neural networks (ANNs), and the tree ensembles are compared against their performance towards the generation of high-quality virtual data. Both boosting and bagging algorithms, such as, the Gradient boosting trees (XGBoost), the AdaBoost and the Random Forests (RFs) were trained on the augmented data to evaluate the performance improvement for lymphoma classification and HCM risk stratification. Our results revealed the favorable performance of the tree ensemble generators, in both domains, yielding virtual data with goodness-of-fit 0.021 and KL-divergence 0.029 in pSS and 0.029, 0.027 in HCM, respectively. The application of the XGBoost on the augmented data revealed an increase by 10.9% in accuracy, 10.7% in sensitivity, 11.5% in specificity for lymphoma classification and 16.1% in accuracy, 16.9% in sensitivity, 13.7% in specificity in HCM risk stratification.

Conditional distribution modeling as an alternative method for covariates simulation: Comparison with joint multivariate normal and bootstrap techniques

Clinical trial simulation (CTS) is a valuable tool in drug development. To obtain realistic scenarios, the subjects included in the CTS must be representative of the target population. Common ways of generating virtual subjects are based upon bootstrap (BS) procedures or multivariate normal distributions (MVNDs). Here, we investigated the performance of an alternative method based on conditional distributions (CDs). Covariate data from a hypertension drug development program were used. The methods were evaluated based on the original data set (internal evaluation) and on their ability to reproduce an older, unobserved population (extrapolation). Similar results were obtained in the internal evaluation for summary statistics, yet BS was able to preserve the correlation structure of the empirical distribution, which was not adequately reproduced by MVND; CD was in between BS and MVND. BS does not allow to extrapolate to an unobserved population. When the data set used to inform the extrapolation was well approximated by an MVND, the results from CD and MVND were comparable. However, improved extrapolation performance was observed for CD when deviations from normality assumptions occurred. If CTS is used to simulate within the observed distribution, BS is the preferred method. When extrapolating to new populations, a parametric method like CD/MVND is needed. In case the empirical multivariate distribution is characterized by linearly related covariates and unimodal marginal distributions, MVND can be used because of the simpler statistical framework and well‐established use; however, if uncertainty about the MVND assumptions exists, CD will increase the confidence in the simulations compared to MVND.

Complete populations of virtual patients for in silico clinical trials

MOTIVATION

Model-based approaches to safety and efficacy assessment of pharmacological drugs, treatment strategies, or medical devices (In Silico Clinical Trial, ISCT) aim to decrease time and cost for the needed experimentations, reduce animal and human testing, and enable precision medicine. Unfortunately, in presence of non-identifiable models (e.g., reaction networks), parameter estimation is not enough to generate complete populations of Virtual Patient (VPs), i.e., populations guaranteed to show the entire spectrum of model behaviours (phenotypes), thus ensuring representativeness of the trial.

RESULTS

We present methods and software based on global search driven by statistical model checking that, starting from a (non-identifiable) quantitative model of the human physiology (plus drugs PK/PD) and suitable biological and medical knowledge elicited from experts, compute a population of VPs whose behaviours are representative of the whole spectrum of phenotypes entailed by the model (completeness) and pairwise distinguishable according to user-provided criteria. This enables full granularity control on the size of the population to employ in an ISCT, guaranteeing representativeness while avoiding over-representation of behaviours.We proved the effectiveness of our algorithm on a non-identifiable ODE-based model of the female Hypothalamic-Pituitary-Gonadal axis, by generating a population of 4 830 264 VPs stratified into 7 levels (at different granularity of behaviours), and assessed its representativeness against 86 retrospective health records from Pfizer, Hannover Medical School and University Hospital of Lausanne. The datasets are respectively covered by our VPs within Average Normalised Mean Absolute Error of 15%, 20%, and 35% (90% of the latter dataset is covered within 20% error).

Translational approaches to treating dynamical diseases through in silico clinical trials

The primary goal of drug developers is to establish efficient and effective therapeutic protocols. Multifactorial pathologies, including dynamical diseases and complex disorders, can be difficult to treat, given the high degree of inter- and intra-patient variability and nonlinear physiological relationships. Quantitative approaches combining mechanistic disease modeling and computational strategies are increasingly leveraged to rationalize pre-clinical and clinical studies and to establish effective treatment strategies. The development of clinical trials has led to new computational methods that allow for large clinical data sets to be combined with pharmacokinetic and pharmacodynamic models of diseases. Here, we discuss recent progress using in silico clinical trials to explore treatments for a variety of complex diseases, ultimately demonstrating the immense utility of quantitative methods in drug development and medicine.

In-Silico Trials for Treatment of Acute Ischemic Stroke

Despite improved treatment, a large portion of patients with acute ischemic stroke due to a large vessel occlusion have poor functional outcome. Further research exploring novel treatments and better patient selection has therefore been initiated. The feasibility of new treatments and optimized patient selection are commonly tested in extensive and expensive randomized clinical trials. in-silico trials, computer-based simulation of randomized clinical trials, have been proposed to aid clinical trials. In this white paper, we present our vision and approach to set up in-silico trials focusing on treatment and selection of patients with an acute ischemic stroke. The INSIST project (IN-Silico trials for treatment of acute Ischemic STroke, www.insist-h2020.eu) is a collaboration of multiple experts in computational science, cardiovascular biology, biophysics, biomedical engineering, epidemiology, radiology, and neurology. INSIST will generate virtual populations of acute ischemic stroke patients based on anonymized data from the recent stroke trials and registry, and build on the existing and emerging in-silico models for acute ischemic stroke, its treatment (thrombolysis and thrombectomy) and the resulting perfusion changes. These models will be used to design a platform for in-silico trials that will be validated with existing data and be used to provide a proof of concept of the potential efficacy of this emerging technology. The platform will be used for preliminary evaluation of the potential suitability and safety of medication, new thrombectomy device configurations and methods to select patient subpopulations for better treatment outcome. This could allow generating, exploring and refining relavant hypotheses on potential causal pathways (which may follow from the evidence obtained from clinical trials) and improving clinical trial design. Importantly, the findings of the in-silico trials will require validation under the controlled settings of randomized clinical trials.

In silico clinical trials: concepts and early adoptions

Innovations in information and communication technology infuse all branches of science, including life sciences. Nevertheless, healthcare is historically slow in adopting technological innovation, compared with other industrial sectors. In recent years, new approaches in modelling and simulation have started to provide important insights in biomedicine, opening the way for their potential use in the reduction, refinement and partial substitution of both animal and human experimentation. In light of this evidence, the European Parliament and the United States Congress made similar recommendations to their respective regulators to allow wider use of modelling and simulation within the regulatory process. In the context of in silico medicine, the term 'in silico clinical trials' refers to the development of patient-specific models to form virtual cohorts for testing the safety and/or efficacy of new drugs and of new medical devices. Moreover, it could be envisaged that a virtual set of patients could complement a clinical trial (reducing the number of enrolled patients and improving statistical significance), and/or advise clinical decisions. This article will review the current state of in silico clinical trials and outline directions for a full-scale adoption of patient-specific modelling and simulation in the regulatory evaluation of biomedical products. In particular, we will focus on the development of vaccine therapies, which represents, in our opinion, an ideal target for this innovative approach.

Generation of virtual patient data for in-silico cardiomyopathies drug development using tree ensembles: a comparative study

In-silico clinical platforms have been recently used as a new revolutionary path for virtual patients (VP) generation and further analysis, such as, drug development. Advanced individualized models have been developed to enhance flexibility and reliability of the virtual patient cohorts. This study focuses on the implementation and comparison of three different methodologies for generating virtual data for in-silico clinical trials. Towards this direction, three computational methods, namely: (i) the multivariate log-normal distribution (log- MVND), (ii) the supervised tree ensembles, and (iii) the unsupervised tree ensembles are deployed and evaluated against their performance towards the generation of high-quality virtual data using the goodness of fit (gof) and the dataset correlation matrix as performance evaluation measures. Our results reveal the dominance of the tree ensembles towards the generation of virtual data with similar distributions (gof values less than 0.2) and correlation patterns (average difference less than 0.03).

A methodological framework for drug development in rare diseases. The CRESim program: Epilogue and perspectives

Based on the 'European Child-Rare-Euro-Simulation' (CRESim) project, this article proposes a generalizable strategy utilizing datasets analysis in combination with modeling and simulation, in order to optimize the clinical drug development applied in the field of rare diseases. The global process includes: (i) the simulation of a realistic virtual population of patients (modeled from a real dataset of patients), (ii) the modeling of disease pathophysiological components and of pharmacokinetic-pharmacodynamic relations of the drug(s) of interest, (iii) the modeling of several randomized controlled clinical trials (RCTs) designs and (iv) the analysis of the results (multi-dimensional approach for RCTs durations and precision of the estimation of the treatment effect). However, whereas modeling and numerical simulation may provide supplementary tools for drug development, they cannot be considered as a substitute for RCTs performed in 'real' patients.

Safety of sacubitril/valsartan initiated during hospitalization: data from a non-selected cohort

Aims

Sacubitril/valsartan is safe when initiated during hospitalization in a clinical trial setting. Its safety in real‐life population is not stablished. We compared the initiation of sacubitril/valsartan during hospitalization in a non‐selected population, in the PIONEER‐HF trial, and in non‐selected outpatients.

Methods and results

Multicentre registry included 527 patients: 100 were started on sacubitril/valsartan during hospitalization (19.0%) and 427 as outpatients (81.0%). Compared with those in the pivotal trial, inpatients in our cohort were older (71 ± 12 vs. 61 ± 14 years; P < 0.001); had more frequently Functional Class II (41 [41.0%] vs. 100 [22.7%]; P < 0.001), higher levels of N‐terminal pro‐B type natriuretic peptide (4044 [1630–8680] vs. 2013 [1002–4132] pg/mL; P < 0.001), better glomerular filtration rate (63.5 [51.0–80.0] vs. 58.4 [47.5–71.5] mL/min; P = 0.01), and higher systolic blood pressure (121 [110–136] vs. 118 [110–133] mmHg; P = 0.03); and received angiotensin‐converting enzyme inhibitors/angiotensin receptor blockers more frequently (92 [92.0%] vs. 208 [52.7%]; P < 0.001). Compared with non‐selected outpatients, inpatients were older (71 ± 12 vs. 68 ± 12 years, P = 0.02), had more frequent Functional Class III–IV (58 [58.0%] vs. 129 [30.3%], P < 0.001), had higher levels of N‐terminal pro‐B type natriuretic peptide (4044 [1630–8680] vs. 2182 [1134–4172]; P < 0.001), and were receiving angiotensin‐converting enzyme inhibitors/angiotensin receptor blockers target dose less frequently (55 [55.0%] vs. 335 [78.5%]; P < 0.001). They also started sacubitril/valsartan with a low dose (50 mg/12 h) more frequently (80 [80.0%] vs. 209 [48.8%], P < 0.001). The initiation of sacubitril/valsartan in outpatients was an independent predictor of high‐dose use (OR 3.1; 95% confidence interval 1.7–5.6, P < 0.001). The follow‐up time in both cohorts, including all patients enrolled, was similar (7.0 ± 0.1 vs. 7.2 ± 2.6 months, P = 0.72). All‐cause admissions during follow‐up were more frequent in inpatients (30 [30.0%] vs. 68 outpatients [15.9%], P = 0.001), with no relevant differences in all‐cause mortality. There was no significant difference in sacubitril/valsartan withdrawal rate (17 inpatients [17.0%] vs. 49 outpatients [11.5%], P = 0.13). The incidence of adverse effects was also similar: hypotension (16 inpatients [16.0%] vs. 71 outpatients [16.7%], P = 0.88), worsening renal function (7 inpatients [7.0%] vs. 29 outpatients [6.8%], P = 0.94), and hyperkalaemia (1 inpatient [1.0%] vs. 21 outpatients [4.9%], P = 0.09). We did not register any case of angioedema.

Conclusions

It is safe to initiate sacubitril/valsartan during hospitalization in daily clinical practice. Inpatients have a higher risk profile and receive low starting doses more frequently than outpatients.

Optimizing drug development in oncology by clinical trial simulation: Why and how?

In therapeutic research, the safety and efficacy of pharmaceutical products are necessarily tested on humans via clinical trials after an extensive and expensive preclinical development period. Methodologies such as computer modeling and clinical trial simulation (CTS) might represent a valuable option to reduce animal and human assays. The relevance of these methods is well recognized in pharmacokinetics and pharmacodynamics from the preclinical phase to postmarketing. However, they are barely used and are poorly regarded for drug approval, despite Food and Drug Administration and European Medicines Agency recommendations. The generalization of CTS could be greatly facilitated by the availability of software for modeling biological systems, by clinical trial studies and hospital databases. Data sharing and data merging raise legal, policy and technical issues that will need to be addressed. Development of future molecules will have to use CTS for faster development and thus enable better patient management. Drug activity modeling coupled with disease modeling, optimal use of medical data and increased computing speed should allow this leap forward. The realization of CTS requires not only bioinformatics tools to allow interconnection and global integration of all clinical data but also a universal legal framework to protect the privacy of every patient. While recognizing that CTS can never replace 'real-life' trials, they should be implemented in future drug development schemes to provide quantitative support for decision-making. This in silico medicine opens the way to the P4 medicine: predictive, preventive, personalized and participatory.

The Hepato-Hypothalamic-Pituitary-Adrenal-Renal Axis: Mathematical Modeling of Cortisol’s Production, Metabolism, and Seasonal Variation

Cortisol dynamics are governed by the integration of influences from the suprachiasmatic nucleus (SCN), the hypothalamic-pituitary-adrenal (HPA) axis, and metabolic enzymes, such as the 11β–hydroxysteroid dehydrogenase (HSD) family, which are highly expressed in hepatic and renal tissue. The coordinated regulation of cortisol dynamics is essential for the maintenance of a healthy state, and aberrant cortisol circadian rhythms are associated with various pathophysiological conditions. The duration of the light-dark cycle, or photoperiod, which regulates SCN activity, varies seasonally, and the shorter photoperiod winter season is associated with elevated cortisol levels, peak inflammatory disease incidence, and symptom exacerbation. Elevated expression and activity of 11β-HSD1 protein, assumed to also occur during the winter, have been allied with numerous inflammatory conditions. A comprehensive understanding of the communication between the underlying regulatory mechanisms of cortisol as well as how changes in their activity could lead to the development of disease is yet to be elucidated. In this work, we propose the use of a semimechanistic mathematical model to explore the impact of the hepato-hypothalamic-pituitary-adrenal-renal axis in modulating neuroendocrine-immune system dynamics. Our model predicts the predominance of a winter proinflammatory state and that genetic variations could alter 11β-HSD enzyme functionality, rendering certain subpopulations more susceptible to disease as a consequence of HPA axis dysregulation.

Evaluating optimal therapy robustness by virtual expansion of a sample population, with a case study in cancer immunotherapy

Cancer is a highly heterogeneous disease, exhibiting spatial and temporal variations that pose challenges for designing robust therapies. Here, we propose the VEPART (Virtual Expansion of Populations for Analyzing Robustness of Therapies) technique as a platform that integrates experimental data, mathematical modeling, and statistical analyses for identifying robust optimal treatment protocols. VEPART begins with time course experimental data for a sample population, and a mathematical model fit to aggregate data from that sample population. Using nonparametric statistics, the sample population is amplified and used to create a large number of virtual populations. At the final step of VEPART, robustness is assessed by identifying and analyzing the optimal therapy (perhaps restricted to a set of clinically realizable protocols) across each virtual population. As proof of concept, we have applied the VEPART method to study the robustness of treatment response in a mouse model of melanoma subject to treatment with immunostimulatory oncolytic viruses and dendritic cell vaccines. Our analysis (i) showed that every scheduling variant of the experimentally used treatment protocol is fragile (nonrobust) and (ii) discovered an alternative region of dosing space (lower oncolytic virus dose, higher dendritic cell dose) for which a robust optimal protocol exists.