A General Network Pharmacodynamic Model-Based Design Pipeline for Customized Cancer Therapy Applied to the VEGFR Pathway

Cell state-directed therapy - epigenetic modulation of gene transcription demonstrated with a quantitative systems pharmacology model of temozolomide

Cancer therapy continues to be plagued by modest therapeutic advances. This is particularly evident in glioblastoma multiforme (GBM) wherein treatment failures are attributed to intratumoral heterogeneity (ITH), a dynamic process of cell state transitions or plasticity. To address ITH, we introduce the concept of cell state-directed (CSD) therapy through a quantitative systems pharmacology model of temozolomide (TMZ), a cornerstone of GBM drug therapy. The model consisting of multiple modules incorporated an epigenetic-based gene transcription-translation module that enabled CSD therapy. Numerous model simulations were conducted to demonstrate the potential impact of CSD therapy on TMZ activity. The simulations included those based on global sensitivity analyses to identify fragile nodes - MDM2 and XIAP - in the network, and also how an epigenetic modifier (birabresib) could overcome a mechanism of TMZ resistance. The positive results of CSD therapy on TMZ activity supports continued efforts to develop CSD therapy as a new anticancer approach.

Regulation of CSF and Brain Tissue Sodium Levels by the Blood-CSF and Blood-Brain Barriers During Migraine

Cerebrospinal fluid (CSF) and brain tissue sodium levels increase during migraine. However, little is known regarding the underlying mechanisms of sodium homeostasis disturbance in the brain during the onset and propagation of migraine. Exploring the cause of sodium dysregulation in the brain is important, since correction of the altered sodium homeostasis could potentially treat migraine. Under the hypothesis that disturbances in sodium transport mechanisms at the blood-CSF barrier (BCSFB) and/or the blood-brain barrier (BBB) are the underlying cause of the elevated CSF and brain tissue sodium levels during migraines, we developed a mechanistic, differential equation model of a rat's brain to compare the significance of the BCSFB and the BBB in controlling CSF and brain tissue sodium levels. The model includes the ventricular system, subarachnoid space, brain tissue and blood. Sodium transport from blood to CSF across the BCSFB, and from blood to brain tissue across the BBB were modeled by influx permeability coefficients P_BCSFB and P_BBB, respectively, while sodium movement from CSF into blood across the BCSFB, and from brain tissue to blood across the BBB were modeled by efflux permeability coefficients PBCSFB′ and PBBB′, respectively. We then performed a global sensitivity analysis to investigate the sensitivity of the ventricular CSF, subarachnoid CSF and brain tissue sodium concentrations to pathophysiological variations in P_BCSFB, P_BBB, PBCSFB′ and PBBB′. Our results show that the ventricular CSF sodium concentration is highly influenced by perturbations of P_BCSFB, and to a much lesser extent by perturbations of PBCSFB′. Brain tissue and subarachnoid CSF sodium concentrations are more sensitive to pathophysiological variations of P_BBB and PBBB′ than variations of P_BCSFB and PBCSFB′ within 30 min of the onset of the perturbations. However, P_BCSFB is the most sensitive model parameter, followed by P_BBB and PBBB′, in controlling brain tissue and subarachnoid CSF sodium levels within 3 h of the perturbation onset.

Applications of Quantitative Systems Pharmacology in Model-Informed Drug Discovery: Perspective on Impact and Opportunities

Quantitative systems pharmacology (QSP) approaches have been increasingly applied in the pharmaceutical since the landmark white paper published in 2011 by a National Institutes of Health working group brought attention to the discipline. In this perspective, we discuss QSP in the context of other modeling approaches and highlight the impact of QSP across various stages of drug development and therapeutic areas. We discuss challenges to the field as well as future opportunities.

Mathematical model of hypoxia and tumor signaling interplay reveals the importance of hypoxia and cell-to-cell variability in tumor growth inhibition

Background

