A minimal PBPK model to accelerate preclinical development of drugs against tuberculosis

Introduction: Understanding drug exposure at disease target sites is pivotal to profiling new drug candidates in terms of tolerability and efficacy. Such quantification is particularly tedious for anti-tuberculosis (TB) compounds as the heterogeneous pulmonary microenvironment due to the infection may alter lung permeability and affect drug disposition. Murine models have been a longstanding support in TB research so far and are here used as human surrogates to unveil the distribution of several anti-TB compounds at the site-of-action via a novel and centralized PBPK design framework. Methods: As an intermediate approach between data-driven pharmacokinetic (PK) models and whole-body physiologically based (PB) PK models, we propose a parsimonious framework for PK investigation (minimal PBPK approach) that retains key physiological processes involved in TB disease, while reducing computational costs and prior knowledge requirements. By lumping together pulmonary TB-unessential organs, our minimal PBPK model counts 9 equations compared to the 36 of published full models, accelerating the simulation more than 3-folds in Matlab 2022b. Results: The model has been successfully tested and validated against 11 anti-TB compounds-rifampicin, rifapentine, pyrazinamide, ethambutol, isoniazid, moxifloxacin, delamanid, pretomanid, bedaquiline, OPC-167832, GSK2556286 - showing robust predictability power in recapitulating PK dynamics in mice. Structural inspections on the proposed design have ensured global identifiability and listed free fraction in plasma and blood-to-plasma ratio as top sensitive parameters for PK metrics. The platform-oriented implementation allows fast comparison of the compounds in terms of exposure and target attainment. Discrepancies in plasma and lung levels for the latest BPaMZ and HPMZ regimens have been analyzed in terms of their impact on preclinical experiment design and on PK/PD indices. Conclusion: The framework we developed requires limited drug- and species-specific information to reconstruct accurate PK dynamics, delivering a unified viewpoint on anti-TB drug distribution at the site-of-action and a flexible fit-for-purpose tool to accelerate model-informed drug design pipelines and facilitate translation into the clinic.

Rhizosphere engineering for sustainable crop production: entropy-based insights

There is a growing interest in exploring interactions at root-soil interface in natural and agricultural ecosystems, but an entropy-based understanding of these dynamic rhizosphere processes is lacking. We have developed a new conceptual model of rhizosphere regulation by localized nutrient supply using thermodynamic entropy. Increased nutrient-use efficiency is achieved by rhizosphere management based on self-organization and minimized entropy via equilibrium attractors comprising (i) optimized root strategies for nutrient acquisition and (ii) improved information exchange related to root-soil-microbe interactions. The cascading effects through different hierarchical levels amplify the underlying processes in plant-soil system. We propose a strategy for manipulating rhizosphere dynamics and improving nutrient-use efficiency by localized nutrient supply with minimization of entropy to underpin sustainable food/feed/fiber production.

An age‐dependent mathematical model of neurofilament trafficking in healthy conditions

Neurofilaments (Nfs) are the major structural component of neurons. Their role as a potential biomarker of several neurodegenerative diseases has been investigated in past years with promising results. However, even under physiological conditions, little is known about the leaking of Nfs from the neuronal system and their detection in the cerebrospinal fluid (CSF) and blood. This study aimed at developing a mathematical model of Nf transport in healthy subjects in the 20–90 age range. The model was implemented as a set of ordinary differential equations describing the trafficking of Nfs from the nervous system to the periphery. Model parameters were calibrated on typical Nf levels obtained from the literature. An age‐dependent function modeled on CSF data was also included and validated on data measured in serum. We computed a global sensitivity analysis of model rates and volumes to identify the most sensitive parameters affecting the model’s steady state. Age, Nf synthesis, and degradation rates proved to be relevant for all model variables. Nf levels in the CSF and in blood were observed to be sensitive to the Nf leakage rates from neurons and to the blood clearance rate, and CSF levels were also sensitive to rates representing CSF turnover. An additional parameter perturbation analysis was also performed to investigate possible transient effects on the model variables not captured by the sensitivity analysis. The model provides useful insights into Nf transport and constitutes the basis for implementing quantitative system pharmacology extensions to investigate Nf trafficking in neurodegenerative diseases.

