In silico design of clinical trials: A method coming of age

Theoretical evaluation of the impact of diverse treatment conditions by calculation of the tumor control probability (TCP) of simulated cervical cancer Hyperthermia-Radiotherapy (HT-RT) treatments in-silico

INTRODUCTION

Hyperthermia (HT) induces various cellular biological processes, such as repair impairment and direct HT cell killing. In this context, in-silico biophysical models that translate deviations in the treatment conditions into clinical outcome variations may be used to study the extent of such processes and their influence on combined hyperthermia plus radiotherapy (HT + RT) treatments under varying conditions.

METHODS

An extended linear-quadratic model calibrated for SiHa and HeLa cell lines (cervical cancer) was used to theoretically study the impact of varying HT treatment conditions on radiosensitization and direct HT cell killing effect. Simulated patients were generated to compute the Tumor Control Probability (TCP) under different HT conditions (number of HT sessions, temperature and time interval), which were randomly selected within margins based on reported patient data.

RESULTS

Under the studied conditions, model-based simulations suggested a treatment improvement with a total CEM43 thermal dose of approximately 10 min. Additionally, for a given thermal dose, TCP increased with the number of HT sessions. Furthermore, in the simulations, we showed that the TCP dependence on the temperature/time interval is more correlated with the mean value than with the minimum/maximum value and that comparing the treatment outcome with the mean temperature can be an excellent strategy for studying the time interval effect.

CONCLUSION

The use of thermoradiobiological models allows us to theoretically study the impact of varying thermal conditions on HT + RT treatment outcomes. This approach can be used to optimize HT treatments, design clinical trials, and interpret patient data.

On the role of benchmarking data sets and simulations in method comparison studies

Method comparisons are essential to provide recommendations and guidance for applied researchers, who often have to choose from a plethora of available approaches. While many comparisons exist in the literature, these are often not neutral but favor a novel method. Apart from the choice of design and a proper reporting of the findings, there are different approaches concerning the underlying data for such method comparison studies. Most manuscripts on statistical methodology rely on simulation studies and provide a single real-world data set as an example to motivate and illustrate the methodology investigated. In the context of supervised learning, in contrast, methods are often evaluated using so-called benchmarking data sets, that is, real-world data that serve as gold standard in the community. Simulation studies, on the other hand, are much less common in this context. The aim of this paper is to investigate differences and similarities between these approaches, to discuss their advantages and disadvantages, and ultimately to develop new approaches to the evaluation of methods picking the best of both worlds. To this aim, we borrow ideas from different contexts such as mixed methods research and Clinical Scenario Evaluation.

Towards systems immunology of critical illness at scale: from single cell 'omics to digital twins

Single-cell 'omics methodology has yielded unprecedented insights based largely on data-centric informatics for reducing, and thus interpreting, massive datasets. In parallel, parsimonious mathematical modeling based on abstractions of pathobiology has also yielded major insights into inflammation and immunity, with these models being extended to describe multi-organ disease pathophysiology as the basis of 'digital twins' and in silico clinical trials. The integration of these distinct methods at scale can drive both basic and translational advances, especially in the context of critical illness, including diseases such as COVID-19. Here, I explore achievements and argue the challenges that are inherent to the integration of data-driven and mechanistic modeling approaches, highlighting the potential of modeling-based strategies for rational immune system reprogramming.

Embracing complexity in sepsis

Sepsis involves the dynamic interplay between a pathogen, the host response, the failure of organ systems, medical interventions and a myriad of other factors. This together results in a complex, dynamic and dysregulated state that has remained ungovernable thus far. While it is generally accepted that sepsis is very complex indeed, the concepts, approaches and methods that are necessary to understand this complexity remain underappreciated. In this perspective we view sepsis through the lens of complexity theory. We describe the concepts that support viewing sepsis as a state of a highly complex, non-linear and spatio-dynamic system. We argue that methods from the field of complex systems are pivotal for a fuller understanding of sepsis, and we highlight the progress that has been made over the last decades in this respect. Still, despite these considerable advancements, methods like computational modelling and network-based analyses continue to fly under the general scientific radar. We discuss what barriers contribute to this disconnect, and what we can do to embrace complexity with regards to measurements, research approaches and clinical applications. Specifically, we advocate a focus on longitudinal, more continuous biological data collection in sepsis. Understanding the complexity of sepsis will require a huge multidisciplinary effort, in which computational approaches derived from complex systems science must be supported by, and integrated with, biological data. Such integration could finetune computational models, guide validation experiments, and identify key pathways that could be targeted to modulate the system to the benefit of the host. We offer an example for immunological predictive modelling, which may inform agile trials that could be adjusted throughout the trajectory of disease. Overall, we argue that we should expand our current mental frameworks of sepsis, and embrace nonlinear, system-based thinking in order to move the field forward.

An immune dysfunction score for stratification of patients with acute infection based on whole-blood gene expression

Dysregulated host responses to infection can lead to organ dysfunction and sepsis, causing millions of global deaths each year. To alleviate this burden, improved prognostication and biomarkers of response are urgently needed. We investigated the use of whole-blood transcriptomics for stratification of patients with severe infection by integrating data from 3149 samples from patients with sepsis due to community-acquired pneumonia or fecal peritonitis admitted to intensive care and healthy individuals into a gene expression reference map. We used this map to derive a quantitative sepsis response signature (SRSq) score reflective of immune dysfunction and predictive of clinical outcomes, which can be estimated using a 7- or 12-gene signature. Last, we built a machine learning framework, SepstratifieR, to deploy SRSq in adult and pediatric bacterial and viral sepsis, H1N1 influenza, and COVID-19, demonstrating clinically relevant stratification across diseases and revealing some of the physiological alterations linking immune dysregulation to mortality. Our method enables early identification of individuals with dysfunctional immune profiles, bringing us closer to precision medicine in infection.

Integrated computational and in vivo models reveal Key Insights into macrophage behavior during bone healing