Human tumor is a complex tissue with multiple heterogeneous hypoxic regions and significant cell-to-cell variability. Due to the complexity of the disease, the explanation of why anticancer therapies fail cannot be attributed to intrinsic or acquired drug resistance alone. Furthermore, there are inconsistent reports of hypoxia-induced kinase activities in different cancer cell-lines, where increase, decreases, or no change has been observed. Thus, we asked, why are there widely contrasting results in kinase activity under hypoxia in different cancer cell-lines and how does hypoxia play a role in anti-cancer drug sensitivity?

Results

We took a modeling approach to address these questions by analyzing the model simulation to explain why hypoxia driven signals can have dissimilar impact on tumor growth and alter the efficacy of anti-cancer drugs. Repeated simulations with varying concentrations of biomolecules followed by decision tree analysis reveal that the highly differential effects among heterogeneous subpopulation of tumor cells could be governed by varying concentrations of just a few key biomolecules. These biomolecules include activated serine/threonine-specific protein kinases (pRAF), mitogen-activated protein kinase kinase (pMEK), protein kinase B (pAkt), or phosphoinositide-4,5-bisphosphate 3-kinase (pPI3K). Additionally, the ratio of activated extracellular signal-regulated kinases (pERK) or pAkt to its respective total was a key factor in determining the sensitivity of pERK or pAkt to hypoxia.

Conclusion

This work offers a mechanistic insight into how hypoxia can affect the efficacy of anti-cancer drug that targets tumor signaling and provides a framework to identify the types of tumor cells that are either sensitive or resistant to anti-cancer therapy.

Translational Quantitative Systems Pharmacology in Drug Development: from Current Landscape to Good Practices

Systems pharmacology approaches have the capability of quantitatively linking the key biological molecules relevant to a drug candidate's mechanism of action (drug-induced signaling pathways) to the clinical biomarkers associated with the proposed target disease, thereby quantitatively facilitating its development and life cycle management. In this review, the model attributes of published quantitative systems pharmacology (QSP) modeling for lowering cholesterol, treating salt-sensitive hypertension, and treating rare diseases as well as describing bone homeostasis and related pharmacological effects are critically reviewed with respect to model quality, calibration, validation, and performance. We further reviewed the common practices in optimizing QSP modeling and prediction. Notably, leveraging genetics and genomic studies for model calibration and validation is common. Statistical and quantitative assessment of QSP prediction and handling of model uncertainty are, however, mostly lacking as are the quantitative and statistical criteria for assessing QSP predictions and the covariance matrix of coefficients between the parameters in a validated virtual population. To accelerate advances and application of QSP with consistent quality, a list of key questions is proposed to be addressed when assessing the quality of a QSP model in hopes of stimulating the scientific community to set common expectations. The common expectations as to what constitutes the best QSP modeling practices, which the scientific community supports, will advance QSP modeling in the realm of informed drug development. In the long run, good practices will extend the life cycles of QSP models beyond the life cycles of individual drugs.

Pharmacodynamic Models of Differential Bortezomib Signaling Across Several Cell Lines of Multiple Myeloma

The heterogeneous polyclonal nature of multiple myeloma complicates the identification of protein biomarkers predictive of drug response. In this study, a pharmacodynamic systems modeling approach was used to link in vitro bortezomib exposure and myeloma cell death. The exposure‐response was integrated through a network of important protein biomarker dynamics activated by bortezomib in four myeloma cell lines. The pharmacodynamic models reasonably characterized the protein and myeloma cell dynamics simultaneously following bortezomib (20 nM) treatment. The models were used to identify differences in pathway dynamics across cell lines from model‐estimated protein biomarker turnover parameters and global sensitivity analyses. Additionally, a statistical correlation analysis between drug sensitivity and model‐fitted protein activation profiles (i.e., cumulative area under the protein expression‐time curves) supported the identification of shared biomarkers associated with sensitivity differences among the cell lines. Both types of analysis identified similar important proteins associated with bortezomib pharmacodynamics, such as phosphorylated Nuclear Factor kappa‐light‐chain‐enhancer of activated B cells (pNFkappaB), phosphorylated protein kinase B (pAKT), and caspase‐8 (Cas 8).

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.

Multi-Scale Network Model Supported by Proteomics for Analysis of Combined Gemcitabine and Birinapant Effects in Pancreatic Cancer Cells

