Requirements for multi-level systems pharmacology models to reach end-usage: the case of type 2 diabetes

Applying quantitative and systems pharmacology to drug development and beyond: An introduction to clinical pharmacologists

ABSTRACT

Quantitative and systems pharmacology (QSP) is an innovative and integrative approach combining physiology and pharmacology to accelerate medical research. This review focuses on QSP's pivotal role in drug development and its broader applications, introducing clinical pharmacologists/researchers to QSP's quantitative approach and the potential to enhance their practice and decision-making. The history of QSP adoption reveals its impact in diverse areas, including glucose regulation, oncology, autoimmune disease, and HIV treatment. By considering receptor-ligand interactions of various cell types, metabolic pathways, signaling networks, and disease biomarkers simultaneously, QSP provides a holistic understanding of interactions between the human body, diseases, and drugs. Integrating knowledge across multiple time and space scales enhances versatility, enabling insights into personalized responses and general trends. QSP consolidates vast data into robust mathematical models, predicting clinical trial outcomes and optimizing dosing based on preclinical data. QSP operates under a "learn and confirm paradigm," integrating experimental findings to generate testable hypotheses and refine them through precise experimental designs. An interdisciplinary collaboration involving expertise in pharmacology, biochemistry, genetics, mathematics, and medicine is vital. QSP's utility in drug development is demonstrated through integration in various stages, predicting drug responses, optimizing dosing, and evaluating combination therapies. Challenges exist in model complexity, communication, and peer review. Standardized workflows and evaluation methods ensure reliability and transparency.

Practical QSP application from the preclinical phase to enhance the probability of clinical success: Insights from case studies in oncology

Quantitative Systems Pharmacology (QSP) has emerged as a promising modeling and simulation (M&S) approach in drug development, with potential to improve clinical success rates. While conventional M&S has significantly contributed to quantitative understanding in late preclinical and clinical phases, it falls short in explaining unexpected phenomena and testing hypotheses in the early research phase. QSP presents a solution to these limitations. To harness the full potential of QSP in early preclinical stages, preclinical modelers who are familiar with conventional M&S need to update their understanding of the differences between conventional M&S and QSP. This review focuses on QSP applications during the preclinical stage, citing case examples and sharing our experiences in oncology. We emphasize the critical role of QSP in increasing the probability of success for clinical proof of concept (PoC) when applied from the early preclinical stage. Enhancing the quality of both hypotheses and QSP models from early preclinical stage is of critical importance. Once a QSP model achieves credibility, it facilitates predictions of clinical responses and potential biomarkers. We propose that sequential QSP applications from preclinical stages can improve success rates of clinical PoC, and emphasize the importance of refining both hypotheses and QSP models throughout the process.

Integrated experimental-computational analysis of a HepaRG liver-islet microphysiological system for human-centric diabetes research

Microphysiological systems (MPS) are powerful tools for emulating human physiology and replicating disease progression in vitro. MPS could be better predictors of human outcome than current animal models, but mechanistic interpretation and in vivo extrapolation of the experimental results remain significant challenges. Here, we address these challenges using an integrated experimental-computational approach. This approach allows for in silico representation and predictions of glucose metabolism in a previously reported MPS with two organ compartments (liver and pancreas) connected in a closed loop with circulating medium. We developed a computational model describing glucose metabolism over 15 days of culture in the MPS. The model was calibrated on an experiment-specific basis using data from seven experiments, where HepaRG single-liver or liver-islet cultures were exposed to both normal and hyperglycemic conditions resembling high blood glucose levels in diabetes. The calibrated models reproduced the fast (i.e. hourly) variations in glucose and insulin observed in the MPS experiments, as well as the long-term (i.e. over weeks) decline in both glucose tolerance and insulin secretion. We also investigated the behaviour of the system under hypoglycemia by simulating this condition in silico, and the model could correctly predict the glucose and insulin responses measured in new MPS experiments. Last, we used the computational model to translate the experimental results to humans, showing good agreement with published data of the glucose response to a meal in healthy subjects. The integrated experimental-computational framework opens new avenues for future investigations toward disease mechanisms and the development of new therapies for metabolic disorders.