Myeloid-derived monocyte and macrophages are key cells in the bone that contribute to remodeling and injury repair. However, their temporal polarization status and control of bone-resorbing osteoclasts and bone-forming osteoblasts responses is largely unknown. In this study, we focused on two aspects of monocyte/macrophage dynamics and polarization states over time: 1) the injury-triggered pro- and anti-inflammatory monocytes/macrophages temporal profiles, 2) the contributions of pro- versus anti-inflammatory monocytes/macrophages in coordinating healing response. Bone healing is a complex multicellular dynamic process. While traditional in vitro and in vivo experimentation may capture the behavior of select populations with high resolution, they cannot simultaneously track the behavior of multiple populations. To address this, we have used an integrated coupled ordinary differential equations (ODEs)-based framework describing multiple cellular species to in vivo bone injury data in order to identify and test various hypotheses regarding bone cell populations dynamics. Our approach allowed us to infer several biological insights including, but not limited to,: 1) anti-inflammatory macrophages are key for early osteoclast inhibition and pro-inflammatory macrophage suppression, 2) pro-inflammatory macrophages are involved in osteoclast bone resorptive activity, whereas osteoblasts promote osteoclast differentiation, 3) Pro-inflammatory monocytes/macrophages rise during two expansion waves, which can be explained by the anti-inflammatory macrophages-mediated inhibition phase between the two waves. In addition, we further tested the robustness of the mathematical model by comparing simulation results to an independent experimental dataset. Taken together, this novel comprehensive mathematical framework allowed us to identify biological mechanisms that best recapitulate bone injury data and that explain the coupled cellular population dynamics involved in the process. Furthermore, our hypothesis testing methodology could be used in other contexts to decipher mechanisms in complex multicellular processes.

Myeloid-derived monocytes/macrophages are key cells for bone remodeling and injury repair. However, their temporal polarization status and control of bone-resorbing osteoclasts and bone-forming osteoblasts responses is largely unknown. In this study, we focused on two aspects of monocyte/macrophage population dynamics: 1) the injury-triggered pro- and anti-inflammatory monocytes/macrophages temporal profiles, 2) the contributions of pro- versus anti-inflammatory monocytes/macrophages in coordinating healing response. In order to test various hypotheses regarding bone cell populations dynamics, we have integrated a coupled ordinary differential equations-based framework describing multiple cellular species to in vivo bone injury data. Our approach allowed us to infer several biological insights including: 1) anti-inflammatory macrophages are key for early osteoclast inhibition and pro-inflammatory macrophage suppression, 2) pro-inflammatory macrophages are involved in osteoclast bone resorptive activity, whereas osteoblasts promote osteoclast differentiation, 3) Pro-inflammatory monocytes/macrophages rise during two expansion waves, which can be explained by the anti-inflammatory macrophages-mediated inhibition phase between the two waves. Taken together, this mathematical framework allowed us to identify biological mechanisms that recapitulate bone injury data and that explain the coupled cellular population dynamics involved in the process.

Unifying cardiovascular modelling with deep reinforcement learning for uncertainty aware control of sepsis treatment

Sepsis is a potentially life-threatening inflammatory response to infection or severe tissue damage. It has a highly variable clinical course, requiring constant monitoring of the patient's state to guide the management of intravenous fluids and vasopressors, among other interventions. Despite decades of research, there's still debate among experts on optimal treatment. Here, we combine for the first time, distributional deep reinforcement learning with mechanistic physiological models to find personalized sepsis treatment strategies. Our method handles partial observability by leveraging known cardiovascular physiology, introducing a novel physiology-driven recurrent autoencoder, and quantifies the uncertainty of its own results. Moreover, we introduce a framework for uncertainty-aware decision support with humans in the loop. We show that our method learns physiologically explainable, robust policies, that are consistent with clinical knowledge. Further our method consistently identifies high-risk states that lead to death, which could potentially benefit from more frequent vasopressor administration, providing valuable guidance for future research.

An artificial intelligence framework integrating longitudinal electronic health records with real-world data enables continuous pan-cancer prognostication

Despite widespread adoption of electronic health records (EHRs), most hospitals are not ready to implement data science research in the clinical pipelines. Here, we develop MEDomics, a continuously learning infrastructure through which multimodal health data are systematically organized and data quality is assessed with the goal of applying artificial intelligence for individual prognosis. Using this framework, currently composed of thousands of individuals with cancer and millions of data points over a decade of data recording, we demonstrate prognostic utility of this framework in oncology. As proof of concept, we report an analysis using this infrastructure, which identified the Framingham risk score to be robustly associated with mortality among individuals with early-stage and advanced-stage cancer, a potentially actionable finding from a real-world cohort of individuals with cancer. Finally, we show how natural language processing (NLP) of medical notes could be used to continuously update estimates of prognosis as a given individual's disease course unfolds.

A data-driven model of the role of energy in sepsis

Exposure to pathogens elicits a complex immune response involving multiple interdependent pathways. This response may mitigate detrimental effects and restore health but, if imbalanced, can lead to negative outcomes including sepsis. This complexity and need for balance pose a challenge for clinicians and have attracted attention from modelers seeking to apply computational tools to guide therapeutic approaches. In this work, we address a shortcoming of such past efforts by incorporating the dynamics of energy production and consumption into a computational model of the acute immune response. With this addition, we performed fits of model dynamics to data obtained from non-human primates exposed to Escherichia coli. Our analysis identifies parameters that may be crucial in determining survival outcomes and also highlights energy-related factors that modulate the immune response across baseline and altered glucose conditions.

Divergent COVID-19 Disease Trajectories Predicted by a DAMP-Centered Immune Network Model

COVID-19 presentations range from mild to moderate through severe disease but also manifest with persistent illness or viral recrudescence. We hypothesized that the spectrum of COVID-19 disease manifestations was a consequence of SARS-CoV-2-mediated delay in the pathogen-associated molecular pattern (PAMP) response, including dampened type I interferon signaling, thereby shifting the balance of the immune response to be dominated by damage-associated molecular pattern (DAMP) signaling. To test the hypothesis, we constructed a parsimonious mechanistic mathematical model. After calibration of the model for initial viral load and then by varying a few key parameters, we show that the core model generates four distinct viral load, immune response and associated disease trajectories termed “patient archetypes”, whose temporal dynamics are reflected in clinical data from hospitalized COVID-19 patients. The model also accounts for responses to corticosteroid therapy and predicts that vaccine-induced neutralizing antibodies and cellular memory will be protective, including from severe COVID-19 disease. This generalizable modeling framework could be used to analyze protective and pathogenic immune responses to diverse viral infections.

