Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology.

Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.

Decoding the Foodome: Molecular Networks Connecting Diet and Health.

Diet, a modifiable risk factor, plays a pivotal role in most diseases, from cardiovascular disease to type 2 diabetes mellitus, cancer, and obesity. However, our understanding of the mechanistic role of the chemical compounds found in food remains incomplete. In this review, we explore the "dark matter" of nutrition, going beyond the macro- and micronutrients documented by national databases to unveil the exceptional chemical diversity of food composition. We also discuss the need to explore the impact of each compound in the presence of associated chemicals and relevant food sources and describe the tools that will allow us to do so. Finally, we discuss the role of network medicine in understanding the mechanism of action of each food molecule. Overall, we illustrate the important role of network science and artificial intelligence in our ability to reveal nutrition's multifaceted role in health and disease.

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.

Model-Informed Precision Dosing: Background, Requirements, Validation, Implementation, and Forward Trajectory of Individualizing Drug Therapy.

Model-informed precision dosing (MIPD) has become synonymous with modern approaches for individualizing drug therapy, in which the characteristics of each patient are considered as opposed to applying a one-size-fits-all alternative. This review provides a brief account of the current knowledge, practices, and opinions on MIPD while defining an achievable vision for MIPD in clinical care based on available evidence. We begin with a historical perspective on variability in dose requirements and then discuss technical aspects of MIPD, including the need for clinical decision support tools, practical validation, and implementation of MIPD in health care. We also discuss novel ways to characterize patient variability beyond the common perceptions of genetic control. Finally, we address current debates on MIPD from the perspectives of the new drug development, health economics, and drug regulations.

Efficacy difference of antipsychotics in Alzheimer's disease and schizophrenia: explained with network efficiency and pathway analysis methods.

Approximately 50% of Alzheimer's disease (AD) patients will develop psychotic symptoms and these patients will experience severe rapid cognitive decline compared with those without psychosis (AD-P). Currently, no medication has been approved by the Food and Drug Administration for AD with psychosis (AD+P) specifically, although atypical antipsychotics are widely used in clinical practice. These drugs have demonstrated modest efficacy in managing psychosis in individuals with AD, with an increased frequency of adverse events, including excess mortality. We compared the differences between the genetic variations/genes associated with AD+P and schizophrenia from existing Genome-Wide Association Study and differentially expressed genes (DEGs). We also constructed disease-specific protein-protein interaction networks for AD+P and schizophrenia. Network efficiency was then calculated to characterize the topological structures of these two networks. The efficiency of antipsychotics in these two networks was calculated. A weight adjustment based on binding affinity to drug targets was later applied to refine our results, and 2013 and 2123 genes were identified as related to AD+P and schizophrenia, respectively, with only 115 genes shared. Antipsychotics showed a significantly lower efficiency in the AD+P network than in the schizophrenia network (P < 0.001) indicating that antipsychotics may have less impact in AD+P than in schizophrenia. AD+P may be caused by mechanisms distinct from those in schizophrenia which result in a decreased efficacy of antipsychotics in AD+P. In addition, the network analysis methods provided quantitative explanations of the lower efficacy of antipsychotics in AD+P.

High dimensionality reduction by matrix factorization for systems pharmacology.

The extraction of predictive features from the complex high-dimensional multi-omic data is necessary for decoding and overcoming the therapeutic responses in systems pharmacology. Developing computational methods to reduce high-dimensional space of features in in vitro, in vivo and clinical data is essential to discover the evolution and mechanisms of the drug responses and drug resistance. In this paper, we have utilized the matrix factorization (MF) as a modality for high dimensionality reduction in systems pharmacology. In this respect, we have proposed three novel feature selection methods using the mathematical conception of a basis for features. We have applied these techniques as well as three other MF methods to analyze eight different gene expression datasets to investigate and compare their performance for feature selection. Our results show that these methods are capable of reducing the feature spaces and find predictive features in terms of phenotype determination. The three proposed techniques outperform the other methods used and can extract a 2-gene signature predictive of a tyrosine kinase inhibitor treatment response in the Cancer Cell Line Encyclopedia.

Quantitative Proteomics for Translational Pharmacology and Precision Medicine: State of The Art and Future Outlook.

