Zinnia P. Parra-Guillen

Pharmacometrics Markup Language (PharmML): Opening New Perspectives for Model Exchange in Drug Development

The lack of a common exchange format for mathematical models in pharmacometrics has been a long‐standing problem. Such a format has the potential to increase productivity and analysis quality, simplify the handling of complex workflows, ensure reproducibility of research, and facilitate the reuse of existing model resources. Pharmacometrics Markup Language (PharmML), currently under development by the Drug Disease Model Resources (DDMoRe) consortium, is intended to become an exchange standard in pharmacometrics by providing means to encode models, trial designs, and modeling steps.

Mathematical Model Approach to Describe Tumour Response in Mice After Vaccine Administration and its Applicability to Immune-Stimulatory Cytokine-Based Strategies

Immunotherapy is a growing therapeutic strategy in oncology based on the stimulation of innate and adaptive immune systems to induce the death of tumour cells. In this paper, we have developed a population semi-mechanistic model able to characterize the mechanisms implied in tumour growth dynamic after the administration of CyaA-E7, a vaccine able to target antigen to dendritic cells, thus triggering a potent immune response. The mathematical model developed presented the following main components: (1) tumour progression in the animals without treatment was described with a linear model, (2) vaccine effects were modelled assuming that vaccine triggers a non-instantaneous immune response inducing cell death. Delayed response was described with a series of two transit compartments, (3) a resistance effect decreasing vaccine efficiency was also incorporated through a regulator compartment dependent upon tumour size, and (4) a mixture model at the level of the elimination of the induced signal vaccine (k2) to model tumour relapse after treatment, observed in a small percentage of animals (15.6%). The proposed model structure was successfully applied to describe antitumor effect of IL-12, suggesting its applicability to different immune-stimulatory therapies. In addition, a simulation exercise to evaluate in silico the impact on tumour size of possible combination therapies has been shown. This type of mathematical approaches may be helpful to maximize the information obtained from experiments in mice, reducing the number of animals and the cost of developing new antitumor immunotherapies.

Pharmacokinetics and Pharmacokinetic–Pharmacodynamic Relationships of Monoclonal Antibodies in Children

Monoclonal antibodies (mAbs) constitute a therapeutically and economically important drug class with increasing use in both adult and paediatric patients. The rather complex pharmacokinetic and pharmacodynamic properties of mAbs have been extensively reviewed in adults. In children, however, limited information is currently available. This paper aims to comprehensively review published pharmacokinetic and pharmacokinetic–pharmacodynamic studies of mAbs in children. The current status of mAbs in the USA and in Europe is outlined, including a critical discussion of the dosing strategies of approved mAbs. The pharmacokinetic properties of mAbs in children are exhaustively summarised along with comparisons to reports in adults: for each pharmacokinetic process, we discuss the general principles and mechanisms of the pharmacokinetic/pharmacodynamic characteristics of mAbs, as well as key growth and maturational processes in children that might impact these characteristics. Throughout this review, considerable knowledge gaps are identified, especially regarding children-specific properties that influence pharmacokinetics, pharmacodynamics and immunogenicity. Furthermore, the large heterogeneity in the presentation of pharmacokinetic/pharmacodynamic data limited clinical inferences in many aspects of paediatric mAb therapy. Overall, further studies are needed to fully understand the impact of body size and maturational changes on drug exposure and response. To maximise future knowledge gain, we propose a ‘Guideline for Best Practice’ on how to report pharmacokinetic and pharmacokinetic–pharmacodynamic results from mAb studies in children which also facilitates comparisons. Finally, we advocate the use of more sophisticated modelling strategies (population analysis, physiology-based approaches) to appropriately characterise pharmacokinetic–pharmacodynamic relationships of mAbs and, thus, allow for a more rational use of mAb in the paediatric population.

Modeling Tumor Response after Combined Administration of Different Immune-Stimulatory Agents