Strategies to identify candidate repurposable drugs: COVID-19 treatment as a case example

Drug repurposing is an invaluable strategy to identify new uses for existing drug therapies that overcome many of the time and financial costs associated with novel drug development. The COVID-19 pandemic has driven an unprecedented surge in the development and use of bioinformatic tools to identify candidate repurposable drugs. Using COVID-19 as a case study, we discuss examples of machine-learning and signature-based approaches that have been adapted to rapidly identify candidate drugs. The Library of Integrated Network-based Signatures (LINCS) and Connectivity Map (CMap) are commonly used repositories and have the advantage of being amenable to use by scientists with limited bioinformatic training. Next, we discuss how these recent advances in bioinformatic drug repurposing approaches might be adapted to identify repurposable drugs for CNS disorders. As the development of novel therapies that successfully target the cause of neuropsychiatric and neurological disorders has stalled, there is a pressing need for innovative strategies to treat these complex brain disorders. Bioinformatic approaches to identify repurposable drugs provide an exciting avenue of research that offer promise for improved treatments for CNS disorders.

Methods to Develop an in silico Clinical Trial: Computational Head-to-Head Comparison of Lisdexamfetamine and Methylphenidate

Regulatory agencies encourage computer modeling and simulation to reduce the time and cost of clinical trials. Although still not classified in formal guidelines, system biology-based models represent a powerful tool for generating hypotheses with great molecular detail. Herein, we have applied a mechanistic head-to-head in silico clinical trial (ISCT) between two treatments for attention-deficit/hyperactivity disorder, to wit lisdexamfetamine (LDX) and methylphenidate (MPH). The ISCT was generated through three phases comprising (i) the molecular characterization of drugs and pathologies, (ii) the generation of adult and children virtual populations (vPOPs) totaling 2,600 individuals and the creation of physiologically based pharmacokinetic (PBPK) and quantitative systems pharmacology (QSP) models, and (iii) data analysis with artificial intelligence methods. The characteristics of our vPOPs were in close agreement with real reference populations extracted from clinical trials, as did our PBPK models with in vivo parameters. The mechanisms of action of LDX and MPH were obtained from QSP models combining PBPK modeling of dosing schemes and systems biology-based modeling technology, i.e., therapeutic performance mapping system. The step-by-step process described here to undertake a head-to-head ISCT would allow obtaining mechanistic conclusions that could be extrapolated or used for predictions to a certain extent at the clinical level. Altogether, these computational techniques are proven an excellent tool for hypothesis-generation and would help reach a personalized medicine.

Alternative strategies in cardiac preclinical research and new clinical trial formats

An efficient and safe drug development process is crucial for the establishment of new drugs on the market aiming to increase quality of life and life-span of our patients. Despite technological advances in the past decade, successful launches of drug candidates per year remain low. We here give an overview about some of these advances and suggest improvements for implementation to boost preclinical and clinical drug development with a focus on the cardiovascular field. We highlight advantages and disadvantages of animal experimentation and thoroughly review alternatives in the field of three-dimensional cell culture as well as preclinical use of spheroids and organoids. Microfluidic devices and their potential as organ-on-a-chip systems, as well as the use of living animal and human cardiac tissues are additionally introduced. In the second part, we examine recent gold standard randomized clinical trials and present possible modifications to increase lead candidate throughput: adaptive designs, master protocols, and drug repurposing. In silico and N-of-1 trials have the potential to redefine clinical drug candidate evaluation. Finally, we briefly discuss clinical trial designs during pandemic times.

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.

A Dynamic Quantitative Systems Pharmacology Model of Inflammatory Bowel Disease: Part 1 - Model Framework

A mechanistic, multistate, mathematical model of inflammatory bowel disease (IBD) was developed by including key biological mechanisms in blood and gut, including cell differentiation, cytokine production, and clinical biomarkers. The model structure is consistent between healthy volunteers and IBD disease phenotype, with 24 parameters changed between diseases. Modular nature of the model allows for easy incorporation of new mechanisms or modification of existing interactions. Model simulations for steady-state levels of proteins and cells in the blood and gut using a population approach are consistent with published data. By simulating the response of two clinical biomarkers, C-reactive protein and fecal calprotectin, to parameter perturbations, the model explores hypotheses for possible treatment mechanisms. With additional experimental validation and addition of drug treatments, the model provides a platform to test hypothesis on treatment effects in IBD.

Dynamical modeling of pro- and anti-inflammatory cytokines in the early stage of septic shock

A dynamical model of the pathophysiological behaviors of IL18 and IL10 cytokines with their receptors is tested against data for the case of early sepsis. The proposed approach considers the surroundings (organs and bone marrow) and the different subsystems (cells and cyctokines). The interactions between blood cells, cytokines and the surroundings are described via mass balances. Cytokines are adsorbed onto associated receptors at the cell surface. The adsorption is described by the Langmuir model and gives rise to the production of more cytokines and associated receptors inside the cell. The quantities of pro and anti-inflammatory cytokines present in the body are combined to give global information via an inflammation level function which describes the patient’s state. Data for parameter estimation comes from the Sepsis 48 H database. Comparisons between patient data and simulations are presented and are in good agreement. For the IL18/IL10 cytokine pair, 5 key parameters have been found. They are linked to pro-inflammatory IL18 cytokine and show that the early sepsis is driven by components of inflammatory character.

Creation and application of virtual patient cohorts of heart models

Patient-specific cardiac models are now being used to guide therapies. The increased use of patient-specific cardiac simulations in clinical care will give rise to the development of virtual cohorts of cardiac models. These cohorts will allow cardiac simulations to capture and quantify inter-patient variability. However, the development of virtual cohorts of cardiac models will require the transformation of cardiac modelling from small numbers of bespoke models to robust and rapid workflows that can create large numbers of models. In this review, we describe the state of the art in virtual cohorts of cardiac models, the process of creating virtual cohorts of cardiac models, and how to generate the individual cohort member models, followed by a discussion of the potential and future applications of virtual cohorts of cardiac models. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Refining rodent models of spinal cord injury

