The relevance and potential roles of microphysiological systems in biology and medicine

Microphysiological systems (MPS), consisting of interacting organs-on-chips or tissue-engineered, 3D organ constructs that use human cells, present an opportunity to bring new tools to biology, medicine, pharmacology, physiology, and toxicology. This issue of Experimental Biology and Medicine describes the ongoing development of MPS that can serve as in-vitro models for bone and cartilage, brain, gastrointestinal tract, lung, liver, microvasculature, reproductive tract, skeletal muscle, and skin. Related topics addressed here are the interconnection of organs-on-chips to support physiologically based pharmacokinetics and drug discovery and screening, and the microscale technologies that regulate stem cell differentiation. The initial motivation for creating MPS was to increase the speed, efficiency, and safety of pharmaceutical development and testing, paying particular regard to the fact that neither monolayer monocultures of immortal or primary cell lines nor animal studies can adequately recapitulate the dynamics of drug-organ, drug-drug, and drug-organ-organ interactions in humans. Other applications include studies of the effect of environmental toxins on humans, identification, characterization, and neutralization of chemical and biological weapons, controlled studies of the microbiome and infectious disease that cannot be conducted in humans, controlled differentiation of induced pluripotent stem cells into specific adult cellular phenotypes, and studies of the dynamics of metabolism and signaling within and between human organs. The technical challenges are being addressed by many investigators, and in the process, it seems highly likely that significant progress will be made toward providing more physiologically realistic alternatives to monolayer monocultures or whole animal studies. The effectiveness of this effort will be determined in part by how easy the constructs are to use, how well they function, how accurately they recapitulate and report human pharmacology and toxicology, whether they can be generated in large numbers to enable parallel studies, and if their use can be standardized consistent with the practices of regulatory science.

Pharmacodynamic Drug-Drug Interactions

Pharmacodynamic drug-drug interactions (DDIs) occur when the pharmacological effect of one drug is altered by that of another drug in a combination regimen. DDIs often are classified as synergistic, additive, or antagonistic in nature, albeit these terms are frequently misused. Within a complex pathophysiological system, the mechanism of interaction may occur at the same target or through alternate pathways. Quantitative evaluation of pharmacodynamic DDIs by employing modeling and simulation approaches is needed to identify and optimize safe and effective combination therapy regimens. This review investigates the opportunities and challenges in pharmacodynamic DDI studies and highlights examples of quantitative methods for evaluating pharmacodynamic DDIs, with a particular emphasis on the use of mechanism-based modeling and simulation in DDI studies. Advancements in both experimental and computational techniques will enable the application of better, model-informed assessments of pharmacodynamic DDIs in drug discovery, development, and therapeutics.

Biologically Relevant Heterogeneity: Metrics and Practical Insights

Heterogeneity is a fundamental property of biological systems at all scales that must be addressed in a wide range of biomedical applications, including basic biomedical research, drug discovery, diagnostics, and the implementation of precision medicine. There are a number of published approaches to characterizing heterogeneity in cells in vitro and in tissue sections. However, there are no generally accepted approaches for the detection and quantitation of heterogeneity that can be applied in a relatively high-throughput workflow. This review and perspective emphasizes the experimental methods that capture multiplexed cell-level data, as well as the need for standard metrics of the spatial, temporal, and population components of heterogeneity. A recommendation is made for the adoption of a set of three heterogeneity indices that can be implemented in any high-throughput workflow to optimize the decision-making process. In addition, a pairwise mutual information method is suggested as an approach to characterizing the spatial features of heterogeneity, especially in tissue-based imaging. Furthermore, metrics for temporal heterogeneity are in the early stages of development. Example studies indicate that the analysis of functional phenotypic heterogeneity can be exploited to guide decisions in the interpretation of biomedical experiments, drug discovery, diagnostics, and the design of optimal therapeutic strategies for individual patients.

A Computational Model of Neoadjuvant PD-1 Inhibition in Non-Small Cell Lung Cancer

