Emulation of Quantitative Systems Pharmacology models to accelerate virtual population inference in immuno-oncology.

Quantitative Systems Pharmacology (QSP) models are increasingly being applied for target discovery and dose selection in immuno-oncology (IO). Typical application involves virtual trial, a simulation of a virtual population of hundreds of model instances with model inputs reflecting individual variability. While the structure of the model and initial parameterisation are based on literature describing the underlying biology, calibration of the virtual population by existing clinical data is frequently required to create tumour and patient population specific model instances. Since comparison of a virtual trial with clinical output requires hundreds of large-scale, non-linear model evaluations, the inference of a virtual population is computationally expensive, frequently becoming a bottleneck. Here, we present novel approach to virtual population inference in IO using emulation of the QSP model and an objective function based on Kolmogorov-Smirnov statistics to maximise congruence of simulated and observed clinical tumour size distributions. We sample the parameter space of a QSP IO model to collect a set of tumour growth time profiles. We evaluate performance of several machine learning approaches in interpolating these time profiles and create a surrogate model, which computes tumor growth profiles faster than the original model and allows examination of tens of millions of virtual patients. We use the surrogate model to infer a virtual population maximising congruence with the waterfall plot of a pembrolizumab clinical trial. We believe that our approach is applicable not only in QSP IO, but also in other applications where virtual populations need to be inferred for computationally expensive mechanistic models.

Mitochondrial Metabolism in the Spotlight: Maintaining Balanced RNAP III Activity Ensures Cellular Homeostasis.

RNA polymerase III (RNAP III) holoenzyme activity and the processing of its products have been linked to several metabolic dysfunctions in lower and higher eukaryotes. Alterations in the activity of RNAP III-driven synthesis of non-coding RNA cause extensive changes in glucose metabolism. Increased RNAP III activity in the S. cerevisiae maf1Δ strain is lethal when grown on a non-fermentable carbon source. This lethal phenotype is suppressed by reducing tRNA synthesis. Neither the cause of the lack of growth nor the underlying molecular mechanism have been deciphered, and this area has been awaiting scientific explanation for a decade. Our previous proteomics data suggested mitochondrial dysfunction in the strain. Using model mutant strains maf1Δ (with increased tRNA abundance) and rpc128-1007 (with reduced tRNA abundance), we collected data showing major changes in the TCA cycle metabolism of the mutants that explain the phenotypic observations. Based on 13C flux data and analysis of TCA enzyme activities, the present study identifies the flux constraints in the mitochondrial metabolic network. The lack of growth is associated with a decrease in TCA cycle activity and downregulation of the flux towards glutamate, aspartate and phosphoenolpyruvate (PEP), the metabolic intermediate feeding the gluconeogenic pathway. rpc128-1007, the strain that is unable to increase tRNA synthesis due to a mutation in the C128 subunit, has increased TCA cycle activity under non-fermentable conditions. To summarize, cells with non-optimal activity of RNAP III undergo substantial adaptation to a new metabolic state, which makes them vulnerable under specific growth conditions. Our results strongly suggest that balanced, non-coding RNA synthesis that is coupled to glucose signaling is a fundamental requirement to sustain a cell's intracellular homeostasis and flexibility under changing growth conditions. The presented results provide insight into the possible role of RNAP III in the mitochondrial metabolism of other cell types.

QSP Designer: Quantitative systems pharmacology modeling with modular biological process map notation and multiple language code generation.

Typical Quantitative Systems Pharmacology (QSP) workflows involve discussion of biology, supported by graphical diagrams, followed by construction of large Ordinary Differential Equation models. QSP Designer facilitates this process by providing enhanced graphical notation, which enables hierarchical presentation with modules and handling of combinatorial complexity with diagram node arrays. Whereas the software includes a simulation engine, a major feature is full model code generation in MATLAB, R, C, and Julia to support multiple modeling communities.

Learn-confirm in model-informed drug development: Assessing an immunogenicity quantitative systems pharmacology platform.