This report was produced by an Expert Working Group (EWG) consisting of UK-based researchers, veterinarians and regulators of animal experiments with specialist knowledge of the use of animal models of spinal cord injury (SCI). It aims to facilitate the implementation of the Three Rs (Replacement, Reduction and Refinement), with an emphasis on refinement. Specific animal welfare issues were identified and discussed, and practical measures proposed, with the aim of reducing animal use and suffering, reducing experimental variability, and increasing translatability within this critically important research field.

Complex Stability and an Irrevertible Transition Reverted by Peptide and Fibroblasts in a Dynamic Model of Innate Immunity

We here apply a control analysis and various types of stability analysis to an in silico model of innate immunity that addresses the management of inflammation by a therapeutic peptide. Motivation is the observation, both in silico and in experiments, that this therapy is not robust. Our modeling results demonstrate how (1) the biological phenomena of acute and chronic modes of inflammation may reflect an inherently complex bistability with an irrevertible flip between the two modes, (2) the chronic mode of the model has stable, sometimes unique, steady states, while its acute-mode steady states are stable but not unique, (3) as witnessed by TNF levels, acute inflammation is controlled by multiple processes, whereas its chronic-mode inflammation is only controlled by TNF synthesis and washout, (4) only when the antigen load is close to the acute mode's flipping point, many processes impact very strongly on cells and cytokines, (5) there is no antigen exposure level below which reduction of the antigen load alone initiates a flip back to the acute mode, and (6) adding healthy fibroblasts makes the transition from acute to chronic inflammation revertible, although (7) there is a window of antigen load where such a therapy cannot be effective. This suggests that triple therapies may be essential to overcome chronic inflammation. These may comprise (1) anti-immunoglobulin light chain peptides, (2) a temporarily reduced antigen load, and (3a) fibroblast repopulation or (3b) stem cell strategies.

Towards a Physiological Scale of Vocal Fold Agent-Based Models of Surgical Injury and Repair: Sensitivity Analysis, Calibration and Verification

Agent based models (ABM) were developed to numerically simulate the biological response to surgical vocal fold injury and repair at the physiological level. This study aimed to improve the representation of existing ABM through a combination of empirical and computational experiments. Empirical data of vocal fold cell populations including neutrophils, macrophages and fibroblasts were obtained using flow cytometry up to four weeks following surgical injury. Random Forests were used as a sensitivity analysis method to identify model parameters that were most influential to ABM outputs. Statistical Parameter Optimization Tool for Python was used to calibrate those parameter values to match the ABM-simulation data with the corresponding empirical data from Day 1 to Day 5 following surgery. Model performance was evaluated by verifying if the empirical data fell within the 95% confidence intervals of ABM outputs of cell quantities at Day 7, Week 2 and Week 4. For Day 7, all empirical data were within the ABM output ranges. The trends of ABM-simulated cell populations were also qualitatively comparable to those of the empirical data beyond Day 7. Exact values, however, fell outside of the 95% statistical confidence intervals. Parameters related to fibroblast proliferation were indicative to the ABM-simulation of fibroblast dynamics in final stages of wound healing.

A conceptual time window-based model for the early stratification of trauma patients

Progress in the testing of therapies targeting the immune response following trauma, a leading cause of morbidity and mortality worldwide, has been slow. We propose that the design of interventional trials in trauma would benefit from a scheme or platform that could support the identification and implementation of prognostic strategies for patient stratification. Here, we propose a stratification scheme based on defined time periods or windows following the traumatic event. This 'time-window' model allows for the incorporation of prognostic variables ranging from circulating biomarkers and clinical data to patient-specific information such as gene variants to predict adverse short- or long-term outcomes. A number of circulating biomarkers, including cell injury markers and damage-associated molecular patterns (DAMPs), and inflammatory mediators have been shown to correlate with adverse outcomes after trauma. Likewise, several single nucleotide polymorphisms (SNPs) associate with complications or death in trauma patients. This review summarizes the status of our understanding of the prognostic value of these classes of variables in predicting outcomes in trauma patients. Strategies for the incorporation of these prognostic variables into schemes designed to stratify trauma patients, such as our time-window model, are also discussed.

Inflammation and Disease: Modelling and Modulation of the Inflammatory Response to Alleviate Critical Illness

Critical illness, a constellation of interrelated inflammatory and physiological derangements occurring subsequent to severe infection or injury, affects a large number of individuals in both developed and developing countries. The prototypical complex system embodied in critical illness has largely defied therapy beyond supportive care. We have focused on the utility of data-driven and mechanistic computational modelling to help address the complexity of critical illness and provide pathways towards discovering potential therapeutic options and combinations. Herein, we review recent progress in this field, with a focus on both animal and computational models of critical illness. We suggest that therapy for critical illness can be posed as a model-based dynamic control problem, and discuss novel theoretical and experimental approaches involving biohybrid devices aimed at reprogramming inflammation dynamically. Together, these advances offer the potential for Model-based Precision Medicine for critical illness.

Arguing for Adaptive Clinical Trials in Sepsis

Sepsis is life-threatening organ dysfunction due to dysregulated response to infection. Patients with sepsis exhibit wide heterogeneity stemming from genetic, molecular, and clinical factors as well as differences in pathogens, creating challenges for the development of effective treatments. Several gaps in knowledge also contribute: (i) biomarkers that identify patients likely to benefit from specific treatments are unknown; (ii) therapeutic dose and duration is often poorly understood; and (iii) short-term mortality, a common outcome measure, is frequently criticized for being insensitive. To date, the majority of sepsis trials use traditional design features, and have largely failed to identify new treatments with incremental benefit over standard of care. Traditional trials are also frequently conducted as part of a drug evaluation process that is segmented into several phases, each requiring separate trials, with a long time delay from inception through design and execution to incorporation of results into clinical practice. By contrast, adaptive clinical trial designs facilitate the evaluation of several candidate treatments simultaneously, learn from emergent discoveries during the course of the trial, and can be structured efficiently to lead to more timely conclusions compared to traditional trial designs. Adoption of new treatments in clinical practice can be accelerated if these trials are incorporated in electronic health records as part of a learning health system. In this review, we discuss challenges in the evaluation of treatments for sepsis, and explore potential benefits and weaknesses of recent advances in adaptive trial methodologies to address these challenges.

