Effects of IL-1β-Blocking Therapies in Type 2 Diabetes Mellitus: A Quantitative Systems Pharmacology Modeling Approach to Explore Underlying Mechanisms

A data-driven computational model for obesity-driven diabetes onset and remission through weight loss

Obesity is a major risk factor for the development of type 2 diabetes (T2D), where a sustained weight loss may result in T2D remission in individuals with obesity. To design effective and feasible intervention strategies to prevent or reverse T2D, it is imperative to study the progression of T2D and remission together. Unfortunately, this is not possible through experimental and observational studies. To address this issue, we introduce a data-driven computational model and use human data to investigate the progression of T2D with obesity and remission through weight loss on the same timeline. We identify thresholds for the emergence of T2D and necessary conditions for remission. We explain why remission is only possible within a window of opportunity and the way that window depends on the progression history of T2D, individual's metabolic state, and calorie restrictions. These findings can help to optimize therapeutic intervention strategies for T2D prevention or treatment.

Applications of Quantitative Systems Pharmacology in Model-Informed Drug Discovery: Perspective on Impact and Opportunities

Quantitative systems pharmacology (QSP) approaches have been increasingly applied in the pharmaceutical since the landmark white paper published in 2011 by a National Institutes of Health working group brought attention to the discipline. In this perspective, we discuss QSP in the context of other modeling approaches and highlight the impact of QSP across various stages of drug development and therapeutic areas. We discuss challenges to the field as well as future opportunities.

A novel fast-slow model of diabetes progression: Insights into mechanisms of response to the interventions in the Diabetes Prevention Program

Several models for the long-term development of T2DM already exist, focusing on the dynamics of the interaction between glycemia, insulinemia and β-cell mass. Current models consider representative (fasting or daily average) glycemia and insulinemia as characterizing the compensation state of the subject at some instant in slow time. This implies that only these representative levels can be followed through time and that the role of fast glycemic oscillations is neglected. An improved model (DPM15) for the long-term progression of T2DM is proposed, introducing separate peripheral and hepatic (liver and kidney) insulin actions. The DPM15 model no longer uses near-equilibrium approximation to separate fast and slow time scales, but rather describes, at each step in slow time, a complete day in the life of the virtual subject in fast time. The model can thus represent both fasting and postprandial glycemic levels and describe the effect of interventions acting on insulin-enhanced tissue glucose disposal or on insulin-inhibited hepatic glucose output, as well as on insulin secretion and β-cell replicating ability. The model can simulate long-term variations of commonly used clinical indices (HOMA-B, HOMA-IR, insulinogenic index) as well as of Oral Glucose Tolerance or Euglycemic Hyperinsulinemic Clamp test results. The model has been calibrated against observational data from the Diabetes Prevention Program study: it shows good adaptation to observations as a function of very plausible values of the parameters describing the effect of such interventions as Placebo, Intensive LifeStyle and Metformin administration.

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.

A SIRT1 agonist reduces cognitive decline in type 2 diabetic rats through antioxidative and anti‑inflammatory mechanisms

Sirtuin 1 (SIRT1) is an NAD+‑dependent protein deacetylase that is involved in cell differentiation, aging, apoptosis, physiological rhythms, metabolic regulation, oxidative stress and numerous other important biological processes. In the present study, the ability of a sirtuin‑1 (SIRT1) agonist, SRT1720, to reduce cognitive decline in type 2 diabetes mellitus (T2DM) was investigated. Streptozotocin‑induced male Sprague‑Dawley rats were used to establish a T2DM model and the protective effect of SRT1720 and its underlying mechanisms were investigated. Body weight and fasting blood glucose (FBG) were recorded and cognitive function was measured with the Morris water maze. Levels of oxidative stress, inflammation, caspase‑3 activity and nuclear factor κB (NF‑κB) mRNA expression were detected with a series of commercial assay kits and reverse transcription‑quantitative polymerase chain reaction, respectively. Western blot analysis was performed to determine the protein expression of NF‑κB, endothelial nitric oxide synthase (eNOS), peroxisome proliferator‑activated receptor γ (PPARγ), AMP‑activated protein kinase (AMPK), heat shock 70 kDa protein (HSP70), SIRT1, nuclear factor erythroid 2‑related factor 2 (Nrf2) and heme oxygenase 1 (HO‑1). The results revealed that SRT1720 significantly increased body weight, decreased FBG, improved cognitive function and reduced the levels of proteins associated with oxidative stress and inflammation damage in T2DM rats. Additionally, SRT1720 significantly decreased NF‑κB p65 mRNA expression and increased eNOS and PPARγ expression. SRT1720 significantly reduced caspase‑3 activity and HSP70 protein expression, and increased p‑AMPK, SIRT1, Nrf2 and HO‑1 protein expression. Collectively, the results indicate that SRT1720 may reduce cognitive decline in T2DM rats through antioxidative and anti‑inflammatory action via NF‑κB and AMPK‑dependent mechanisms.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