Microphysiological systems (MPS) are powerful tools to unravel biological knowledge underlying disease. MPS provide a physiologically relevant, human-based in vitro setting, which can potentially yield better translatability to humans than current animal models and traditional cell cultures. However, mechanistic interpretation and extrapolation of the experimental results to human outcome remain significant challenges. In this study, we confront these challenges using an integrated experimental-computational approach. We present a computational model describing glucose metabolism in a previously reported MPS integrating liver and pancreas. This MPS supports a homeostatic feedback loop between HepaRG/HHSteC spheroids and pancreatic islets, and allows for detailed investigations of mechanisms underlying type 2 diabetes in humans. We show that the computational model captures the complex dynamics of glucose-insulin regulation observed in the system, and can provide mechanistic insight into disease progression features, such as insulin resistance and β-cell dynamics. Furthermore, the computational model can explain key differences in temporal dynamics between MPS and human responses, and thus provides a tool for translating experimental insights into human outcome. The integrated experimental-computational framework opens new avenues for future investigations toward disease mechanisms and the development of new therapies for metabolic disorders.

Digital twin predicting diet response before and after long-term fasting

Today, there is great interest in diets proposing new combinations of macronutrient compositions and fasting schedules. Unfortunately, there is little consensus regarding the impact of these different diets, since available studies measure different sets of variables in different populations, thus only providing partial, non-connected insights. We lack an approach for integrating all such partial insights into a useful and interconnected big picture. Herein, we present such an integrating tool. The tool uses a novel mathematical model that describes mechanisms regulating diet response and fasting metabolic fluxes, both for organ-organ crosstalk, and inside the liver. The tool can mechanistically explain and integrate data from several clinical studies, and correctly predict new independent data, including data from a new study. Using this model, we can predict non-measured variables, e.g. hepatic glycogen and gluconeogenesis, in response to fasting and different diets. Furthermore, we exemplify how such metabolic responses can be successfully adapted to a specific individual’s sex, weight, height, as well as to the individual’s historical data on metabolite dynamics. This tool enables an offline digital twin technology.

Fasting and diet are central components of prevention against cardiovascular disease. Unfortunately, there is little consensus regarding which diet schemes are optimal. This is partially because different clinical studies contribute with different non-connected pieces of knowledge, which have not been fully integrated into a useful and interconnected big picture. In principle, mathematical models describing meal responses could be used for such an integration. However, today’s models still lack critical mechanisms, such as protein metabolism and a dynamic glycogen regulation. Herein, we present a) a new expanded model structure including these mechanisms; b) a set of parameters which can simultaneously describe a wide array of complementary estimation data, in both healthy and diabetic populations; c) a personalisation-script, which allows these generic parameters to be tuned to an individual/sub-population, using demographics (age, weight, height, diabetes status) and historic metabolic data. We exemplify how this personalisation can be used to predict new independent data, including a new clinical study, where a qualitatively new prediction is validated: that an oral protein tolerance test gives a clear response in plasma glucose, after, but not before, a 48h fasting period. Our combined model, parameters, and fitting script lay the foundation for an offline digital twin.

A systems biology analysis of lipolysis and fatty acid release from adipocytes in vitro and from adipose tissue in vivo

Lipolysis and the release of fatty acids to supply energy fuel to other organs, such as between meals, during exercise, and starvation, are fundamental functions of the adipose tissue. The intracellular lipolytic pathway in adipocytes is activated by adrenaline and noradrenaline, and inhibited by insulin. Circulating fatty acids are elevated in type 2 diabetic individuals. The mechanisms behind this elevation are not fully known, and to increase the knowledge a link between the systemic circulation and intracellular lipolysis is key. However, data on lipolysis and knowledge from in vitro systems have not been linked to corresponding in vivo data and knowledge in vivo. Here, we use mathematical modelling to provide such a link. We examine mechanisms of insulin action by combining in vivo and in vitro data into an integrated mathematical model that can explain all data. Furthermore, the model can describe independent data not used for training the model. We show the usefulness of the model by simulating new and more challenging experimental setups in silico, e.g. the extracellular concentration of fatty acids during an insulin clamp, and the difference in such simulations between individuals with and without type 2 diabetes. Our work provides a new platform for model-based analysis of adipose tissue lipolysis, under both non-diabetic and type 2 diabetic conditions.