Gemcitabine combined with birinapant, an inhibitor of apoptosis protein antagonist, acts synergistically to reduce pancreatic cancer cell proliferation. A large‐scale proteomics dataset provided rich time‐series data on proteome‐level changes that reflect the underlying biological system and mechanisms of action of these drugs. A multiscale network model was developed to link the signaling pathways of cell cycle regulation, DNA damage response, DNA repair, apoptosis, nuclear factor‐kappa β (NF‐κβ), and mitogen‐activated protein kinase (MAPK)‐p38 to cell cycle progression, proliferation, and death. After validating the network model under different conditions, the Sobol Sensitivity Analysis was applied to identify promising targets to enhance gemcitabine efficacy. The effects of p53 silencing and combining curcumin with gemcitabine were also tested with the developed model. Merging proteomics analysis with systems modeling facilitates the characterization of quantitative relations among relevant signaling pathways in drug action and resistance, and such multiscale network models could be applied for prediction of combination efficacy and target selection.

Computer Simulation of TSP1 Inhibition of VEGF-Akt-eNOS: An Angiogenesis Triple Threat

The matricellular protein thrombospondin-1 (TSP1) is a potent inhibitor of angiogenesis. Specifically, TSP1 has been experimentally shown to inhibit signaling downstream of vascular endothelial growth factor (VEGF). The molecular mechanism of this inhibition is not entirely clear. We developed a detailed computational model of VEGF signaling to Akt-endothelial nitric oxide synthase (eNOS) to investigate the quantitative molecular mechanism of TSP1 inhibition. The model demonstrated that TSP1 acceleration of VEGFR2 degradation is sufficient to explain the inhibition of VEGFR2 and eNOS phosphorylation. However, Akt inhibition requires TSP1-induced phosphatase recruitment to VEGFR2. The model was then utilized to test various strategies for the rescue of VEGF signaling to Akt and eNOS. Inhibiting TSP1 was predicted to be not as effective as CD47 depletion in rescuing signaling to Akt. The model further predicts that combination strategy involving depletion of CD47 and inhibition of TSP1 binding to CD47 is necessary for effective recovery of signaling to eNOS. In all, computational modeling offers insight to molecular mechanisms involving TSP1 interaction with VEGF signaling and provides strategies for rescuing angiogenesis by targeting TSP1-CD47 axis.

Multiscale systems pharmacological analysis of everolimus action in hepatocellular carcinoma

Dysregulation of mTOR pathway is common in hepatocellular carcinoma (HCC). A translational quantitative systems pharmacology (QSP), pharmacokinetic (PK), and pharmacodynamic (PD) model dissecting the circuitry of this pathway was developed to predict HCC patients' response to everolimus, an mTOR inhibitor. The time course of key signaling proteins in the mTOR pathway, HCC cells viability, tumor volume (TV) and everolimus plasma and tumor concentrations in xenograft mice, clinical PK of everolimus and progression free survival (PFS) in placebo and everolimus-treated patients were extracted from literature. A comprehensive and multiscale QSP/PK/PD model was developed, qualified, and translated to clinical settings. Model fittings and simulations were performed using Monolix software. The S6-kinase protein was identified as critical in the mTOR signaling pathway for describing everolimus lack of efficacy in HCC patients. The net growth rate constant (kg) of HCC cells was estimated at 0.02 h^-1 (2.88%RSE). The partition coefficient of everolimus into the tumor (kp) was determined at 0.06 (12.98%RSE). The kg in patients was calculated from the doubling time of TV in naturally progressing HCC patients, and was determined at 0.004 day^-1. Model-predicted and observed PFS were in good agreement for placebo and everolimus-treated patients. In conclusion, a multiscale QSP/PK/PD model elucidating everolimus lack of efficacy in HCC patients was successfully developed and predicted PFS reasonably well compared to observed clinical findings. This model may provide insights into clinical response to everolimus-based therapy and serve as a valuable tool for the clinical translation of efficacy for novel mTOR inhibitors.

How to mathematically optimize drug regimens using optimal control