There are excellent methods to mathematically model time-resolved biological data on a small scale using accurate mechanistic models. Despite the rapidly increasing availability of such data, mechanistic models have not been applied on a genome-wide level due to excessive runtimes and the non-identifiability of model parameters. However, genome-wide, mechanistic models could potentially answer key clinical questions, such as finding the best drug combinations to induce an expression change from a disease to a healthy state. We present LASSIM, which is a toolbox built to infer parameters within mechanistic models on a genomic scale. This is made possible due to a property shared across biological systems, namely the existence of a subset of master regulators, here denoted the core system. The introduction of a core system of genes simplifies the network inference into small solvable sub-problems, and implies that all main regulatory actions on peripheral genes come from a small set of regulator genes. This separation allows substantial parts of computations to be solved in parallel, i.e. permitting the use of a computer cluster, which substantially reduces computation time.

Advanced Boolean modeling of biological networks applied to systems pharmacology

Motivation

Literature on complex diseases is abundant but not always quantitative. Many molecular pathways are qualitatively well described but this information cannot be used in traditional quantitative mathematical models employed in drug development. Tools for analysis of discrete networks are useful to capture the available information in the literature but have not been efficiently integrated by the pharmaceutical industry. We propose an expansion of the usual analysis of discrete networks that facilitates the identification/validation of therapeutic targets.

Results

In this article, we propose a methodology to perform Boolean modeling of Systems Biology/Pharmacology networks by using SPIDDOR (Systems Pharmacology for effIcient Drug Development On R) R package. The resulting models can be used to analyze the dynamics of signaling networks associated to diseases to predict the pathogenesis mechanisms and identify potential therapeutic targets.

Availability and Implementation

The source code is available at https://github.com/SPIDDOR/SPIDDOR .

Contact

itzirurzun@alumni.unav.es , itroconiz@unav.es.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

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

We are currently in the middle of a major shift in biomedical research: unprecedented and rapidly growing amounts of data may be obtained today, from in vitro, in vivo and clinical studies, at molecular, physiological and clinical levels. To make use of these large-scale, multi-level datasets, corresponding multi-level mathematical models are needed, i.e. models that simultaneously capture multiple layers of the biological, physiological and disease-level organization (also referred to as quantitative systems pharmacology-QSP-models). However, today's multi-level models are not yet embedded in end-usage applications, neither in drug research and development nor in the clinic. Given the expectations and claims made historically, this seemingly slow adoption may seem surprising. Therefore, we herein consider a specific example-type 2 diabetes-and critically review the current status and identify key remaining steps for these models to become mainstream in the future. This overview reveals how, today, we may use models to ask scientific questions concerning, e.g., the cellular origin of insulin resistance, and how this translates to the whole-body level and short-term meal responses. However, before these multi-level models can become truly useful, they need to be linked with the capabilities of other important existing models, in order to make them 'personalized' (e.g. specific to certain patient phenotypes) and capable of describing long-term disease progression. To be useful in drug development, it is also critical that the developed models and their underlying data and assumptions are easily accessible. For clinical end-usage, in addition, model links to decision-support systems combined with the engagement of other disciplines are needed to create user-friendly and cost-efficient software packages.

BioModels: Content, Features, Functionality, and Use

BioModels is a reference repository hosting mathematical models that describe the dynamic interactions of biological components at various scales. The resource provides access to over 1,200 models described in literature and over 140,000 models automatically generated from pathway resources. Most model components are cross-linked with external resources to facilitate interoperability. A large proportion of models are manually curated to ensure reproducibility of simulation results. This tutorial presents BioModels' content, features, functionality, and usage.

Anti-inflammatory agents to treat or prevent type 2 diabetes, metabolic syndrome and cardiovascular disease

INTRODUCTION

There is a growing body of evidence to suggest that chronic silent inflammation is a key feature in abdominal obesity, metabolic syndrome, type 2 diabetes (T2DM) and cardiovascular disease (CVD). These observations suggest that pharmacological strategies, which reduce inflammation, may be therapeutically useful in treating obesity, type 2 diabetes and associated CVD.

AREA COVERED

The article covers novel strategies, using either small molecules or monoclonal antibodies. These strategies include: approaches targeting IKK-b-NF-kB (salicylates, salsalate), TNF-α (etanercept, infliximab, adalimumab), IL-1β (anakinra, canakinumab) and IL-6 (tocilizumab), AMP-activated protein kinase activators, sirtuin-1 activators, mammalian target of rapamycin inhibitors and C-C motif chemokine receptor 2 antagonists.

EXPERT OPINION

The available data supports the concept that targeting inflammation improves insulin sensitivity and β-cell function; it also ameliorates glucose control in insulin-resistant patients with inflammatory rheumatoid diseases as well in patients with metabolic syndrome or T2DM. Although promising, the observed metabolic effects remain rather modest in most clinical trials. The potential use of combined anti-inflammatory agents targeting both insulin resistance and insulin secretion appears appealing but remains unexplored. Large-scale prospective clinical trials are underway to investigate the safety and efficacy of different anti-inflammatory drugs. Further evidence is needed to support the concept that targeting inflammation pathways may represent a valuable option to tackle the cardiometabolic complications of obesity.