Over the past 20 years, quantitative proteomics has contributed a wealth of protein expression data, which are currently used for a variety of systems pharmacology applications, as a complement or a surrogate for activity of the corresponding proteins. A symposium at the 25th North American ISSX meeting, in Boston, in September 2023, was held to explore current and emerging applications of quantitative proteomics in translational pharmacology and strategies for improved integration into model-informed drug development based on practical experience of each of the presenters. A summary of the talks and discussions is presented in this perspective alongside future outlooks that were outlined for future meetings. Significance Statement This perspective explores current and emerging applications of quantitative proteomics in translational pharmacology and precision medicine, and outlines outlooks for improved integration into model-informed drug development.

An In Silico Platform to Predict Cardiotoxicity Risk of Anti-tumor Drug Combination with hiPSC-CMs Based In Vitro Study.

OBJECTIVE: Antineoplastic agent-induced systolic dysfunction is a major reason for interruption of anticancer treatment. Although targeted anticancer agents infrequently cause systolic dysfunction, their combinations with chemotherapies remarkably increase the incidence. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) provide a potent in vitro model to assess cardiovascular safety. However, quantitatively predicting the reduction of ejection fraction based on hiPSC-CMs is challenging due to the absence of the body's regulatory response to cardiomyocyte injury. METHODS: Here, we developed and validated an in vitro-in vivo translational platform to assess the reduction of ejection fraction induced by antineoplastic drugs based on hiPSC-CMs. The translational platform integrates drug exposure, drug-cardiomyocyte interaction, and systemic response. The drug-cardiomyocyte interaction was implemented as a mechanism-based toxicodynamic (TD) model, which was then integrated into a quantitative system pharmacology-physiological-based pharmacokinetics (QSP-PBPK) model to form a complete translational platform. The platform was validated by comparing the model-predicted and clinically observed incidence of doxorubicin and trastuzumab-induced systolic dysfunction. RESULTS: A total of 33,418 virtual patients were incorporated to receive doxorubicin and trastuzumab alone or in combination. For doxorubicin, the QSP-PBPK-TD model successfully captured the overall trend of systolic dysfunction incidences against the cumulative doses. For trastuzumab, the predicted incidence interval was 0.31-2.7% for single-agent treatment and 0.15-10% for trastuzumab-doxorubicin sequential treatment, covering the observations in clinical reports (0.50-1.0% and 1.5-8.3%, respectively). CONCLUSIONS: In conclusion, the in vitro-in vivo translational platform is capable of predicting systolic dysfunction incidence almost merely depend on hiPSC-CMs, which could facilitate optimizing the treatment protocol of antineoplastic agents.

Repurposing ibudilast to mitigate Alzheimer's disease by targeting inflammation.

Alzheimer's disease is a multifactorial disease that exhibits cognitive deficits, neuronal loss, amyloid plaques, neurofibrillary tangles and neuroinflammation in the brain. Hence, a multi-target drug would improve treatment efficacy. We applied a new multi-scale predictive modelling framework that integrates machine learning with biophysics and systems pharmacology to screen drugs for Alzheimer's disease using patients' tissue samples. Our predictive modelling framework identified ibudilast as a drug with repurposing potential to treat Alzheimer's disease. Ibudilast is a multi-target drug, as it is a phosphodiesterase inhibitor and toll-like receptor 4 (TLR4) antagonist. In addition, we predict that ibudilast inhibits off-target kinases (e.g. IRAK1 and GSG2). In Japan and other Asian countries, ibudilast is approved for treating asthma and stroke due to its anti-inflammatory potential. Based on these previous studies and on our predictions, we tested for the first time the efficacy of ibudilast in Fisher transgenic 344-AD rats. This transgenic rat model is unique as it exhibits hippocampal-dependent spatial learning and memory deficits and Alzheimer's disease pathology, including hippocampal amyloid plaques, tau paired-helical filaments, neuronal loss and microgliosis, in a progressive age-dependent manner that mimics the pathology observed in Alzheimer's disease patients. Following long-term treatment with ibudilast, transgenic rats were evaluated at 11 months of age for spatial memory performance and Alzheimer's disease pathology. We demonstrate that ibudilast-treatment of transgenic rats mitigated hippocampal-dependent spatial memory deficits, as well as hippocampal (hilar subregion) amyloid plaque and tau paired-helical filament load, and microgliosis compared to untreated transgenic rat. Neuronal density analysed across all hippocampal regions was similar in ibudilast-treated transgenic compared to untreated transgenic rats. Interestingly, RNA sequencing analysis of hippocampal tissue showed that ibudilast-treatment affects gene expression levels of the TLR and ubiquitin-proteasome pathways differentially in male and female transgenic rats. Based on the TLR4 signalling pathway, our RNA sequencing data suggest that ibudilast-treatment inhibits IRAK1 activity by increasing expression of its negative regulator IRAK3, and/or by altering TRAF6 and other TLR-related ubiquitin ligase and conjugase levels. Our results support that ibudilast can serve as a repurposed drug that targets multiple pathways including TLR signalling and the ubiquitin/proteasome pathway to reduce cognitive deficits and pathology relevant to Alzheimer's disease.