Immunotherapy and immune checkpoint blocking antibodies such as anti-PD-1 are approved and significantly improve the survival of advanced non-small cell lung cancer (NSCLC) patients, but there has been little success in identifying biomarkers capable of separating the responders from non-responders before the onset of the therapy. In this study, we developed a quantitative system pharmacology (QSP) model to represent the anti-tumor immune response in human NSCLC that integrated our knowledge of tumor growth, antigen processing and presentation, T cell activation and distribution, antibody pharmacokinetics, and immune checkpoint dynamics. The model was calibrated with the available data and was used to identify potential biomarkers as well as patient-specific response based on the patient parameters. The model predicted that in addition to tumor mutational burden (TMB), a known biomarker for anti-PD-1 therapy in NSCLC, the number of effector T cells and regulatory T cells in the tumor and blood is a predictor of the responders. Furthermore, the model simulated a set of 12 patients with known TMB and MHC/antigen-binding affinity from a recent clinical trial (ClinicalTrials.gov number, NCT02259621) on neoadjuvant nivolumab therapy in resectable lung cancer and predicted an augmented durable response in patients with adjuvant nivolumab treatment in addition to the clinical trial protocol of neoadjuvant nivolumab treatment followed by resection. Overall, the model provides a valuable framework to model tumor immunity and response to immune checkpoint blockers to enhance biomarker discovery and performing virtual clinical trials to aid in design and interpretation of the current trials with fewer patients.

QSP Toolbox: Computational Implementation of Integrated Workflow Components for Deploying Multi-Scale Mechanistic Models

Quantitative systems pharmacology (QSP) modeling has become increasingly important in pharmaceutical research and development, and is a powerful tool to gain mechanistic insights into the complex dynamics of biological systems in response to drug treatment. However, even once a suitable mathematical framework to describe the pathophysiology and mechanisms of interest is established, final model calibration and the exploration of variability can be challenging and time consuming. QSP models are often formulated as multi-scale, multi-compartment nonlinear systems of ordinary differential equations. Commonly accepted modeling strategies, workflows, and tools have promise to greatly improve the efficiency of QSP methods and improve productivity. In this paper, we present the QSP Toolbox, a set of functions, structure array conventions, and class definitions that computationally implement critical elements of QSP workflows including data integration, model calibration, and variability exploration. We present the application of the toolbox to an ordinary differential equations-based model for antibody drug conjugates. As opposed to a single stepwise reference model calibration, the toolbox also facilitates simultaneous parameter optimization and variation across multiple in vitro, in vivo, and clinical assays to more comprehensively generate alternate mechanistic hypotheses that are in quantitative agreement with available data. The toolbox also includes scripts for developing and applying virtual populations to mechanistic exploration of biomarkers and efficacy. We anticipate that the QSP Toolbox will be a useful resource that will facilitate implementation, evaluation, and sharing of new methodologies in a common framework that will greatly benefit the community.

Systems Pharmacology: Defining the Interactions of Drug Combinations

The majority of diseases are associated with alterations in multiple molecular pathways and complex interactions at the cellular and organ levels. Single-target monotherapies therefore have intrinsic limitations with respect to their maximum therapeutic benefits. The potential of combination drug therapies has received interest for the treatment of many diseases and is well established in some areas, such as oncology. Combination drug treatments may allow us to identify synergistic drug effects, reduce adverse drug reactions, and address variability in disease characteristics between patients. Identification of combination therapies remains challenging. We discuss current state-of-the-art systems pharmacology approaches to enable rational identification of combination therapies. These approaches, which include characterization of mechanisms of disease and drug action at a systems level, can enable understanding of drug interactions at the molecular, cellular, physiological, and organismal levels. Such multiscale understanding can enable precision medicine by promoting the rational development of combination therapy at the level of individual patients for many diseases.

A Translational Quantitative Systems Pharmacology Model for CD3 Bispecific Molecules: Application to Quantify T Cell-Mediated Tumor Cell Killing by P-Cadherin LP DART®

CD3 bispecific antibody constructs recruit cytolytic T cells to kill tumor cells, offering a potent approach to treat cancer. T cell activation is driven by the formation of a trimolecular complex (trimer) between drugs, T cells, and tumor cells, mimicking an immune synapse. A translational quantitative systems pharmacology (QSP) model is proposed for CD3 bispecific molecules capable of predicting trimer concentration and linking it to tumor cell killing. The model was used to quantify the pharmacokinetic (PK)/pharmacodynamic (PD) relationship of a CD3 bispecific targeting P-cadherin (PF-06671008). It describes the disposition of PF-06671008 in the central compartment and tumor in mouse xenograft models, including binding to target and T cells in the tumor to form the trimer. The model incorporates T cell distribution to the tumor, proliferation, and contraction. PK/PD parameters were estimated for PF-06671008 and a tumor stasis concentration (TSC) was calculated as an estimate of minimum efficacious trimer concentration. TSC values ranged from 0.0092 to 0.064 pM across mouse tumor models. The model was translated to the clinic and used to predict the disposition of PF-06671008 in patients, including the impact of binding to soluble P-cadherin. The predicted terminal half-life of PF-06671008 in the clinic was approximately 1 day, and P-cadherin expression and number of T cells in the tumor were shown to be sensitive parameters impacting clinical efficacy. A translational QSP model is presented for CD3 bispecific molecules, which integrates in silico, in vitro and in vivo data in a mechanistic framework, to quantify and predict efficacy across species.