This article gives an overview of a technique called optimal control, which is used to optimize real-world quantities represented by mathematical models. I include background information about the historical development of the technique and applications in a variety of fields. The main focus here is the application to diseases and therapies, particularly the optimization of combination therapies, and I highlight several such examples. I also describe the basic theory of optimal control, and illustrate each of the steps with an example that optimizes the doses in a combination regimen for leukemia. References are provided for more complex cases. The article is aimed at modelers working in drug development, who have not used optimal control previously. My goal is to make this technique more accessible in the biopharma community.

Sobol Sensitivity Analysis: A Tool to Guide the Development and Evaluation of Systems Pharmacology Models

A systems pharmacology model typically integrates pharmacokinetic, biochemical network, and systems biology concepts into a unifying approach. It typically consists of a large number of parameters and reaction species that are interlinked based upon the underlying (patho)physiology and the mechanism of drug action. The more complex these models are, the greater the challenge of reliably identifying and estimating respective model parameters. Global sensitivity analysis provides an innovative tool that can meet this challenge. CPT Pharmacometrics Syst. Pharmacol. (2015) 4, 69-79; doi:10.1002/psp4.6; published online 25 February 2015.

Integrating Transcriptomic Data with Mechanistic Systems Pharmacology Models for Virtual Drug Combination Trials

Monotherapy clinical trials with mutation-targeted kinase inhibitors, despite some success in other cancers, have yet to impact glioblastoma (GBM). Besides insufficient blood-brain barrier penetration, combinations are key to overcoming obstacles such as intratumoral heterogeneity, adaptive resistance, and the epistatic nature of tumor genomics that cause mutation-targeted therapies to fail. With now hundreds of potential drugs, exploring the combination space clinically and preclinically is daunting. We are building a simulation-based approach that integrates patient-specific data with a mechanistic computational model of pan-cancer driver pathways (receptor tyrosine kinases, RAS/RAF/ERK, PI3K/AKT/mTOR, cell cycle, apoptosis, and DNA damage) to prioritize drug combinations by their simulated effects on tumor cell proliferation and death. Here we illustrate a first step, tailoring the model to 14 GBM patients from The Cancer Genome Atlas defined by an mRNA-seq transcriptome, and then simulating responses to three promiscuous FDA-approved kinase inhibitors (bosutinib, ibrutinib, and cabozantinib) with evidence for blood-brain barrier penetration. The model captures binding of the drug to primary targets and off-targets based on published affinity data and simulates responses of 100 heterogeneous tumor cells within a patient. Single drugs are marginally effective or even counterproductive. Common copy number alterations (PTEN loss, EGFR amplification, and NF1 loss) have a negligible correlation with single-drug or combination efficacy, reinforcing the importance of postgenetic approaches that account for kinase inhibitor promiscuity to match drugs to patients. Drug combinations tend to be either cytostatic or cytotoxic, but seldom both, highlighting the need for considering targeted and nontargeted therapy. Although we focus on GBM, the approach is generally applicable.

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.

Computational investigation of sphingosine kinase 1 (SphK1) and calcium dependent ERK1/2 activation downstream of VEGFR2 in endothelial cells

Vascular endothelial growth factor (VEGF) is a powerful regulator of neovascularization. VEGF binding to its cognate receptor, VEGFR2, activates a number of signaling pathways including ERK1/2. Activation of ERK1/2 is experimentally shown to involve sphingosine kinase 1 (SphK1) activation and its calcium-dependent translocation downstream of ERK1/2. Here we construct a rule-based computational model of signaling downstream of VEGFR2, by including SphK1 and calcium positive feedback mechanisms, and investigate their consequences on ERK1/2 activation. The model predicts the existence of VEGF threshold in ERK1/2 activation that can be continuously tuned by cellular concentrations of SphK1 and sphingosine 1 phosphate (S1P). The computer model also predicts powerful effects of perturbations in plasma and ER calcium pump rates and the current through the CRAC channels on ERK1/2 activation dynamics, highlighting the critical role of intracellular calcium in shaping the pERK1/2 signal. The model is then utilized to simulate anti-angiogenic therapeutic interventions targeting VEGFR2-ERK1/2 axis. Simulations indicate that monotherapies that exclusively target VEGFR2 phosphorylation, VEGF, or VEGFR2 are ineffective in shutting down signaling to ERK1/2. By simulating therapeutic strategies that target multiple nodes of the pathway such as Raf and SphK1, we conclude that combination therapy should be much more effective in blocking VEGF signaling to EKR1/2. The model has important implications for interventions that target signaling pathways in angiogenesis relevant to cancer, vascular diseases, and wound healing.