Secondary brain injury: Predicting and preventing insults

Mortality or severe disability affects the majority of patients after severe traumatic brain injury (TBI). Adherence to the brain trauma foundation guidelines has overall improved outcomes; however, traditional as well as novel interventions towards intracranial hypertension and secondary brain injury have come under scrutiny after series of negative randomized controlled trials. In fact, it would not be unfair to say there has been no single major breakthrough in the management of severe TBI in the last two decades. One plausible hypothesis for the aforementioned failures is that by the time treatment is initiated for neuroprotection, or physiologic optimization, irreversible brain injury has already set in. We, and others, have recently developed predictive models based on machine learning from continuous time series of intracranial pressure and partial brain tissue oxygenation. These models provide accurate predictions of physiologic crises events in a timely fashion, offering the opportunity for an earlier application of targeted interventions. In this article, we review the rationale for prediction, discuss available predictive models with examples, and offer suggestions for their future prospective testing in conjunction with preventive clinical algorithms. This article is part of the Special Issue entitled "Novel Treatments for Traumatic Brain Injury".

High-Performance Agent-Based Modeling Applied to Vocal Fold Inflammation and Repair

Fast and accurate computational biology models offer the prospect of accelerating the development of personalized medicine. A tool capable of estimating treatment success can help prevent unnecessary and costly treatments and potential harmful side effects. A novel high-performance Agent-Based Model (ABM) was adopted to simulate and visualize multi-scale complex biological processes arising in vocal fold inflammation and repair. The computational scheme was designed to organize the 3D ABM sub-tasks to fully utilize the resources available on current heterogeneous platforms consisting of multi-core CPUs and many-core GPUs. Subtasks are further parallelized and convolution-based diffusion is used to enhance the performance of the ABM simulation. The scheme was implemented using a client-server protocol allowing the results of each iteration to be analyzed and visualized on the server (i.e., in-situ) while the simulation is running on the same server. The resulting simulation and visualization software enables users to interact with and steer the course of the simulation in real-time as needed. This high-resolution 3D ABM framework was used for a case study of surgical vocal fold injury and repair. The new framework is capable of completing the simulation, visualization and remote result delivery in under 7 s per iteration, where each iteration of the simulation represents 30 min in the real world. The case study model was simulated at the physiological scale of a human vocal fold. This simulation tracks 17 million biological cells as well as a total of 1.7 billion signaling chemical and structural protein data points. The visualization component processes and renders all simulated biological cells and 154 million signaling chemical data points. The proposed high-performance 3D ABM was verified through comparisons with empirical vocal fold data. Representative trends of biomarker predictions in surgically injured vocal folds were observed.

Diagnosis and Management of Polytraumatized Patients With Severe Extremity Trauma

Multiply injured patients with severe extremity trauma are at risk of acute systemic complications and are at high risk of developing longer term orthopaedic complications including soft-tissue infection, osteomyelitis, posttraumatic osteoarthritis, and nonunion. It is becoming increasingly recognized that injury magnitude and response to injury have major jurisdiction pertaining to patient outcomes and complications. The complexities of injury and injury response that affect outcomes present opportunities to apply precision approaches to understand and quantify injury magnitude and injury response on a patient-specific basis. Here, we present novel approaches to measure injury magnitude by adopting methods that quantify both mechanical and ischemic tissue injury specific to each patient. We also present evolving computational approaches that have provided new insight into the complexities of inflammation and immunologic response to injury specific to each patient. These precision approaches are on the forefront of understanding how to stratify individualized injury and injury response in an effort to optimize titrated orthopaedic surgical interventions, which invariably involve most of the multiply injured patients. Finally, we present novel methods directed at mangled limbs with severe soft-tissue injury that comprise severely injured patients. Specifically, methods being developed to treat mangled limbs with volumetric muscle loss have the potential to improve limb outcomes and also mitigate uncompensated inflammation that occurs in these patients.

In silico clinical trials for pediatric orphan diseases

To date poor treatment options are available for patients with congenital pseudarthrosis of the tibia (CPT), a pediatric orphan disease. In this study we have performed an in silico clinical trial on 200 virtual subjects, generated from a previously established model of murine bone regeneration, to tackle the challenges associated with the small, pediatric patient population. Each virtual subject was simulated to receive no treatment and bone morphogenetic protein (BMP) treatment. We have shown that the degree of severity of CPT is significantly reduced with BMP treatment, although the effect is highly subject-specific. Using machine learning techniques we were also able to stratify the virtual subject population in adverse responders, non-responders, responders and asymptomatic. In summary, this study shows the potential of in silico medicine technologies as well as their implications for other orphan diseases.

The Virtual Anemia Trial: An Assessment of Model-Based In Silico Clinical Trials of Anemia Treatment Algorithms in Patients With Hemodialysis

In silico approaches have been proposed as a novel strategy to increase the repertoire of clinical trial designs. Realistic simulations of clinical trials can provide valuable information regarding safety and limitations of treatment protocols and have been shown to assist in the cost‐effective planning of clinical studies. In this report, we present a blueprint for the stepwise integration of internal, external, and ecological validity considerations in virtual clinical trials (VCTs). We exemplify this approach in the context of a model‐based in silico clinical trial aimed at anemia treatment in patients undergoing hemodialysis (HD). Hemoglobin levels and subsequent anemia treatment were simulated on a per patient level over the course of a year and compared to real‐life clinical data of 79,426 patients undergoing HD. The novel strategies presented here, aimed to improve external and ecological validity of a VCT, significantly increased the predictive power of the discussed in silico trial.

Mathematical models used to inform study design or surveillance systems in infectious diseases: a systematic review

BACKGROUND

Mathematical models offer the possibility to investigate the infectious disease dynamics over time and may help in informing design of studies. A systematic review was performed in order to determine to what extent mathematical models have been incorporated into the process of planning studies and hence inform study design for infectious diseases transmitted between humans and/or animals.