A Perspective on Implementing a Quantitative Systems Pharmacology Platform for Drug Discovery and the Advancement of Personalized Medicine

Drug candidates exhibiting well-defined pharmacokinetic and pharmacodynamic profiles that are otherwise safe often fail to demonstrate proof-of-concept in phase II and III trials. Innovation in drug discovery and development has been identified as a critical need for improving the efficiency of drug discovery, especially through collaborations between academia, government agencies, and industry. To address the innovation challenge, we describe a comprehensive, unbiased, integrated, and iterative quantitative systems pharmacology (QSP)-driven drug discovery and development strategy and platform that we have implemented at the University of Pittsburgh Drug Discovery Institute. Intrinsic to QSP is its integrated use of multiscale experimental and computational methods to identify mechanisms of disease progression and to test predicted therapeutic strategies likely to achieve clinical validation for appropriate subpopulations of patients. The QSP platform can address biological heterogeneity and anticipate the evolution of resistance mechanisms, which are major challenges for drug development. The implementation of this platform is dedicated to gaining an understanding of mechanism(s) of disease progression to enable the identification of novel therapeutic strategies as well as repurposing drugs. The QSP platform will help promote the paradigm shift from reactive population-based medicine to proactive personalized medicine by focusing on the patient as the starting and the end point.

Conducting a Virtual Clinical Trial in HER2-Negative Breast Cancer Using a Quantitative Systems Pharmacology Model With an Epigenetic Modulator and Immune Checkpoint Inhibitors

The survival rate of patients with breast cancer has been improved by immune checkpoint blockade therapies, and the efficacy of their combinations with epigenetic modulators has shown promising results in preclinical studies. In this prospective study, we propose an ordinary differential equation (ODE)-based quantitative systems pharmacology (QSP) model to conduct an in silico virtual clinical trial and analyze potential predictive biomarkers to improve the anti-tumor response in HER2-negative breast cancer. The model is comprised of four compartments: central, peripheral, tumor, and tumor-draining lymph node, and describes immune activation, suppression, T cell trafficking, and pharmacokinetics and pharmacodynamics (PK/PD) of the therapeutic agents. We implement theoretical mechanisms of action for checkpoint inhibitors and the epigenetic modulator based on preclinical studies to investigate their effects on anti-tumor response. According to model-based simulations, we confirm the synergistic effect of the epigenetic modulator and that pre-treatment tumor mutational burden, tumor-infiltrating effector T cell (Teff) density, and Teff to regulatory T cell (Treg) ratio are significantly higher in responders, which can be potential biomarkers to be considered in clinical trials. Overall, we present a readily reproducible modular model to conduct in silico virtual clinical trials on patient cohorts of interest, which is a step toward personalized medicine in cancer immunotherapy.

Paths to Successful Translation of New Therapies for Severe Traumatic Brain Injury in the Golden Age of Traumatic Brain Injury Research: A Pittsburgh Vision