Immunogenicity against therapeutic proteins frequently causes attrition owing to its potential impact on pharmacokinetics, pharmacodynamics, efficacy, and safety. Predicting immunogenicity is complex because of its multifactorial drivers, including compound properties, subject characteristics, and treatment parameters. To integrate these, the Immunogenicity Simulator was developed using published, predominantly late-stage trial data from 15 therapeutic proteins. This single-blinded evaluation with subject-level data from 10 further monoclonals assesses the Immunogenicity Simulator's credibility for application during the drug development process.

Mechanistic Models of COVID-19: Insights into Disease Progression, Vaccines, and Therapeutics.

The COVID-19 pandemic has severely impacted health systems and economies worldwide. Significant global efforts are therefore ongoing to improve vaccine efficacies, optimize vaccine deployment, and develop new antiviral therapies to combat the pandemic. Mechanistic viral dynamics and quantitative systems pharmacology models of SARS-CoV-2 infection, vaccines, immunomodulatory agents, and antiviral therapeutics have played a key role in advancing our understanding of SARS-CoV-2 pathogenesis and transmission, the interplay between innate and adaptive immunity to influence the outcomes of infection, effectiveness of treatments, mechanisms and performance of COVID-19 vaccines, and the impact of emerging SARS-CoV-2 variants. Here, we review some of the critical insights provided by these models and discuss the challenges ahead.

Application of quantitative systems pharmacology to guide the optimal dosing of COVID-19 vaccines.

Optimal use and distribution of coronavirus disease 2019 (COVID-19) vaccines involves adjustments of dosing. Due to the rapidly evolving pandemic, such adjustments often need to be introduced before full efficacy data are available. As demonstrated in other areas of drug development, quantitative systems pharmacology (QSP) is well placed to guide such extrapolation in a rational and timely manner. Here, we propose for the first time how QSP can be applied in the context of COVID-19 vaccine development.

Quantitative Systems Pharmacology Approaches for Immuno-Oncology: Adding Virtual Patients to the Development Paradigm.

Drug development in oncology commonly exploits the tools of molecular biology to gain therapeutic benefit through reprograming of cellular responses. In immuno-oncology (IO) the aim is to direct the patient's own immune system to fight cancer. After remarkable successes of antibodies targeting PD1/PD-L1 and CTLA4 receptors in targeted patient populations, the focus of further development has shifted toward combination therapies. However, the current drug-development approach of exploiting a vast number of possible combination targets and dosing regimens has proven to be challenging and is arguably inefficient. In particular, the unprecedented number of clinical trials testing different combinations may no longer be sustainable by the population of available patients. Further development in IO requires a step change in selection and validation of candidate therapies to decrease development attrition rate and limit the number of clinical trials. Quantitative systems pharmacology (QSP) proposes to tackle this challenge through mechanistic modeling and simulation. Compounds' pharmacokinetics, target binding, and mechanisms of action as well as existing knowledge on the underlying tumor and immune system biology are described by quantitative, dynamic models aiming to predict clinical results for novel combinations. Here, we review the current QSP approaches, the legacy of mathematical models available to quantitative clinical pharmacologists describing interaction between tumor and immune system, and the recent development of IO QSP platform models. We argue that QSP and virtual patients can be integrated as a new tool in existing IO drug development approaches to increase the efficiency and effectiveness of the search for novel combination therapies.

Integration of Omics Data Sources to Inform Mechanistic Modeling of Immune-Oncology Therapies: A Tutorial for Clinical Pharmacologists.

Application of contemporary molecular biology techniques to clinical samples in oncology resulted in the accumulation of unprecedented experimental data. These "omics" data are mined for discovery of therapeutic target combinations and diagnostic biomarkers. It is less appreciated that omics resources could also revolutionize development of the mechanistic models informing clinical pharmacology quantitative decisions about dose amount, timing, and sequence. We discuss the integration of omics data to inform mechanistic models supporting drug development in immuno-oncology. To illustrate our arguments, we present a minimal clinical model of the Cancer Immunity Cycle (CIC), calibrated for non-small cell lung carcinoma using tumor microenvironment composition inferred from transcriptomics of clinical samples. We review omics data resources, which can be integrated to parameterize mechanistic models of the CIC. We propose that virtual trial simulations with clinical Quantitative Systems Pharmacology platforms informed by omics data will be making increasing impact in the development of cancer immunotherapies.