A systems biology analysis of adrenergically stimulated adiponectin exocytosis in white adipocytes

Circulating levels of the adipocyte hormone adiponectin are typically reduced in obesity, and this deficiency has been linked to metabolic diseases. It is thus important to understand the mechanisms controlling adiponectin exocytosis. This understanding is hindered by the high complexity of both the available data and the underlying signaling network. To deal with this complexity, we have previously investigated how different intracellular concentrations of Ca²⁺, cAMP, and ATP affect adiponectin exocytosis, using both patch-clamp recordings and systems biology mathematical modeling. Recent work has shown that adiponectin exocytosis is physiologically triggered via signaling pathways involving adrenergic β₃ receptors (β₃ARs). Therefore, we developed a mathematical model that also includes adiponectin exocytosis stimulated by extracellular epinephrine or the β₃AR agonist CL 316243. Our new model is consistent with all previous patch-clamp data as well as new data (collected from stimulations with a combination of the intracellular mediators and extracellular adrenergic stimuli) and can predict independent validation data. We used this model to perform new in silico experiments where corresponding wet lab experiments would be difficult to perform. We simulated adiponectin exocytosis in single cells in response to the reduction of β₃ARs that is observed in adipocytes from animals with obesity-induced diabetes. Finally, we used our model to investigate intracellular dynamics and to predict both cAMP levels and adiponectin release by scaling the model from single-cell to a population of cells-predictions corroborated by experimental data. Our work brings us one step closer to understanding the intricate regulation of adiponectin exocytosis.

An Updated Organ-Based Multi-Level Model for Glucose Homeostasis: Organ Distributions, Timing, and Impact of Blood Flow

Glucose homeostasis is the tight control of glucose in the blood. This complex control is important, due to its malfunction in serious diseases like diabetes, and not yet sufficiently understood. Due to the involvement of numerous organs and sub-systems, each with their own intra-cellular control, we have developed a multi-level mathematical model, for glucose homeostasis, which integrates a variety of data. Over the last 10 years, this model has been used to insert new insights from the intra-cellular level into the larger whole-body perspective. However, the original cell-organ-body translation has during these years never been updated, despite several critical shortcomings, which also have not been resolved by other modeling efforts. For this reason, we here present an updated multi-level model. This model provides a more accurate sub-division of how much glucose is being taken up by the different organs. Unlike the original model, we now also account for the different dynamics seen in the different organs. The new model also incorporates the central impact of blood flow on insulin-stimulated glucose uptake. Each new improvement is clear upon visual inspection, and they are also supported by statistical tests. The final multi-level model describes >300 data points in >40 time-series and dose-response curves, resulting from a large variety of perturbations, describing both intra-cellular processes, organ fluxes, and whole-body meal responses. We hope that this model will serve as an improved basis for future data integration, useful for research and drug developments within diabetes.

Type 2 diabetes progression: A regulatory network approach

In order to elucidate central elements underlying type 2 diabetes, we constructed a regulatory network model involving 37 components (molecules, receptors, processes, etc.) associated to signaling pathways of pancreatic beta-cells. In a first approximation, the network topology was described by Boolean rules whose interacting dynamics predicted stationary patterns broadly classified as health, metabolic syndrome, and diabetes stages. A subsequent approximation based on a continuous logic analysis allowed us to characterize the progression of the disease as transitions between these states associated to alterations of cell homeostasis due to exhaustion or exacerbation of specific regulatory signals. The method allowed the identification of key transcription factors involved in metabolic stress as essential for the progression of the disease. Integration of the present analysis with existent mathematical models designed to yield accurate account of experimental data in human or animal essays leads to reliable predictions for beta-cell mass, insulinemia, glycemia, and glycosylated hemoglobin in diabetic fatty rats.

Type 2 diabetes: one disease, many pathways