New neuroprotective therapies for severe traumatic brain injury (TBI) have not translated from pre-clinical to clinical success. Numerous explanations have been suggested in both the pre-clinical and clinical arenas. Coverage of TBI in the lay press has reinvigorated interest, creating a golden age of TBI research with innovative strategies to circumvent roadblocks. We discuss the need for more robust therapies. We present concepts for traditional and novel approaches to defining therapeutic targets. We review lessons learned from the ongoing work of the pre-clinical drug and biomarker screening consortium Operation Brain Trauma Therapy and suggest ways to further enhance pre-clinical consortia. Biomarkers have emerged that empower choice and assessment of target engagement by candidate therapies. Drug combinations may be needed, and it may require moving beyond conventional drug therapies. Precision medicine may also link the right therapy to the right patient, including new approaches to TBI classification beyond the Glasgow Coma Scale or anatomical phenotyping-incorporating new genetic and physiologic approaches. Therapeutic breakthroughs may also come from alternative approaches in clinical investigation (comparative effectiveness, adaptive trial design, use of the electronic medical record, and big data). The full continuum of care must also be represented in translational studies, given the important clinical role of pre-hospital events, extracerebral insults in the intensive care unit, and rehabilitation. TBI research from concussion to coma can cross-pollinate and further advancement of new therapies. Misconceptions can stifle/misdirect TBI research and deserve special attention. Finally, we synthesize an approach to deliver therapeutic breakthroughs in this golden age of TBI research.

Strategies for clinical dose optimization of T cell-engaging therapies in oncology

Innovative approaches in the design of T cell-engaging (TCE) molecules are ushering in a new wave of promising immunotherapies for the treatment of cancer. Their mechanism of action, which generates an in trans interaction to create a synthetic immune synapse, leads to complex and interconnected relationships between the exposure, efficacy, and toxicity of these drugs. Challenges thus arise when designing optimal clinical dose regimens for TCEs with narrow therapeutic windows, with a variety of dosing strategies being evaluated to mitigate key side effects such as cytokine release syndrome, neurotoxicity, and on-target off-tumor toxicities. This review evaluates the current approaches to dose optimization throughout the preclinical and clinical development of TCEs, along with perspectives for improvement of these strategies. Quantitative approaches used to aid the understanding of dose-exposure-response relationships are highlighted, along with opportunities to guide the rational design of next-generation TCE molecules, and optimize their dose regimens in patients.

Mechanistic Systems Modeling to Improve Understanding and Prediction of Cardiotoxicity Caused by Targeted Cancer Therapeutics

Tyrosine kinase inhibitors (TKIs) are highly potent cancer therapeutics that have been linked with serious cardiotoxicity, including left ventricular dysfunction, heart failure, and QT prolongation. TKI-induced cardiotoxicity is thought to result from interference with tyrosine kinase activity in cardiomyocytes, where these signaling pathways help to control critical processes such as survival signaling, energy homeostasis, and excitation–contraction coupling. However, mechanistic understanding is limited at present due to the complexities of tyrosine kinase signaling, and the wide range of targets inhibited by TKIs. Here, we review the use of TKIs in cancer and the cardiotoxicities that have been reported, discuss potential mechanisms underlying cardiotoxicity, and describe recent progress in achieving a more systematic understanding of cardiotoxicity via the use of mechanistic models. In particular, we argue that future advances are likely to be enabled by studies that combine large-scale experimental measurements with Quantitative Systems Pharmacology (QSP) models describing biological mechanisms and dynamics. As such approaches have proven extremely valuable for understanding and predicting other drug toxicities, it is likely that QSP modeling can be successfully applied to cardiotoxicity induced by TKIs. We conclude by discussing a potential strategy for integrating genome-wide expression measurements with models, illustrate initial advances in applying this approach to cardiotoxicity, and describe challenges that must be overcome to truly develop a mechanistic and systematic understanding of cardiotoxicity caused by TKIs.

FDA-Industry Scientific Exchange on assessing quantitative systems pharmacology models in clinical drug development: a meeting report, summary of challenges/gaps, and future perspective

The pharmaceutical industry is actively applying quantitative systems pharmacology (QSP) to make internal decisions and guide drug development. To facilitate the eventual development of a common framework for assessing the credibility of QSP models for clinical drug development, scientists from US Food and Drug Administration and the pharmaceutical industry organized a full-day virtual Scientific Exchange on July 1, 2020. An assessment form was used to ensure consistency in the evaluation process. Among the cases presented, QSP was applied to various therapeutic areas. Applications mostly focused on phase 2 dose selection. Model transparency, including details on expert knowledge and data used for model development, was identified as a major factor for robust model assessment. The case studies demonstrated some commonalities in the workflow of QSP model development, calibration, and validation but differ in the size, scope, and complexity of QSP models, in the acceptance criteria for model calibration and validation, and in the algorithms/approaches used for creating virtual patient populations. Though efforts are being made to build the credibility of QSP models and the confidence is increasing in applying QSP for internal decisions at the clinical stages of drug development, there are still many challenges facing QSP application to late stage drug development. The QSP community needs a strategic plan that includes the ability and flexibility to Adapt, to establish Common expectations for model Credibility needed to inform drug Labeling and patient care, and to AIM to achieve the goal (ACCLAIM).