Genome-wide analysis of the role of the antibiotic biosynthesis regulator AbsA2 in Streptomyces coelicolor A3(2).

The AbsA1-AbsA2 two component signalling system of Streptomyces coelicolor has long been known to exert a powerful negative influence on the production of the antibiotics actinorhodin, undecylprodiginine and the Calcium-Dependent Antibiotic (CDA). Here we report the analysis of a ΔabsA2 deletion strain, which exhibits the classic precocious antibiotic hyper-production phenotype, and its complementation by an N-terminal triple-FLAG-tagged version of AbsA2. The complemented and non-complemented ΔabsA2 mutant strains were used in large-scale microarray-based time-course experiments to investigate the effect of deleting absA2 on gene expression and to identify the in vivo AbsA2 DNA-binding target sites using ChIP-on chip. We show that in addition to binding to the promoter regions of redZ and actII-orfIV AbsA2 binds to several previously unidentified sites within the cda biosynthetic gene cluster within and/or upstream of SCO3215-SCO3216, SCO3217, SCO3229-SCO3230, and SCO3226, and we relate the pattern of AbsA2 binding to the results of the transcriptomic study and antibiotic phenotypic assays. Interestingly, dual 'biphasic' ChIP peaks were observed with AbsA2 binding across the regulatory genes actII-orfIV and redZ and the absA2 gene itself, while more conventional single promoter-proximal peaks were seen at the CDA biosynthetic genes suggesting a different mechanism of regulation of the former loci. Taken together the results shed light on the complex mechanism of regulation of antibiotic biosynthesis in Streptomyces coelicolor and the important role of AbsA2 in controlling the expression of three antibiotic biosynthetic gene clusters.

Multi-scale, whole-system models of liver metabolic adaptation to fat and sugar in non-alcoholic fatty liver disease.

Non-alcoholic fatty liver disease (NAFLD) is a serious public health issue associated with high fat, high sugar diets. However, the molecular mechanisms mediating NAFLD pathogenesis are only partially understood. Here we adopt an iterative multi-scale, systems biology approach coupled to in vitro experimentation to investigate the roles of sugar and fat metabolism in NAFLD pathogenesis. The use of fructose as a sweetening agent is controversial; to explore this, we developed a predictive model of human monosaccharide transport, signalling and metabolism. The resulting quantitative model comprising a kinetic model describing monosaccharide transport and insulin signalling integrated with a hepatocyte-specific genome-scale metabolic network (GSMN). Differential kinetics for the utilisation of glucose and fructose were predicted, but the resultant triacylglycerol production was predicted to be similar for monosaccharides; these predictions were verified by in vitro data. The role of physiological adaptation to lipid overload was explored through the comprehensive reconstruction of the peroxisome proliferator activated receptor alpha (PPARα) regulome integrated with a hepatocyte-specific GSMN. The resulting qualitative model reproduced metabolic responses to increased fatty acid levels and mimicked lipid loading in vitro. The model predicted that activation of PPARα by lipids produces a biphasic response, which initially exacerbates steatosis. Our data support the evidence that it is the quantity of sugar rather than the type that is critical in driving the steatotic response. Furthermore, we predict PPARα-mediated adaptations to hepatic lipid overload, shedding light on potential challenges for the use of PPARα agonists to treat NAFLD.

Integration of Genome Scale Metabolic Networks and Gene Regulation of Metabolic Enzymes With Physiologically Based Pharmacokinetics.

The scope of physiologically based pharmacokinetic (PBPK) modeling can be expanded by assimilation of the mechanistic models of intracellular processes from systems biology field. The genome scale metabolic networks (GSMNs) represent a whole set of metabolic enzymes expressed in human tissues. Dynamic models of the gene regulation of key drug metabolism enzymes are available. Here, we introduce GSMNs and review ongoing work on integration of PBPK, GSMNs, and metabolic gene regulation. We demonstrate example models.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types.

UNLABELLED: Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis. SIGNIFICANCE: We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

The prognostic significance of specific HOX gene expression patterns in ovarian cancer.