Diabetes is a chronic, progressive disease that calls for longitudinal data and analysis. We introduce a longitudinal mathematical model that is capable of representing the metabolic state of an individual at any point in time during their progression from normal glucose tolerance to type 2 diabetes (T2D) over a period of years. As an application of the model, we account for the diversity of pathways typically followed, focusing on two extreme alternatives, one that goes through impaired fasting glucose (IFG) first and one that goes through impaired glucose tolerance (IGT) first. These two pathways are widely recognized to stem from distinct metabolic abnormalities in hepatic glucose production and peripheral glucose uptake, respectively. We confirm this but go beyond to show that IFG and IGT lie on a continuum ranging from high hepatic insulin resistance and low peripheral insulin resistance to low hepatic resistance and high peripheral resistance. We show that IFG generally incurs IGT and IGT generally incurs IFG on the way to T2D, highlighting the difference between innate and acquired defects and the need to assess patients early to determine their underlying primary impairment and appropriately target therapy. We also consider other mechanisms, showing that IFG can result from impaired insulin secretion, that non-insulin-dependent glucose uptake can also mediate or interact with these pathways, and that impaired incretin signaling can accelerate T2D progression. We consider whether hyperinsulinemia can cause insulin resistance in addition to being a response to it and suggest that this is a minor effect.

Pediatric Extrapolation in Type 2 Diabetes: Future Implications of a Workshop

Extrapolation from adults to youth with type 2 diabetes (T2D) is challenged by differences in disease progression and manifestation. This manuscript presents the results of a mock-team workshop focused on examining the typical team-based decision process used to propose a pediatric development plan for T2D addressing the viability of extrapolation. The workshop was held at the American Society for Clinical Pharmacology and Therapeutics (ASCPT) in Orlando, Florida on March 21, 2018.

Model-inferred mechanisms of liver function from magnetic resonance imaging data: Validation and variation across a clinically relevant cohort

Estimation of liver function is important to monitor progression of chronic liver disease (CLD). A promising method is magnetic resonance imaging (MRI) combined with gadoxetate, a liver-specific contrast agent. For this method, we have previously developed a model for an average healthy human. Herein, we extended this model, by combining it with a patient-specific non-linear mixed-effects modeling framework. We validated the model by recruiting 100 patients with CLD of varying severity and etiologies. The model explained all MRI data and adequately predicted both timepoints saved for validation and gadoxetate concentrations in both plasma and biopsies. The validated model provides a new and deeper look into how the mechanisms of liver function vary across a wide variety of liver diseases. The basic mechanisms remain the same, but increasing fibrosis reduces uptake and increases excretion of gadoxetate. These mechanisms are shared across many liver functions and can now be estimated from standard clinical images.

Being able to accurately and reliably estimate liver function is important when monitoring the progression of patients with liver disease, as well as when identifying drug-induced liver injury during drug development. A promising method for quantifying liver function is to use magnetic resonance imaging combined with gadoxetate. Gadoxetate is a liver-specific contrast agent, which is taken up by the hepatocytes and excreted into the bile. We have previously developed a mechanistic model for gadoxetate dynamics using averaged data from healthy volunteers. In this work, we extended our model with a non-linear mixed-effects modeling framework to give patient-specific estimates of the gadoxetate transport-rates. We validated the model by recruiting 100 patients with liver disease, covering a range of severity and etiologies. All patients underwent an MRI-examination and provided both blood and liver biopsies. Our validated model provides a new and deeper look into how the mechanisms of liver function varies across a wide variety of liver diseases. The basic mechanisms remain the same, but increasing fibrosis reduces uptake and increases excretion of gadoxetate.

Quantitative Systems Pharmacology: An Exemplar Model-Building Workflow With Applications in Cardiovascular, Metabolic, and Oncology Drug Development

Quantitative systems pharmacology (QSP), a mechanistically oriented form of drug and disease modeling, seeks to address a diverse set of problems in the discovery and development of therapies. These problems bring a considerable amount of variability and uncertainty inherent in the nonclinical and clinical data. Likewise, the available modeling techniques and related software tools are manifold. Appropriately, the development, qualification, application, and impact of QSP models have been similarly varied. In this review, we describe the progressive maturation of a QSP modeling workflow: a necessary step for the efficient, reproducible development and qualification of QSP models, which themselves are highly iterative and evolutive. Furthermore, we describe three applications of QSP to impact drug development; one supporting new indications for an approved antidiabetic clinical asset through mechanistic hypothesis generation, one highlighting efficacy and safety differentiation within the sodium‐glucose cotransporter‐2 inhibitor drug class, and one enabling rational selection of immuno‐oncology drug combinations.