Quantitative systems pharmacology modeling of HER2-positive metastatic breast cancer for translational efficacy evaluation and combination assessment across therapeutic modalities.

HER2-positive (HER2+) metastatic breast cancer (mBC) is highly aggressive and a major threat to human health. Despite the significant improvement in patients' prognosis given the drug development efforts during the past several decades, many clinical questions still remain to be addressed such as efficacy when combining different therapeutic modalities, best treatment sequences, interindividual variability as well as resistance and potential coping strategies. To better answer these questions, we developed a mechanistic quantitative systems pharmacology model of the pathophysiology of HER2+ mBC that was extensively calibrated and validated against multiscale data to quantitatively predict and characterize the signal transduction and preclinical tumor growth kinetics under different therapeutic interventions. Focusing on the second-line treatment for HER2+ mBC, e.g., antibody-drug conjugates (ADC), small molecule inhibitors/TKI and chemotherapy, the model accurately predicted the efficacy of various drug combinations and dosing regimens at the in vitro and in vivo levels. Sensitivity analyses and subsequent heterogeneous phenotype simulations revealed important insights into the design of new drug combinations to effectively overcome various resistance scenarios in HER2+ mBC treatments. In addition, the model predicted a better efficacy of the new TKI plus ADC combination which can potentially reduce drug dosage and toxicity, while it also shed light on the optimal treatment ordering of ADC versus TKI plus capecitabine regimens, and these findings were validated by new in vivo experiments. Our model is the first that mechanistically integrates multiple key drug modalities in HER2+ mBC research and it can serve as a high-throughput computational platform to guide future model-informed drug development and clinical translation.

Prediction of the liver safety profile of a first-in-class myeloperoxidase inhibitor using quantitative systems toxicology modeling.

The novel myeloperoxidase inhibitor verdiperstat was developed as a treatment for neuroinflammatory and neurodegenerative diseases. During development, a computational prediction of verdiperstat liver safety was performed using DILIsym v8A, a quantitative systems toxicology (QST) model of liver safety.A physiologically-based pharmacokinetic (PBPK) model of verdiperstat was constructed in GastroPlus 9.8, and outputs for liver and plasma time courses of verdiperstat were input into DILIsym. In vitro experiments measured the likelihood that verdiperstat would inhibit mitochondrial function, inhibit bile acid transporters, and generate reactive oxygen species (ROS); these results were used as inputs into DILIsym, with two alternate sets of parameters used in order to fully explore the sensitivity of model predictions. Verdiperstat dosing protocols up to 600 mg BID were simulated for up to 48 weeks using a simulated population (SimPops) in DILIsym.Verdiperstat was predicted to be safe, with only very rare, mild liver enzyme increases as a potential possibility in highly sensitive individuals. Subsequent Phase 3 clinical trials found that ALT elevations in the verdiperstat treatment group were generally similar to those in the placebo group. This validates the DILIsym simulation results and demonstrates the power of QST modelling to predict the liver safety profile of novel therapeutics.

Multiomic Analysis Reveals Disruption of Cholesterol Homeostasis by Cannabidiol in Human Cell Lines.