The aims of this work were as follows: 1) to develop a semimechanistic pharmacodynamic model describing tumor shrinkage after administration of a previously developed antitumor vaccine (CyaA-E7) in combination with CpG (a TLR9 ligand) and/or cyclophosphamide (CTX), and 2) to assess the translational capability of the model to describe tumor effects of different immune-based treatments. Population approach with NONMEM version 7.2 was used to analyze the previously published data. These data were generated by injecting 5 × 105 tumor cells expressing human papillomavirus (HPV)-E7 proteins into C57BL/6 mice. Large and established tumors were treated with CpG and/or CTX administered alone or in combination with CyaA-E7. Applications of the model were assessed by comparing model-based simulations with preclinical and clinical outcomes obtained from literature. CpG effects were modeled: 1) as an amplification of the immune signal triggered by the vaccine and 2) by shortening the delayed response of the vaccine. CTX effects were included through a direct decrease of the tumor-induced inhibition of vaccine efficacy over time, along with a delayed induction of tumor cell death. A pharmacodynamic model, built based on plausible biologic mechanisms known for the coadjuvants, successfully characterized tumor response in all experimental scenarios. The model developed was satisfactory applied to reproduce clinical outcomes when CpG or CTX was used in combination with different vaccines. The results found after simulation exercise indicated that the contribution of the coadjuvants to the tumor response elicited by vaccines can be predicted for other immune-based treatments.

Review on modeling anti-antibody responses to monoclonal antibodies

Monoclonal antibodies (mAbs) represent a therapeutic strategy that has been increasingly used in different diseases. mAbs are highly specific for their targets leading to induce specific effector functions. Despite their therapeutic benefits, the presence of immunogenic reactions is of growing concern. The immunogenicity identified as anti-drug antibodies (ADA) production due to the continuous administration of mAbs may affect the pharmacokinetics (PK) and/or the pharmacodynamics (PD) of mAbs administered to patients. Therefore, the immunogenicity and its clinical impact have been studied by several authors using PK modeling approaches. In this review, the authors try to present all those models under a unique theoretical mechanism-based framework incorporating the main considerations related to ADA formation, and how ADA may affect the efficacy or toxicity profile of some therapeutic biomolecules.

The Fusion Protein of IFN-α and Apolipoprotein A-I Crosses the Blood–Brain Barrier by a Saturable Transport Mechanism

IFN-α is widely used for the treatment of chronic viral hepatitis and malignancies. However, systemic IFN-α treatment causes severe neuropsychiatric complications in humans, including depression, anxiety, and cognitive impairments. We have previously reported that the fusion protein formed by IFN-α and apolipoprotein A-I (IA) circulates bound to high-density lipoproteins (HDLs) and exhibits liver targeting, increased half-life, enhanced immunostimulatory activity, and reduced cytotoxicity. As the transport of HDLs across the blood–brain barrier is a highly complex and regulated process, in this study, we examine the effects of IA on the brain. Determination of IFN-α in brain and serum after hydrodynamic administration of different doses of a plasmid encoding IFN-α or IA showed that IA penetrated into the brain by a saturable transport mechanism. Thus, at high serum levels of the transgenes, the induction of IFN-sensitive genes and the number of phospho-STAT1+ cell nuclei in the brain were substantially higher with IFN-α than with IA. This was associated with attenuation of neurodepression in mice given IA, as manifested by shorter immobility time in the tail suspension test. However, when given low doses of rIFN-α or the same antiviral units of HDLs containing IA, the induction of IFN-stimulated genes in the brain was significantly greater with the latter. In conclusion, IA crosses the blood–brain barrier not by diffusion, as is the case of IFN-α, but by a facilitated saturable transport mechanism. Thus, linkage to apolipoprotein A-I may serve to modulate the effects of IFN-α on the CNS.

Utilising the EGFR interactome to identify mechanisms of drug resistance in non-small cell lung cancer – Proof of concept towards a systems pharmacology approach

Drug treatment of epidermal growth factor receptor (EGFR) positive non-small cell lung cancer has improved substantially by targeting activating mutations within the receptor tyrosine kinase domain. However, the development of drug resistance still limits this approach. As root causes, large heterogeneity between tumour entities but also within tumour cells have been suggested. Therefore, approaches to identify these multitude and complex mechanisms are urgently required. Affinity purification coupled with high resolution mass spectrometry was applied to isolate and characterise the EGFR interactome from HCC4006 non-small cell lung cancer cells and their variant HCC4006rERLO0.5 adapted to grow in the presence of therapeutically relevant concentrations of erlotinib. Bioinformatics analyses were carried out to identify proteins and their related molecular functions that interact differentially with EGFR in the untreated state or when incubated with erlotinib prior to EGFR activation. Across all experimental conditions 375 proteins were detected to participate in the EGFR interactome, 90% of which constituted a complex protein interaction network that was bioinformatically reconstructed from literature data. Treatment of HCC4006rERLO0.5 cells carrying a resistance phenotype to erlotinib was associated with an increase of protein levels of members of the clathrin-associated adaptor protein family AP2 (AP2A1, AP2A2, AP2B1), structural proteins of cytoskeleton rearrangement as well as signalling molecules such as Shc. Validation experiments confirmed activation of the Ras-Raf-Mek-Erk (MAPK)-pathway, of which Shc is an initiating adaptor molecule, in HCC4006rERLO0.5 cells. Taken together, differential proteins in the EGFR interactome of HCC4006rERLO0.5 cells were identified that could be related to multiple resistance mechanisms including alterations in growth factor receptor expression, cellular remodelling processes suggesting epithelial-to-mesenchymal transition as well as alterations in downstream signalling. Knowledge of these mechanisms is a pivotal step to build an integrative model of drug resistance in a systems pharmacology manner and to be able to investigate the interplay of these mechanisms and ultimately recommend combinatorial treatment strategies to overcome drug resistance.