Quantitative Systems Pharmacology: A Framework for Context

Quantitative Systems Pharmacology (QSP) is receiving increased attention. As the momentum builds and the expectations grow it is important to (re)assess and formalize the basic concepts and approaches. In this short review, I argue that QSP, in addition to enabling the rational integration of data and development of complex models, maybe more importantly, provides the foundations for developing an integrated framework for the assessment of drugs and their impact on disease within a broader context expanding the envelope to account in great detail for physiology, environment and prior history. I articulate some of the critical enablers, major obstacles and exciting opportunities manifesting themselves along the way. Charting such overarching themes will enable practitioners to identify major and defining factors as the field progressively moves towards personalized and precision health care delivery.

A mathematical model to identify optimal combinations of drug targets for dupilumab poor responders in atopic dermatitis

BACKGROUND

Several biologics for atopic dermatitis (AD) have demonstrated good efficacy in clinical trials, but with a substantial proportion of patients being identified as poor responders. This study aims to understand the pathophysiological backgrounds of patient variability in drug response, especially for dupilumab, and to identify promising drug targets in dupilumab poor responders.

METHODS

We conducted model-based meta-analysis of recent clinical trials of AD biologics and developed a mathematical model that reproduces reported clinical efficacies for nine biological drugs (dupilumab, lebrikizumab, tralokinumab, secukinumab, fezakinumab, nemolizumab, tezepelumab, GBR 830, and recombinant interferon-gamma) by describing system-level AD pathogenesis. Using this model, we simulated the clinical efficacy of hypothetical therapies on virtual patients.

RESULTS

Our model reproduced reported time courses of %improved EASI and EASI-75 of the nine drugs. The global sensitivity analysis and model simulation indicated the baseline level of IL-13 could stratify dupilumab good responders. Model simulation on the efficacies of hypothetical therapies revealed that simultaneous inhibition of IL-13 and IL-22 was effective, whereas application of the nine biologic drugs was ineffective, for dupilumab poor responders (EASI-75 at 24 weeks: 21.6% vs. max. 1.9%).

CONCLUSION

Our model identified IL-13 as a potential predictive biomarker to stratify dupilumab good responders, and simultaneous inhibition of IL-13 and IL-22 as a promising drug therapy for dupilumab poor responders. This model will serve as a computational platform for model-informed drug development for precision medicine, as it allows evaluation of the effects of new potential drug targets and the mechanisms behind patient variability in drug response.

Development of bispecific T cell engagers: harnessing quantitative systems pharmacology

Bispecific T cell engagers (bsTCEs) have emerged as a promising class of cancer immunotherapy. Several bsTCEs have achieved marketing approval; dozens more are under clinical investigation. However, the clinical development of bsTCEs remains rife with challenges, including nuanced pharmacology, limited translatability of preclinical findings, frequent on-target toxicity, and convoluted dosing regimens. In this opinion article we present a distinct perspective on how quantitative systems pharmacology (QSP) can serve as a powerful tool for overcoming these obstacles. Recent advances in QSP modeling have empowered developers of bsTCEs to gain a deeper understanding of their context-dependent pharmacology, bridge gaps in experimental data, guide first-in-human (FIH) dose selection, design dosing regimens with expanded therapeutic windows, and improve long-term treatment outcomes. We use recent case studies to exemplify the potential of QSP techniques to support future bsTCE development.

Virtual Systems Pharmacology (ViSP) software for simulation from mechanistic systems-level models

Multiple software programs are available for designing and running large scale system-level pharmacology models used in the drug development process. Depending on the problem, scientists may be forced to use several modeling tools that could increase model development time, IT costs and so on. Therefore, it is desirable to have a single platform that allows setting up and running large-scale simulations for the models that have been developed with different modeling tools. We developed a workflow and a software platform in which a model file is compiled into a self-contained executable that is no longer dependent on the software that was used to create the model. At the same time the full model specifics is preserved by presenting all model parameters as input parameters for the executable. This platform was implemented as a model agnostic, therapeutic area agnostic and web-based application with a database back-end that can be used to configure, manage and execute large-scale simulations for multiple models by multiple users. The user interface is designed to be easily configurable to reflect the specifics of the model and the user's particular needs and the back-end database has been implemented to store and manage all aspects of the systems, such as Models, Virtual Patients, User Interface Settings, and Results. The platform can be adapted and deployed on an existing cluster or cloud computing environment. Its use was demonstrated with a metabolic disease systems pharmacology model that simulates the effects of two antidiabetic drugs, metformin and fasiglifam, in type 2 diabetes mellitus patients.

Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis

Autophagy plays an essential role in cell survival/death and functioning. Modulation of autophagy has been recognized as a promising therapeutic strategy against diseases/disorders associated with uncontrolled growth or accumulation of biomolecular aggregates, organelles, or cells including those caused by cancer, aging, neurodegeneration, and liver diseases such as α1-antitrypsin deficiency. Numerous pharmacological agents that enhance or suppress autophagy have been discovered. However, their molecular mechanisms of action are far from clear. Here, we collected a set of 225 autophagy modulators and carried out a comprehensive quantitative systems pharmacology (QSP) analysis of their targets using both existing databases and predictions made by our machine learning algorithm. Autophagy modulators include several highly promiscuous drugs (e.g., artenimol and olanzapine acting as activators, fostamatinib as an inhibitor, or melatonin as a dual-modulator) as well as selected drugs that uniquely target specific proteins (~30% of modulators). They are mediated by three layers of regulation: (i) pathways involving core autophagy-related (ATG) proteins such as mTOR, AKT, and AMPK; (ii) upstream signaling events that regulate the activity of ATG pathways such as calcium-, cAMP-, and MAPK-signaling pathways; and (iii) transcription factors regulating the expression of ATG proteins such as TFEB, TFE3, HIF-1, FoxO, and NF-κB. Our results suggest that PKA serves as a linker, bridging various signal transduction events and autophagy. These new insights contribute to a better assessment of the mechanism of action of autophagy modulators as well as their side effects, development of novel polypharmacological strategies, and identification of drug repurposing opportunities.

Quantitative Systems Pharmacology Model-Based Predictions of Clinical Endpoints to Optimize Warfarin and Rivaroxaban Anti-Thrombosis Therapy

Background

Tight monitoring of efficacy and safety of anticoagulants such as warfarin is imperative to optimize the benefit-risk ratio of anticoagulants in patients. The standard tests used are measurements of prothrombin time (PT), usually expressed as international normalized ratio (INR), and activated partial thromboplastin time (aPTT).

Objective

To leverage a previously developed quantitative systems pharmacology (QSP) model of the human coagulation network to predict INR and aPTT for warfarin and rivaroxaban, respectively.

Methods

A modeling and simulation approach was used to predict INR and aPTT measurements of patients receiving steady-state anticoagulation therapy. A previously developed QSP model was leveraged for the present analysis. The effect of genetic polymorphisms known to influence dose response of warfarin (CYP2C9, VKORC1) were implemented into the model by modifying warfarin clearance (CYP2C9 *1: 0.2 L/h; *2: 0.14 L/h, *3: 0.04 L/h) and the concentration of available vitamin K (VKORC1 GA: −22% from normal vitamin K concentration; AA: −44% from normal vitamin K concentration). Virtual patient populations were used to assess the ability of the model to accurately predict routine INR and aPTT measurements from patients under long-term anticoagulant therapy.

Results

The introduced model accurately described the observed INR of patients receiving long-term warfarin treatment. The model was able to demonstrate the influence of genetic polymorphisms of CYP2C9 and VKORC1 on the INR measurements. Additionally, the model was successfully used to predict aPTT measurements for patients receiving long-term rivaroxaban therapy.

Conclusion

The QSP model accurately predicted INR and aPTT measurements observed during routine therapeutic drug monitoring. This is an exemplar of how a QSP model can be adapted and used as a model-based precision dosing tool during clinical practice and drug development to predict efficacy and safety of anticoagulants to ultimately help optimize anti-thrombotic therapy.

Development of Prodrugs for PDT-Based Combination Therapy Using a Singlet-Oxygen-Sensitive Linker and Quantitative Systems Pharmacology