METHODS

We searched Ovid Medline and two trial registry platforms (Cochrane, WHO) using search terms related to infection, mathematical model, and study design from the earliest dates to October 2016. Eligible publications and registered trials included mathematical models (compartmental, individual-based, or Markov) which were described and used to inform the design of infectious disease studies. We extracted information about the investigated infection, population, model characteristics, and study design.

RESULTS

We identified 28 unique publications but no registered trials. Focusing on compartmental and individual-based models we found 12 observational/surveillance studies and 11 clinical trials. Infections studied were equally animal and human infectious diseases for the observational/surveillance studies, while all but one between humans for clinical trials. The mathematical models were used to inform, amongst other things, the required sample size (n = 16), the statistical power (n = 9), the frequency at which samples should be taken (n = 6), and from whom (n = 6).

CONCLUSIONS

Despite the fact that mathematical models have been advocated to be used at the planning stage of studies or surveillance systems, they are used scarcely. With only one exception, the publications described theoretical studies, hence, not being utilised in real studies.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Solving Immunology?

Emergent responses of the immune system result from the integration of molecular and cellular networks over time and across multiple organs. High-content and high-throughput analysis technologies, concomitantly with data-driven and mechanistic modeling, hold promise for the systematic interrogation of these complex pathways. However, connecting genetic variation and molecular mechanisms to individual phenotypes and health outcomes has proven elusive. Gaps remain in data, and disagreements persist about the value of mechanistic modeling for immunology. Here, we present the perspectives that emerged from the National Institute of Allergy and Infectious Disease (NIAID) workshop 'Complex Systems Science, Modeling and Immunity' and subsequent discussions regarding the potential synergy of high-throughput data acquisition, data-driven modeling, and mechanistic modeling to define new mechanisms of immunological disease and to accelerate the translation of these insights into therapies.

Phase i trials in melanoma: A framework to translate preclinical findings to the clinic

BACKGROUND

One of major issues in clinical trials in oncology is their high failure rate, despite the fact that the trials were designed based on the data from successful equivalent preclinical studies. This is in part due to the intrinsic homogeneity of preclinical model systems and the contrasting heterogeneity of actual patient responses.

METHODS

We present a mathematical model-driven framework, phase i (virtual/imaginary) trials, that integrates the heterogeneity of actual patient responses and preclinical studies through a cohort of virtual patients. The framework includes an experimentally calibrated mathematical model, a cohort of heterogeneous virtual patients, an assessment of stratification factors, and treatment optimisation. We show the detailed process through the lens of melanoma combination therapy (chemotherapy and an AKT inhibitor), using both preclinical and clinical data.

RESULTS

The mathematical model predicts melanoma treatment response and resistance to mono and combination therapies and was calibrated and then validated with in vitro experimental data. The validated model and a genetic algorithm were used to generate virtual patients whose tumour volume responses to the combination therapy matched statistically the actual heterogeneous patient responses in the clinical trial. Analyses on simulated cohorts revealed key model parameters such as a tumour volume doubling rate and a therapy-induced phenotypic switch rate that may have clinical correlates. Finally, our approach predicts optimal AKT inhibitor scheduling suggesting more effective but less toxic treatment strategies.

CONCLUSION

Our proposed computational framework to implement phase i trials in cancer can readily capture observed heterogeneous clinical outcomes and predict patient survival. Importantly, phase i trials can be used to optimise future clinical trial design.

Multiscale Modeling in the Clinic: Drug Design and Development

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multiscale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multiscale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions to guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multiscale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical and computational techniques employed for multiscale modeling approaches used in pharmacometric and systems pharmacology models in drug development and present several examples illustrating the current state-of-the-art models for (1) excitable systems and applications in cardiac disease; (2) stem cell driven complex biosystems; (3) nanoparticle delivery, with applications to angiogenesis and cancer therapy; (4) host-pathogen interactions and their use in metabolic disorders, inflammation and sepsis; and (5) computer-aided design of nanomedical systems. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multiscale models.

Systems Biology of Immunomodulation for Post-Stroke Neuroplasticity: Multimodal Implications of Pharmacotherapy and Neurorehabilitation

AIMS

Recent studies indicate that anti-inflammatory drugs, act as a double-edged sword, not only exacerbating secondary brain injury but also contributing to neurological recovery after stroke. Our aim is to explore whether there is a beneficial role for neuroprotection and functional recovery using anti-inflammatory drug along with neurorehabilitation therapy using transcranial direct current stimulation (tDCS) and repetitive transcranial magnetic stimulation (rTMS), so as to improve functional recovery after ischemic stroke.

METHODS

We develop a computational systems biology approach from preclinical data, using ordinary differential equations, to study the behavior of both phenotypes of microglia, such as M1 type (pro-inflammatory) vis-à-vis M2 type (anti-inflammatory) under anti-inflammatory drug action (minocycline). We explore whether pharmacological treatment along with cerebral stimulation using tDCS and rTMS is beneficial or not. We utilize the systems pathway analysis of minocycline in nuclear factor kappa beta (NF-κB) signaling and neurorehabilitation therapy using tDCS and rTMS that act through brain-derived neurotrophic factor (BDNF) and tropomyosin-related kinase B (TrkB) signaling pathways.

RESULTS

We demarcate the role of neuroinflammation and immunomodulation in post-stroke recovery, under minocycline activated-microglia and neuroprotection together with improved neurogenesis, synaptogenesis, and functional recovery under the action of rTMS or tDCS. We elucidate the feasibility of utilizing rTMS/tDCS to increase neuroprotection across the reperfusion stage during minocycline administration. We delineate that the signaling pathways of minocycline by modulation of inflammatory genes in NF-κB and proteins activated by tDCS and rTMS through BDNF, TrkB, and calmodulin kinase (CaMK) signaling. Utilizing systems biology approach, we show that the activation pathways for pharmacotherapy (minocycline) and neurorehabilitation (rTMS applied to ipsilesional cortex and tDCS) results into increased neuronal and synaptic activity that commonly occur through activation of N-methyl-d-aspartate receptors. We construe that considerable additive neuroprotection effect would be obtained and delayed reperfusion injury can be remedied, if one uses multimodal intervention of minocycline together with tDCS and rTMS.