Vascular endothelial growth factor (VEGF) signaling is a potent regulator of angiogenesis, the growth and development of new vessels out of a preexisting vascular network. Angiogenesis requires enhanced survival, proliferation, and motility of the vascular endothelial cells. Crucial signaling endpoints in VEGF-mediated angiogenic response include elevation in intracellular calcium and the activation of the proteins ERK1 and 2 (ERK1/2). In this study, we have developed a novel computer model for the activation of ERK1/2 and calcium downstream of VEGF receptor type 2 (VEGFR2). Our model is the first of its kind to incorporate and investigate the consequences of calcium elevation and the role of a cellular lipid modifier known as sphingosine kinase 1 (SphK1). We also utilize the model to simulate therapeutic strategies targeting VEGF signaling to ERK1/2 indicating inefficiency of single therapies known as tyrosine kinase inhibitors (TKI) that target receptor phosphorylation. Computer simulations indicate that combination therapy is essential for effective blockade of this important pathway. Our results have important implications for human diseases such as cancer where plethora of anti-VEGF therapies are currently employed. Overall, our computer model sheds new light on a complex feedback involving SphK1 and calcium that radically alters the response of cells to VEGF.

Precision or Personalized Medicine for Cancer Chemotherapy: Is there a Role for Herbal Medicine

Although over 100 chemotherapeutic agents are currently available for the treatment of cancer patients, the overall long term clinical benefit is disappointing due to the lack of effectiveness or severe side effects from these agents. In order to improve the therapeutic outcome, a new approach called precision medicine or personalized medicine has been proposed and initiated by the U.S. National Institutes of Health. However, the limited availability of effective medications and the high cost are still the major barriers for many cancer patients. Thus alternative approaches such as herbal medicines could be a feasible and less costly option. Unfortunately, scientific evidence for the efficacy of a majority of herbal medicines is still lacking and their development to meet FDA approval or other regulatory agencies is a big challenge. However, herbal medicines may be able to play an important role in precision medicine or personalized medicine. This review will focus on the existing and future technologies that could speed the development of herbal products for treatment of resistant cancer in individual patients. Specifically, it will concentrate on reviewing the phenotypic (activity based) rather than genotypic (mechanism based) approach to develop herbal medicine useful for personalized cancer chemotherapy.

CNS Anticancer Drug Discovery and Development Conference White Paper

Following the first CNS Anticancer Drug Discovery and Development Conference, the speakers from the first 4 sessions and organizers of the conference created this White Paper hoping to stimulate more and better CNS anticancer drug discovery and development. The first part of the White Paper reviews, comments, and, in some cases, expands on the 4 session areas critical to new drug development: pharmacological challenges, recent drug approaches, drug targets and discovery, and clinical paths. Following this concise review of the science and clinical aspects of new CNS anticancer drug discovery and development, we discuss, under the rubric "Accelerating Drug Discovery and Development for Brain Tumors," further reasons why the pharmaceutical industry and academia have failed to develop new anticancer drugs for CNS malignancies and what it will take to change the current status quo and develop the drugs so desperately needed by our patients with malignant CNS tumors. While this White Paper is not a formal roadmap to that end, it should be an educational guide to clinicians and scientists to help move a stagnant field forward.

Optimization of human dose prediction by using quantitative and translational pharmacology in drug discovery

