Aims Pharmacokinetics of tacrolimus after sublingual administration is not characterized in paediatric liver transplant patients. Therefore, we aimed to develop a population pharmacokinetic model of sublingually administered tacrolimus in patients who cannot swallow the capsules due to their age, sedation status and/or mechanical ventilation during the first weeks post-transplantation. Methods Demographic, clinical and pharmacological variables, including tacrolimus whole blood concentrations obtained from therapeutic drug monitoring and data from dense-sampling pharmacokinetic profiles, were recorded in 26 paediatric patients with biliary atresia who underwent liver transplantation between 2016 and 2021. Population pharmacokinetic analysis was performed with NONMEM v7.4. Results Disposition of tacrolimus was best characterized by a 2-compartment model with clearance achieving half of the maximum elimination capacity (CLMAX = 4.1 L/h) at 4.6 days post-transplantation (T-50). Compared to sedated patients, nonsedated status showed an increased first-order absorption rate constant (1.1 vs. 0.1 h(-1)) and a 24% reduction in bioavailability (F-NS) at 14 days post-transplant. The model was able to explain the oral absorption pattern in nonsedated patients as the result of gut bioavailability (0.9) and hepatic extraction ratio, with the latter being responsible for first-pass effects. Estimates of interindividual variability remained moderate (25.9% for the gut bioavailability) to high (79.8% for the apparent volume of distribution of the central compartment, and 101% for T-50). Conclusion A population pharmacokinetic model of sublingually administered tacrolimus in paediatric patients was developed to characterize different absorption mechanisms. Once the model is externally validated, the effect of post-transplant time on clearance and the sedation status may be considered in routine dosing management.

Background and Purpose Acute intermittent porphyria (AIP) is a rare disease caused by a genetic mutation in the hepatic activity of the porphobilinogen-deaminase. We aimed to develop a mechanistic model of the enzymatic restoration effects of a novel therapy based on the administration of different formulations of recombinant human-PBGD (rhPBGD) linked to the ApoAI lipoprotein. This fusion protein circulates in blood, incorporating into HDL and penetrating hepatocytes. Experimental Approach Single i.v. dose of different formulations of rhPBGD linked to ApoAI were administered to AIP mice in which a porphyric attack was triggered by i.p. phenobarbital. Data consist on 24 h urine excreted amounts of heme precursors, 5-aminolevulinic acid (ALA), PBG and total porphyrins that were analysed using non-linear mixed-effects analysis. Key Results The mechanistic model successfully characterized over time the amounts excreted in urine of the three heme precursors for different formulations of rhPBGD and unravelled several mechanisms in the heme pathway, such as the regulation in ALA synthesis by heme. Treatment with rhPBGD formulations restored PBGD activity, increasing up to 51 times the value of the rate of tPOR formation estimated from baseline. Model-based simulations showed that several formulation prototypes provided efficient protective effects when administered up to 1 week prior to the occurrence of the AIP attack. Conclusion and Implications The model developed had excellent performance over a range of doses and formulation type. This mechanistic model warrants use beyond ApoAI-conjugates and represents a useful tool towards more efficient drug treatments of other enzymopenias as well as for acute intermittent porphyria.

Aims The aims of this work were to build a semi-mechanistic tumour growth inhibition (TGI) model for metastatic colorectal cancer (mCRC) patients receiving either cetuximab + chemotherapy or chemotherapy alone and to identify early predictors of overall survival (OS). Methods A total of 1716 patients from 4 mCRC clinical studies were included in the analysis. The TGI model was built with 8973 tumour size measurements where the probability of drop-out was also included and modelled as a time-to-event variable using parametric survival models, as it was the case in the OS analysis. The effects of patient- and tumour-related covariates on model parameters were explored. Results Chemotherapy and cetuximab effects were included in an additive form in the TGI model. Development of resistance was found to be faster for chemotherapy (drug effect halved at wk 8) compared to cetuximab (drug effect halved at wk 12). KRAS wild-type status and presenting a right-sided primary lesion were related to a 3.5-fold increase in cetuximab drug effect and a 4.7x larger cetuximab resistance, respectively. The early appearance of a new lesion (HR = 4.14), a large tumour size at baseline (HR = 1.62) and tumour heterogeneity (HR = 1.36) were the main predictors of OS. Conclusions Semi-mechanistic TGI and OS models have been developed in a large population of mCRC patients receiving chemotherapy in combination or not with cetuximab. Tumour-related predictors, including a machine learning derived-index of tumour heterogeneity, were linked to changes in drug effect, resistance to treatment or OS, contributing to the understanding of the variability in clinical response.