The nonpsychoactive cannabinoid, cannabidiol (CBD), is Food and Dug Administration approved for treatment of two drug-resistant epileptic disorders and is seeing increased use among the general public, yet the mechanisms that underlie its therapeutic effects and side-effect profiles remain unclear. Here, we report a systems-level analysis of CBD action in human cell lines using temporal multiomic profiling. FRET-based biosensor screening revealed that CBD elicits a sharp rise in cytosolic calcium, and activation of AMP-activated protein kinase in human keratinocyte and neuroblastoma cell lines. CBD treatment leads to alterations in the abundance of metabolites, mRNA transcripts, and proteins associated with activation of cholesterol biosynthesis, transport, and storage. We found that CBD rapidly incorporates into cellular membranes, alters cholesterol accessibility, and disrupts cholesterol-dependent membrane properties. Sustained treatment with high concentrations of CBD induces apoptosis in a dose-dependent manner. CBD-induced apoptosis is rescued by inhibition of cholesterol synthesis and potentiated by compounds that disrupt cholesterol trafficking and storage. Our data point to a pharmacological interaction of CBD with cholesterol homeostasis pathways, with potential implications in its therapeutic use.

An industry perspective on current QSP trends in drug development.

2023 marks the 10th anniversary of Natpara's submission to the US FDA, which led to the first recorded regulatory interaction where a decision was supported by Quantitative and Systems Pharmacology (QSP) simulations. It had taken about 5 years for the timid QSP discipline to emerge as an effective Model-Informed Drug Development (MIDD) tool with visible impact in the pharmaceutical industry. Since then, the presence of QSP in the regulatory environment has continued to increase, to the point that the Agency reported 60 QSP submissions in 2020 alone, representing ~ 4% of their annual IND submissions [1]. What sort of industry mindset has enabled QSP to reach this level of success? How does QSP fit within the MIDD paradigm? Does QSP mean the same to Discovery and to Clinical Development projects? How do 'platforms' compare to 'fit-for-purpose' QSP models in an industrial setting? Can QSP and empirical Pharmacokinetic-Pharmacodynamic (PKPD) modelling be complementary? What level of validation is required to inform drug development decisions? This article reflects on all these questions, in particular addressing those audiences with limited line-of-sight into the drug industry decision-making machinery.

A quantitative systems pharmacology model of plasma kallikrein-kinin system dysregulation in hereditary angioedema.

Hereditary angioedema (HAE) due to C1-inhibitor deficiency is a rare, debilitating, genetic disorder characterized by recurrent, unpredictable, attacks of edema. The clinical symptoms of HAE arise from excess bradykinin generation due to dysregulation of the plasma kallikrein-kinin system (KKS). A quantitative systems pharmacology (QSP) model that mechanistically describes the KKS and its role in HAE pathophysiology was developed based on HAE attacks being triggered by autoactivation of factor XII (FXII) to activated FXII (FXIIa), resulting in kallikrein production from prekallikrein. A base pharmacodynamic model was constructed and parameterized from literature data and ex vivo assays measuring inhibition of kallikrein activity in plasma of HAE patients or healthy volunteers who received lanadelumab. HAE attacks were simulated using a virtual patient population, with attacks recorded when systemic bradykinin levels exceeded 20 pM. The model was validated by comparing the simulations to observations from lanadelumab and plasma-derived C1-inhibitor clinical trials. The model was then applied to analyze the impact of nonadherence to a daily oral preventive therapy; simulations showed a correlation between the number of missed doses per month and reduced drug effectiveness. The impact of reducing lanadelumab dosing frequency from 300 mg every 2 weeks (Q2W) to every 4 weeks (Q4W) was also examined and showed that while attack rates with Q4W dosing were substantially reduced, the extent of reduction was greater with Q2W dosing. Overall, the QSP model showed good agreement with clinical data and could be used for hypothesis testing and outcome predictions.

Virtual clinical trials via a QSP immuno-oncology model to simulate the response to a conditionally activated PD-L1 targeting antibody in NSCLC.

