PKPD modeling of acquired resistance to anti-cancer drug treatment

Dosage optimization for reducing tumor burden using a phenotype-structured population model with a drug-resistance continuum

Drug resistance is a significant obstacle to effective cancer treatment. To gain insights into how drug resistance develops, we adopted a concept called fitness landscape and employed a phenotype-structured population model by fitting to a set of experimental data on a drug used for ovarian cancer, olaparib. Our modeling approach allowed us to understand how a drug affects the fitness landscape and track the evolution of a population of cancer cells structured with a spectrum of drug resistance. We also incorporated pharmacokinetic (PK) modeling to identify the optimal dosages of the drug that could lead to long-term tumor reduction. We derived a formula that indicates that maximizing variation in plasma drug concentration over a dosing interval could be important in reducing drug resistance. Our findings suggest that it may be possible to achieve better treatment outcomes with a drug dose lower than the levels recommended by the drug label. Acknowledging the current limitations of our work, we believe that our approach, which combines modeling of both PK and drug resistance evolution, could contribute to a new direction for better designing drug treatment regimens to improve cancer treatment.

Comparison of Tumor Binding Across Tumor Types and Cell Lines to Support Free Drug Considerations for Oncology Drug Discovery

Tumor binding is an important parameter to derive unbound tumor concentration to explore pharmacokinetics (PK) and pharmacodynamics (PD) relationships for oncology disease targets. Tumor binding was evaluated using eleven matrices, including various commonly used ex vivo human and mouse xenograft and syngeneic tumors, tumor cell lines and liver as a surrogate tissue. The results showed that tumor binding is highly correlated among the different tumors and tumor cell lines except for the mouse melanoma (B16F10) tumor type. Liver fraction unbound (fu) has a good correlation with B16F10 tumor binding. Liver also demonstrates a two-fold equivalency, on average, with binding of other tumor types when a scaling factor is applied. Predictive models were developed for tumor binding, with correlations established with LogD (acids), predicted muscle fu (neutrals) and measured plasma protein binding (bases) to estimate tumor fu when experimental data are not available. Many approaches can be applied to obtain and estimate tumor binding values. One strategy proposed is to use a surrogate tumor tissue, such as mouse xenograft ovarian cancer (OVCAR3) tumor, as a surrogate for tumor binding (except for B16F10) to provide an early assessment of unbound tumor concentrations for development of PK/PD relationships.

Investigating brain uptake of a non-targeting monoclonal antibody after intravenous and intracerebroventricular administration

Monoclonal antibodies play an important role in the treatment of various diseases. However, the development of these drugs against neurological disorders where the drug target is located in the brain is challenging and requires a good understanding of the local drug concentration in the brain. In this original research, we investigated the systemic and local pharmacokinetics in the brain of healthy rats after either intravenous (IV) or intracerebroventricular (ICV) administration of EGFRvIII-T-Cell bispecific (TCB), a bispecific monoclonal antibody. We established an experimental protocol that allows serial sampling in serum, cerebrospinal fluid (CSF) and interstitial fluid (ISF) of the prefrontal cortex in freely moving rats. For detection of drug concentration in ISF, a push-pull microdialysis technique with large pore membranes was applied. Brain uptake into CSF and ISF was characterized and quantified with a reduced brain physiologically-based pharmacokinetic model. The model allowed us to interpret the pharmacokinetic processes of brain uptake after different routes of administration. The proposed model capturing the pharmacokinetics in serum, CSF and ISF of the prefrontal cortex suggests a barrier function between the CSF and ISF that impedes free antibody transfer. This finding suggests that ICV administration may not be better suited to reach higher local drug exposure as compared to IV administration. The model enabled us to quantify the relative contribution of the blood-brain barrier (BBB) and Blood-CSF-Barrier to the uptake into the interstitial fluid of the brain. In addition, we compared the brain uptake of three monoclonal antibodies after IV dosing. In summary, the presented approach can be applied to profile compounds based on their relative uptake in the brain and provides quantitative insights into which pathways are contributing to the net exposure in the brain.

Characterization of the Safety Profile of Sweet Chestnut Wood Distillate Employed in Agriculture

In organic agriculture, synthetic pesticides and treatments are substituted by natural remedies with interesting success for product yield and environmental outcomes, but the safety of these bio-based products needs to be assessed in vertebrate and human models. Therefore, in this paper we assessed the safety profile of sweet chestnut (Castanea sativa) wood distillate (WD) on the different cellular components of tissues implied in transcutaneous absorption. We investigated the viability of different cell lines mimicking the skin (HaCaT keratinocytes), mucosa (A431), connective (normal human dermal fibroblasts, NHDF) and vascular (human umbilical vein endothelial cells, HUVEC) tissues after exposure to increasing concentrations (0.04–0.5%, v/v, corresponding to 1:2800–1:200 dilutions) of WD. A short exposure to increasing doses of WD was well tolerated up to the highest concentration. Instead, following a prolonged treatment, a concentration dependent cytotoxic effect was observed. Notably, a different behavior was found with the various cell lines, with higher sensitivity to cytotoxicity by the cells with higher proliferation rate and reduced doubling time (human keratinocytes). Moreover, to exclude an inflammatory effect at the not cytotoxic WD concentrations, the expression of the main inducible markers of inflammation, cyclooxygenase-2 (COX-2) and microsomal prostaglandin E synthase-1 (mPGES-1), were assessed, and no improvement was found both after brief and prolonged exposure. In conclusion, our data exclude any inflammatory and cytotoxic effect at the lowest WD concentrations, namely 0.07% and 0.04%, mimicking some recommended dilutions of the product and the potential exposure doses for the operators in agriculture. Nevertheless, higher concentrations showed a safe profile for short time usage, but caution should be used by farmers following persistent product exposure.