HOX genes are vital for all aspects of mammalian growth and differentiation, and their dysregulated expression is related to ovarian carcinogenesis. The aim of the current study was to establish the prognostic value of HOX dysregulation as well as its role in platinum resistance. The potential to target HOX proteins through the HOX/PBX interaction was also explored in the context of platinum resistance. HOX gene expression was determined in ovarian cancer cell lines and primary EOCs by QPCR, and compared to expression in normal ovarian epithelium and fallopian tube tissue samples. Statistical analysis included one-way ANOVA and t-tests, using statistical software R and GraphPad. The analysis identified 36 of the 39 HOX genes as being overexpressed in high grade serous EOC compared to normal tissue. We detected a molecular HOX gene-signature that predicted poor outcome. Overexpression of HOXB4 and HOXB9 was identified in high grade serous cell lines after platinum resistance developed. Targeting the HOX/PBX dimer with the HXR9 peptide enhanced the cytotoxicity of cisplatin in platinum-resistant ovarian cancer. In conclusion, this study has shown the HOX genes are highly dysregulated in ovarian cancer with high expression of HOXA13, B6, C13, D1 and D13 being predictive of poor clinical outcome. Targeting the HOX/PBX dimer in platinum-resistant cancer represents a potentially new therapeutic option that should be further developed and tested in clinical trials.

Generation of 2,000 breast cancer metabolic landscapes reveals a poor prognosis group with active serotonin production.

A major roadblock in the effective treatment of cancers is their heterogeneity, whereby multiple molecular landscapes are classified as a single disease. To explore the contribution of cellular metabolism to cancer heterogeneity, we analyse the Metabric dataset, a landmark genomic and transcriptomic study of 2,000 individual breast tumours, in the context of the human genome-scale metabolic network. We create personalized metabolic landscapes for each tumour by exploring sets of active reactions that satisfy constraints derived from human biochemistry and maximize congruency with the Metabric transcriptome data. Classification of the personalized landscapes derived from 997 tumour samples within the Metabric discovery dataset reveals a novel poor prognosis cluster, reproducible in the 995-sample validation dataset. We experimentally follow mechanistic hypotheses resulting from the computational study and establish that active serotonin production is a major metabolic feature of the poor prognosis group. These data support the reconsideration of concomitant serotonin-specific uptake inhibitors treatment during breast cancer chemotherapy.

MUFINS: multi-formalism interaction network simulator.

Systems Biology has established numerous approaches for mechanistic modeling of molecular networks in the cell and a legacy of models. The current frontier is the integration of models expressed in different formalisms to address the multi-scale biological system organization challenge. We present MUFINS (MUlti-Formalism Interaction Network Simulator) software, implementing a unique set of approaches for multi-formalism simulation of interaction networks. We extend the constraint-based modeling (CBM) framework by incorporation of linear inhibition constraints, enabling for the first time linear modeling of networks simultaneously describing gene regulation, signaling and whole-cell metabolism at steady state. We present a use case where a logical hypergraph model of a regulatory network is expressed by linear constraints and integrated with a Genome-Scale Metabolic Network (GSMN) of mouse macrophage. We experimentally validate predictions, demonstrating application of our software in an iterative cycle of hypothesis generation, validation and model refinement. MUFINS incorporates an extended version of our Quasi-Steady State Petri Net approach to integrate dynamic models with CBM, which we demonstrate through a dynamic model of cortisol signaling integrated with the human Recon2 GSMN and a model of nutrient dynamics in physiological compartments. Finally, we implement a number of methods for deriving metabolic states from ~omics data, including our new variant of the iMAT congruency approach. We compare our approach with iMAT through the analysis of 262 individual tumor transcriptomes, recovering features of metabolic reprogramming in cancer. The software provides graphics user interface with network visualization, which facilitates use by researchers who are not experienced in coding and mathematical modeling environments.

Integrating Kinetic Model of E. coli with Genome Scale Metabolic Fluxes Overcomes Its Open System Problem and Reveals Bistability in Central Metabolism.