Recently, immunotherapies for antitumoral response have adopted conditionally activated molecules with the objective of reducing systemic toxicity. Amongst these are conditionally activated antibodies, such as PROBODY® activatable therapeutics (Pb-Tx), engineered to be proteolytically activated by proteases found locally in the tumor microenvironment (TME). These PROBODY® therapeutics molecules have shown potential as PD-L1 checkpoint inhibitors in several cancer types, including both effectiveness and locality of action of the molecule as shown by several clinical trials and imaging studies. Here, we perform an exploratory study using our recently published quantitative systems pharmacology model, previously validated for triple-negative breast cancer (TNBC), to computationally predict the effectiveness and targeting specificity of a PROBODY® therapeutics drug compared to the non-modified antibody. We begin with the analysis of anti-PD-L1 immunotherapy in non-small cell lung cancer (NSCLC). As a first contribution, we have improved previous virtual patient selection methods using the omics data provided by the iAtlas database portal compared to methods previously published in literature. Furthermore, our results suggest that masking an antibody maintains its efficacy while improving the localization of active therapeutic in the TME. Additionally, we generalize the model by evaluating the dependence of the response to the tumor mutational burden, independently of cancer type, as well as to other key biomarkers, such as CD8/Treg Tcell and M1/M2 macrophage ratio. While our results are obtained from simulations on NSCLC, our findings are generalizable to other cancer types and suggest that an effective and highly selective conditionally activated PROBODY® therapeutics molecule is a feasible option.

Penicillin-Binding Protein 5/6 Acting as a Decoy Target in Pseudomonas aeruginosa Identified by Whole-Cell Receptor Binding and Quantitative Systems Pharmacology.

The ß-lactam antibiotics have been successfully used for decades to combat susceptible Pseudomonas aeruginosa, which has a notoriously difficult to penetrate outer membrane (OM). However, there is a dearth of data on target site penetration and covalent binding of penicillin-binding proteins (PBP) for ß-lactams and ß-lactamase inhibitors in intact bacteria. We aimed to determine the time course of PBP binding in intact and lysed cells and estimate the target site penetration and PBP access for 15 compounds in P. aeruginosa PAO1. All ß-lactams (at 2 × MIC) considerably bound PBPs 1 to 4 in lysed bacteria. However, PBP binding in intact bacteria was substantially attenuated for slow but not for rapid penetrating ß-lactams. Imipenem yielded 1.5 ± 0.11 log10 killing at 1h compared to <0.5 log10 killing for all other drugs. Relative to imipenem, the rate of net influx and PBP access was ~ 2-fold slower for doripenem and meropenem, 7.6-fold for avibactam, 14-fold for ceftazidime, 45-fold for cefepime, 50-fold for sulbactam, 72-fold for ertapenem, ~ 249-fold for piperacillin and aztreonam, 358-fold for tazobactam, ~547-fold for carbenicillin and ticarcillin, and 1,019-fold for cefoxitin. At 2 × MIC, the extent of PBP5/6 binding was highly correlated (r2 = 0.96) with the rate of net influx and PBP access, suggesting that PBP5/6 acted as a decoy target that should be avoided by slowly penetrating, future ß-lactams. This first comprehensive assessment of the time course of PBP binding in intact and lysed P. aeruginosa explained why only imipenem killed rapidly. The developed novel covalent binding assay in intact bacteria accounts for all expressed resistance mechanisms.

Virus-CKB: an integrated bioinformatics platform and analysis resource for COVID-19 research.

Given the scale and rapid spread of the coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), there is an urgent need for medicines that can help before vaccines are available. In this study, we present a viral-associated disease-specific chemogenomics knowledgebase (Virus-CKB) and apply our computational systems pharmacology-target mapping to rapidly predict the FDA-approved drugs which can quickly progress into clinical trials to meet the urgent demand of the COVID-19 outbreak. Virus-CKB reuses the underlying platform of our DAKB-GPCRs but adds new features like multiple-compound support, multi-cavity protein support and customizable symbol display. Our one-stop computing platform describes the chemical molecules, genes and proteins involved in viral-associated diseases regulation. To date, Virus-CKB archived 65 antiviral drugs in the market, 107 viral-related targets with 189 available 3D crystal or cryo-EM structures and 2698 chemical agents reported for these target proteins. Moreover, Virus-CKB is implemented with web applications for the prediction of the relevant protein targets and analysis and visualization of the outputs, including HTDocking, TargetHunter, BBB predictor, NGL Viewer, Spider Plot, etc. The Virus-CKB server is accessible at https://www.cbligand.org/g/virus-ckb.