In this perspective article, we explain how quantitative and translational pharmacology, when well-implemented, is believed to lead to improved clinical candidates and drug targets that are differentiated from current treatment options. Quantitative and translational pharmacology aims to build and continuously improve the quantitative relationship between drug exposure, target engagement, efficacy, safety and its interspecies relationship at every phase of drug discovery. Drug hunters should consider and apply these concepts to develop compounds with a higher probability of interrogating the clinical biological hypothesis. We offer different approaches to set an initial effective concentration or pharmacokinetic-pharmacodynamic target in man and to predict human pharmacokinetics that determine together the predicted human dose and dose schedule. All concepts are illustrated with ample literature examples.

Molecular mechanism matters: Benefits of mechanistic computational models for drug development

Making drug development a more efficient and cost-effective process will have a transformative effect on human health. A key, yet underutilized, tool to aid in this transformation is mechanistic computational modeling. By incorporating decades of hard-won prior knowledge of molecular interactions, cellular signaling, and cellular behavior, mechanistic models can achieve a level of predictiveness that is not feasible using solely empirical characterization of drug pharmacodynamics. These models can integrate diverse types of data from cell culture and animal experiments, including high-throughput systems biology experiments, and translate the results into the context of human disease. This provides a framework for identification of new drug targets, measurable biomarkers for drug action in target tissues, and patient populations for which a drug is likely to be effective or ineffective. Additionally, mechanistic models are valuable in virtual screening of new therapeutic strategies, such as gene or cell therapy and tissue regeneration, identifying the key requirements for these approaches to succeed in a heterogeneous patient population. These capabilities, which are distinct from and complementary to those of existing drug development strategies, demonstrate the opportunity to improve success rates in the drug development pipeline through the use of mechanistic computational models.

Systems biology of the microvasculature

The vascular network carries blood throughout the body, delivering oxygen to tissues and providing a pathway for communication between distant organs. The network is hierarchical and structured, but also dynamic, especially at the smaller scales. Remodeling of the microvasculature occurs in response to local changes in oxygen, gene expression, cell-cell communication, and chemical and mechanical stimuli from the microenvironment. These local changes occur as a result of physiological processes such as growth and exercise, as well as acute and chronic diseases including stroke, cancer, and diabetes, and pharmacological intervention. While the vasculature is an important therapeutic target in many diseases, drugs designed to inhibit vascular growth have achieved only limited success, and no drug has yet been approved to promote therapeutic vascular remodeling. This highlights the challenges involved in identifying appropriate therapeutic targets in a system as complex as the vasculature. Systems biology approaches provide a means to bridge current understanding of the vascular system, from detailed signaling dynamics measured in vitro and pre-clinical animal models of vascular disease, to a more complete picture of vascular regulation in vivo. This will translate to an improved ability to identify multi-component biomarkers for diagnosis, prognosis, and monitoring of therapy that are easy to measure in vivo, as well as better drug targets for specific disease states. In this review, we summarize systems biology approaches that have advanced our understanding of vascular function and dysfunction in vivo, with a focus on computational modeling.

Network pharmacodynamic models for customized cancer therapy

Pharmacokinetics (PKs) and pharmacodynamics (PDs) have always been integral to the design of rational drug dosing regimens. Early on PK-driven approaches came under the auspices of therapeutic drug monitoring that progressed into population-based PK and PK/PD modeling analyses. As the availability of tissue samples for measurement of drug concentrations is limited in patients, the bulk of such model-based methods relied on plasma drug concentrations to both build models and monitor therapy. The continued advances in systems biology and the spawning of systems pharmacology propelled the creation of enhanced PD (ePD) models. One of the main characteristic of ePD models is that they are derived from mechanistically grounded biochemical reaction networks. These models are commonly represented as systems of coupled ordinary differential equations with the ability to tailor each reaction and protein concentration to an individual's genomic/proteomic profile. As patient genomic analyses become more common, many genetic and protein abnormalities can be represented in the ePD models, and thus offer a path toward personalized anticancer therapies. By linking PK models to ePD models, a full spectrum of pharmacological simulation tools is available to design sophisticated multidrug regimens. However, ePD models are not a panacea and face challenges in model identifiability, scaling and parameter estimation. Nonetheless, as new technologies evolve and are coupled with fresh ideas on model implementation, it is likely that ePD and PK/ePD models will be considered a viable enterprise to customize anticancer drug therapy.