Can Systems Biology Advance Clinical Precision Oncology?

Precision oncology is perceived as a way forward to treat individual cancer patients. However, knowing particular cancer mutations is not enough for optimal therapeutic treatment, because cancer genotype-phenotype relationships are nonlinear and dynamic. Systems biology studies the biological processes at the systems' level, using an array of techniques, ranging from statistical methods to network reconstruction and analysis, to mathematical modeling. Its goal is to reconstruct the complex and often counterintuitive dynamic behavior of biological systems and quantitatively predict their responses to environmental perturbations. In this paper, we review the impact of systems biology on precision oncology. We show examples of how the analysis of signal transduction networks allows to dissect resistance to targeted therapies and inform the choice of combinations of targeted drugs based on tumor molecular alterations. Patient-specific biomarkers based on dynamical models of signaling networks can have a greater prognostic value than conventional biomarkers. These examples support systems biology models as valuable tools to advance clinical and translational oncological research.

A robust computational pipeline for model-based and data-driven phenotype clustering

MOTIVATION

Precision medicine is a promising field that proposes, in contrast to a one-size-fits-all approach, the tailoring of medical decisions, treatments or products. In this context, it is crucial to introduce innovative methods to stratify a population of patients on the basis of an accurate system-level knowledge of the disease. This is particularly important in very challenging conditions, where the use of standard statistical methods can be prevented by poor data availability or by the need of oversimplifying the processes regulating a complex disease.

RESULTS

We define an innovative method for phenotype classification that combines experimental data and a mathematical description of the disease biology. The methodology exploits the mathematical model for inferring additional subject features relevant for the classification. Finally, the algorithm identifies the optimal number of clusters and classifies the samples on the basis of a subset of the features estimated during the model fit. We tested the algorithm in two test cases: an in silico case in the context of dyslipidemia, a complex disease for which a large population of patients has been generated, and a clinical test case, in the context of a lysosomal rare disorder, for which the amount of available data was limited. In both the scenarios, our methodology proved to be accurate and robust, and allowed the inference of an additional phenotype division that the experimental data did not show.

AVAILABILITY AND IMPLEMENTATION

The code to reproduce the in silico results has been implemented in MATLAB v.2017b and it is available in the Supplementary Material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

History and Future Perspectives on the Discipline of Quantitative Systems Pharmacology Modeling and Its Applications

Mathematical biology and pharmacology models have a long and rich history in the fields of medicine and physiology, impacting our understanding of disease mechanisms and the development of novel therapeutics. With an increased focus on the pharmacology application of system models and the advances in data science spanning mechanistic and empirical approaches, there is a significant opportunity and promise to leverage these advancements to enhance the development and application of the systems pharmacology field. In this paper, we will review milestones in the evolution of mathematical biology and pharmacology models, highlight some of the gaps and challenges in developing and applying systems pharmacology models, and provide a vision for an integrated strategy that leverages advances in adjacent fields to overcome these challenges.

A QSP model of prostate cancer immunotherapy to identify effective combination therapies

Immunotherapy, by enhancing the endogenous anti-tumor immune responses, is showing promising results for the treatment of numerous cancers refractory to conventional therapies. However, its effectiveness for advanced castration-resistant prostate cancer remains unsatisfactory and new therapeutic strategies need to be developed. To this end, systems pharmacology modeling provides a quantitative framework to test in silico the efficacy of new treatments and combination therapies. In this paper we present a new Quantitative Systems Pharmacology (QSP) model of prostate cancer immunotherapy, calibrated using data from pre-clinical experiments in prostate cancer mouse models. We developed the model by using Ordinary Differential Equations (ODEs) describing the tumor, key components of the immune system, and seven treatments. Numerous combination therapies were evaluated considering both the degree of tumor inhibition and the predicted synergistic effects, integrated into a decision tree. Our simulations predicted cancer vaccine combined with immune checkpoint blockade as the most effective dual-drug combination immunotherapy for subjects treated with androgen-deprivation therapy that developed resistance. Overall, the model presented here serves as a computational framework to support drug development, by generating hypotheses that can be tested experimentally in pre-clinical models.