Systems pharmacology: a combination strategy for improving efficacy of PD-1/PD-L1 blockade.

Targeting tumor microenvironment (TME), such as immune checkpoint blockade (ICB), has achieved increased overall response rates in many advanced cancers, such as non-small cell lung cancer (NSCLC), however, only in a fraction of patients. To improve the overall and durable response rates, combining other therapeutics, such as natural products, with ICB therapy is under investigation. Unfortunately, due to the lack of systematic methods to characterize the relationship between TME and ICB, development of rational immune-combination therapy is a critical challenge. Here, we proposed a systems pharmacology strategy to identify resistance regulators of PD-1/PD-L1 blockade and develop its combinatorial drug by integrating multidimensional omics and pharmacological methods. First, a high-resolution TME cell atlas was inferred from bulk sequencing data by referring to a high-resolution single-cell data and was used to predict potential resistance regulators of PD-1/PD-L1 blockade through TME stratification analysis. Second, to explore the drug targeting the resistance regulator, we carried out the large-scale target fishing and the network analysis between multi-target drug and the resistance regulator. Finally, we predicted and verified that oxymatrine significantly enhances the infiltration of CD8+ T cells into TME and is a powerful combination agent to enhance the therapeutic effect of anti-PD-L1 in a mouse model of lung adenocarcinoma. Overall, the systems pharmacology strategy offers a paradigm to identify combinatorial drugs for ICB therapy with a systems biology perspective of drug-target-pathway-TME phenotype-ICB combination.

Clarifying the mechanism of apigenin against blood-brain barrier disruption in ischemic stroke using systems pharmacology.

Currently, recombinant tissue plasminogen activator (rtPA) is an effective therapy for ischemic stroke (IS). However, blood-brain barrier (BBB) disruption is a serious side effect of rtPA therapy and may lead to patients' death. The natural polyphenol apigenin has a good therapeutic effect on IS. Apigenin has potential BBB protection, but the mechanism by which it protects the BBB integrity is not clear. In this study, we used network pharmacology, bioinformatics, molecular docking and molecular dynamics simulation to reveal the mechanisms by which apigenin protects the BBB. Among the 146 targets of apigenin for the treatment of IS, 20 proteins were identified as core targets (e.g., MMP-9, TLR4, STAT3). Apigenin protects BBB integrity by inhibiting the activity of MMPs through anti-inflammation and anti-oxidative stress. These mechanisms included JAK/STAT, the toll-like receptor signaling pathway, and Nitrogen metabolism signaling pathways. The findings of this study contribute to a more comprehensive understanding of the mechanism of apigenin in the treatment of BBB disruption and provide ideas for the development of drugs to treat IS.

Towards a platform quantitative systems pharmacology (QSP) model for preclinical to clinical translation of antibody drug conjugates (ADCs).

A next generation multiscale quantitative systems pharmacology (QSP) model for antibody drug conjugates (ADCs) is presented, for preclinical to clinical translation of ADC efficacy. Two HER2 ADCs (trastuzumab-DM1 and trastuzumab-DXd) were used for model development, calibration, and validation. The model integrates drug specific experimental data including in vitro cellular disposition data, pharmacokinetic (PK) and tumor growth inhibition (TGI) data for T-DM1 and T-DXd, as well as system specific data such as properties of HER2, tumor growth rates, and volumes. The model incorporates mechanistic detail at the intracellular level, to account for different mechanisms of ADC processing and payload release. It describes the disposition of the ADC, antibody, and payload inside and outside of the tumor, including binding to off-tumor, on-target sinks. The resulting multiscale PK model predicts plasma and tumor concentrations of ADC and payload. Tumor payload concentrations predicted by the model were linked to a TGI model and used to describe responses following ADC administration to xenograft mice. The model was translated to humans and virtual clinical trial simulations were performed that successfully predicted progression free survival response for T-DM1 and T-DXd for the treatment of HER2+ metastatic breast cancer, including differential efficacy based upon HER2 expression status. In conclusion, the presented model is a step toward a platform QSP model and strategy for ADCs, integrating multiple types of data and knowledge to predict ADC efficacy. The model has potential application to facilitate ADC design, lead candidate selection, and clinical dosing schedule optimization.