Systems pharmacology approaches for optimization of antiangiogenic therapies: challenges and opportunities

Targeted therapies have become an important therapeutic paradigm for multiple malignancies. The rapid development of resistance to these therapies impedes the successful management of advanced cancer. Due to the redundancy in angiogenic signaling, alternative proangiogenic factors are activated upon treatment with anti-VEGF agents. Higher doses of the agents lead to greater stimulation of compensatory proangiogenic pathways that limit the therapeutic efficacy of VEGF-targeted drugs and produce escape mechanisms for tumor. Evidence suggests that dose intensity and schedules affect the dynamics of the development of this resistance. Thus, an optimal dosing regimen is crucial to maximizing the therapeutic benefit of antiangiogenic agents and limiting treatment resistance. A systems pharmacology approach using multiscale computational modeling can facilitate a mechanistic understanding of these dynamics of angiogenic biomarkers and their impacts on tumor reduction and resistance. Herein, we discuss a systems pharmacology approach integrating the biology of VEGF-targeted therapy resistance, including circulating biomarkers, and pharmacodynamics to enable the optimization of antiangiogenic therapy for therapeutic gains.

Optimizing oncology therapeutics through quantitative translational and clinical pharmacology: challenges and opportunities

Despite advances in biomedical research that have deepened our understanding of cancer hallmarks, resulting in the discovery and development of targeted therapies, the success rates of oncology drug development remain low. Opportunities remain for objective dose selection informed by exposure-response understanding to optimize the benefit-risk balance of novel therapies for cancer patients. This review article discusses the principles and applications of modeling and simulation approaches across the lifecycle of development of oncology therapeutics. Illustrative examples are used to convey the value gained from integration of quantitative clinical pharmacology strategies from the preclinical-translational phase through confirmatory clinical evaluation of efficacy and safety.

Quantitative insight in utilizing circulating angiogenic factors as biomarkers for antiangiogenic therapy: systems pharmacology approach

Circulating angiogenic factors (CAF) like vascular endothelial growth factor (VEGF), placental growth factor (PlGF), and sVEGFR2 have potential as biomarkers for antiangiogenic therapy. The interpretation of changes in CAF is complicated by the dynamic nature of the tumor and host cells emanating CAF in response to VEGF pathway inhibition. We developed a systems pharmacology model of anti-VEGF agents to investigate CAF modulation by tumor and host cells, and the relationship between overall CAF changes in response to sunitinib and antitumor efficacy. This model distinguishes between the tumor cells' contributions from tumor-independent response to therapy and total plasma CAF correlating with antitumor activity. Altered VEGF is more likely to serve as a useful biomarker reflecting tumor responses in cancer patients whose pretreatment VEGF is higher than baseline VEGF in healthy subjects. Our findings provide a mechanistic insight into tumor modulation of angiogenic molecules, and may explain the inconsistent results found in previous biomarker studies.

Multiscale design of cell-type-specific pharmacokinetic/pharmacodynamic models for personalized medicine: application to temozolomide in brain tumors

Optimizing anticancer therapeutics needs to account for variable drug responses in heterogeneous cell populations within the tumor as well as in organs of toxicity. To address cell heterogeneity, we propose a multiscale modeling approach—from in vitro to preclinical and clinical studies—to develop cell-type–specific pharmacokinetic–pharmacodynamic (PK-PD) models. A physiologically based mechanistic modeling approach integrating data from aqueous solutions, U87 glioma cells, mice, and cancer patients was utilized to characterize the brain disposition of temozolomide (TMZ), the cornerstone of chemotherapy against glioblastoma multiforme. The final model represented intracellular normal brain and brain tumor compartments in which TMZ pH-dependent conversion to the DNA-alkylating species leads to the formation of DNA adducts that serve as an entry point for a PD model. This multiscale protocol can be extended to account for TMZ PK-PD in different cell populations, thus providing a critical tool to personalize TMZ-based chemotherapy on a cell-type–specific basis.