Photodynamic therapy (PDT) has become an effective treatment for certain types of solid tumors. The combination of PDT with other therapies has been extensively investigated in recent years to improve its effectiveness and expand its applications. This focused review summarizes the development of a prodrug system in which anticancer drugs are activated locally at tumor sites during PDT treatment. The development of a singlet-oxygen-sensitive linker that can be conveniently conjugated to various drugs and efficiently cleaved to release intact drugs is recapitulated. The initial design of prodrugs, preliminary efficacy evaluation, pharmacokinetics study, and optimization using quantitative systems pharmacology is discussed. Current treatment optimization in animal models using physiologically based a pharmacokinetic (PBPK) modeling approach is also explored.

Inflammation as a Risk Factor in Cardiotoxicity: An Important Consideration for Screening During Drug Development

Numerous commonly prescribed drugs, including antiarrhythmics, antihistamines, and antibiotics, carry a proarrhythmic risk and may induce dangerous arrhythmias, including the potentially fatal Torsades de Pointes. For this reason, cardiotoxicity testing has become essential in drug development and a required step in the approval of any medication for use in humans. Blockade of the hERG K⁺ channel and the consequent prolongation of the QT interval on the ECG have been considered the gold standard to predict the arrhythmogenic risk of drugs. In recent years, however, preclinical safety pharmacology has begun to adopt a more integrative approach that incorporates mathematical modeling and considers the effects of drugs on multiple ion channels. Despite these advances, early stage drug screening research only evaluates QT prolongation in experimental and computational models that represent healthy individuals. We suggest here that integrating disease modeling with cardiotoxicity testing can improve drug risk stratification by predicting how disease processes and additional comorbidities may influence the risks posed by specific drugs. In particular, chronic systemic inflammation, a condition associated with many diseases, affects heart function and can exacerbate medications’ cardiotoxic effects. We discuss emerging research implicating the role of inflammation in cardiac electrophysiology, and we offer a perspective on how in silico modeling of inflammation may lead to improved evaluation of the proarrhythmic risk of drugs at their early stage of development.

The timing of cyclic cytotoxic chemotherapy can worsen neutropenia and neutrophilia

Despite recent advances in immunotherapies, cytotoxic chemotherapy continues to be a first-line treatment option for the majority of cancers. Unfortunately, a common side effect in patients undergoing chemotherapy treatment is neutropenia. To mitigate the risk of neutropenia and febrile neutropenia, prophylactic treatment with granulocyte-colony stimulating factor (G-CSF) is administered. Extensive pharmacokinetic/pharmacodynamic modelling of myelosuppression during chemotherapy has suggested avenues for therapy optimization to mitigate this neutropenia. However, the issue of resonance, whereby neutrophil oscillations are induced by the periodic administration of cytotoxic chemotherapy and the coadministration of G-CSF, potentially aggravating a patient's neutropenic/neutrophilic status, is not well-characterized in the clinical literature. Here, through analysis of neutrophil data from young acute lymphoblastic leukaemia patients, we find that resonance is occurring during cyclic chemotherapy treatment in 26% of these patients. Motivated by these data and our previous modelling studies on adult lymphoma patients, we examined resonance during treatment with or without G-CSF. Using our quantitative systems pharmacology model of granulopoiesis, we show that the timing of cyclic chemotherapy can worsen neutropenia or neutrophilia, and suggest clinically-actionable schedules to reduce the resonant effect. We emphasize that delaying supportive G-CSF therapy to 6-7 days after chemotherapy can mitigate myelosuppressive effects. This study therefore highlights the importance of quantitative systems pharmacology for the clinical practice for developing rational therapeutic strategies.

Quantitative Systems Pharmacological Analysis of Drugs of Abuse Reveals the Pleiotropy of Their Targets and the Effector Role of mTORC1

Existing treatments against drug addiction are often ineffective due to the complexity of the networks of protein-drug and protein-protein interactions (PPIs) that mediate the development of drug addiction and related neurobiological disorders. There is an urgent need for understanding the molecular mechanisms that underlie drug addiction toward designing novel preventive or therapeutic strategies. The rapidly accumulating data on addictive drugs and their targets as well as advances in machine learning methods and computing technology now present an opportunity to systematically mine existing data and draw inferences on potential new strategies. To this aim, we carried out a comprehensive analysis of cellular pathways implicated in a diverse set of 50 drugs of abuse using quantitative systems pharmacology methods. The analysis of the drug/ligand-target interactions compiled in DrugBank and STITCH databases revealed 142 known and 48 newly predicted targets, which have been further analyzed to identify the KEGG pathways enriched at different stages of drug addiction cycle, as well as those implicated in cell signaling and regulation events associated with drug abuse. Apart from synaptic neurotransmission pathways detected as upstream signaling modules that “sense” the early effects of drugs of abuse, pathways involved in neuroplasticity are distinguished as determinants of neuronal morphological changes. Notably, many signaling pathways converge on important targets such as mTORC1. The latter emerges as a universal effector of the persistent restructuring of neurons in response to continued use of drugs of abuse.