CONCLUSION

Additive beneficial effect is, thus, noticed for pharmacotherapy along with neurorehabilitation therapy, by maneuvering the dynamics of immunomodulation using anti-inflammatory drug and cerebral stimulation for augmenting the functional recovery after stroke, which may engender clinical applicability for enhancing plasticity, rehabilitation, and neurorestoration.

Deciphering the Counterplay of Aspergillus fumigatus Infection and Host Inflammation by Evolutionary Games on Graphs

Microbial invaders are ubiquitously present and pose the constant risk of infections that are opposed by various defence mechanisms of the human immune system. A tight regulation of the immune response ensures clearance of microbial invaders and concomitantly limits host damage that is crucial for host viability. To investigate the counterplay of infection and inflammation, we simulated the invasion of the human-pathogenic fungus Aspergillus fumigatus in lung alveoli by evolutionary games on graphs. The layered structure of the innate immune system is represented by a sequence of games in the virtual model. We show that the inflammatory cascade of the immune response is essential for microbial clearance and that the inflammation level correlates with the infection-dose. At low infection-doses, corresponding to daily inhalation of conidia, the resident alveolar macrophages may be sufficient to clear infections, however, at higher infection-doses their primary task shifts towards recruitment of neutrophils to infection sites.

Implications of Heterogeneity of Treatment Effect for Reporting and Analysis of Randomized Trials in Critical Care

Randomized clinical trials (RCTs) are conducted to guide clinicians' selection of therapies for individual patients. Currently, RCTs in critical care often report an overall mean effect and selected individual subgroups. Yet work in other fields suggests that such reporting practices can be improved. Specifically, this Critical Care Perspective reviews recent work on so-called "heterogeneity of treatment effect" (HTE) by baseline risk and extends that work to examine its applicability to trials of acute respiratory failure and severe sepsis. Because patients in RCTs in critical care medicine-and patients in intensive care units-have wide variability in their risk of death, these patients will have wide variability in the absolute benefit that they can derive from a given therapy. If the side effects of the therapy are not perfectly collinear with the treatment benefits, this will result in HTE, where different patients experience quite different expected benefits of a therapy. We use simulations of RCTs to demonstrate that such HTE could result in apparent paradoxes, including: (1) positive trials of therapies that are beneficial overall but consistently harm or have little benefit to low-risk patients who met enrollment criteria, and (2) overall negative trials of therapies that still consistently benefit high-risk patients. We further show that these results persist even in the presence of causes of death unmodified by the treatment under study. These results have implications for reporting and analyzing RCT data, both to better understand how our therapies work and to improve the bedside applicability of RCTs. We suggest a plan for measurement in future RCTs in the critically ill.

Trauma in silico: Individual-specific mathematical models and virtual clinical populations

Trauma-induced critical illness is driven by acute inflammation, and elevated systemic interleukin-6 (IL-6) after trauma is a biomarker of adverse outcomes. We constructed a multicompartment, ordinary differential equation model that represents a virtual trauma patient. Individual-specific variants of this model reproduced both systemic inflammation and outcomes of 33 blunt trauma survivors, from which a cohort of 10,000 virtual trauma patients was generated. Model-predicted length of stay in the intensive care unit, degree of multiple organ dysfunction, and IL-6 area under the curve as a function of injury severity were in concordance with the results from a validation cohort of 147 blunt trauma patients. In a subcohort of 98 trauma patients, those with high-IL-6 single-nucleotide polymorphisms (SNPs) exhibited higher plasma IL-6 levels than those with low IL-6 SNPs, matching model predictions. Although IL-6 could drive mortality in individual virtual patients, simulated outcomes in the overall cohort were independent of the propensity to produce IL-6, a prediction verified in the 98-patient subcohort. In silico randomized clinical trials suggested a small survival benefit of IL-6 inhibition, little benefit of IL-1β inhibition, and worse survival after tumor necrosis factor-α inhibition. This study demonstrates the limitations of extrapolating from reductionist mechanisms to outcomes in individuals and populations and demonstrates the use of mechanistic simulation in complex diseases.

Integration of advances in social media and mHealth technology are pivotal to successful cancer prevention and control

The successful prevention and treatment of cancer is dependent upon efficient and reliable communication between healthcare workers and patients. Advances in social media and mHealth platforms have provided new ways in which to enhance the sharing of cancer related information. Other benefits of embracing this technology include utilising its analytic capabilities which can process the vast quantity of information generated from genome exploration in a highly efficient manner. The aim of this review is to provide an overview of the rapidly evolving areas through which digital engagement is proving useful in the prevention and control of cancer.

Reverse Engineering the Inflammatory "Clock": From Computational Modeling to Rational Resetting

Properly-regulated inflammation is central to homeostasis. Traumatic injury, hemorrhagic shock, septic shock, and other injury-related processes such as wound healing are associated with dysregulated inflammation. Like many biological processes, inflammation is a dynamic, complex system whose function, like that of an analog clock, cannot be discerned simply from a laundry list of its parts (data). The advent of multiplexed platforms for gathering biological data, while providing an unprecedented level of detailed information about the inflammatory response, has paradoxically also proven to be overwhelming. This problem is especially acute when the datasets involve time courses, since typical statistical analyses and data-driven modeling are geared towards single time points. Various groups have addressed this problem using dynamic approaches to data-driven and mechanistic computational modeling. These modeling tools can be thought of as the "gears" and "hands" of the "clock," and have led to insights regarding principal drivers, dynamic networks, feedbacks, and regulatory switches that characterize and perhaps regulate the inflammatory response. In parallel, mechanistic computational models have given an abstracted sense of how the inflammatory "clock" works, leading to in silico models of critically ill individuals and populations. Integrating data-driven and mechanistic modeling may point the way to a rational "resetting" of inflammation via model-driven precision medicine.

Insights into the Role of Chemokines, Damage-Associated Molecular Patterns, and Lymphocyte-Derived Mediators from Computational Models of Trauma-Induced Inflammation

SIGNIFICANCE