Model Description Language (MDL): A Standard for Modeling and Simulation

Recent work on Model Informed Drug Discovery and Development (MID3) has noted the need for clarity in model description used in quantitative disciplines such as pharmacology and statistics. 1-3 Cur ...

Semi-mechanistic bone marrow exhaustion pharmacokinetic/pharmacodynamic model for chemotherapy-induced cumulative neutropenia

Paclitaxel is a commonly used cytotoxic anticancer drug with potentially life-threatening toxicity at therapeutic doses and high interindividual pharmacokinetic variability. Thus, drug and effect monitoring is indicated to control dose-limiting neutropenia. Joerger et al. (2016) developed a dose individualization algorithm based on a pharmacokinetic (PK)/pharmacodynamic (PD) model describing paclitaxel and neutrophil concentrations. Furthermore, the algorithm was prospectively compared in a clinical trial against standard dosing (Central European Society for Anticancer Drug Research Study of Paclitaxel Therapeutic Drug Monitoring; 365 patients, 720 cycles) but did not substantially improve neutropenia. This might be caused by misspecifications in the PK/PD model underlying the algorithm, especially without consideration of the observed cumulative pattern of neutropenia or the platinum-based combination therapy, both impacting neutropenia. This work aimed to externally evaluate the original PK/PD model for potential misspecifications and to refine the PK/PD model while considering the cumulative neutropenia pattern and the combination therapy. An underprediction was observed for the PK (658 samples), the PK parameters, and these parameters were re-estimated using the original estimates as prior information. Neutrophil concentrations (3274 samples) were overpredicted by the PK/PD model, especially for later treatment cycles when the cumulative pattern aggravated neutropenia. Three different modeling approaches (two from the literature and one newly developed) were investigated. The newly developed model, which implemented the bone marrow hypothesis semiphysiologically, was superior. This model further included an additive effect for toxicity of carboplatin combination therapy. Overall, a physiologically plausible PK/PD model was developed that can be used for dose adaptation simulations and prospective studies to further improve paclitaxel/carboplatin combination therapy.

Role of Cytochrome P450 3A4 and 1A2 Phenotyping in Patients with Advanced Non-Small-Cell Lung Cancer (NSCLC) Receiving Erlotinib Treatment.

Erlotinib is metabolized by cytochrome p450 (CYP) 3A and CYP1A. This study assessed CYP3A4 (midazolam) and CYP1A2 (caffeine) phenotyping in plasma and dried blood spots (DBS) for predicting the pharmacokinetics and toxicity of erlotinib in 36 patients with advanced NSCLC. On day 1, erlotinib 150 mg OD was initiated, and the two oral probe drugs midazolam (2 mg) and caffeine (100 mg) were added on day 1. Plasma and DBS were collected for erlotinib, OSI‐420 and probe drugs for up to 6 hr on day 1 and 2‐weekly up to week 10. Probe drugs, erlotinib and OSI‐420 were analysed using LC‐MS‐MS, and PK data were processed using population modelling. A high correlation was found between plasma and DBS concentrations for erlotinib (R2 = 0.960, p < 0.0001), OSI‐420 (R2 = 0.971, p < 0.0001), midazolam (R2 = 0.995, p < 0.0001) and caffeine (R2 = 0.968, p < 0.0001). Apparent oral caffeine clearance was significantly correlated with erlotinib clearance (R2 = 0.33, p = 0.048), while midazolam clearance was not (R2 = −0.09, p = 0.596). Erlotinib clearance was lower in patients experiencing grade 2 or 3 rash as compared to patients experiencing grade 0 or 1 rash (3.15 versus 3.93 L/hr, p = 0.086 for Students t‐test). The results suggest that probe drug phenotyping is unlikely to substitute therapeutic drug monitoring of erlotinib in patients with advanced NSCLC, but erlotinib PK sampling from DBS may replace more invasive venous sampling and facilitate TDM in patients with cancer.