Toward Better Understanding of Insulin Therapy by Translation of a PK-PD Model to Visualize Insulin and Glucose Action Profiles

Insulin replacement therapy is a fundamental treatment for glycemic control for managing diabetes. The engineering of insulin analogues has focused on providing formulations with action profiles that mimic as closely as possible the pattern of physiological insulin secretion that normally occurs in healthy individuals without diabetes. Hence, it may be helpful to practitioners to visualize insulin concentration profiles and associated glucose action profiles. Expanding on a previous analysis that established a pharmacokinetic (PK) model to describe typical profiles of insulin concentration over time following subcutaneous administration of various insulin formulations, the goal of the current analysis was to link the PK model to an integrated glucose‐insulin (IGI) systems pharmacology model. After the pharmacokinetic‐pharmacodynamic (PK‐PD) model was qualified by comparing model predictions with clinical observations, it was used to project insulin (PK) and glucose (PD) profiles of common insulin regimens and dosing scenarios. The application of the PK‐PD model to clinical scenarios was further explored by incorporating the impact of several hypothetical factors together, such as changing the timing or frequency of administration in a multiple‐dosing regimen over the course of a day, administration of more than 1 insulin formulation, or insulin dosing adjusted for carbohydrates in meals. Visualizations of insulin and glucose profiles for commonly prescribed regimens could be rapidly generated by implementing the linked subcutaneous insulin PK‐IGI model using the R statistical program (version 3.4.4) and a contemporary web‐based interface, which could enhance clinical education on glycemic control with insulin therapy.

In vivo and in silico dynamics of the development of Metabolic Syndrome

The Metabolic Syndrome (MetS) is a complex, multifactorial disorder that develops slowly over time presenting itself with large differences among MetS patients. We applied a systems biology approach to describe and predict the onset and progressive development of MetS, in a study that combined in vivo and in silico models. A new data-driven, physiological model (MINGLeD: Model INtegrating Glucose and Lipid Dynamics) was developed, describing glucose, lipid and cholesterol metabolism. Since classic kinetic models cannot describe slowly progressing disorders, a simulation method (ADAPT) was used to describe longitudinal dynamics and to predict metabolic concentrations and fluxes. This approach yielded a novel model that can describe long-term MetS development and progression. This model was integrated with longitudinal in vivo data that was obtained from male APOE*3-Leiden.CETP mice fed a high-fat, high-cholesterol diet for three months and that developed MetS as reflected by classical symptoms including obesity and glucose intolerance. Two distinct subgroups were identified: those who developed dyslipidemia, and those who did not. The combination of MINGLeD with ADAPT could correctly predict both phenotypes, without making any prior assumptions about changes in kinetic rates or metabolic regulation. Modeling and flux trajectory analysis revealed that differences in liver fluxes and dietary cholesterol absorption could explain this occurrence of the two different phenotypes. In individual mice with dyslipidemia dietary cholesterol absorption and hepatic turnover of metabolites, including lipid fluxes, were higher compared to those without dyslipidemia. Predicted differences were also observed in gene expression data, and consistent with the emergence of insulin resistance and hepatic steatosis, two well-known MetS co-morbidities. Whereas MINGLeD specifically models the metabolic derangements underlying MetS, the simulation method ADAPT is generic and can be applied to other diseases where dynamic modeling and longitudinal data are available.

Impact of amyloid-beta changes on cognitive outcomes in Alzheimer's disease: analysis of clinical trials using a quantitative systems pharmacology model

BACKGROUND

Despite a tremendous amount of information on the role of amyloid in Alzheimer's disease (AD), almost all clinical trials testing this hypothesis have failed to generate clinically relevant cognitive effects.

METHODS

We present an advanced mechanism-based and biophysically realistic quantitative systems pharmacology computer model of an Alzheimer-type neuronal cortical network that has been calibrated with Alzheimer Disease Assessment Scale, cognitive subscale (ADAS-Cog) readouts from historical clinical trials and simulated the differential impact of amyloid-beta (Aβ40 and Aβ42) oligomers on glutamate and nicotinic neurotransmission.

RESULTS