Mechanistic evaluation of the effect of sodium-dependent glucose transporter 2 inhibitors on delayed glucose absorption in patients with type 2 diabetes mellitus using a quantitative systems pharmacology model of human systemic glucose dynamics

Sodium-dependent glucose transporter (SGLT) 2 is specifically expressed in the kidney, while SGLT1 is present in the kidneys and small intestine. SGLT2 inhibitors are a class of oral antidiabetic drugs that lower elevated plasma glucose levels by promoting the urinary excretion of excess glucose through the inhibition of renal glucose reuptake. The inhibition selectivity for SGLT2 over SGLT1 (SGLT2/1 selectivity) of marketed SGLT2 inhibitors is diverse, while SGLT2/1 selectivity of canagliflozin is relatively low. Although canagliflozin suppresses postprandial glucose levels, the degree of contribution for SGLT1 inhibition to this effect remains unproven. To analyze the effect of SGLT2 inhibitors on postprandial glucose level, we constructed a novel quantitative systems pharmacology (QSP) model, called human systemic glucose dynamics (HSGD) model, integrating intestinal absorption, metabolism, and renal reabsorption of glucose. This HSGD model reproduced the postprandial plasma glucose concentration-time profiles during a meal tolerance test under different clinical trial conditions. Simulations after canagliflozin administration showed a dose-dependent delay of time (T_max,glc ) to reach maximum concentration of glucose (C_max,glc ), and the delay of T_max,glc disappeared when inhibition of SGLT1 was negated. In addition, contribution ratio of intestinal SGLT1 inhibition to the decrease in C_max,glc was estimated to be 23%-28%, when 100 and 300 mg of canagliflozin are administered. This HSGD model enabled us to provide the partial contribution of intestinal SGLT1 inhibition to the improvement of postprandial hyperglycemia as well as to quantitatively describe the plasma glucose dynamics following SGLT2 inhibitors.

Evaluation of therapeutic strategies targeting BCAA catabolism using a systems pharmacology model

Background: Abnormal branched-chained amino acids (BCAA) accumulation in cardiomyocytes is associated with cardiac remodeling in heart failure. Administration of branched-chain α-keto acid dehydrogenase (BCKD) kinase inhibitor BT2 has been shown to reduce cardiac BCAA levels and demonstrated positive effects on cardiac function in a preclinical setting. The current study is focused on evaluating the impact of BT2 on the systemic and cardiac levels of BCAA and their metabolites as well as activities of BCAA catabolic enzymes using a quantitative systems pharmacology model. Methods: The model is composed of an ordinary differential equation system characterizing BCAA consumption with food, disposal in the proteins, reversible branched-chain-amino-acid aminotransferase (BCAT)-mediated transamination to branched-chain keto-acids (BCKA), followed by BCKD-mediated oxidation. Activity of BCKD is regulated by the balance of BCKDK and protein phosphatase 2Cm (PP2Cm) activities, affected by BT2 treatment. Cardiac BCAA levels are assumed to directly affect left ventricular ejection fraction (LVEF). Biochemical characteristics of the enzymes are taken from the public domains, while plasma and cardiac BCAA and BCKA levels in BT2 treated mice are used to inform the model parameters. Results: The model provides adequate reproduction of the experimental data and predicts synchronous BCAA responses in the systemic and cardiac space, dictated by rapid BCAA equilibration between the tissues. The model-based simulations indicate maximum possible effect of BT2 treatment on BCAA reduction to be 40% corresponding to 12% increase in LVEF. Model sensitivity analysis demonstrates strong impact of BCKDK and PP2Cm activities as well as total BCKD and co-substrate levels (glutamate, ketoglutarate and ATP) on BCAA and BCKA levels. Conclusion: Model based simulations confirms using of plasma measurements as a marker of cardiac BCAA changes under BCKDK inhibition. The proposed model can be used for optimization of preclinical study design for novel compounds targeting BCAA catabolism.