An understanding of the dynamics of the metabolic profile of a bacterial cell is sought from a dynamical systems analysis of kinetic models. This modelling formalism relies on a deterministic mathematical description of enzyme kinetics and their metabolite regulation. However, it is severely impeded by the lack of available kinetic information, limiting the size of the system that can be modelled. Furthermore, the subsystem of the metabolic network whose dynamics can be modelled is faced with three problems: how to parameterize the model with mostly incomplete steady state data, how to close what is now an inherently open system, and how to account for the impact on growth. In this study we address these challenges of kinetic modelling by capitalizing on multi-'omics' steady state data and a genome-scale metabolic network model. We use these to generate parameters that integrate knowledge embedded in the genome-scale metabolic network model, into the most comprehensive kinetic model of the central carbon metabolism of E. coli realized to date. As an application, we performed a dynamical systems analysis of the resulting enriched model. This revealed bistability of the central carbon metabolism and thus its potential to express two distinct metabolic states. Furthermore, since our model-informing technique ensures both stable states are constrained by the same thermodynamically feasible steady state growth rate, the ensuing bistability represents a temporal coexistence of the two states, and by extension, reveals the emergence of a phenotypically heterogeneous population.

Lipid metabolism and Type VII secretion systems dominate the genome scale virulence profile of Mycobacterium tuberculosis in human dendritic cells.

BACKGROUND: Mycobacterium tuberculosis continues to kill more people than any other bacterium. Although its archetypal host cell is the macrophage, it also enters, and survives within, dendritic cells (DCs). By modulating the behaviour of the DC, M. tuberculosis is able to manipulate the host's immune response and establish an infection. To identify the M. tuberculosis genes required for survival within DCs we infected primary human DCs with an M. tuberculosis transposon library and identified mutations with a reduced ability to survive. RESULTS: Parallel sequencing of the transposon inserts of the surviving mutants identified a large number of genes as being required for optimal intracellular fitness in DCs. Loci whose mutation attenuated intracellular survival included those involved in synthesising cell wall lipids, not only the well-established virulence factors, pDIM and cord factor, but also sulfolipids and PGL, which have not previously been identified as having a direct virulence role in cells. Other attenuated loci included the secretion systems ESX-1, ESX-2 and ESX-4, alongside many PPE genes, implicating a role for ESX-5. In contrast the canonical ESAT-6 family of ESX substrates did not have intra-DC fitness costs suggesting an alternative ESX-1 associated virulence mechanism. With the aid of a gene-nutrient interaction model, metabolic processes such as cholesterol side chain catabolism, nitrate reductase and cysteine-methionine metabolism were also identified as important for survival in DCs. CONCLUSION: We conclude that many of the virulence factors required for survival in DC are shared with macrophages, but that survival in DCs also requires several additional functions, such as cysteine-methionine metabolism, PGLs, sulfolipids, ESX systems and PPE genes.

[Significance of biomechanical criteria in proper treatment planning in knee osteoarthritis].

Introduction: Following the established and unequivocal criteria when choosing the treatment of knee osteoarthritis, despite possibility of precise imaging, is still very problematic. This is partly because doctors of different disciplines: general practitioners, orthopaedic surgeons, rheumatologists, physiotherapists, are involved in the treatment of this disease. For most of them the basic criteria to implement the treatment are: pain and assessment of X -ray in the supine position. As a result of that, despite slight and doubtful improvement, treatment is improperly targeted and extended. The aim of the study was to prove that the correct diagnosis, including biomechanical characteristics of the affected limb, choice of correct treatment, and observance of indications, can shorten treatment, make it more efficient and less expensive. Material and methods: We analyzed of 103 patients qualified for total knee arthroplasty between 2006 and 2011. The indication was primary and advanced knee osteoarthritis in phases III and IV according to Ahlbäcka scale medial displacement of mechanical axis deviation in the lower limb with club- -foot (I to IV degrees) and with centre of rotation of angulation (CORA) in the proximal part of the tibia. Only subjects with body mass index under 25 were included in the study. Surveys were used to assess the pre -operative duration of illness, non-surgical methods of treatment and their effects. Final evaluation was made using the Visual Analogue Scale (10 -degrees) of pain (during the effort and rest) and gait efficiency 30 days after completed treatment. Results: The 66.5% of patients showed no significant improvement in all assessed parameters, which would justify to continuation of the therapy as described above.