Preclinical data suggest a beneficial effect of shorter Aβ forms within a limited dose range. Such a beneficial effect of Aβ40 on glutamate neurotransmission in human patients is absolutely necessary to reproduce clinical data on the ADAS-Cog in minimal cognitive impairment (MCI) patients with and without amyloid load, the effect of APOE genotype effect on the slope of the cognitive trajectory over time in placebo AD patients and higher sensitivity to cholinergic manipulation with scopolamine associated with higher Aβ in MCI subjects. We further derive a relationship between units of Aβ load in our model and the standard uptake value ratio from amyloid imaging. When introducing the documented clinical pharmacodynamic effects on Aβ levels for various amyloid-related clinical interventions in patients with low Aβ baseline, the platform predicts an overall significant worsening for passive vaccination with solanezumab, beta-secretase inhibitor verubecestat and gamma-secretase inhibitor semagacestat. In contrast, all three interventions improved cognition in subjects with moderate to high baseline Aβ levels, with verubecestat anticipated to have the greatest effect (around ADAS-Cog value 1.5 points), solanezumab the lowest (0.8 ADAS-Cog value points) and semagacestat in between. This could explain the success of many amyloid interventions in transgene animals with an artificial high level of Aβ, but not in AD patients with a large variability of amyloid loads.

CONCLUSIONS

If these predictions are confirmed in post-hoc analyses of failed clinical amyloid-modulating trials, one should question the rationale behind testing these interventions in early and prodromal subjects with low or zero amyloid load.

Cross-talks via mTORC2 can explain enhanced activation in response to insulin in diabetic patients

The molecular mechanisms of insulin resistance in Type 2 diabetes have been extensively studied in primary human adipocytes, and mathematical modelling has clarified the central role of attenuation of mammalian target of rapamycin (mTOR) complex 1 (mTORC1) activity in the diabetic state. Attenuation of mTORC1 in diabetes quells insulin-signalling network-wide, except for the mTOR in complex 2 (mTORC2)-catalysed phosphorylation of protein kinase B (PKB) at Ser⁴⁷³ (PKB-S473P), which is increased. This unique increase could potentially be explained by feedback and interbranch cross-talk signals. To examine if such mechanisms operate in adipocytes, we herein analysed data from an unbiased phosphoproteomic screen in 3T3-L1 adipocytes. Using a mathematical modelling approach, we showed that a negative signal from mTORC1-p70 S6 kinase (S6K) to rictor–mTORC2 in combination with a positive signal from PKB to SIN1–mTORC2 are compatible with the experimental data. This combined cross-branch signalling predicted an increased PKB-S473P in response to attenuation of mTORC1 – a distinguishing feature of the insulin resistant state in human adipocytes. This aspect of insulin signalling was then verified for our comprehensive model of insulin signalling in human adipocytes. Introduction of the cross-branch signals was compatible with all data for insulin signalling in human adipocytes, and the resulting model can explain all data network-wide, including the increased PKB-S473P in the diabetic state. Our approach was to first identify potential mechanisms in data from a phosphoproteomic screen in a cell line, and then verify such mechanisms in primary human cells, which demonstrates how an unbiased approach can support a direct knowledge-based study.

The promises of quantitative systems pharmacology modelling for drug development

Recent growth in annual new therapeutic entity (NTE) approvals by the U.S. Food and Drug Administration (FDA) suggests a positive trend in current research and development (R&D) output. Prior to this, the cost of each NTE was considered to be rising exponentially, with compound failure occurring mainly in clinical phases. Quantitative systems pharmacology (QSP) modelling, as an additional tool in the drug discovery arsenal, aims to further reduce NTE costs and improve drug development success. Through in silico mathematical modelling, QSP can simulate drug activity as perturbations in biological systems and thus understand the fundamental interactions which drive disease pathology, compound pharmacology and patient response. Here we review QSP, pharmacometrics and systems biology models with respect to the diseases covered as well as their clinical relevance and applications. Overall, the majority of modelling focus was aligned with the priority of drug-discovery and clinical trials. However, a few clinically important disease categories, such as Immune System Diseases and Respiratory Tract Diseases, were poorly covered by computational models. This suggests a possible disconnect between clinical and modelling agendas. As a standard element of the drug discovery pipeline the uptake of QSP might help to increase the efficiency of drug development across all therapeutic indications.