Hepatitis B liver infection is caused by hepatitis B virus (HBV) and represents a major global disease problem when it becomes chronic, as is the case for 80-90% of vertical or early life infections. However, in the vast majority (>95%) of adult exposures, the infected individuals are capable of mounting an effective immune response leading to infection resolution. A good understanding of HBV dynamics and the interaction between the virus and immune system during acute infection represents an essential step to characterize and understand the key biological processes involved in disease resolution, which may help to identify potential interventions to prevent chronic hepatitis B. In this work, a quantitative systems pharmacology model for acute hepatitis B characterizing viral dynamics and the main components of the innate, adaptive, and tolerant immune response has been successfully developed. To do so, information from multiple sources and across different organization levels has been integrated in a common mechanistic framework. The final model adequately describes the chronology and plausibility of an HBV-triggered immune response, as well as clinical data from acute patients reported in the literature. Given the holistic nature of the framework, the model can be used to illustrate the relevance of the different immune pathways and biological processes to ultimate response, observing the negligible contribution of the innate response and the key contribution of the cellular response on viral clearance. More specifically, moderate reductions of the proliferation of activated cytotoxic CD8+ lymphocytes or increased immunoregulatory effects can drive the system towards chronicity. (c) 2021 The Author(s). Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.

V937 is an investigational novel oncolytic non-genetically modified Kuykendall strain of Coxsackievirus A21 which is in clinical development for the treatment of advanced solid tumor malignancies. V937 infects and lyses tumor cells expressing the intercellular adhesion molecule I (ICAM-I) receptor. We integrated in vitro and in vivo data from six different preclinical studies to build a mechanistic model that allowed a quantitative analysis of the biological processes of V937 viral kinetics and dynamics, viral distribution to tumor, and anti-tumor response elicited by V937 in human xenograft models in immunodeficient mice following intratumoral and intravenous administration. Estimates of viral infection and replication which were calculated from in vitro experiments were successfully used to describe the tumor response in vivo under various experimental conditions. Despite the predicted high clearance rate of V937 in systemic circulation (t(1/2) = 4.3 min), high viral replication was observed in immunodeficient mice which resulted in tumor shrinkage with both intratumoral and intravenous administration. The described framework represents a step towards the quantitative characterization of viral distribution, replication, and oncolytic effect of a novel oncolytic virus following intratumoral and intravenous administrations in the absence of an immune response. This model may further be expanded to integrate the role of the immune system on viral and tumor dynamics to support the clinical development of oncolytic viruses.

Simple Summary: The clinical efficacy of immunotherapies when treating cold tumors is still low, and different treatment combinations are needed when dealing with this challenging scenario. In this work, a middle-out strategy was followed to develop a model describing the antitumor efficacy of different immune-modulator combinations, including an antigen, a toll-like receptor-3 agonist, and an immune checkpoint inhibitor in mice treated with non-inflamed tumor cells. Our results support that clinical response requires antigen-presenting cell activation and also relies on the amount of CD8 T cells and tumor resistance mechanisms present. This mathematical model is a very useful platform to evaluate different immuno-oncology combinations in both preclinical and clinical settings. Immune checkpoint inhibitors, administered as single agents, have demonstrated clinical efficacy. However, when treating cold tumors, different combination strategies are needed. This work aims to develop a semi-mechanistic model describing the antitumor efficacy of immunotherapy combinations in cold tumors. Tumor size of mice treated with TC-1/A9 non-inflamed tumors and the drug effects of an antigen, a toll-like receptor-3 agonist (PIC), and an immune checkpoint inhibitor (anti-programmed cell death 1 antibody) were modeled using Monolix and following a middle-out strategy. Tumor growth was best characterized by an exponential model with an estimated initial tumor size of 19.5 mm(3) and a doubling time of 3.6 days. In the treatment groups, contrary to the lack of response observed in monotherapy, combinations including the antigen were able to induce an antitumor response. The final model successfully captured the 23% increase in the probability of cure from bi-therapy to triple-therapy. Moreover, our work supports that CD8(+) T lymphocytes and resistance mechanisms are strongly related to the clinical outcome. The activation of antigen-presenting cells might be needed to achieve an antitumor response in reduced immunogenic tumors when combined with other immunotherapies. These models can be used as a platform to evaluate different immuno-oncology combinations in preclinical and clinical scenarios.