Optimum multi-drug regime for compartment model of tumour: cell-cycle-specific dynamics in the presence of resistance

This work is focused on multi-objective optimisation of a multi-drug chemotherapy schedule for cell-cycle-specific cancer treatment under the influence of drug resistance. The acquired drug resistance to chemotherapeutic agents is incorporated into the existing compartmental model of breast cancer. Furthermore, the toxic effect of drugs on healthy cells and overall drug concentration in the patient body are also constrained in the proposed model. The objective is to determine the optimal drug schedule according to the patient's physiological condition so that the tumour burden is minimised. A multi-objective optimisation algorithm, non-dominated sorting genetic algorithm-II (NSGA-II) is utilised to solve the problem. The obtained results are thoroughly analysed to illustrate the impact of drug resistance on the treatment. The capability of optimised schedules to deal with parametric uncertainty is also analysed. The drug schedules obtained in this work align well with the clinical standards. It is also revealed that the NSGA-II optimised drug schedule with proper rest period between successive dosages yields the minimum cancer load at the end of the treatment.

Acquired Resistance in Lung Cancer

The last decade has witnessed a transformation in the treatment of advanced-stage lung cancer from a largely palliative approach to one where long-term durable remissions and even cures might be within reach. In this review, we discuss the current state of oncogene-directed precision medicine therapies in lung cancer and focus on the major cause of mortality for lung cancer patients: acquired resistance. We consider the multifaceted resistance mechanisms tumors utilize, often simultaneously. We then present areas for future scientific and clinical investigation with an emphasis on population dynamics, early detection, combinatorial therapies targeting resistance mechanisms, and understanding the drug-tolerant persister state.

Optimization of Cancer Treatment in the Frequency Domain

Thorough exploration of alternative dosing frequencies is often not performed in conventional pharmacometrics approaches. Quantitative systems pharmacology (QSP) can provide novel insights into optimal dosing regimen and drug behaviors which could add a new dimension to the design of novel treatments. However, methods for such an approach are currently lacking. Recently, we illustrated the utility of frequency-domain response analysis (FdRA), an analytical method used in control engineering, using several generic pharmacokinetic-pharmacodynamic case studies. While FdRA is not applicable to models harboring ever increasing variables such as those describing tumor growth, studying such models in the frequency domain provides valuable insight into optimal dosing frequencies. Through the analysis of three distinct tumor growth models (cell cycle-specific, metronomic, and acquired resistance), we demonstrate the application of a simulation-based analysis in the frequency domain to optimize cancer treatments. We study the response of tumor growth to dosing frequencies while simultaneously examining treatment safety, and found for all three models that above a certain dosing frequency, tumor size is insensitive to an increase in dosing frequency, e.g., for the cell cycle-specific model, one dose per 3 days, and an hourly dose yield the same reduction of tumor size to 3% of the initial size after 1 year of treatment. Additionally, we explore the effect of drug elimination rate changes on the tumor growth response. In summary, we show that the frequency-domain view of three models of tumor growth dynamics can help in optimizing drug dosing regimen to improve treatment success.

Systematic Review of Patient-Derived Xenograft Models for Preclinical Studies of Anti-Cancer Drugs in Solid Tumors

Patient-derived xenograft (PDX) models are used as powerful tools for understanding cancer biology in PDX clinical trials and co-clinical trials. In this systematic review, we focus on PDX clinical trials or co-clinical trials for drug development in solid tumors and summarize the utility of PDX models in the development of anti-cancer drugs, as well as the challenges involved in this approach, following the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines. Recently, the assessment of drug efficacy by PDX clinical and co-clinical trials has become an important method. PDX clinical trials can be used for the development of anti-cancer drugs before clinical trials, with their efficacy assessed by the modified response evaluation criteria in solid tumors (mRECIST). A few dozen cases of PDX models have completed enrollment, and the efficacy of the drugs is assessed by 1 × 1 × 1 or 3 × 1 × 1 approaches in the PDX clinical trials. Furthermore, co-clinical trials can be used for personalized care or precision medicine with the evaluation of a new drug or a novel combination. Several PDX models from patients in clinical trials have been used to assess the efficacy of individual drugs or drug combinations in co-clinical trials.

Pharmacodynamic modelling of resistance to epidermal growth factor receptor inhibition in brain metastasis mouse models

Purpose

Epidermal growth factor receptor (EGFR) is thought to play a role in the regulation of cell proliferation; with its activation stimulating tumour growth. EGFR inhibitors have shown promise in the treatment of cancer, particularly in non-small cell lung cancer, however, resistance is observed in the majority of patients. A tumour growth model was developed aiming to explain this resistance.

Methods

The model incorporating populations of both sensitive and resistant cells were fitted to data from a study of EGFR inhibitor AZD3759 in brain metastasis mouse models. The observed regrowth of tumours in higher dose groups suggested the development of resistance to treatment. The bioluminescence observations were highly variable, covering many orders of magnitude, so to assess how reliable the model was, the parameter estimates were compared to those found in less noisy subcutaneous mouse models.

Results

The fitted model suggested that resistance was mainly due to a proportion of cells being resistant at baseline, and the contribution of mutations occurring during the study leading to resistance was negligible. Estimated growth rate and dose–response was found to be comparable between brain metastasis and subcutaneous mouse models.

Conclusions

The developed model to describe resistance suggests that the resistance to EGFR-inhibition seen in these xenografts is best described by assuming a small percentage of cells are resistant to treatment at baseline. This model suggests changes to dosing and dosing schedule may not prevent resistance to treatment developing, and that additional treatments would need to be used in combination to overcome resistance.