Traumatic injury elicits a complex, dynamic, multidimensional inflammatory response that is intertwined with complications such as multiple organ dysfunction and nosocomial infection. The complex interplay between inflammation and physiology in critical illness remains a challenge for translational research, including the extrapolation to human disease from animal models.

RECENT ADVANCES

Over the past decade, we and others have attempted to decipher the biocomplexity of inflammation in these settings of acute illness, using computational models to improve clinical translation. In silico modeling has been suggested as a computationally based framework for integrating data derived from basic biology experiments as well as preclinical and clinical studies.

CRITICAL ISSUES

Extensive studies in cells, mice, and human blunt trauma patients have led us to suggest (i) that while an adequate level of inflammation is required for healing post-trauma, inflammation can be harmful when it becomes self-sustaining via a damage-associated molecular pattern/Toll-like receptor-driven feed-forward circuit; (ii) that chemokines play a central regulatory role in driving either self-resolving or self-maintaining inflammation that drives the early activation of both classical innate and more recently recognized lymphoid pathways; and (iii) the presence of multiple thresholds and feedback loops, which could significantly affect the propagation of inflammation across multiple body compartments.

FUTURE DIRECTIONS

These insights from data-driven models into the primary drivers and interconnected networks of inflammation have been used to generate mechanistic computational models. Together, these models may be used to gain basic insights as well as serving to help define novel biomarkers and therapeutic targets.

Mathematical Modeling of Early Cellular Innate and Adaptive Immune Responses to Ischemia/Reperfusion Injury and Solid Organ Allotransplantation

A mathematical model of the early inflammatory response in transplantation is formulated with ordinary differential equations. We first consider the inflammatory events associated only with the initial surgical procedure and the subsequent ischemia/reperfusion (I/R) events that cause tissue damage to the host as well as the donor graft. These events release damage-associated molecular pattern molecules (DAMPs), thereby initiating an acute inflammatory response. In simulations of this model, resolution of inflammation depends on the severity of the tissue damage caused by these events and the patient's (co)-morbidities. We augment a portion of a previously published mathematical model of acute inflammation with the inflammatory effects of T cells in the absence of antigenic allograft mismatch (but with DAMP release proportional to the degree of graft damage prior to transplant). Finally, we include the antigenic mismatch of the graft, which leads to the stimulation of potent memory T cell responses, leading to further DAMP release from the graft and concomitant increase in allograft damage. Regulatory mechanisms are also included at the final stage. Our simulations suggest that surgical injury and I/R-induced graft damage can be well-tolerated by the recipient when each is present alone, but that their combination (along with antigenic mismatch) may lead to acute rejection, as seen clinically in a subset of patients. An emergent phenomenon from our simulations is that low-level DAMP release can tolerize the recipient to a mismatched allograft, whereas different restimulation regimens resulted in an exaggerated rejection response, in agreement with published studies. We suggest that mechanistic mathematical models might serve as an adjunct for patient- or sub-group-specific predictions, simulated clinical studies, and rational design of immunosuppression.

A Computational, Tissue-Realistic Model of Pressure Ulcer Formation in Individuals with Spinal Cord Injury

People with spinal cord injury (SCI) are predisposed to pressure ulcers (PU). PU remain a significant burden in cost of care and quality of life despite improved mechanistic understanding and advanced interventions. An agent-based model (ABM) of ischemia/reperfusion-induced inflammation and PU (the PUABM) was created, calibrated to serial images of post-SCI PU, and used to investigate potential treatments in silico. Tissue-level features of the PUABM recapitulated visual patterns of ulcer formation in individuals with SCI. These morphological features, along with simulated cell counts and mediator concentrations, suggested that the influence of inflammatory dynamics caused simulations to be committed to "better" vs. "worse" outcomes by 4 days of simulated time and prior to ulcer formation. Sensitivity analysis of model parameters suggested that increasing oxygen availability would reduce PU incidence. Using the PUABM, in silico trials of anti-inflammatory treatments such as corticosteroids and a neutralizing antibody targeted at Damage-Associated Molecular Pattern molecules (DAMPs) suggested that, at best, early application at a sufficiently high dose could attenuate local inflammation and reduce pressure-associated tissue damage, but could not reduce PU incidence. The PUABM thus shows promise as an adjunct for mechanistic understanding, diagnosis, and design of therapies in the setting of PU.

The inflammatory response to influenza A virus (H1N1): An experimental and mathematical study

Mortality from influenza infections continues as a global public health issue, with the host inflammatory response contributing to fatalities related to the primary infection. Based on Ordinary Differential Equation (ODE) formalism, a computational model was developed for the in-host response to influenza A virus, merging inflammatory, innate, adaptive and humoral responses to virus and linking severity of infection, the inflammatory response, and mortality. The model was calibrated using dense cytokine and cell data from adult BALB/c mice infected with the H1N1 influenza strain A/PR/8/34 in sublethal and lethal doses. Uncertainty in model parameters and disease mechanisms was quantified using Bayesian inference and ensemble model methodology that generates probabilistic predictions of survival, defined as viral clearance and recovery of the respiratory epithelium. The ensemble recovers the expected relationship between magnitude of viral exposure and the duration of survival, and suggests mechanisms primarily responsible for survival, which could guide the development of immuno-modulatory interventions as adjuncts to current anti-viral treatments. The model is employed to extrapolate from available data survival curves for the population and their dependence on initial viral aliquot. In addition, the model allows us to illustrate the positive effect of controlled inflammation on influenza survival.

Computational Studies of the Intestinal Host-Microbiota Interactome

A large and growing body of research implicates aberrant immune response and compositional shifts of the intestinal microbiota in the pathogenesis of many intestinal disorders. The molecular and physical interaction between the host and the microbiota, known as the host-microbiota interactome, is one of the key drivers in the pathophysiology of many of these disorders. This host-microbiota interactome is a set of dynamic and complex processes, and needs to be treated as a distinct entity and subject for study. Disentangling this complex web of interactions will require novel approaches, using a combination of data-driven bioinformatics with knowledge-driven computational modeling. This review describes the computational approaches for investigating the host-microbiota interactome, with emphasis on the human intestinal tract and innate immunity, and highlights open challenges and existing gaps in the computation methodology for advancing our knowledge about this important facet of human health.