Crohn's disease (CD) is a complex inflammatory bowel disease whose pathogenesis appears to involve several immunologic defects causing functional impairment of the gut. Its complexity and the reported loss of effectiveness over time of standard of care together with the increase in its worldwide incidence require the application of techniques aiming to find new therapeutic strategies. Currently, systems pharmacology modeling has been gaining importance as it integrates the available knowledge of the system into a single computational model. In this work, the following workflow for robust application of systems pharmacology modeling was followed: 1) scope definition; 2) species selection and circulating plasma levels based on a search in the literature; 3) representation of model topology and parametrization of the interactions, after literature data extraction and curation, and the implementation of ordinary differential equations in SimBiology (MATLAB version R2018b); and 4) model curation and evaluation by visual comparison of simulated interleukin (IL) concentrations with the reported levels in plasma, and sensitivity analysis performed to confirm model robustness and identify the most influential parameters. Finally, 5) exposure to two dose levels of recombinant human IL10 was evaluated by simulation and comparison with reported clinical study results. In summary, we present a quantitative systems pharmacology model for the main ILs involved in CD developed using a standardized methodology and supported by a comprehensive repository summarizing the most relevant literature in the field. However, it has to be taken into account that external validation is still pending as available clinical data were primarily used for model training. SIGNIFICANCE STATEMENT Crohn's disease (CD) is a complex heterogeneous inflammatory bowel disorder. Systems pharmacology modeling offers a great opportunity for integration of the available knowledge on the disease using a computational framework. As a result of this work, a comprehensive repository along with a quantitative systems pharmacology model for the main interleukins involved in CD is provided. This model is useful for the in silico evaluation of biomarkers and potential therapeutic targets and can be adapted to address research gaps regarding CD.

Background and Purpose Acute intermittent porphyria (AIP) results from haplo-insufficiency of the porphobilinogen deaminase (PBGD) gene encoding the third enzyme in the haem biosynthesis pathway. As liver is the main organ of pathology for AIP, emerging therapies that restore enzyme hepatic levels are appealing. The objective of this work was to develop a mechanistic-based computational framework to describe the effects of novel PBGD mRNA therapy on the accumulation of neurotoxic haem precursors in small and large animal models. Experimental Approach Liver PBGD activity data and/or 24-hr urinary haem precursors were obtained from genetic AIP mice and wild-type mice, rats, rabbits, and macaques. To mimic acute attacks, porphyrogenic drugs were administered over one or multiple challenges, and animals were used as controls or treated with different PBGD mRNA products. Available experimental data were sequentially used to build and validate a semi-mechanistic mathematical model using non-linear mixed-effects approach. Key Results The developed framework accounts for the different biological processes involved (i.e., mRNA sequence, release from lipid nanoparticle and degradation, mRNA translation, increased PBGD activity in liver, and haem precursor metabolism) in a simplified mechanistic fashion. The model, validated using external data, shows robustness in the extrapolation of PBGD activity data in rat, rabbit, and non-human primate species. Conclusion and Implications This quantitative framework provides a valuable tool to compare PBGD mRNA drug products during early preclinical stages, optimize the amount of experimental data required, and project results to humans, thus supporting drug development and clinical dose and dosing regimen selection.

Total tumor size (TS) metrics used in TS models in oncology do not consider tumor heterogeneity, which could help to better predict drug efficacy. We analyzed individual target lesions (iTLs) of patients with metastatic colorectal carcinoma (mCRC) to determine differences in TS dynamics by using the ClassIfication Clustering of Individual Lesions (CICIL) methodology. Results from subgroup analyses comparing genetic mutations and TS metrics were assessed and applied to survival analyses. Data from four mCRC clinical studies were analyzed (1781 patients, 6369 iTLs). CICIL was used to assess differences in lesion TS dynamics within a tissue (intra-class) or across different tissues (inter-class). First, lesions were automatically classified based on their location. Cross-correlation coefficients (CCs) determined if each pair of lesions followed similar or opposite dynamics. Finally, CCs were grouped by using the K-means clustering method. Heterogeneity in tumor dynamics was lower in the intra-class analysis than in the inter-class analysis for patients receiving cetuximab. More tumor heterogeneity was found in KRAS mutated patients compared to KRAS wild-type (KRASwt) patients and when using sum of longest diameters versus sum of products of diameters. Tumor heterogeneity quantified as the median patient's CC was found to be a predictor of overall survival (OS) (HR = 1.44, 95% CI 1.08-1.92), especially in KRASwt patients. Intra- and inter-tumor tissue heterogeneities were assessed with CICIL. Derived metrics of heterogeneity were found to be a predictor of OS time. Considering differences between lesions' TS dynamics could improve oncology models in favor of a better prediction of OS.

Over the past decade, the insulin-like growth factor (IGF)-signaling pathway has gained substantial interest as potential therapeutic target in oncology. Xentuzumab, a humanized IgG1 monoclonal antibody, binds to IGF-I and IGF-II thereby inhibiting the downstream signaling essential for survival and tumor growth. This pathway is further regulated by circulating IGF binding proteins (IGFBPs). In this work, a mechanistic model characterizing the dynamics and interactions of IGFs, IGFBPs, and Xentuzumab has been developed to guide dose selection. Therefore, in vitro and in vivo literature information was combined with temporal IGF-I, IGF-II, and IGFBP-3 total plasma concentrations from two phase I studies. Based on the established quantitative framework, the time-course of free IGFs as ultimate drug targets not measured in clinics was predicted. Finally, a dose of 1000 mg/week-predicted to reduce free IGF-I and free IGF-II at steady-state by at least 90% and 64%, respectively-was suggested for phase II.